CN104006781B - 曲面法矢测量精度的计算方法 - Google Patents

曲面法矢测量精度的计算方法 Download PDF

Info

Publication number
CN104006781B
CN104006781B CN201410272692.5A CN201410272692A CN104006781B CN 104006781 B CN104006781 B CN 104006781B CN 201410272692 A CN201410272692 A CN 201410272692A CN 104006781 B CN104006781 B CN 104006781B
Authority
CN
China
Prior art keywords
error
curve
cross
section
max
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.)
Active
Application number
CN201410272692.5A
Other languages
English (en)
Other versions
CN104006781A (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.)
Tsinghua University
Chengdu Aircraft Industrial Group Co Ltd
Original Assignee
Tsinghua University
Chengdu Aircraft Industrial Group Co Ltd
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 Tsinghua University, Chengdu Aircraft Industrial Group Co Ltd filed Critical Tsinghua University
Priority to CN201410272692.5A priority Critical patent/CN104006781B/zh
Publication of CN104006781A publication Critical patent/CN104006781A/zh
Application granted granted Critical
Publication of CN104006781B publication Critical patent/CN104006781B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Length Measuring Devices With Unspecified Measuring Means (AREA)

Abstract

本发明提供了一种曲面法矢测量精度的计算方法,其包括步骤:S1,计算菱形布局距离传感器在截面YOZ中的误差:用截面YOZ截取待测曲面获得待测曲线,建立二维法矢测量误差模型,并将二维法矢测量误差δX分为由曲线曲率变化带来的误差δXP和由距离传感器测距带来的误差δXI;计算δXI的范围为假定距离传感器测距没有误差,获得tanδXP的范围为根据得出的δXI及δXP的范围,得到YOZ截面内测得的曲线法矢绕X转轴的最大误差δX=|δXI|max+|δXP|max;S2,计算菱形布局距离传感器在截面XOZ中的误差:在XOZ截面中采用与步骤S1相同的计算过程,得到XOZ截面内测得的曲线法矢绕Y转轴的最大误差δY=|δYI|max+|δYP|max;S3,计算三维法矢测量的精度:三维法矢测量的精度为Δ,则三维法矢测量的精度Δ计算式为(tanΔ)2=(tanδX)2+(tanδY)2

Description

曲面法矢测量精度的计算方法
技术领域
本发明涉及自动化装备及制造领域,尤其涉及一种曲面法矢测量精度的计算方法。
背景技术
在现代的自动化装备及制造中需要用到曲面法矢测量装置,其中常用的一种是利用四个距离传感器的曲面法矢测量装置。参照图1,该曲面法矢测量装置利用四个组成菱形的距离传感器3测量工件曲面4上待测点O’的法矢。该曲面法矢测量装置包括:法兰基准面1,安装于移动驱动装置(未示出)(例如经由连接孔11);四个支架2,安装在法兰基准面1上;四个距离传感器3,分别安装于四个支架2,使各距离传感器3的轴线垂直于法兰基准面1且顶表面设置有测头点31,用于检测工件曲面4的距离。其中,各距离传感器3的测头点31到法兰基准面1的垂直投影点的连线构成菱形ABCD的四个顶点,四个距离传感器3的测头点31到工件曲面的垂直投影点分别对应为A’、B’、C’、D’,且使菱形ABCD的中心点O为工件曲面上待测点O’在法兰基准面1上的投影;以菱形ABCD的中心点O为原点、菱形ABCD的对角线AC为Y轴及Y转轴、菱形ABCD的对角线BD为X轴及X转轴、菱形ABCD的中心点O与待测点O’之间的连线为Z轴建立测量坐标系;通过距离传感器3测得菱形ABCD的四个顶点到工件曲面4的垂直距离分别为AA’=hA、BB’=hB、CC’=hC、DD’=hD,且AO=CO=RY、BO=DO=RX;工件曲面4上待测点O’的法矢在测量坐标系中绕X转轴和Y转轴的偏角计算公式为:
α X = a r c t a n ( h D - h B 2 R X )
α Y = a r c t a n ( h C - h A 2 R Y )
精度是各类传感器最重要的指标。这种利用四个距离传感器3测量工件的曲面法矢测量装置,其工作对象是复杂曲面,而工作原理是将待测点O’附近区域看做小平面进行测量,因此会带来待测点O’的法矢测量误差;同时距离传感器3本身的测距误差,也会给待测点O’的法矢测量带来误差。
发明内容
鉴于背景技术中存在的问题,本发明的目的在于提供一种曲面法矢测量精度的计算方法,其能够量化地求出曲面法矢测量装置的精度。
为了实现上述目的,本发明提供了一种曲面法矢测量精度的计算方法,其包括步骤:S1,计算菱形布局距离传感器在截面YOZ中的误差;S2,计算菱形布局距离传感器在截面XOZ中的误差;以及S3,计算三维法矢测量的精度。
计算菱形布局距离传感器在截面YOZ中的误差的步骤S1包括子步骤:
S11,用截面YOZ截取待测曲面获得待测曲线,建立二维法矢测量误差模型,并将二维法矢测量误差δX分为由曲线曲率变化带来的误差δXP和由距离传感器测距带来的误差δXI:假定通过测量和调整使距离传感器A点和C点到工件曲面的垂直距离hA和hC读数相等,由于测距传感器误差,hA不严格等于hC,因此带来误差δXI;由于曲率变化,O’法矢并不垂直于A’C’,因此带来误差δXP
S12,计算δXI:假定测距传感器到中心的距离为RY,其测量最大误差为±δh,则δXI的范围为
S13,假定距离传感器测距没有误差,计算δXP
建立二维数学模型,假定待测曲线函数式为y=f(x),误差tanδXP=|f'(0)|;
针对由曲线曲率变化带来的误差δXI的部分,假定传感器测距无误差,即有f(-R)=f(R);通过微分中值定理可得f'(0)=ζf”(η),其中ζ∈(-R,R)使f'(ζ)=0,η∈(0,ζ);
曲率范围为|K|min≤|K|≤|K|max,为了保证测量高精度,设计距离传感器布局时保证R≤0.1×|r|min=(10×|K|max)-1,通过计算可求得|f”(η)|≤1.02|K|max从而获得tanδXP的范围为
S14,计算曲线法矢测量在截面YOZ中的最大误差:
根据得出的δXI及δXP的范围,得到YOZ截面内测得的曲线法矢绕X转轴的最大误差δX=|δXI|max+|δXP|max
计算菱形布局距离传感器在截面XOZ中的误差的步骤S2为:在XOZ截面中采用与步骤S1的子步骤S11-S14相同的计算过程,得到XOZ截面内测得的曲线法矢绕Y转轴的最大误差δY=|δYI|max+|δYP|max
计算三维法矢测量的精度的步骤S3为:三维法矢测量的精度为Δ,则三维法矢测量的精度Δ计算式为(tanΔ)2=(tanδX)2+(tanδY)2
本发明的有益效果如下:
本发明的曲面法矢测量精度的计算方法能够通过数学模型建立量化地求出上述曲面法矢测量装置的精度,而且所得精度数学表达式明确,可以帮助量化分析法矢测量精度的影响因素并对法矢测量装置的参数设计进行指导。
附图说明
图1为曲面法矢测量装置的立体图;
图2为曲面法矢测量装置测量工件曲面上待测点的法矢的示意图;
图3为根据本发明的曲面法矢测量精度的计算方法在YOZ平面中的二维法矢测量精度的示意图;
图4为根据本发明的曲面法矢测量精度的计算方法在YOZ平面中的二维法矢测量精度中针对δXP计算的示意图;
图5为根据本发明的曲面法矢测量精度的计算方法在YOZ平面中的二维法矢测量精度中针对δXP计算时求|f'(x)|max的示意图;
图6为根据本发明的曲面法矢测量精度的计算方法在YOZ平面中的二维法矢测量精度中针对δXP计算时y=f'(x)的函数曲线范围图;
图7为根据本发明的曲面法矢测量精度的计算方法在YOZ平面中的利用二维法矢测量精度计算三维法矢测量精度的示意图。
其中,附图标记说明如下:
1 法兰基准面 3 距离传感器
11 连接孔 31 测头点
2 支架 4 工件曲面
具体实施方式
下面参照附图来详细说明根据本发明的曲面法矢测量精度的计算方法。
参照图2至图7,根据本发明的曲面法矢测量精度的计算方法包括步骤:S1,计算菱形布局距离传感器在截面YOZ中的误差;S2,计算菱形布局距离传感器在截面XOZ中的误差;以及S3,计算三维法矢测量的精度。
计算菱形布局距离传感器在截面YOZ中的误差的步骤S1包括子步骤:
S11,用截面YOZ截取待测曲面获得待测曲线,建立二维法矢测量误差模型,并将二维法矢测量误差δX分为由曲线曲率变化带来的误差δXP和由距离传感器测距带来的误差δXI:假定通过测量和调整使距离传感器A点和C点到工件曲面的垂直距离hA和hC读数相等,由于测距传感器误差,hA不严格等于hC,因此带来误差δXI;由于曲率变化,O’法矢并不垂直于A’C’,因此带来误差δXP
S12,计算δXI:假定测距传感器到中心的距离为RY,其测量最大误差为±δh,则δXI的范围为
S13,假定距离传感器测距没有误差,计算δXP
建立二维数学模型,假定待测曲线函数式为y=f(x),误差tanδXP=|f'(0)|;
针对由曲线曲率变化带来的误差δXI的部分,假定传感器测距无误差,即有f(-R)=f(R);通过微分中值定理可得f'(0)=ζf”(η),其中ζ∈(-R,R)使f'(ζ)=0,η∈(0,ζ);
曲率范围为|K|min≤|K|≤|K|max,为了保证测量高精度,设计距离传感器布局时保证R≤0.1×|r|min=(10×|K|max)-1,通过计算可求得|f”(η)|≤1.02|K|max从而获得tanδXP的范围为
S14,计算曲线法矢测量在截面YOZ中的最大误差:
根据得出的δXI及δXP的范围,得到YOZ截面内测得的曲线法矢绕X转轴的最大误差δX=|δXI|max+|δXP|max
计算菱形布局距离传感器在截面XOZ中的误差的步骤S2为:在XOZ截面中采用与步骤S1的子步骤S11-S14相同的计算过程,得到XOZ截面内测得的曲线法矢绕Y转轴的最大误差δY=|δYI|max+|δYP|max
计算三维法矢测量的精度的步骤S3为:三维法矢测量的精度为Δ,则三维法矢测量的精度Δ计算式为(tanΔ)2=(tanδX)2+(tanδY)2
以下将详细说明上述曲面法矢测量精度的计算方法的过程。
参照图2,本发明所述的曲面法矢测量精度的计算方法针对一种曲面法矢测量装置,该曲面法矢测量装置利用四个组成菱形的距离传感器3测量工件曲面4上待测点O’的法矢。该曲面法矢测量装置测量工件曲面上待测点O’的法矢在测量坐标系中绕X转轴和Y转轴的偏角计算公式为:
用上述曲面法矢测量装置测量工件曲面4的法矢,其算法原理是将工件曲面4在待测点O’附近看做小平面,且距离传感器3本身具有误差,由此带来法矢测量的误差。因此,本发明的曲面法矢测量精度的计算方法针对以上问题,通过数学模型建立,量化地计算出法矢检测精度,可以计算出曲面法矢测量装置的精度指标,并能对精度影响因素的分析进行指导。
根据本发明的曲面法矢测量精度的计算原理为:用截面YOZ、XOZ截取待测曲面获得待测曲线,建立二维法矢测量误差模型,并将二维法矢测量误差δX、δY分为由曲线曲率变化带来的误差δXP、δYP和由距离传感器测距带来的误差δXI、δYI;各截面中的计算δXI、δYI;假定距离传感器测距没有误差,计算各截面中的δXP、δYP;基于两个截面的计算结果来计算三维法矢测量的精度。
以下将主要以截面YOZ为例进行详细说明。
参照图3,用截面YOZ截取待测曲面获得待测曲线;假定通过测量和调整使A点和C点距离传感器读数相等,此时认为测得法矢为Z轴方向,而实际法矢为O’J方向,即法矢测量的误差。由于测距传感器误差,hA不严格等于hC,因此带来误差δXI;由于曲率变化,O’法矢并不垂直于A’C’,因此带来误差δXP。其中,测距传感器到中心的距离为RY,其测量最大误差为±δh;因此计算可得δXI的范围为:
参照图4,建立二维数学模型,假定待测曲线函数式为y=f(x),所以误差tanδXP=|f'(0)|。针对由曲线曲率变化带来的误差δXI的部分,假定传感器测距无误差,即有f(-R)=f(R);通过微分中值定理可得,存在ζ∈(-R,R)使
且存在η∈(0,ζ)使
由式(3)及式(5)可得:
tanδXP=f'(0)=ζf”(η)……(式6)
工件曲面的曲率K范围为|K|min≤|K|≤|K|max,为了保证测量精度(±0.5°内),设计距离传感器布局时保证R≤0.1×|r|min=(10×|K|max)-1,其中r为曲率半径。
参照图5,当曲线为曲率恒为|K|max且R=(10×|K|max)-1时|f'(x)|能取到最大值,所以有:
根据曲率公式及曲率的范围并结合式(7)以及曲率和曲率半径之间的关系(即K×r=1),可得:
|f”(x)|≤1.02|K|max……(式8)
参照图6,由于式(8)所得的结论,所以曲线y=f'(x)在两条过点(ζ,0)、斜率为±1.02|K|max的直线之间,所以有:
又因为所以可得:
并代入式(11)可得:
又根据式(9)及式(10)可以求得
0<|K|min≤K1≤1.02|K|max……(式13)
0<|K|min≤K2≤1.02|K|max……(式14)
联合式(12)、式(13)及式(14)可得:
综上,即可求得tanδXP的范围为:
根据式(3)及式(16)中分别求出的δXI及δXP的范围,可得该YOZ截面内测得的曲线法矢绕X转轴的最大误差δX=|δXI|max+|δXP|max
对于截面XOZ中的误差计算,截取XOZ截面采用与在XOZ截面内的上述计算过程相同,从而同样可得到XOZ截面内测得的曲线法矢绕Y转轴的最大误差δY=|δYI|max+|δYP|max
在YOZ截面、XOZ截面的计算完成之后,参照图7,利用勾股定理,根据各三角表达式及|FO|2+|GO|2=|HO|2可得:
(tanΔ)2=(tanδX)2+(tanδY)2……(式17)
即可计算出三维法矢测量的精度Δ。
本发明的曲面法矢测量精度的计算方法能够通过数学模型建立量化地求出上述曲面法矢测量装置的精度,而且所得精度数学表达式明确,可以帮助分析法矢测量精度的影响因素并对法矢测量装置的参数设计进行指导。

Claims (1)

1.一种曲面法矢测量精度的计算方法,其中,曲面法矢测量装置包括:法兰基准面,安装于移动驱动装置;四个支架,安装在法兰基准面上;四个距离传感器,分别安装于四个支架,使各距离传感器的轴线垂直于法兰基准面且顶表面设置有测头点,用于检测工件曲面的距离,其中,各距离传感器的测头点到法兰基准面的垂直投影点的连线构成菱形ABCD的四个顶点,四个距离传感器的测头点到工件曲面的垂直投影点分别对应为A’、B’、C’、D’,且使菱形ABCD的中心点O为工件曲面上待测点O’在法兰基准面上的投影;以菱形ABCD的中心点O为原点、菱形ABCD的对角线AC为Y轴及Y转轴、菱形ABCD的对角线BD为X轴及X转轴、菱形ABCD的中心点O与待测点O’之间的连线为Z轴建立测量坐标系,通过距离传感器测得菱形ABCD的四个顶点到工件曲面的垂直距离分别为AA’=hA、BB’=hB、CC’=hC、DD’=hD,且AO=CO=RY、BO=DO=RX,其特征在于,所述曲面法矢测量精度的计算方法包括步骤:
S1,计算菱形布局距离传感器在截面YOZ中的误差,包括子步骤:
S11,用截面YOZ截取待测曲面获得待测曲线,建立二维法矢测量误差模型,并将二维法矢测量误差δX分为由曲线曲率变化带来的误差δXP和由距离传感器测距带来的误差δXI:假定通过测量和调整使距离传感器A点和C点到工件曲面的垂直距离hA和hC读数相等,由于距离传感器误差,hA不严格等于hC,因此带来误差δXI;由于曲率变化,O’法矢并不垂直于A’C’,因此带来误差δXP
S12,计算δXI:假定距离传感器到中心的距离为RY,其测量最大误差为±δh,则δXI的范围为
S13,假定距离传感器测距没有误差,计算δXP
建立二维数学模型,假定待测曲线函数式为y=f(x),误差tanδXP=|f′(0)|;
针对由曲线曲率变化带来的误差δXI的部分,假定传感器测距无误差,即有f(-R)=f(R);通过微分中值定理可得f′(0)=ζf″(η),其中ζ∈(-R,R)使f′(ζ)=0,η∈(0,ζ);
曲率范围为|K|min≤|K|≤|K|max,为了保证测量高精度,设计距离传感器布局时保证R≤0.1×|r|min=(10×|K|max)-1,其中r为曲率半径,通过计算可求得|f″(η)|≤1.02|K|max从而获得tanδXP的范围为
S14,计算曲线法矢测量在截面YOZ中的最大误差:
根据得出的δXI及δXP的范围,得到YOZ截面内测得的曲线法矢绕X转轴的最大误差δX=|δXI|max+|δXP|max
S2,计算菱形布局距离传感器在截面XOZ中的误差
在XOZ截面中采用与步骤S1的子步骤S11-S14相同的计算过程,得到XOZ截面内测得的曲线法矢绕Y转轴的最大误差δY=|δYI|max+|δYP|max
S3,计算三维法矢测量的精度
三维法矢测量的精度为Δ,则三维法矢测量的精度Δ计算式为(tanΔ)2=(tanδX)2+(tanδY)2
CN201410272692.5A 2014-06-18 2014-06-18 曲面法矢测量精度的计算方法 Active CN104006781B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410272692.5A CN104006781B (zh) 2014-06-18 2014-06-18 曲面法矢测量精度的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410272692.5A CN104006781B (zh) 2014-06-18 2014-06-18 曲面法矢测量精度的计算方法

Publications (2)

Publication Number Publication Date
CN104006781A CN104006781A (zh) 2014-08-27
CN104006781B true CN104006781B (zh) 2017-01-04

Family

ID=51367556

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410272692.5A Active CN104006781B (zh) 2014-06-18 2014-06-18 曲面法矢测量精度的计算方法

Country Status (1)

Country Link
CN (1) CN104006781B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112033307B (zh) * 2020-07-15 2021-08-03 成都飞机工业(集团)有限责任公司 一种法矢测量装置
CN112525131B (zh) * 2020-10-23 2021-09-14 清华大学 非接触式距离传感器阵列测量工件曲率的方法及***

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102430779A (zh) * 2011-09-22 2012-05-02 西北工业大学 一种测量自由曲面任意点处法向矢量的装置及其测量方法
CN102445171A (zh) * 2011-09-22 2012-05-09 西北工业大学 一种确定自由曲面上任意点处法向矢量的方法
CN102607484A (zh) * 2012-03-27 2012-07-25 西北工业大学 一种变间距自适应搜寻测量点的飞机壁板法矢确定方法
CN103433810A (zh) * 2013-07-19 2013-12-11 华中科技大学 一种复杂曲面法矢在机检测装置及方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5936039B2 (ja) * 2012-02-29 2016-06-15 三菱重工業株式会社 法線検出方法、法線検出装置および法線検出機能を備えた加工機

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102430779A (zh) * 2011-09-22 2012-05-02 西北工业大学 一种测量自由曲面任意点处法向矢量的装置及其测量方法
CN102445171A (zh) * 2011-09-22 2012-05-09 西北工业大学 一种确定自由曲面上任意点处法向矢量的方法
CN102607484A (zh) * 2012-03-27 2012-07-25 西北工业大学 一种变间距自适应搜寻测量点的飞机壁板法矢确定方法
CN103433810A (zh) * 2013-07-19 2013-12-11 华中科技大学 一种复杂曲面法矢在机检测装置及方法

Also Published As

Publication number Publication date
CN104006781A (zh) 2014-08-27

Similar Documents

Publication Publication Date Title
Zhong et al. Dynamic accuracy evaluation for five-axis machine tools using S trajectory deviation based on R-test measurement
CN109032069B (zh) 一种采用电涡流位移传感器的非接触式R-test测量仪球心坐标计算方法
WO2021128424A1 (zh) 绝对电角度检测方法、***及计算机可读存储介质
JP6149337B1 (ja) 表面形状測定装置
US8290740B2 (en) Form measuring instrument, form measuring method, and program
CN104142209B (zh) 质量重心测量装置及方法
CN105473981A (zh) 接触式探针的校准
Zakharov et al. Minimization of the systematic error in centerless measurement of the roundness of parts
CN107063060A (zh) 一种确定表面平面度的方法及装置
KR20190137131A (ko) 머신 비전 시스템
CN104897105B (zh) 直线形状测定方法及直线形状测定装置
JP5652631B2 (ja) 真円度測定装置における心ずれ量算出方法
CN106767443A (zh) 一种新型全自动二次元影像检测仪及测量方法
CN103363923A (zh) 一种激光视觉测距的非接触式齿轮齿向测量方法
Janecki et al. Problems of measurement of barrel-and saddle-shaped elements using the radial method
CN104567773B (zh) 臂型三维测量机和支撑该测量机的基部的倾斜校正方法
Yalovoy et al. The centerless measurement of roundness with optimal adjustment
CN104006781B (zh) 曲面法矢测量精度的计算方法
JP2017161252A (ja) 表面形状測定方法及び表面形状測定装置
CN105094051B (zh) 运动平台***的平面定位补偿方法
JP4890188B2 (ja) 運動誤差測定基準体及び運動誤差測定装置
JP2008268118A (ja) 形状測定方法及び装置
CN106247997B (zh) 用于正交导轨平台的垂直度误差测量方法
CN104930953B (zh) 机械零件圆柱度误差的四点法测量装置及测量方法
JP2010256277A (ja) 被測定物の外形測定方法及び外形測定装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C41 Transfer of patent application or patent right or utility model
CB03 Change of inventor or designer information

Inventor after: Wu Dan

Inventor after: Yang Xiangdong

Inventor after: Fu Chenglong

Inventor after: Liu Li

Inventor after: Xu Jing

Inventor after: Xie Ying

Inventor after: Chen Ken

Inventor after: Gao Yuhao

Inventor after: Zheng Linbin

Inventor after: He Fengtao

Inventor after: Xu Jing

Inventor after: Wang Guolei

Inventor after: Guo Xifeng

Inventor after: Song Libin

Inventor before: Chen Ken

Inventor before: Yang Dongchao

Inventor before: Gao Yuhao

Inventor before: Wu Dan

Inventor before: Wang Guolei

Inventor before: Song Libin

Inventor before: Yang Xiangdong

Inventor before: Fu Chenglong

Inventor before: Xu Jing

Inventor before: Liu Li

COR Change of bibliographic data
TA01 Transfer of patent application right

Effective date of registration: 20160108

Address after: 100084 Beijing City, Haidian District Tsinghua Yuan

Applicant after: Tsinghua University

Applicant after: Chengdu Airplane Industry (Group) Corp., Ltd.

Address before: 100084 Beijing City, Haidian District Tsinghua Yuan

Applicant before: Tsinghua University

C14 Grant of patent or utility model
GR01 Patent grant