如何使用 python 坐标计算机翼厚度

how to calculate airfoil thickness using coordinate with python

到目前为止,我有一个机翼文件 S823,它具有以下坐标。

  1.000000  0.000000
  0.996260  0.000570
  0.985390  0.002590
  0.968200  0.006440
  0.945590  0.012020
  0.918230  0.018900
  0.886490  0.026690
  0.850640  0.035270
  0.811260  0.044730
  0.769210  0.054910
  0.725320  0.065490
  0.680480  0.075890
  0.635300  0.085260
  0.589910  0.092770
  0.544140  0.097990
  0.497850  0.100890
  0.451090  0.101740
  0.404250  0.100920
  0.357860  0.098640
  0.312450  0.095050
  0.268520  0.090260
  0.226560  0.084410
  0.187020  0.077610
  0.150340  0.070010
  0.116930  0.061730
  0.087140  0.052940
  0.061300  0.043770
  0.039710  0.034410
  0.022590  0.025040
  0.010170  0.015870
  0.002610  0.007130
  0.000000  0.000000
  0.000020 -0.000570
  0.003550 -0.007470
  0.013320 -0.014500
  0.028210 -0.021620
  0.047970 -0.028480
  0.072370 -0.034920
  0.101170 -0.040770
  0.134060 -0.045950
  0.170710 -0.050360
  0.210750 -0.053940
  0.253750 -0.056640
  0.299260 -0.058430
  0.346770 -0.059270
  0.395780 -0.059140
  0.445720 -0.058010
  0.496040 -0.055800
  0.546170 -0.052370
  0.595730 -0.047450
  0.644760 -0.040840
  0.693550 -0.033050
  0.741760 -0.025080
  0.788480 -0.017710
  0.832770 -0.011380
  0.873640 -0.006380
  0.910130 -0.002830
  0.941340 -0.000660
  0.966490  0.000320
  0.984930  0.000470
  0.996200  0.000190
  1.000000  0.000000

坐标从后缘开始到前缘然后再回到后缘。

由于x坐标与x轴上下的点不相同(坐标从前缘移回后缘时没有x=0.002610)。我很难直接减去坐标来获得我的最大厚度。

我想知道是否有一种方法可以平滑机翼坐标并获得最大厚度?

是的,您通常希望使用诸如三次 B 样条之类的工具来平滑曲线。请注意,B 样条曲线通常不会精确地通过您指定的每个点。沿着 top/bottom 边,你通常有一个足够大的曲线,这不是什么问题。

在前沿附近,您经常需要填充一些点才能使事情正常进行。不幸的是,在大多数情况下(可能包括这个,只有 62 分)很难确定您要填写的那些分应该放在哪里。

当我不得不这样做时,我使用 xfoil 来模拟机翼,并使用 LE 直到我得到拖动 bucket/LD/Cp 曲线来匹配(至少漂亮密切)与那些给设计。但这并不总是那么容易。