CN108955981A - Suitable for rotation boundary layer wall surface shear stress measurement method and device - Google Patents

Suitable for rotation boundary layer wall surface shear stress measurement method and device Download PDF

Info

Publication number
CN108955981A
CN108955981A CN201810923270.8A CN201810923270A CN108955981A CN 108955981 A CN108955981 A CN 108955981A CN 201810923270 A CN201810923270 A CN 201810923270A CN 108955981 A CN108955981 A CN 108955981A
Authority
CN
China
Prior art keywords
boundary layer
shear stress
velocity
wall
wall surface
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
CN201810923270.8A
Other languages
Chinese (zh)
Other versions
CN108955981B (en
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.)
Beijing Feiqing Technology Co ltd
Original Assignee
Beijing Ruiao Technology Co Ltd
Beihang University
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 Beijing Ruiao Technology Co Ltd, Beihang University filed Critical Beijing Ruiao Technology Co Ltd
Priority to CN201810923270.8A priority Critical patent/CN108955981B/en
Publication of CN108955981A publication Critical patent/CN108955981A/en
Application granted granted Critical
Publication of CN108955981B publication Critical patent/CN108955981B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01LMEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L5/00Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

The invention discloses one kind to be suitable for rotation boundary layer wall surface shear stress measurement method, presets to probe;The speed at the rotation multiple and different pre-determined distances of boundary layer wall surface is measured by pre-set probe, the multiple measurement results that will acquire are transmitted;The velocity function that one section of linear distribution is selected in the nearest region of distance rotation boundary layer wall surface, to velocity function linear fit;Measurement is completed by translation, hypothesis estimation and nondimensionalization operation to the velocity gradient that linear fit operation obtains.This method do not need to carry out any processing to wall surface, does not influence near wall flowing, it can be achieved that cut off the measurement of stress to wall surface, simple and convenient;It can be realized simultaneously wall surface and cut off stress τwAnd the determination of wall surface distance y;The miniaturization test macro of use may be implemented the measurement of rotating condition lower wall surface shear stress, and have the high efficiency and ease for use of application.The invention also discloses be suitable for rotation boundary layer wall surface shear stress measuring device.

Description

适用于旋转边界层壁面剪切应力测量方法和装置Method and device for measuring wall shear stress of rotating boundary layer

技术领域technical field

本发明涉及测试测量技术领域,特别是涉及一种适用于旋转边界层壁面剪切应力测量方法和装置。The invention relates to the technical field of test and measurement, in particular to a method and device suitable for measuring shear stress of a rotating boundary layer wall.

背景技术Background technique

在旋转机械的研究中,壁面截切应力是一个非常重要的量。常规的静止条件下,由于安装空间大,可以采用多种方法实现剪切应力的测量。但是在旋转条件下,由于空间受限,因此壁面剪切应力的测量一直以来都是一个难题。In the study of rotating machinery, wall shear stress is a very important quantity. Under normal static conditions, due to the large installation space, various methods can be used to measure the shear stress. However, under rotating conditions, the measurement of wall shear stress has always been a difficult problem due to the limited space.

常规测量壁面剪切应力的方法一般采用薄膜传感器对其壁面的剪切应力进行测量,但是为了防止薄膜传感器对壁面流动的影响,一般在壁面开一个槽,将薄膜传感器埋入槽中,并保证薄膜传感器的外表面与壁面一平,进而消除对壁面附近流动的影响。但是这种方法需要对通道壁面进行特殊的开槽处理,较为复杂;尤其在壁面加热时,由于还需要布置加热装置,该方法则更难实现。The conventional method of measuring wall shear stress generally uses a thin film sensor to measure the shear stress on the wall, but in order to prevent the influence of the thin film sensor on the flow of the wall, a groove is generally opened on the wall, and the thin film sensor is buried in the groove, and ensure The outer surface of the thin film sensor is flush with the wall surface, thereby eliminating the influence on the flow near the wall surface. However, this method requires special groove treatment on the wall surface of the channel, which is relatively complicated; especially when the wall surface is heated, because a heating device needs to be arranged, this method is more difficult to implement.

发明内容Contents of the invention

基于此,有必要针对传统技术存在的问题,提供一种适用于旋转边界层壁面剪切应力测量方法和装置。具体的,在本公开中,为了解决旋转机械中的壁面截切应力测量的问题,本公开提出了适用于旋转边界层壁面剪切应力测量方法,具体为,基于线性底层假设的旋转边界层壁面剪切应力测量方法。该方法利用边界层热线探针、位移机构、小型CTA模块、小型数模转换模块、滑环引电器以及计算机,对旋转条件下边界层速度进行测量,最终通过线性底层假设,计算出壁面剪切应力。通过该方法的提出,即不需要对壁面进行任何处理、不影响壁面附近流动,即可实现对壁面截切应力的测量,简单方便;同时可以实现壁面截切应力τw以及壁面距离y的确定;进一步地,采用的小型化测试系统,可以实现旋转条件下壁面剪切应力的测量。Based on this, it is necessary to provide a method and device suitable for measuring the wall shear stress of the rotating boundary layer for the problems existing in the traditional technology. Specifically, in this disclosure, in order to solve the problem of wall shear stress measurement in rotating machinery, this disclosure proposes a method for measuring wall shear stress applicable to rotating boundary layers, specifically, the rotating boundary layer wall surface based on the assumption of a linear bottom layer Shear stress measurement method. This method uses a boundary layer hot wire probe, a displacement mechanism, a small CTA module, a small digital-to-analog conversion module, a slip ring lead, and a computer to measure the boundary layer velocity under rotating conditions, and finally calculates the wall shear through the assumption of a linear bottom layer. stress. Through the proposal of this method, the measurement of the wall shear stress can be realized without any treatment on the wall and without affecting the flow near the wall, which is simple and convenient; at the same time, the determination of the wall shear stress τ w and the wall distance y can be realized ; Furthermore, the miniaturized test system adopted can realize the measurement of the wall shear stress under the rotating condition.

第一方面,本发明实施例提供了一种适用于旋转边界层壁面剪切应力测量方法,所述方法包括:对探针进行预先设置;通过预先设置的所述探针对旋转边界层壁面多个不同预设距离处的速度进行测量,并将获取的多个测量结果进行传输;在距离旋转边界层壁面最近的区域内选择一段线性分布的速度函数,对所述速度函数进行线性拟合操作;对线性拟合操作获取的速度梯度依次通过平移、假设估计以及无量纲化操作完成测量。In the first aspect, the embodiment of the present invention provides a method suitable for measuring the shear stress of the wall surface of the rotating boundary layer. The method includes: pre-setting the probe; measure the velocity at different preset distances, and transmit the obtained multiple measurement results; select a linearly distributed velocity function in the area closest to the wall of the rotating boundary layer, and perform a linear fitting operation on the velocity function. ; The velocity gradient obtained by the linear fitting operation is sequentially measured through translation, hypothesis estimation and non-dimensionalization operations.

在其中一个实施例中,所述对探针进行预先设置包括:将边界层热线探针通过支杆固定设置于位移机构上方,并将所述边界层热线探针固定设置于旋转边界层壁面预设距离处。In one of the embodiments, the preset setting of the probe includes: fixing the boundary layer hot-wire probe above the displacement mechanism through a support rod, and fixing the boundary layer hot-wire probe on the rotating boundary layer wall preset Set the distance.

在其中一个实施例中,所述通过预先设置的所述探针对旋转边界层壁面多个不同预设距离处的速度进行测量包括:通过预先设置的边界层热线探针对距离旋转边界层壁面最近的一点的速度进行测量,以及通过预先设置的边界层热线探针对距离旋转边界层壁面最远的一点的速度进行测量。In one of the embodiments, the measuring the velocity at multiple preset distances from the wall of the rotating boundary layer through the preset probe includes: measuring the distance from the wall of the rotating boundary layer through the preset hot wire probe Velocity is measured at the nearest point and at the point furthest from the rotating boundary layer wall by a pre-set boundary layer hot wire probe.

在其中一个实施例中,还包括:通过预先设置的所述边界层热线探针对距离旋转边界层壁面所述最近的一点与所述最远的一点中间的多个点的速度进行测量。In one of the embodiments, the method further includes: measuring the velocity of multiple points between the closest point and the farthest point from the wall surface of the rotating boundary layer by using the pre-set boundary layer hot wire probe.

在其中一个实施例中,所述将获取的多个测量结果进行传输包括:依次通过CTA模块、数模转化模块以及滑环引电器将获取的所述多个测量结果传输至计算机上。In one embodiment, the transmitting the multiple acquired measurement results includes: sequentially transmitting the multiple acquired measurement results to the computer through the CTA module, the digital-to-analog conversion module, and the slip ring lead.

在其中一个实施例中,还包括:检验用于进行所述线性拟合操作所选取的所述速度函数上的坐标点是否与所有处在预设范围内的坐标点全部重合;若用于进行所述线性拟合操作所选取的所述速度函数上的坐标点与所有处在预设范围内的坐标点全部重合,则对线性拟合操作获取的所述速度梯度依次通过平移、假设估计获取的旋转边界层壁面剪切应力为真实的旋转边界层壁面剪切应力。In one of the embodiments, it also includes: checking whether the coordinate points on the velocity function selected for the linear fitting operation coincide with all the coordinate points within the preset range; The coordinate points on the velocity function selected by the linear fitting operation coincide with all the coordinate points within the preset range, and then the velocity gradient obtained by the linear fitting operation is sequentially acquired through translation and hypothetical estimation The wall shear stress of the rotating boundary layer is the real wall shear stress of the rotating boundary layer.

在其中一个实施例中,还包括:若用于进行所述线性拟合操作所选取的所述速度函数上的坐标点与所有处在预设范围内的坐标点未全部重合,选择所述预设范围内的坐标点的所述速度函数依次通过平移、假设估计以及无量纲化循环操作直至获取所述真实的所述旋转边界层壁面剪切应力以及与其对应的纵坐标。In one of the embodiments, it also includes: if the coordinate points on the speed function selected for the linear fitting operation do not coincide with all the coordinate points within the preset range, selecting the preset The velocity function of the coordinate points within the set range is sequentially operated through translation, hypothesis estimation and non-dimensionalization until the real wall shear stress of the rotating boundary layer and its corresponding ordinate are obtained.

第二方面,本发明实施例提供了一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现上述第一方面的适用于旋转边界层壁面剪切应力测量方法。In the second aspect, an embodiment of the present invention provides a computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the above-mentioned first aspect applicable to the rotation boundary can be realized. Method for measuring shear stress on layer wall.

第三方面,本发明实施例提供了一种包含指令的计算机程序产品,当该计算机程序产品在计算机上运行时,使得计算机执行上述第一方面所述的方法。In a third aspect, an embodiment of the present invention provides a computer program product including instructions, and when the computer program product is run on a computer, it causes the computer to execute the method described in the first aspect above.

第四方面,本发明实施例还提供了一种适用于旋转边界层壁面剪切应力测量装置,所述装置包括:设置模块,用于对探针进行预先设置;测量与传输模块,用于通过预先设置的所述探针对旋转边界层壁面多个不同预设距离处的速度进行测量,并将获取的多个测量结果进行传输;拟合处理模块,用于在距离旋转边界层壁面最近的区域内选择一段线性分布的速度函数,对所述速度函数进行线性拟合操作;测量模块,用于对线性拟合操作获取的速度梯度依次通过平移、假设估计以及无量纲化操作完成测量。In the fourth aspect, the embodiment of the present invention also provides a device suitable for measuring the shear stress of the wall surface of the rotating boundary layer, the device includes: a setting module, used to pre-set the probe; a measurement and transmission module, used to pass The preset probe measures the velocity at multiple preset distances from the wall of the rotating boundary layer, and transmits the obtained multiple measurement results; the fitting processing module is used for Select a section of linearly distributed velocity function in the region, and perform a linear fitting operation on the velocity function; the measurement module is used to measure the velocity gradient obtained by the linear fitting operation through translation, hypothesis estimation, and dimensionless operation in sequence.

本发明提供的一种适用于旋转边界层壁面剪切应力测量方法和装置,对探针进行预先设置;通过预先设置的探针对旋转边界层壁面多个不同预设距离处的速度进行测量,并将获取的多个测量结果进行传输;在距离旋转边界层壁面最近的区域内选择一段线性分布的速度函数,对速度函数进行线性拟合操作;对线性拟合操作获取的速度梯度依次通过平移、假设估计以及无量纲化操作完成测量。该方法利用边界层热线探针、位移机构、小型CTA模块、小型数模转换模块、滑环引电器以及计算机,对旋转条件下边界层速度进行测量,最终通过线性底层假设,计算出壁面剪切应力。通过该方法的提出,即不需要对壁面进行任何处理、不影响壁面附近流动,即可实现对壁面截切应力的测量,简单方便;同时可以实现壁面截切应力τw以及壁面距离y的确定;进一步地,采用的小型化测试系统,可以实现旋转条件下壁面剪切应力的测量,且具有应用的高效性与易用性。The present invention provides a method and device suitable for measuring the shear stress of the wall surface of the rotating boundary layer. The probe is preset; the velocity at multiple preset distances on the wall surface of the rotating boundary layer is measured by the preset probe, And transmit the multiple measurement results obtained; select a linearly distributed velocity function in the area closest to the wall of the rotating boundary layer, and perform a linear fitting operation on the velocity function; the velocity gradient obtained by the linear fitting operation is sequentially passed through translation , hypothetical estimation, and dimensionless operations to complete the measurement. This method uses a boundary layer hot wire probe, a displacement mechanism, a small CTA module, a small digital-to-analog conversion module, a slip ring lead, and a computer to measure the boundary layer velocity under rotating conditions, and finally calculates the wall shear through the assumption of a linear bottom layer. stress. Through the proposal of this method, the measurement of the wall shear stress can be realized without any treatment on the wall and without affecting the flow near the wall, which is simple and convenient; at the same time, the determination of the wall shear stress τ w and the wall distance y can be realized ; Further, the adopted miniaturized test system can realize the measurement of the wall shear stress under the rotating condition, and has high efficiency and ease of use.

附图说明Description of drawings

图1为本发明一个实施例中的一种适用于旋转边界层壁面剪切应力测量方法的步骤流程示意图;Fig. 1 is a schematic flow chart of a method for measuring wall shear stress of a rotating boundary layer in one embodiment of the present invention;

图2为本发明一个实施例中的一种适用于旋转边界层壁面剪切应力测量装置的结构示意图。Fig. 2 is a schematic structural diagram of a device suitable for measuring wall shear stress of a rotating boundary layer in an embodiment of the present invention.

具体实施方式Detailed ways

为了使本发明的目的、技术方案及优点更加清楚明白,以下通过实施例,并结合附图,对本发明适用于旋转边界层壁面剪切应力测量方法和装置的具体实施方式进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。In order to make the purpose, technical solution and advantages of the present invention clearer, the specific implementation of the present invention, which is applicable to the method and device for measuring the wall shear stress of the rotating boundary layer, will be further described in detail through the following examples and in conjunction with the accompanying drawings. It should be understood that the specific embodiments described here are only used to explain the present invention, not to limit the present invention.

如图1所示,为一个实施例中的一种适用于旋转边界层壁面剪切应力测量方法的流程示意图。具体包括以下步骤:As shown in FIG. 1 , it is a schematic flowchart of a method for measuring wall shear stress of a rotating boundary layer in an embodiment. Specifically include the following steps:

步骤102,对探针进行预先设置。需要说明的是,对探针进行预先设置包括:将边界层热线探针通过支杆固定设置于位移机构上方,并将边界层热线探针固定设置于旋转边界层壁面预设距离处。Step 102, presetting the probe. It should be noted that the presetting of the probe includes: fixing the boundary layer hot-wire probe above the displacement mechanism through the support rod, and fixing the boundary layer hot-wire probe at a preset distance from the rotating boundary layer wall.

步骤104,通过预先设置的探针对旋转边界层壁面多个不同预设距离处的速度进行测量,并将获取的多个测量结果进行传输。Step 104: Measure the velocity at multiple preset distances from the wall surface of the rotating boundary layer by using a preset probe, and transmit the obtained multiple measurement results.

在一个实施例中,通过预先设置的探针对旋转边界层壁面多个不同预设距离处的速度进行测量包括:通过预先设置的边界层热线探针对距离旋转边界层壁面最近的一点的速度进行测量,以及通过预先设置的边界层热线探针对距离旋转边界层壁面最远的一点的速度进行测量。In one embodiment, using a preset probe to measure the velocity at multiple preset distances from the wall of the rotating boundary layer includes: measuring the velocity of a point closest to the wall of the rotating boundary layer using a preset hot wire probe of the boundary layer Measurements are taken, as well as the velocity at the point farthest from the rotating boundary layer wall by a pre-set boundary layer hot wire probe.

进一步地,本公开涉及的一种适用于旋转边界层壁面剪切应力测量方法还包括:通过预先设置的边界层热线探针对距离旋转边界层壁面最近的一点与最远的一点中间的多个点的速度进行测量。Further, the method for measuring the shear stress of the wall surface of the rotating boundary layer according to the present disclosure further includes: using a pre-set boundary layer hot wire probe to pair a plurality of The velocity of the point is measured.

此外,还需要说明的是,将获取的多个测量结果进行传输包括:依次通过CTA(Constant Temperature Anemometer,恒温式热线风速仪)模块、数模转化模块以及滑环引电器将获取的多个测量结果传输至计算机上。由此,提高了计算机实时获取测量结果的高效性与易用性。In addition, it should also be noted that the transmission of multiple measurement results obtained includes: multiple measurements that will be obtained through the CTA (Constant Temperature Anemometer, constant temperature hot wire anemometer) module, digital-to-analog conversion module, and slip ring lead-in device in turn. The results are transferred to a computer. As a result, the efficiency and ease of use of the computer for real-time acquisition of measurement results are improved.

步骤106,在距离旋转边界层壁面最近的区域内选择一段线性分布的速度函数,对速度函数进行线性拟合操作。Step 106, select a section of linearly distributed velocity function in the area closest to the wall of the rotating boundary layer, and perform a linear fitting operation on the velocity function.

步骤108,对线性拟合操作获取的速度梯度依次通过平移、假设估计以及无量纲化操作完成测量。In step 108, the velocity gradient acquired by the linear fitting operation is sequentially measured through translation, hypothesis estimation, and dimensionless operations.

在一个实施例中,本公开涉及的一种适用于旋转边界层壁面剪切应力测量方法还包括:检验用于进行线性拟合操作所选取的速度函数上的坐标点是否与所有处在预设范围内的坐标点全部重合;若用于进行线性拟合操作所选取的速度函数上的坐标点与所有处在预设范围内的坐标点全部重合,则对线性拟合操作获取的速度梯度依次通过平移、假设估计获取的旋转边界层壁面剪切应力为真实的旋转边界层壁面剪切应力。In one embodiment, the method for measuring the wall shear stress of the rotating boundary layer according to the present disclosure further includes: checking whether the coordinate points on the velocity function selected for the linear fitting operation are consistent with all the preset The coordinate points within the range are all coincident; if the coordinate points on the velocity function selected for the linear fitting operation coincide with all the coordinate points within the preset range, then the velocity gradient obtained by the linear fitting operation is sequentially The wall shear stress of the rotating boundary layer obtained by translation and assumption estimation is the real wall shear stress of the rotating boundary layer.

此外,需要说明的是,本公开涉及的一种适用于旋转边界层壁面剪切应力测量方法还包括:若用于进行线性拟合操作所选取的速度函数上的坐标点与所有处在预设范围内的坐标点未全部重合,选择预设范围内的坐标点的速度函数依次通过平移、假设估计以及无量纲化循环操作直至获取真实的旋转边界层壁面剪切应力以及与其对应的纵坐标。In addition, it should be noted that a method for measuring the wall shear stress of the rotating boundary layer involved in the present disclosure also includes: if the coordinate points on the velocity function selected for the linear fitting operation are all in the preset The coordinate points in the range are not all coincident, and the velocity function of the coordinate points in the preset range is selected through translation, assumption estimation and non-dimensionalization cycle operation until the real wall shear stress of the rotating boundary layer and its corresponding ordinate are obtained.

为了进一步地理解并运用本公开提出的一种适用于旋转边界层壁面剪切应力测量方法,进行以下示例。需要说明的是,本公开所保护的范围不限于以下示例。In order to further understand and apply a method for measuring wall shear stress in a rotating boundary layer proposed by the present disclosure, the following example is performed. It should be noted that the protection scope of the present disclosure is not limited to the following examples.

具体的,本公开涉及旋转涡轮机械边界层速度测量领域,尤其涉及一种基于线性底层假设的旋转边界层壁面剪切应力测量方法。基于线性底层假设的旋转边界层壁面剪切应力测量方法,包括以下步骤:Specifically, the present disclosure relates to the field of velocity measurement of the boundary layer of a rotating turbomachinery, and in particular to a method for measuring wall shear stress of a rotating boundary layer based on a linear bottom layer assumption. The method for measuring the wall shear stress of the rotating boundary layer based on the assumption of a linear bottom layer includes the following steps:

首先,将边界层探针通过支杆固定在位移机构上,并将探针置于靠近壁面位置处,假设距离壁面距离为y0;其次,探针测量的靠近壁面处的速度,通过小型CTA(ConstantTemperature Anemometer,恒温式热线风速仪)模块、小型数模转化模块以及滑环引电器传递到静止的计算机上,实现对壁面附近第一个点的速度(v0)的测量,假设第一点探针靠近壁面的距离为y0First, the boundary layer probe is fixed on the displacement mechanism through the support rod, and the probe is placed close to the wall, assuming that the distance from the wall is y 0 ; secondly, the velocity near the wall measured by the probe is measured by a small CTA (ConstantTemperature Anemometer, constant temperature hot-wire anemometer) module, small digital-to-analog conversion module and slip ring leads are transmitted to the static computer to realize the measurement of the velocity (v 0 ) of the first point near the wall, assuming the first point The distance between the probe and the wall is y 0 .

进一步地,在位移机构的控制下,将热线探针向远离壁面的位置移动△y的距离(y1),再次测量速度v1;更进一步地,重复上述的过程,直至测量到足够原理壁面的距离的速度yn,vnFurther, under the control of the displacement mechanism, move the hot wire probe away from the wall by a distance of △y (y 1 ), and measure the velocity v 1 again; furthermore, repeat the above process until a sufficient wall surface is measured. The speed y n ,v n of the distance.

需要说明的是,通过以上获得的速度(v0,v1···vn)和距离(y0,y1···yn),在靠近壁面的区域内选择一段接近线性分布的速度(vm~vt),对其进行线性拟合获得速度梯度;平移y坐标系(δy),使得拟合直线能够通过原点(y’0,y’1···y’n,其中y’0=δy+y0)。It should be noted that, through the velocity (v 0 , v 1 ···v n ) and distance (y 0 , y 1 ···y n ) obtained above, select a section of velocity close to the linear distribution in the area close to the wall (v m ~v t ), linearly fit it to obtain the velocity gradient; translate the y coordinate system (δy), so that the fitted line can pass through the origin (y' 0 ,y' 1 ···y' n , where y ' 0 =δy+y 0 ).

可以理解的是,假设这个速度梯度就是线性底层速度梯度,并且根据这个速度梯度初步估计壁面剪切应力(vm<v<vt,y’m<y’<y’t)。需要说明的是,τw是壁面剪切应力,其中右上角标记符号代表的是第一次迭代的结果,且其中μ是流体的粘性系数。根据初步估计的壁面剪切应力对平移后的y坐标(y’0,y’1···y’n)以及速度v(v0,v1···vn)进行无量纲化,即获得无量纲的需要说明的是,在获得无量纲的第一个结果中,其中y+代表的是无量纲的测点距离壁面距离,ρ代表密度,此处分母的v代表流体动力粘度;在获得无量纲的第二个结果中,U+代表无量纲的速度,其中,公式中的分子v代表速度。It is understandable that the velocity gradient is assumed to be a linear bottom velocity gradient, and the wall shear stress is preliminarily estimated based on this velocity gradient ( v m <v<v t ,y' m <y'<y' t ). It should be noted that τ w is the wall shear stress, where the symbol in the upper right corner represents the result of the first iteration, and where μ is the viscosity coefficient of the fluid. According to the preliminary estimated wall shear stress, the translated y-coordinate (y' 0 ,y' 1 ···y' n ) and velocity v(v 0 ,v 1 ···v n ) are dimensionless, namely get dimensionless It should be noted that, in the first dimensionless result, y+ represents the dimensionless distance from the measuring point to the wall, ρ represents the density, and v in the denominator here represents the hydrodynamic viscosity; In the two results, U+ represents the dimensionless velocity, and the numerator v in the formula represents the velocity.

最终,检查用于拟合直线所选取的速度上的坐标点(vm<v<vt,y’m<y’<y’t)是否与所有处在3.5<y+<5的范围内的坐标点全部重合。如果是,则之前获得的壁面剪切应力就是真实的壁面剪切应力,平移后的y坐标就是真实的y坐标。如果不是,则选择所有处在3.5<y+<5的范围内的坐标点的速度重复平移、假设估计以及计算获取无量纲的步骤,直至得到真实的壁面剪切应力和y坐标为止。Finally, check whether the coordinate points on the speed selected for fitting the straight line (v m <v<v t ,y' m <y'<y' t ) are within the range of 3.5<y + <5 The coordinate points of all coincide. If yes, the wall shear stress obtained before is the real wall shear stress, and the translated y-coordinate is the real y-coordinate. If not, select the velocities of all coordinate points within the range of 3.5<y + <5 and repeat the steps of translation, hypothesis estimation and calculation to obtain dimensionless until the real wall shear stress and y coordinates are obtained.

本发明提供的一种适用于旋转边界层壁面剪切应力测量方法,对探针进行预先设置;通过预先设置的探针对旋转边界层壁面多个不同预设距离处的速度进行测量,并将获取的多个测量结果进行传输;在距离旋转边界层壁面最近的区域内选择一段线性分布的速度函数,对速度函数进行线性拟合操作;对线性拟合操作获取的速度梯度依次通过平移、假设估计以及无量纲化操作完成测量。该方法利用边界层热线探针、位移机构、小型CTA模块、小型数模转换模块、滑环引电器以及计算机,对旋转条件下边界层速度进行测量,最终通过线性底层假设,计算出壁面剪切应力。该方法即不需要对壁面进行任何处理、不影响壁面附近流动,即可实现对壁面截切应力的测量,简单方便;同时可以实现壁面截切应力τw以及壁面距离y的确定;进一步地,采用的小型化测试系统,可以实现旋转条件下壁面剪切应力的测量,且具有应用的高效性与易用性。The present invention provides a method suitable for measuring the shear stress of the wall surface of the rotating boundary layer. The probe is preset; the velocity at multiple preset distances on the wall surface of the rotating boundary layer is measured through the preset probe, and the The multiple measurement results obtained are transmitted; a velocity function with linear distribution is selected in the area closest to the wall of the rotating boundary layer, and a linear fitting operation is performed on the velocity function; the velocity gradient obtained by the linear fitting operation is sequentially passed through translation, assumption Estimation as well as dimensionless operations are performed on measurements. This method uses a boundary layer hot wire probe, a displacement mechanism, a small CTA module, a small digital-to-analog conversion module, a slip ring lead, and a computer to measure the boundary layer velocity under rotating conditions, and finally calculates the wall shear through the assumption of a linear bottom layer. stress. This method does not require any treatment on the wall surface and does not affect the flow near the wall surface, and can realize the measurement of the wall shear stress, which is simple and convenient; at the same time, it can realize the determination of the wall shear stress τ w and the wall distance y; further, The miniaturized test system adopted can realize the measurement of wall shear stress under rotating conditions, and has high efficiency and ease of use.

基于同一发明构思,还提供了一种适用于旋转边界层壁面剪切应力测量装置。由于此装置解决问题的原理与前述一种适用于旋转边界层壁面剪切应力测量方法相似,因此,该装置的实施可以按照前述方法的具体步骤实现,重复之处不再赘述。Based on the same inventive concept, a device suitable for measuring wall shear stress of a rotating boundary layer is also provided. Since the problem-solving principle of this device is similar to the aforementioned method for measuring the wall shear stress of the rotating boundary layer, the implementation of the device can be realized according to the specific steps of the aforementioned method, and the repetition will not be repeated.

如图2所示,为一个实施例中的一种适用于旋转边界层壁面剪切应力测量装置的结构示意图。该适用于旋转边界层壁面剪切应力测量装置10包括:设置模块100、测量与传输模块200、拟合处理模块300和测量模块400。As shown in FIG. 2 , it is a schematic structural diagram of a device suitable for measuring wall shear stress of a rotating boundary layer in an embodiment. The device 10 suitable for measuring wall shear stress of a rotating boundary layer includes: a setting module 100 , a measurement and transmission module 200 , a fitting processing module 300 and a measurement module 400 .

其中,设置模块100用于对探针进行预先设置;测量与传输模块200用于通过预先设置的探针对旋转边界层壁面多个不同预设距离处的速度进行测量,并将获取的多个测量结果进行传输;拟合处理模块300用于在距离旋转边界层壁面最近的区域内选择一段线性分布的速度函数,对速度函数进行线性拟合操作;测量模块400用于对线性拟合操作获取的速度梯度依次通过平移、假设估计以及无量纲化操作完成测量。Among them, the setting module 100 is used to pre-set the probe; the measurement and transmission module 200 is used to measure the velocity at multiple preset distances on the wall surface of the rotating boundary layer through the preset probe, and the obtained multiple The measurement results are transmitted; the fitting processing module 300 is used to select a section of linearly distributed velocity function in the area closest to the wall of the rotating boundary layer, and performs a linear fitting operation on the velocity function; the measurement module 400 is used to obtain the linear fitting operation The velocity gradient of is measured sequentially through translation, assumption estimation, and dimensionless operations.

本发明提供的一种适用于旋转边界层壁面剪切应力测量装置,首先通过设置模块对探针进行预先设置;再通过测量与传输模块通过预先设置的探针对旋转边界层壁面多个不同预设距离处的速度进行测量,并将获取的多个测量结果进行传输;再次通过拟合处理模块在距离旋转边界层壁面最近的区域内选择一段线性分布的速度函数,对速度函数进行线性拟合操作;最终通过测量模块对线性拟合操作获取的速度梯度依次通过平移、假设估计以及无量纲化操作完成测量。该装置利用边界层热线探针、位移机构、小型CTA模块、小型数模转换模块、滑环引电器以及计算机,对旋转条件下边界层速度进行测量,最终通过线性底层假设,计算出壁面剪切应力。该方法即不需要对壁面进行任何处理、不影响壁面附近流动,即可实现对壁面截切应力的测量,简单方便;同时可以实现壁面截切应力τw以及壁面距离y的确定;进一步地,采用的小型化测试系统,可以实现旋转条件下壁面剪切应力的测量,且具有应用的高效性与易用性。The present invention provides a device suitable for measuring the shear stress of the wall surface of the rotating boundary layer. Firstly, the probe is preset through the setting module; Set the velocity at the distance to measure, and transmit the obtained multiple measurement results; again, through the fitting processing module, select a linearly distributed velocity function in the area closest to the wall of the rotating boundary layer, and perform linear fitting on the velocity function Operation; finally, the velocity gradient obtained by the linear fitting operation is finally measured by the measurement module through translation, hypothesis estimation and non-dimensionalization operations. The device uses a boundary layer hot wire probe, a displacement mechanism, a small CTA module, a small digital-to-analog conversion module, a slip ring lead, and a computer to measure the boundary layer velocity under rotating conditions, and finally calculates the wall shear through the assumption of a linear bottom layer. stress. This method does not require any treatment on the wall surface and does not affect the flow near the wall surface, and can realize the measurement of the wall shear stress, which is simple and convenient; at the same time, it can realize the determination of the wall shear stress τ w and the wall distance y; further, The miniaturized test system adopted can realize the measurement of wall shear stress under rotating conditions, and has high efficiency and ease of use.

本发明实施例还提供了一种计算机可读存储介质。该计算机可读存储介质上存储有计算机程序,该程序被图1中处理器执行。The embodiment of the present invention also provides a computer-readable storage medium. A computer program is stored on the computer-readable storage medium, and the program is executed by the processor in FIG. 1 .

本发明实施例还提供了一种包含指令的计算机程序产品。当该计算机程序产品在计算机上运行时,使得计算机执行上述图1的方法。The embodiment of the present invention also provides a computer program product including instructions. When the computer program product is run on the computer, the computer is made to execute the above-mentioned method in FIG. 1 .

本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于一计算机可读取存储介质中,该程序在执行时,可包括如上述各方法的实施例的流程。其中,所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory,ROM)或随机存储记忆体(Random AccessMemory,RAM)等。Those of ordinary skill in the art can understand that all or part of the processes in the methods of the above embodiments can be implemented through computer programs to instruct related hardware, and the programs can be stored in a computer-readable storage medium. During execution, it may include the processes of the embodiments of the above-mentioned methods. Wherein, the storage medium may be a magnetic disk, an optical disk, a read-only memory (Read-Only Memory, ROM) or a random access memory (Random Access Memory, RAM) and the like.

以上所述实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。The technical features of the above-mentioned embodiments can be combined arbitrarily. To make the description concise, all possible combinations of the technical features in the above-mentioned embodiments are not described. However, as long as there is no contradiction in the combination of these technical features, should be considered as within the scope of this specification.

以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。The above-mentioned embodiments only express several implementation modes of the present invention, and the description thereof is relatively specific and detailed, but should not be construed as limiting the patent scope of the present invention. It should be pointed out that those skilled in the art can make several modifications and improvements without departing from the concept of the present invention, and these all belong to the protection scope of the present invention. Therefore, the protection scope of the patent for the present invention should be based on the appended claims.

Claims (10)

1.一种适用于旋转边界层壁面剪切应力测量方法,其特征在于,所述方法包括:1. A method for measuring wall shear stress applicable to rotating boundary layers, is characterized in that, the method comprises: 对探针进行预先设置;Preset the probe; 通过预先设置的所述探针对旋转边界层壁面多个不同预设距离处的速度进行测量,并将获取的多个测量结果进行传输;Using the preset probe to measure the velocity at multiple preset distances on the wall surface of the rotating boundary layer, and transmit the obtained multiple measurement results; 在距离旋转边界层壁面最近的区域内选择一段线性分布的速度函数,对所述速度函数进行线性拟合操作;Select a section of velocity function with linear distribution in the area nearest to the wall of the rotating boundary layer, and perform a linear fitting operation on the velocity function; 对线性拟合操作获取的速度梯度依次通过平移、假设估计以及无量纲化操作完成测量。The velocity gradient obtained by the linear fitting operation is sequentially measured through translation, hypothesis estimation, and dimensionless operation. 2.根据权利要求1所述的适用于旋转边界层壁面剪切应力测量方法,其特征在于,所述对探针进行预先设置包括:将边界层热线探针通过支杆固定设置于位移机构上方,并将所述边界层热线探针固定设置于旋转边界层壁面预设距离处。2. The method for measuring wall shear stress suitable for rotating boundary layers according to claim 1, wherein said presetting the probe comprises: fixing the boundary layer hot wire probe above the displacement mechanism through a support rod , and the boundary layer hot wire probe is fixedly arranged at a preset distance from the wall surface of the rotating boundary layer. 3.根据权利要求1所述的适用于旋转边界层壁面剪切应力测量方法,其特征在于,所述通过预先设置的所述探针对旋转边界层壁面多个不同预设距离处的速度进行测量包括:通过预先设置的边界层热线探针对距离旋转边界层壁面最近的一点的速度进行测量,以及通过预先设置的边界层热线探针对距离旋转边界层壁面最远的一点的速度进行测量。3. according to claim 1, be applicable to rotating boundary layer wall shear stress measuring method, it is characterized in that, described by the described probe that presets to the velocity at a plurality of different preset distances of rotating boundary layer wall Measurements include: measurement of the velocity at the point closest to the wall of the rotating boundary layer by a pre-set boundary layer hot-wire probe, and measurement of the velocity of a point farthest from the wall of the rotating boundary layer by a pre-set hot-wire probe of the boundary layer . 4.根据权利要求3所述的适用于旋转边界层壁面剪切应力测量方法,其特征在于,还包括:通过预先设置的所述边界层热线探针对距离旋转边界层壁面所述最近的一点与所述最远的一点中间的多个点的速度进行测量。4. according to claim 3, be applicable to rotating boundary layer wall shear stress measurement method, it is characterized in that, also comprise: through described boundary layer hot-wire probe set in advance to distance the described nearest point of rotating boundary layer wall The velocities of multiple points intermediate to the furthest point are measured. 5.根据权利要求1所述的适用于旋转边界层壁面剪切应力测量方法,其特征在于,所述将获取的多个测量结果进行传输包括:5. The method for measuring wall shear stress applicable to rotating boundary layers according to claim 1, wherein said transmitting a plurality of measurement results obtained comprises: 依次通过CTA模块、数模转化模块以及滑环引电器将获取的所述多个测量结果传输至计算机上。The multiple measurement results obtained are transmitted to the computer through the CTA module, the digital-to-analog conversion module and the slip ring lead in sequence. 6.根据权利要求1所述的适用于旋转边界层壁面剪切应力测量方法,其特征在于,还包括:检验用于进行所述线性拟合操作所选取的所述速度函数上的坐标点是否与所有处在预设范围内的坐标点全部重合;6. The method for measuring wall shear stress applicable to rotating boundary layers according to claim 1, further comprising: checking whether the coordinate point on the selected velocity function for carrying out the linear fitting operation is It coincides with all the coordinate points within the preset range; 若用于进行所述线性拟合操作所选取的所述速度函数上的坐标点与所有处在预设范围内的坐标点全部重合,则对线性拟合操作获取的所述速度梯度依次通过平移、假设估计获取的旋转边界层壁面剪切应力为真实的旋转边界层壁面剪切应力。If the coordinate points on the velocity function selected for the linear fitting operation coincide with all the coordinate points within the preset range, the velocity gradient obtained by the linear fitting operation is sequentially translated , assuming that the estimated wall shear stress of the rotating boundary layer is the real wall shear stress of the rotating boundary layer. 7.根据权利要求6所述的适用于旋转边界层壁面剪切应力测量方法,其特征在于,还包括:若用于进行所述线性拟合操作所选取的所述速度函数上的坐标点与所有处在预设范围内的坐标点未全部重合,选择所述预设范围内的坐标点的所述速度函数依次通过平移、假设估计以及无量纲化循环操作直至获取所述真实的所述旋转边界层壁面剪切应力以及与其对应的纵坐标。7. the method for measuring wall shear stress applicable to rotating boundary layers according to claim 6, further comprising: if the coordinate point on the selected velocity function is used for carrying out the linear fitting operation and All the coordinate points within the preset range are not all coincident, and the velocity function for selecting the coordinate points within the preset range is sequentially operated through translation, hypothesis estimation, and dimensionless loop operations until the real rotation is obtained. Boundary layer wall shear stress and its corresponding ordinate. 8.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该程序被处理器执行时实现所述权利要求1-7中任一项所述方法的步骤。8. A computer-readable storage medium, on which a computer program is stored, wherein when the program is executed by a processor, the steps of the method according to any one of claims 1-7 are implemented. 9.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现所述权利要求1-7中任一项所述方法的步骤。9. A computer device, comprising a memory, a processor, and a computer program stored on the memory and operable on the processor, characterized in that, when the processor executes the program, it implements claims 1-7 The steps of any one of the methods. 10.一种适用于旋转边界层壁面剪切应力测量装置,其特征在于,所述装置包括:10. A device for measuring wall shear stress of a rotating boundary layer, characterized in that the device comprises: 设置模块,用于对探针进行预先设置;A setting module is used to pre-set the probe; 测量与传输模块,用于通过预先设置的所述探针对旋转边界层壁面多个不同预设距离处的速度进行测量,并将获取的多个测量结果进行传输;The measurement and transmission module is used to measure the velocity at multiple preset distances on the wall surface of the rotating boundary layer through the preset probe, and transmit the obtained multiple measurement results; 拟合处理模块,用于在距离旋转边界层壁面最近的区域内选择一段线性分布的速度函数,对所述速度函数进行线性拟合操作;The fitting processing module is used to select a section of linearly distributed velocity function in the area closest to the wall of the rotating boundary layer, and perform a linear fitting operation on the velocity function; 测量模块,用于对线性拟合操作获取的速度梯度依次通过平移、假设估计以及无量纲化操作完成测量。The measurement module is used to measure the velocity gradient obtained by the linear fitting operation through translation, assumption estimation and dimensionless operation in sequence.
CN201810923270.8A 2018-08-14 2018-08-14 Applicable to the method and device for measuring the shear stress on the wall of the rotating boundary layer Active CN108955981B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810923270.8A CN108955981B (en) 2018-08-14 2018-08-14 Applicable to the method and device for measuring the shear stress on the wall of the rotating boundary layer

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810923270.8A CN108955981B (en) 2018-08-14 2018-08-14 Applicable to the method and device for measuring the shear stress on the wall of the rotating boundary layer

Publications (2)

Publication Number Publication Date
CN108955981A true CN108955981A (en) 2018-12-07
CN108955981B CN108955981B (en) 2020-11-24

Family

ID=64469139

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810923270.8A Active CN108955981B (en) 2018-08-14 2018-08-14 Applicable to the method and device for measuring the shear stress on the wall of the rotating boundary layer

Country Status (1)

Country Link
CN (1) CN108955981B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109871601A (en) * 2019-01-30 2019-06-11 天津大学 Method of Equivalently Simulating Pipeline Water Flow State by Rotating Disk Electrode with Equal Shear Stress
CN110261009A (en) * 2019-07-09 2019-09-20 北京航空航天大学 The wall heat flux density measuring method and device assumed based on linear bottom

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4879899A (en) * 1988-04-27 1989-11-14 Massachusetts Institute Of Technology Shear stress gauge
CN105022919A (en) * 2015-07-06 2015-11-04 北京航空航天大学 Method for predicting size of low-speed impact dent of circular metal sheet
CN106441793A (en) * 2016-11-14 2017-02-22 北京航空航天大学 Adjusting method and apparatus of air-flow temperature in calibration wind tunnel
CN106643879A (en) * 2016-11-14 2017-05-10 北京航空航天大学 Temperature and speed synchronous measurement method and device
CN107561310A (en) * 2017-08-31 2018-01-09 北京航空航天大学 Rotate PIV method of testings and system
CN108279187A (en) * 2018-03-27 2018-07-13 苏州科技大学 The test device and test method of fluid viscosity

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4879899A (en) * 1988-04-27 1989-11-14 Massachusetts Institute Of Technology Shear stress gauge
CN105022919A (en) * 2015-07-06 2015-11-04 北京航空航天大学 Method for predicting size of low-speed impact dent of circular metal sheet
CN106441793A (en) * 2016-11-14 2017-02-22 北京航空航天大学 Adjusting method and apparatus of air-flow temperature in calibration wind tunnel
CN106643879A (en) * 2016-11-14 2017-05-10 北京航空航天大学 Temperature and speed synchronous measurement method and device
CN107561310A (en) * 2017-08-31 2018-01-09 北京航空航天大学 Rotate PIV method of testings and system
CN108279187A (en) * 2018-03-27 2018-07-13 苏州科技大学 The test device and test method of fluid viscosity

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
北京长城计量测试技术研究所: "《航空计量技术》", 31 December 2013, 航空工业出版社 *
魏宽等: "旋转状态下方形通道内部流场特性热线实验", 《航空动力学报》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109871601A (en) * 2019-01-30 2019-06-11 天津大学 Method of Equivalently Simulating Pipeline Water Flow State by Rotating Disk Electrode with Equal Shear Stress
CN109871601B (en) * 2019-01-30 2023-05-23 天津大学 A method for equivalently simulating water flow in pipelines with rotating disc electrodes with equal shear stress
CN110261009A (en) * 2019-07-09 2019-09-20 北京航空航天大学 The wall heat flux density measuring method and device assumed based on linear bottom

Also Published As

Publication number Publication date
CN108955981B (en) 2020-11-24

Similar Documents

Publication Publication Date Title
US10072961B2 (en) Level measuring instrument for determining moisture content
Schubring et al. A model for pressure loss, film thickness, and entrained fraction for gas–liquid annular flow
CN103808381B (en) A kind of temperature influence eliminating method of transit-time ultrasonic flow meter
CN108955981A (en) Suitable for rotation boundary layer wall surface shear stress measurement method and device
Djenidi et al. A spectral chart method for estimating the mean turbulent kinetic energy dissipation rate
CN103049670A (en) Pipe excitation source identification and prediction method of vibration response to pipe excitation source
CN106526225A (en) Single-frame particle imaging tachymeter suitable for high-speed flow filed and speed measurement method
CN111060043B (en) Ultrasonic thickness measurement method based on temperature compensation
Tang et al. Transport equation for the mean turbulent energy dissipation rate on the centreline of a fully developed channel flow
Örlü et al. High-order generalisation of the diagnostic scaling for turbulent boundary layers
JP2020143904A (en) Thermal flow meter and flow rate correction method
CN107561310A (en) Rotate PIV method of testings and system
CN110261009B (en) Wall surface heat flow density measuring method and device based on linear bottom layer hypothesis
JP2009236863A (en) Defect inspection method and device therefor
Ofuchi et al. Void fraction measurement in a Gas-Liquid swirling flow using an ultrasonic sensor
CN109000575B (en) Method and device for measuring distance between hot wire probe and wall surface
CN108511094B (en) Device and method for monitoring reactor core liquid level and cavitation
WO1991014161A1 (en) Thermal sensing
CN106643879B (en) Temperature and speed sync measuring method and device
CN211014087U (en) Simulation detection device for heat dissipation effect of low-thermal-resistance cable filling medium
CN203970408U (en) A kind of X-ray flat panel detector and X-ray machine
US20230105228A1 (en) Estimation device, estimation method, and non-transitory computer-readable recording medium for thickness of deposit
CN207068485U (en) A device for monitoring reactor core liquid level and cavitation
KR101522249B1 (en) a gas mass flow meter program using ultra sonic wave and the measuring device using thereof
CN204439796U (en) A kind of motor of specific radiating condition is to dragging test board

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20220126

Address after: 1701-3, block a, building 1, No. 10, Ronghua Middle Road, Beijing Economic and Technological Development Zone, Daxing District, Beijing 100176

Patentee after: Beijing Feiqing Technology Co.,Ltd.

Address before: 100191 No. 37, Haidian District, Beijing, Xueyuan Road

Patentee before: BEIHANG University

Patentee before: Beijing Ruiao Technology Co., Ltd