CN111539072A - 一种结冰风洞试验冰形冰角量化评估方法和系统 - Google Patents

一种结冰风洞试验冰形冰角量化评估方法和系统 Download PDF

Info

Publication number
CN111539072A
CN111539072A CN202010649538.0A CN202010649538A CN111539072A CN 111539072 A CN111539072 A CN 111539072A CN 202010649538 A CN202010649538 A CN 202010649538A CN 111539072 A CN111539072 A CN 111539072A
Authority
CN
China
Prior art keywords
ice
angle
shaped
airfoil
thickness
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202010649538.0A
Other languages
English (en)
Other versions
CN111539072B (zh
Inventor
郭向东
赵献礼
陈圆圆
白云卯
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Low Speed Aerodynamics Institute of China Aerodynamics Research and Development Center
Original Assignee
Low Speed Aerodynamics Institute of China Aerodynamics Research and Development Center
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Low Speed Aerodynamics Institute of China Aerodynamics Research and Development Center filed Critical Low Speed Aerodynamics Institute of China Aerodynamics Research and Development Center
Priority to CN202010649538.0A priority Critical patent/CN111539072B/zh
Publication of CN111539072A publication Critical patent/CN111539072A/zh
Application granted granted Critical
Publication of CN111539072B publication Critical patent/CN111539072B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/15—Vehicle, aircraft or watercraft design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Automation & Control Theory (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

本发明适用于航空领域，提供了一种结冰风洞试验冰形冰角量化评估方法和系统，首先针对结冰风洞试验中获取的冰形轮廓线，采用冰形几何轮廓线数字化方法，获得冰形轮廓线离散坐标点，进而针对冰形离散点，建立冰形类型判定准则对冰形类型进行判定，并根据不同的冰形类型计算对应的冰形特征参数。本发明的评估方法和系统适用于大规模的试验冰形冰角量化评估任务，计算方法中不涉及复杂的数学公式和耗时的数值迭代过程，可以快速高效的完成冰形冰角量化评估。

Description

一种结冰风洞试验冰形冰角量化评估方法和系统
技术领域
本发明属于航空领域，尤其涉及一种结冰风洞试验冰形冰角量化评估方法和系统。
背景技术
当飞机在云层中飞行时，云层中的过冷水滴（即温度低于冰点的液态水滴）会不断撞击飞机迎风面，导致其表面发生结冰现象。飞机结冰广泛存在于飞行实践中，并严重威胁飞行安全。结冰风洞是开展飞机结冰研究及飞机部件防除冰系统验证的重要地面试验设备，其在飞机结冰适航审定中扮演着重要角色。结冰风洞试验中，积冰外形几何轮廓是最重要的结冰试验结果，是临界冰形选取、防除冰系统性能评估、冰模研制等试验目标的关键参数。
结冰风洞试验冰形主要包括流向冰和角状冰两大类，其中流向冰对应霜冰结冰条件（环境温度较低，液态水含量较低），而角状冰则普遍在明冰和混合冰结冰条件（环境温度较高，液态水含量较高）下形成。这两种冰形都具有明显的角状特征，其中流向冰普遍为单角结构，而角状冰普遍为双角结构。冰角是积冰影响飞机气动特性的重要特征，是冰形定量化评估的关键元素。
目前，结冰风洞普遍采用冰形绘制技术实现积冰二维几何轮廓的测量。该技术利用热刀切割积冰主体，使用铅笔在网格纸上描绘积冰切缝处的二维冰形轮廓，获得二维冰形轮廓线网格图纸。针对绘制的二维冰形网格图纸，通常采用视觉定性与简单测量相结合的评估方法，判断冰形类型和位置，计算冰形冰角的简单特征。这种方法显著受人为主观因素的影响，获得的冰角特征结果可靠性、适用性较差，尤其对于大量冰形试验结果，该方法处理效率较低，无法满足试验周期的要求。
针对结冰风洞试验冰形定量化处理问题，国内外结冰研究机构及学者均开展了一些研究，发展了一些量化评估方法，但是这些方法在实际工程应用中仍存在一些缺陷，缺乏系统、高效的结冰风洞试验冰形量化评估方法。国外，美国NASA Glenn中心为验证结冰计算程序和结冰相似准则，率先对结冰风洞试验冰形定量化处理问题开展了研究，发展了试验冰形与计算冰形定量化比较评估方法，形成了THICK冰形特征参数计算程序，初步实现了冰形特征参数的定量计算，但是该方法无法区分流向冰和角状冰，同时角状冰上下冰角范围仅简单的通过前缘点进行区分，普适性较差。
而在国内，相关研究较少，近年来，随着大型结冰风洞-3米×2米结冰风洞的建成，在结冰试验的急迫需求牵引下，一些学者针对冰形相似度及简化问题开展了初步研究，发展了一些简化的量化评估方法，但是这些研究主要关注冰形特征参数的深度分析，缺乏试验冰形特征参数提取与计算的详细研究，无法完全满足结冰风洞试验中系统、高效、适用性强的冰形定量化评估需求。
发明内容
为解决结冰风洞试验冰形冰角量化评估问题，本发明提供一种结冰风洞试验冰形冰角量化评估方法和系统，针对离散化的冰形和翼型数字点进行计算，计算中不涉及复杂的数学公式和数值迭代过程，可以快速高效的完成冰形特征参数的定量计算，适用于大规模的试验冰形冰角量化评估任务。
本发明是这样实现的， 一种结冰风洞试验冰形冰角量化评估方法，包括如下步骤：
（1）数字化冰形几何轮廓线
对试验中绘制的冰形轮廓线网格图纸进行数字化扫描，定义数字化坐标系，提取翼型和冰形轮廓线离散坐标点，绘制冰形轮廓线图；
（2）计算积冰厚度分布
针对数字化的冰形和翼型离散点，计算积冰厚度分布；
（3）根据翼型试验攻角α，确定冰形分界点I D
3-1. 计算各冰形离散点i对应的积冰厚度线与X轴的夹角θ；
3-2.将各冰形离散点处θ与试验攻角α对比，当角度差绝对值小于设定值时，即选择该冰形离散点为冰形分界点I D ；
（4）判断冰形类型和冰角区域
4-1. 以冰形分界点I D 为界，将冰形离散点分为上冰形区和下冰形区，计算上下冰形区内的最大积冰厚度T u 和T l ；
Figure 201428DEST_PATH_IMAGE001
式中，N IS 表示冰形总离散点个数；同时记录上下冰形区最大厚度对应的冰形离散点序号I u 和I l ，以及翼型投影点序号J u 和J l ；
4-2. 将分界点编号I D 与I u 、I l 进行比较：
如果I D ＝I u 或I l ，则冰形为流向冰形，仅有单冰角结构；
如果I D ≠I u 或I l ，则冰形为角状冰形，包括上下两个冰角，其中I D 同时为上下冰角区分界点；
（5）根据确定的冰形类型，进行冰角特征参数计算。
进一步地，步骤（2）中，计算冰形各离散点与翼型所有离散点间的距离，选取其中最短距离作为冰形离散点到翼型的积冰厚度T i ，同时记录对应的翼型投影点序号J i ，
Figure 56251DEST_PATH_IMAGE003
式中，X IS,i 和Y IS,i 表示冰形离散点i对应的横纵坐标，X A,j 和Y A,j 表示翼型离散点j对应的横纵坐标，N A 表示翼型总离散点个数。
进一步地，步骤3-1中，夹角θ根据下式进行计算：
Figure 835988DEST_PATH_IMAGE004
式中，X IS,i 和Y IS,i 表示冰形离散点i对应的横纵坐标，T i 为积冰厚度，X A,Ji 和Y AJi 表示翼型投影点J i 对应的横纵坐标。
进一步地，步骤5中，计算以下冰角特征参数：
流向冰形的冰角特征参数包括：冰角厚度T s 、冰角角度θ s ；
角状冰形的冰角特征参数包括：上冰角厚度T du 、下冰角厚度T dl 、上冰角角度θ du 、下冰角角度θ dl 。
流向冰形的冰角特征参数采用如下方法计算：
5-1. 计算冰角厚度T s ：
Figure 168881DEST_PATH_IMAGE005
式中，T i 为冰形离散点i对应的积冰厚度；
记录冰角厚度T s 对应的冰形离散点序号I s ；
5-2. 计算冰角角度θ s ：
Figure 983253DEST_PATH_IMAGE006
Figure 743399DEST_PATH_IMAGE007
Figure 10432DEST_PATH_IMAGE008
式中，X IS,Is 和Y IS,Is 表示冰形离散点I s 对应的横纵坐标。
角状冰形的冰角特征参数采用如下方法计算：
5-3. 计算上下冰角厚度T du 和T dl ：
Figure 147015DEST_PATH_IMAGE009
记录上下冰角厚度T du 和T dl 对应的冰形离散点序号I du 和I dl ；
5-4. 计算上下冰角角度θ du 和θ dl ：
Figure 815894DEST_PATH_IMAGE010
Figure 809258DEST_PATH_IMAGE011
Figure 501270DEST_PATH_IMAGE012
Figure 238282DEST_PATH_IMAGE013
Figure 699350DEST_PATH_IMAGE014
Figure 129195DEST_PATH_IMAGE015
式中，X IS,Idl 和Y IS,Idl 表示冰形离散点I dl 对应的横纵坐标; X IS,Idu 和Y IS,Idu 表示冰形离散点I du 对应的横纵坐标。
进一步地，定义数字化坐标系的具体方法为：将翼型前缘点设为坐标原点，X轴沿翼型弦线，指向翼型尾部，Y轴从翼型下翼面指向上翼面。
本发明还提供一种结冰风洞试验冰形冰角量化评估系统，包括扫描仪、坐标系设定模块和冰角量化评估模块，其中，
扫描仪用于扫描试验中绘制的冰形轮廓线网格图纸；
坐标系设定模块用于定义数字化坐标系；
冰角量化评估模块用于执行前述的结冰风洞试验冰形冰角量化评估方法。
本发明还提供一种存储介质，存储了用于执行如上结冰风洞试验冰形冰角量化评估方法的计算机程序。
本发明针对结冰风洞试验中获取的冰形轮廓线，采用冰形几何轮廓线数字化方法，获得冰形轮廓线离散坐标点，进而针对冰形离散点，建立冰形类型判定准则，再根据冰形类型计算冰角特征参数，计算中不涉及复杂的数学公式和数值迭代过程，可以快速高效的完成冰形特征参数的定量计算，适用于大规模的试验冰形冰角量化评估任务。
同时，采用本发明的评估系统，有助于实现参数自动化计算，既能提高计算精度又能提高计算效率。
本发明可以有效解决结冰风洞试验冰形冰角的量化评估问题，实现冰形的定量化后处理，为临界冰形选取、防除冰系统性能评估、冰模研制等多种试验需求提供技术支撑。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对本发明实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面所描述的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明结冰风洞试验冰形冰角量化评估流程图；
图2是本发明实施例的冰形扫描图片；
图3是本发明所定义的数字化坐标系；
图4是本发明实施例的冰形数字化轮廓线图；
图5是本发明实施例的积冰厚度分布图；
图6是本发明实施例的冰形分界点示意图；
图7a是本发明实施例的流向冰形的判定示意图；
图7b是本发明实施例的角状冰形的判定示意图；
图8是流向冰冰角特征参数定义；
图9是角状冰冰角特征参数定义；
图10是另一实施例流向冰冰角特征参数计算结果示意图；
图11是又一实施例角状冰冰角特征参数计算结果示意图。
具体实施方式
下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，旨在用于解释本发明，而不能理解为对本发明的限制。
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。
实施例1
一种结冰风洞试验冰形冰角量化评估方法，如图1所示，包括如下步骤：
（1）数字化冰形几何轮廓线
对试验中绘制的冰形轮廓线网格图纸进行数字化扫描，定义数字化坐标系，提取翼型和冰形轮廓线离散坐标点，绘制冰形轮廓线图；
具体地，
1-1. 利用大幅面扫描仪，扫描冰形轮廓线网格图纸，输出数字化冰形轮廓线图片如图2所示；
1-2. 定义数字化坐标系，如图3所示，将翼型前缘点设为坐标原点，X轴沿翼型弦线，指向翼型尾部，Y轴从翼型下翼面（压力面）指向上翼面（吸力面）；
1-3. 利用图片数字化软件对冰形图片进行数字化处理，提取有序的翼型和冰形轮廓线离散化坐标点；
1-3-1. 按照步骤1-2定义的坐标系，提取翼型轮廓线离散化坐标点，提取顺序从翼型最上部（Y坐标最大值处）沿翼型表面到翼型最下部（Y坐标最小值处）。
1-3-2. 提取冰形轮廓线离散化坐标点，提取顺序从翼型上翼面冰形轮廓线起点出发，沿冰形轮廓线轨迹，到下翼面冰形轮廓线终点结束。对于冰形轮廓线不连续的间断区域，以翼型轮廓线离散点代替。
提取数字点时，应尽量减小数字点间距，以最大程度地保留冰形和翼型细节信息。
1-4.以模型弦长为特征长度，对提取的冰形轮廓线离散化坐标点进行无量纲化处理，绘制无量纲化的冰形轮廓线图，如图4所示。流向冰和角状冰两种典型冰形轮廓线如图7a和7b所示。
（2）计算积冰厚度分布
针对数字化的冰形和翼型离散点，计算积冰厚度分布，具体步骤为：
以翼型前缘点为原点，计算翼型各离散点j沿翼型表面到原点的距离Sj，其中上翼面处Sj为正，下翼面处为负，j=1对应上翼面最大S处。Sj表示为：
Figure 370820DEST_PATH_IMAGE016
j=1
Figure 580697DEST_PATH_IMAGE017
式中，X A,k 和Y A,k 表示翼型离散点k对应的横纵坐标，No表示翼型前缘点对应的翼型离散点数组序号。
计算冰形各离散点与翼型所有离散点间的距离，选取其中最短距离作为冰形离散点到翼型的积冰厚度T i ，同时记录对应的翼型投影点序号J i ，
Figure 224168DEST_PATH_IMAGE018
式中，X IS,i 和Y IS,i 表示冰形离散点i对应的横纵坐标，X A,j 和Y A,j 表示翼型离散点j对应的横纵坐标，N A 表示翼型总离散点个数。
以冰形离散点的积冰厚度T为纵坐标，以对应翼型离散点翼面距离S为横坐标，绘制积冰厚度分布图，如图5所示。
（3）根据翼型试验攻角α，确定冰形分界点I D
3-1. 计算各冰形离散点i对应的积冰厚度线与X轴的夹角θ；
Figure 497018DEST_PATH_IMAGE004
式中，X IS,i 和Y IS,i 表示冰形离散点i对应的横纵坐标，T i 为积冰厚度，X A,Ji 和Y AJi 表示翼型投影点J i 对应的横纵坐标。
3-2.将各冰形离散点处θ与试验攻角α对比，当角度差绝对值小于设定值时，即选择该冰形离散点为冰形分界点I D ，如图6所示。
由于计算得到的θ不一定与α完全相等，一般会存在一个很小的偏差，因此设定一个误差范围，即上述设定值，该设定值的取值根据结冰风洞角度控制精度来确定，此处将设定值取值为0.2°。
试验攻角α为实验时设定的实验参数，不同试验取值不同。
（4）判断冰形类型和冰角区域
4-1. 以冰形分界点I D 为界，将冰形离散点分为上冰形区和下冰形区，计算上下冰形区内的最大积冰厚度T u 和T l ；
Figure 225939DEST_PATH_IMAGE001
式中，N IS 表示冰形总离散点个数；同时记录上下冰形区最大厚度对应的冰形离散点序号I u 和I l ，以及翼型投影点序号J u 和J l ；
4-2. 将分界点编号I D 与I u 、I l 进行比较：
如果I D ＝I u 或I l ，则冰形为流向冰形，仅有单冰角结构，如图7a所示；
如果I D ≠I u 或I l ，则冰形为角状冰形，包括上下两个冰角，其中I D 同时为上下冰角区分界点，如图7b所示；
（5）根据确定的冰形类型，进行冰角特征参数计算。
根据冰形特征，典型流向冰形的冰角特征参数包括：冰角厚度T s 、冰角角度θ s ，如图8所示，分别定义为：冰角厚度T s ，定义为冰角尖部到翼型表面的最短距离，即翼型前缘最大积冰厚度；冰角角度θ s ，定义为冰角尖部与坐标系原点间连线与坐标系X轴的夹角，角度按逆时针方向增大。
流向冰形的冰角特征参数采用如下方法计算：
1.计算冰角厚度T s ：
Figure 242437DEST_PATH_IMAGE005
式中，T i 为冰形离散点i对应的积冰厚度；
记录冰角厚度T s 对应的冰形离散点序号I s ；
2.计算冰角角度θ s ：
Figure 740414DEST_PATH_IMAGE006
Figure 449744DEST_PATH_IMAGE007
Figure 400383DEST_PATH_IMAGE008
式中，X IS,Is 和Y IS,Is 表示冰形离散点I s 对应的横纵坐标。
典型角状冰形的冰角参数，如图9所示，包括上冰角厚度T du 、下冰角厚度T dl 、上冰角角度θ du 、下冰角角度θ dl 四个特征参数：上冰角厚度T du ，定义为上冰角尖部到翼型表面的最短距离；下冰角厚度T dl ，定义为下冰角尖部到翼型表面的最短距离；上冰角角度θ du ，定义为上冰角尖部与坐标系原点间连线与坐标系X轴的夹角，角度按逆时针方向增大；下冰角角度θ dl ，定义为下冰角尖部与坐标系原点间连线与坐标系X轴的夹角，角度按逆时针方向增大。
角状冰形的冰角特征参数包括：上冰角厚度T du 、下冰角厚度T dl 、上冰角角度θ du 、下冰角角度θ dl 。
角状冰形的冰角特征参数采用如下方法计算：
3. 计算上下冰角厚度T du 和T dl ：
Figure 220571DEST_PATH_IMAGE009
记录上下冰角厚度T du 和T dl 对应的冰形离散点序号I du 和I dl ；
4. 计算上下冰角角度θ du 和θ dl ：
Figure 573055DEST_PATH_IMAGE010
Figure 515603DEST_PATH_IMAGE011
Figure 891221DEST_PATH_IMAGE012
Figure 311838DEST_PATH_IMAGE013
Figure 722091DEST_PATH_IMAGE014
Figure 835540DEST_PATH_IMAGE015
式中，X IS,Idl 和Y IS,Idl 表示冰形离散点I dl 对应的横纵坐标; X IS,Idu 和Y IS,Idu 表示冰形离散点I du 对应的横纵坐标。
实施例2
采用本发明的一种结冰风洞试验冰形冰角量化评估方法，计算了流向冰形的冰角特征参数，如图10所示，该冰形冰角厚度T s 和角度θ s 分别为0.033和199°。
实施例3
采用本发明的一种结冰风洞试验冰形冰角量化评估方法，计算了角状冰形的冰角特征参数，如图11所示，中上冰角厚度T du 和角度θ du 分别为0.11和147°，下冰角厚度T dl 和角度θ dl 分别为0.102和236°。
结果表明，本发明可以准确的提取并计算典型冰形冰角特征参数，可以有效解决结冰风洞试验冰形冰角的定量化评估问题，为临界冰形选取、防除冰系统性能评估、冰模研制等试验目标提供技术支撑。
实施例4
一种结冰风洞试验冰形冰角量化评估系统，包括扫描仪、坐标系设定模块和冰角量化评估模块，其中，扫描仪用于扫描试验中绘制的冰形轮廓线网格图纸，输出数字化冰形轮廓线图片；
坐标系设定模块用于定义数字化坐标系；
冰角量化评估模块用于执行本发明的结冰风洞试验冰形冰角量化评估方法。
使用时，先用扫描仪扫描试验中绘制的冰形轮廓线网格图纸，输出数字化冰形轮廓线图片；实验人员手动定义数字化坐标，或者评估系统自动定义数字化坐标系，由实验人员手动确认；系统自动提取翼型和冰形轮廓线离散化坐标点，并按照本发明的结冰风洞试验冰形冰角量化评估方法自动进行冰形类型判定，再根据判定的冰形类型计算冰角特征参数并输出。由此大大减轻了实验人员的工作量，提高了试验效率和计算的准确性。
本发明还提供一种存储介质，存储了用于执行如上结冰风洞试验冰形冰角量化评估方法的计算机程序。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (9)

1.一种结冰风洞试验冰形冰角量化评估方法，其特征在于，包括如下步骤：
（1）数字化冰形几何轮廓线
对试验中绘制的冰形轮廓线网格图纸进行数字化扫描，定义数字化坐标系，提取翼型和冰形轮廓线离散坐标点，绘制冰形轮廓线图；
（2）计算积冰厚度分布
针对数字化的冰形和翼型离散点，计算积冰厚度分布；
（3）根据翼型试验攻角α，确定冰形分界点I D
3-1. 计算各冰形离散点i对应的积冰厚度线与X轴的夹角θ；
3-2.将各冰形离散点处θ与试验攻角α比，当角度差绝对值小于设定值时，即选择该冰形离散点为冰形分界点I D ；
（4）判断冰形类型和冰角区域
4-1. 以冰形分界点I D 为界，将冰形离散点分为上冰形区和下冰形区，计算上下冰形区内的最大积冰厚度T u 和T l ；
Figure 578230DEST_PATH_IMAGE001
式中，N IS 表示冰形总离散点个数；同时记录上下冰形区最大厚度对应的冰形离散点序号I u 和I l ，以及翼型投影点序号J u 和J l ；
4-2. 将分界点I D 与I u 、I l 进行比较：
如果I D ＝I u 或I l ，则冰形为流向冰形；
如果I D ≠I u 或I l ，则冰形为角状冰形，其中I D 同时为上下冰角区分界点；
（5）根据确定的冰形类型，进行冰角特征参数计算。
2.如权利要求1所述的一种结冰风洞试验冰形冰角量化评估方法，其特征在于，所述步骤（2）中，计算冰形各离散点与翼型所有离散点间的距离，选取其中最短距离作为冰形离散点到翼型的积冰厚度T i ，同时记录对应的翼型投影点序号J i ，
Figure 819856DEST_PATH_IMAGE002
式中，X IS,i 和Y IS,i 表示冰形离散点i对应的横纵坐标，X A,j 和Y A,j 表示翼型离散点j对应的横纵坐标，N A 表示翼型总离散点个数。
3.如权利要求1或2所述的一种结冰风洞试验冰形冰角量化评估方法，其特征在于，所述步骤3-1中，夹角θ根据下式进行计算：
Figure 32662DEST_PATH_IMAGE003
式中，X IS,i 和Y IS,i 表示冰形离散点i对应的横纵坐标，T i 为积冰厚度，X A,Ji 和Y AJi 表示翼型投影点J i 对应的横纵坐标。
4.如权利要求1所述的一种结冰风洞试验冰形冰角量化评估方法，其特征在于，所述的步骤（5）中，计算以下冰角特征参数：
流向冰形的冰角特征参数包括：冰角厚度T s 、冰角角度θ s ；
角状冰形的冰角特征参数包括：上冰角厚度T du 、下冰角厚度T dl 、上冰角角度θ du 、下冰角角度θ dl 。
5.如权利要求4所述的一种结冰风洞试验冰形冰角量化评估方法，其特征在于，所述流向冰形的冰角特征参数采用如下方法计算：
5-1. 计算冰角厚度T s ：
Figure 676133DEST_PATH_IMAGE004
式中，T i 为冰形离散点i对应的积冰厚度；
记录冰角厚度T s 对应的冰形离散点序号I s ；
5-2. 计算冰角角度θ s ：
Figure 447518DEST_PATH_IMAGE005
Figure 176440DEST_PATH_IMAGE006
Figure 192937DEST_PATH_IMAGE007
式中，X IS,Is 和Y IS,Is 表示冰形离散点I s 对应的横纵坐标。
6.如权利要求4所述的一种结冰风洞试验冰形冰角量化评估方法，其特征在于，所述角状冰形的冰角特征参数采用如下方法计算：
5-3. 计算上下冰角厚度T du 和T dl ：
Figure 690915DEST_PATH_IMAGE008
记录上下冰角厚度T du 和T dl 对应的冰形离散点I du 和I dl ；
5-4. 计算上下冰角角度θ du 和θ dl ：
Figure 400245DEST_PATH_IMAGE009
Figure 350883DEST_PATH_IMAGE010
Figure 171072DEST_PATH_IMAGE011
Figure 523555DEST_PATH_IMAGE012
Figure 902322DEST_PATH_IMAGE013
Figure 340257DEST_PATH_IMAGE014
式中，X IS,Idl 和Y IS,Idl 表示冰形离散点I dl 对应的横纵坐标; X IS,Idu 和Y IS,Idu 表示冰形离散点I du 对应的横纵坐标。
7.如权利要求1所述的一种结冰风洞试验冰形冰角量化评估方法，其特征在于，所述的定义数字化坐标系的具体方法为：将翼型前缘点设为坐标原点，X轴沿翼型弦线，指向翼型尾部，Y轴从翼型下翼面指向上翼面。
8.一种结冰风洞试验冰形冰角量化评估系统，其特征在于，包括扫描仪、坐标系设定模块和冰角量化评估模块，其中，
扫描仪用于扫描试验中绘制的冰形轮廓线网格图纸；
坐标系设定模块用于定义数字化坐标系；
冰角量化评估模块用于执行权利要求1-7任一所述的结冰风洞试验冰形冰角量化评估方法。
9.一种存储介质，其特征在于，所述存储介质存储了用于执行权利要求1-7任一所述的结冰风洞试验冰形冰角量化评估方法的计算机程序。
CN202010649538.0A 2020-07-08 2020-07-08 一种结冰风洞试验冰形冰角量化评估方法和系统 Active CN111539072B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010649538.0A CN111539072B (zh) 2020-07-08 2020-07-08 一种结冰风洞试验冰形冰角量化评估方法和系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010649538.0A CN111539072B (zh) 2020-07-08 2020-07-08 一种结冰风洞试验冰形冰角量化评估方法和系统

Publications (2)

Publication Number Publication Date
CN111539072A true CN111539072A (zh) 2020-08-14
CN111539072B CN111539072B (zh) 2020-09-29

Family

ID=71978366

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010649538.0A Active CN111539072B (zh) 2020-07-08 2020-07-08 一种结冰风洞试验冰形冰角量化评估方法和系统

Country Status (1)

Country Link
CN (1) CN111539072B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112393873A (zh) * 2021-01-11 2021-02-23 中国空气动力研究与发展中心低速空气动力研究所 一种攻角匹配方法
CN112985753A (zh) * 2021-05-06 2021-06-18 中国空气动力研究与发展中心计算空气动力研究所 一种用于风洞试验的冰形参数化方法
CN114169256A (zh) * 2022-02-14 2022-03-11 中国空气动力研究与发展中心低速空气动力研究所 一种获得结冰风洞缩比试验工况最优风速的方法
CN117191334A (zh) * 2023-11-08 2023-12-08 中国空气动力研究与发展中心低速空气动力研究所 一种复杂冰形测量误差评估方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20130117234A (ko) * 2012-04-18 2013-10-25 한국항공우주산업 주식회사 항공기 엔진의 공기 흡입구 주위에 대한 결빙 증식을 계산하는 방법
CN105783770A (zh) * 2016-01-22 2016-07-20 西南科技大学 一种基于线结构光的冰形轮廓测量的方法
CN105869212A (zh) * 2016-03-11 2016-08-17 空气动力学国家重点实验室 一种冰形简化方法
CN105957055A (zh) * 2016-04-20 2016-09-21 空气动力学国家重点实验室 一种冰形相似度量化评估方法
JP2018146363A (ja) * 2017-03-03 2018-09-20 三菱航空機株式会社 ３次元位置計測システム及び方法
CN109544576A (zh) * 2018-11-23 2019-03-29 中国空气动力研究与发展中心低速空气动力研究所 一种基于机器视觉的飞机模型冰形轮廓提取方法
CN109696289A (zh) * 2019-02-18 2019-04-30 中国空气动力研究与发展中心低速空气动力研究所 一种风力机翼段结冰测试系统及方法
CN110889233A (zh) * 2019-12-04 2020-03-17 中国空气动力研究与发展中心 结冰风洞试验冰形三维实体生成方法及三维实物冰
CN111238759A (zh) * 2020-03-31 2020-06-05 中国空气动力研究与发展中心低速空气动力研究所 一种结冰风洞测压试验方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20130117234A (ko) * 2012-04-18 2013-10-25 한국항공우주산업 주식회사 항공기 엔진의 공기 흡입구 주위에 대한 결빙 증식을 계산하는 방법
CN105783770A (zh) * 2016-01-22 2016-07-20 西南科技大学 一种基于线结构光的冰形轮廓测量的方法
CN105869212A (zh) * 2016-03-11 2016-08-17 空气动力学国家重点实验室 一种冰形简化方法
CN105957055A (zh) * 2016-04-20 2016-09-21 空气动力学国家重点实验室 一种冰形相似度量化评估方法
JP2018146363A (ja) * 2017-03-03 2018-09-20 三菱航空機株式会社 ３次元位置計測システム及び方法
CN109544576A (zh) * 2018-11-23 2019-03-29 中国空气动力研究与发展中心低速空气动力研究所 一种基于机器视觉的飞机模型冰形轮廓提取方法
CN109696289A (zh) * 2019-02-18 2019-04-30 中国空气动力研究与发展中心低速空气动力研究所 一种风力机翼段结冰测试系统及方法
CN110889233A (zh) * 2019-12-04 2020-03-17 中国空气动力研究与发展中心 结冰风洞试验冰形三维实体生成方法及三维实物冰
CN111238759A (zh) * 2020-03-31 2020-06-05 中国空气动力研究与发展中心低速空气动力研究所 一种结冰风洞测压试验方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
郭向东 等: "大型结冰风洞云雾场适航应用符合性验证", 《航空学报》 *
郭向东 等: "结冰风洞中液滴过冷特性数值研究", 《航空学报》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112393873A (zh) * 2021-01-11 2021-02-23 中国空气动力研究与发展中心低速空气动力研究所 一种攻角匹配方法
CN112985753A (zh) * 2021-05-06 2021-06-18 中国空气动力研究与发展中心计算空气动力研究所 一种用于风洞试验的冰形参数化方法
CN114169256A (zh) * 2022-02-14 2022-03-11 中国空气动力研究与发展中心低速空气动力研究所 一种获得结冰风洞缩比试验工况最优风速的方法
CN117191334A (zh) * 2023-11-08 2023-12-08 中国空气动力研究与发展中心低速空气动力研究所 一种复杂冰形测量误差评估方法

Also Published As

Publication number Publication date
CN111539072B (zh) 2020-09-29

Similar Documents

Publication Publication Date Title
CN111539072B (zh) 一种结冰风洞试验冰形冰角量化评估方法和系统
CN111323200B (zh) 一种结冰风洞试验冰形面积计算方法
CN106680798B (zh) 一种机载lidar航带重叠区冗余辨识及消除方法
CN109544576A (zh) 一种基于机器视觉的飞机模型冰形轮廓提取方法
Landwehr et al. Direct flux measurements from mobile platforms at sea: Motion and airflow distortion corrections revisited
Broeren et al. Validation of 3-D ice accretion measurement methodology for experimental aerodynamic simulation
Broeren et al. Aerodynamic simulation of ice accretion on airfoils
CN110641726B (zh) 一种快速确定飞机机翼安装角的方法
CN107357959A (zh) 一种高精度的叶片截面特征参数提取方法
CN108843492B (zh) 通过无人机进行风机偏航角测算方法及系统
CN111310381A (zh) 一种三维水滴收集系数计算方法
CN113420503A (zh) 基于温度序列相似性度量的积冰检测方法、系统、终端及应用
CN113792387B (zh) 飞行器积冰冰形模拟方法、装置、计算机设备及存储介质
CN111400826B (zh) 一种任意时刻冰形预测方法及系统
Zan et al. Analysis of patrol frigate air wakes
CN113642267B (zh) 一种飞行器表面前缘附着线区域提取方法
Broeren et al. Aerodynamic classification of swept-wing ice accretion
Hospers Eulerian method for super-cooled large-droplet ice-accretion on aircraft wings
Broeren et al. Effect of geometric fidelity on the aerodynamics of a swept wing with glaze ice accretion
CN111739048B (zh) 一种结冰风洞试验冰形几何轮廓线数字化方法
KR101409358B1 (ko) 에어포일 임계 결빙형상을 결정하는 방법
CN111964865A (zh) 用于飞机冰风洞试验的砂纸冰型粗糙度的计算方法
CN103390109A (zh) 一种气动力特性快速预测方法
Petterson Low-speed aerodynamic and flowfield characteristics of a UCAV
CN112985753B (zh) 一种用于风洞试验的冰形参数化方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant