CN103940578B - 一种基于应变信号近似熵的风洞试验天平的评估方法 - Google Patents
一种基于应变信号近似熵的风洞试验天平的评估方法 Download PDFInfo
- Publication number
- CN103940578B CN103940578B CN201410164546.0A CN201410164546A CN103940578B CN 103940578 B CN103940578 B CN 103940578B CN 201410164546 A CN201410164546 A CN 201410164546A CN 103940578 B CN103940578 B CN 103940578B
- Authority
- CN
- China
- Prior art keywords
- balance
- strain
- model
- signal
- dimensional
- 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.)
- Expired - Fee Related
Links
Landscapes
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
一种基于应变信号近似熵的风洞试验天平的评估方法,其做法主要是,在飞行器模型顶端内壁、模型框架和天平上分别安装模型、框架、天平三向应变片,对三个片测出的各向应变信号进行傅里叶变换得到频谱信号,进而计算出0~300Hz内六个频带的近似熵特征值,再分别计算模型、框架同向应变信号的近似熵特征值的差值,以及模型与天平同向应变信号的近似熵特征值的差值;当各向的这两种应变近似熵特征值的差值均在规定范围内时,评估结果判定天平的测试数据可信,否则判定不可信。从而保证风洞试验时通过天平测出的模型的力学数据准确、可靠;为航空航天飞行器提供更准确、可靠的试验数据。
Description
技术领域
本发明涉及航空航天飞行器风洞试验天平的校准与评估方法,尤其涉及一种基于应变信号近似熵的风洞试验天平的评估方法。
背景技术
航空航天飞行器的研制过程中,为了了解飞行器性能、降低飞行器研制风险和成本,通常会进行风洞模型试验。在各种试验条件下对飞行器试验模型受到的气动力大小、力矩进行测量,是风洞试验中最基本的实验项目。其测试原理是,在风洞的底部固定一模型支座,模型支座的上端穿入试验模型中与天平的固定端螺纹连接,天平的测量端与试验模型中的模型框架螺纹连接;模型框架又通过螺栓固定连接试验模型。试验时,往风洞通入特定风速的气流,试验模型受到相应的气动力和力矩的作用,该气动力和力矩通过模型框架传递给支座上端的天平,由天平感应和测量出气流作用在模型上的沿垂直坐标分解的三个气动力和力矩分量。由天平感应和测量出六个自由度上的气动力和力矩。将天平(受到)测出的气动力和力矩推断为试验模型受到的气动力和力矩。
由于气动力和力矩从试验模型向模型框架及天平传递时,试验模型与模型框架的连接螺栓,模型框架与天平的连接螺栓,模型框架本身均会对传递的气动力和力矩造成衰减和干扰,导致风洞模型试验中,气动力和力矩的测试误差大,精度低。因此,需要对天平校正或评估。
现有的天平校准技术多为在天平校准设备上,对天平施加标准载荷,并采集天平输出信号,计算天平各分量的输出信号与标准载荷的变化关系。这种校准方法周期长、步骤多,且只能对天平的静态特性进行校准,无法对实际工况下天平的动态特性进行评估,导致实际工况下的风洞试验的力学数据的准确度和可靠性难以确定,不能准确、可靠的评价航空航天飞行器的性能,也不能为航空航天飞行器提供准确、可靠的试验数据。
发明内容
本发明的目的是提供一种基于应变信号近似熵的风洞试验天平的评估方法,该方法能够实现实际工况下的航空航天飞行器风洞模型试验中天平动态输出特性的评估,使航空航天飞行器风洞模型试验的气动力和力矩的测试误差在允许的范围内,风洞实验测出的力学数据更准确、可靠;从而能更准确、可靠的评价航空航天飞行器的性能,也能为航空航天飞行器提供更准确、可靠的试验数据。
本发明为实现其发明目的所采用的技术方案是:一种基于应变信号近似熵的风洞试验天平的评估方法,其步骤如下:
A、试验前,在飞行器模型的顶端内壁、模型框架和天平上分别安装模型三向应变片、框架三向应变片和天平三向应变片,模型三向应变片、框架三向应变片和天平三向应变片的信号输出端均与信号调理仪相连;信号调理仪还通过数据采集设备与计算机相连;
B、模拟实际工况,对飞行器模型进行风洞试验;试验的同时,模型三向应变片、框架三向应变片和天平三向应变片将各自测出的应变信号fpv(t)以采样频率L通过信号调理仪送入计算机;其中p=1,2,3,分别表示模型三向应变片、框架三向应变片和天平三向应变片上采集的应变信号,v=X,Y,Z分别表示三向应变片在X、Y、Z方向上采集到的应变信号,t=1,2,3,…,N,代表应变信号fpv(t)的采样点的序列数,N为总采样点数;
C、计算机通过傅里叶变换将每路应变信号fpv(t)变换成频谱信号φpv(ω);其中ω=1,2,…,L/2为频率;并将每路频谱信号φpv(ω)在0~300Hz频率段均分成六个频带,提取频谱信号φpv(ω)在频带i的应变信号的近似熵特征值i=1,2,…,6,i为频带的序号;
D、将由模型三向应变片得到的X、Y、Z方向上的应变信号在频带i的近似熵特征值与由框架三向应变片得到的X、Y、Z方向上在频带i的应变信号的近似熵特征值分别按下式计算得到X、Y、Z方向上在频带i的模型应变与框架应变的相对差值
式中|·|表示求绝对值;
将由模型三向应变片(1a)得到的X、Y、Z方向上的应变信号在频带i的近似熵特征值与由天平三向应变片(3a)得到的X、Y、Z方向上在频带i的应变信号的近似熵特征值分别按下式计算得到X、Y、Z方向上在频带i的模型应变与天平应变的相对差值
E、若X、Y、Z方向上在频带i的模型应变与框架应变的相对差值中的最大值小于5%,且X、Y、Z方向上在频带i的模型应变与天平应变的相对差值中的最大值小于10%,则判定所有的气动力和力矩从试验模型通过模型框架传递到天平的误差在允许范围内,风洞试验的天平测试数据可信;否则,判定气动力和力矩从试验模型通过模型框架传递到天平的误差超出允许范围,风洞试验的天平测试数据不可信。
与现有技术相比,本发明的有益效果是:
一、在试验模型,模型框架和天平上发生的振动(应变)信号与其受到的气动力和力矩成线性正相关,即在三处安装的三向应变片测出的应变值与各自受到的力成线性正相关。因此,模型应变片与框架应变片测出的同一方向的应变近似熵特征值的差异反映了气动力和力矩从试验模型传递到模型框架上的衰减和干扰;模型应变片与天平应变片测出的同一方向的应变近似熵特征值的差异反映了气动力和力矩从试验模型通过框架传递到天平上的衰减和干扰;当这两种各向应变近似熵特征值差异均在规定范围内时,表明所有的气动力和力矩从试验模型通过模型框架传递到天平的衰减和干扰均在允许范围内,评估结果判定测试数据可信,否则评估判定测试结果不可信,需要对天平和模型框架进行改进、校准。从而保证风洞试验时通过天平测出的模型的力学数据准确、可靠;使风洞试验能更准确、可靠的评价航空航天飞行器的性能,也能为航空航天飞行器提供更准确、可靠的试验数据。
二、本发明的方法不需要对天平结构以及影响天平精度的因素进行严格的分析建模,只需要加装三个三向应变片及其后续处理设备,通过试验与数据处理便可以快速评估天平的测量准确性。设备安装简便,可重复利用,节约了试验成本,缩短试验周期,易于广泛地推广应用于天平性能的检测试验,对航空航天飞行器的研制均具有重要的意义。
上述的步骤C中提取频谱信号φpv(ω)在频带i的应变信号的近似熵特征值的具体做法是:
C1、提取每路频谱信号φpv(ω)在频带i内的频谱信号其中m∈M,M为频谱信号的频率点的集合,M={1,2,…M0},M0为频谱信号的最大频率点;并计算出频带i内的频谱信号的标准差
C2、令j∈M且j≠M0,得到频率点j的二维矢量
C3、计算频率点j的二维矢量与其余频率点b的二维矢量之间的距离djb,2,
式中,b∈M且b≠j,b≠M0;
C4、统计频率点j的二维矢量与其余频率点b的二维矢量之间的距离djb,2小于0.2倍频谱信号的标准差的个数,随后算出该个数与距离djb,2的总个数M0-2的比值,即频率点j与其余频率点b的二维矢量关联度
C5、计算每路频谱信号φpv(ω)在频带i内的频谱信号的二维矢量序列的平均自相关度
C6、令k∈M且k≠M0,k≠M0-1得到频率点k的三维矢量
C7、计算频率点k的三维矢量与其余频率点c的三维矢量之间的距离dkc,3,
其中,c∈M且c≠k,c≠M0c≠M0-1;
C8、统计频率点k的三维矢量与其余频率点c的三维矢量之间的距离dkc,3小于0.2倍频谱信号的标准差的个数,随后算出该个数与距离dkc,3的总个数M0-3的比值,即频率点k的频率点k与其余频率点c的三维矢量关联度
C9、计算每路频谱信号φpv(ω)在频带i内的频谱信号的三维矢量序列的平均自相关度
C10、将C5步得到的二维矢量序列的平均自相关度与C9步得到的三维矢量序列的平均自相关度相减,即得到频谱信号φpv(ω)在频带i的应变信号的近似熵特征值即
采用如上方法计算应变信号近似熵,其计算简便,同时也可以统一应变信号的量纲,便于考察各路信号之间的近似熵差异;客观反映气动力和力矩在传递过程中的衰减和干扰。
下面结合附图和具体实施方式对本发明作进一步详细的描述。
附图说明
图1是用本发明实施例的方法对风洞试验天平进行评估时的相关器件的安装结构示意图。
图1中1b为模型支座,箭头方向为风洞的气流方向。
具体实施方式
实施例
图1示出,本发明的一种具体实施方式是,一种基于应变信号近似熵的风洞试验天平的评估方法,其步骤如下:
A、试验前,在飞行器模型1的顶端内壁、模型框架2和天平3上分别安装模型三向应变片1a、框架三向应变片2a和天平三向应变片3a,模型三向应变片1a、框架三向应变片2a和天平三向应变片3a的信号输出端均与信号调理仪相连;信号调理仪还通过数据采集设备与计算机相连;
B、模拟实际工况,对飞行器模型进行风洞试验;试验的同时,模型三向应变片1a、框架三向应变片2a和天平三向应变片3a将各自测出的应变信号fpv(t)以采样频率L通过信号调理仪送入计算机;其中p=1,2,3,分别表示模型三向应变片1a、框架三向应变片2a和天平三向应变片3a上采集的应变信号,v=X,Y,Z分别表示三向应变片在X、Y、Z方向上采集到的应变信号,t=1,2,3,…,N,代表应变信号fpv(t)的采样点的序列数,N为总采样点数;
C、计算机通过傅里叶变换将每路应变信号fpv(t)变换成频谱信号φpv(ω);其中ω=1,2,…,L/2为频率;并将每路频谱信号φpv(ω)在0~300Hz频率段均分成六个频带,提取频谱信号φpv(ω)在频带i的应变信号的近似熵特征值i=1,2,…,6,i为频带的序号;
D、将由模型三向应变片1a得到的X、Y、Z方向上的应变信号在频带i的近似熵特征值与由框架三向应变片2a得到的X、Y、Z方向上在频带i的应变信号的近似熵特征值分别按下式计算得到X、Y、Z方向上在频带i的模型应变与框架应变的相对差值
式中|·|表示求绝对值;
将由模型三向应变片1a得到的X、Y、Z方向上的应变信号在频带i的近似熵特征值与由天平三向应变片3a得到的X、Y、Z方向上在频带i的应变信号的近似熵特征值分别按下式计算得到X、Y、Z方向上在频带i的模型应变与天平应变的相对差值
E、若X、Y、Z方向上在频带i的模型应变与框架应变的相对差值中的最大值小于5%,且X、Y、Z方向上在频带i的模型应变与天平应变的相对差值中的最大值小于10%,则判定所有的气动力和力矩从试验模型通过模型框架传递到天平的误差在允许范围内,风洞试验的天平测试数据可信;否则,判定气动力和力矩从试验模型通过模型框架传递到天平的误差超出允许范围,风洞试验的天平测试数据不可信。
本例的步骤C中提取频谱信号φpv(ω)在频带i的应变信号的近似熵特征值的具体做法是:
C1、提取每路频谱信号φpv(ω)在频带i内的频谱信号其中m∈M,M为频谱信号的频率点的集合,M={1,2,…M0},M0为频谱信号的最大频率点;并计算出频带i内的频谱信号的标准差
C2、令j∈M且j≠M0,得到频率点j的二维矢量
C3、计算频率点j的二维矢量与其余频率点b的二维矢量之间的距离djb,2,
式中,b∈M且b≠j,b≠M0;
C4、统计频率点j的二维矢量与其余频率点b的二维矢量之间的距离djb,2小于0.2倍频谱信号的标准差的个数,随后算出该个数与距离djb,2的总个数M0-2的比值,即频率点j与其余频率点b的二维矢量关联度
C5、计算每路频谱信号φpv(ω)在频带i内的频谱信号的二维矢量序列的平均自相关度
C6、令k∈M且k≠M0,k≠M0-1得到频率点k的三维矢量
C7、计算频率点k的三维矢量与其余频率点c的三维矢量之间的距离dkc,3,
其中,c∈M且c≠k,c≠M0c≠M0-1;
C8、统计频率点k的三维矢量与其余频率点c的三维矢量之间的距离dkc,3小于0.2倍频谱信号的标准差的个数,随后算出该个数与距离dkc,3的总个数M0-3的比值,即频率点k的频率点k与其余频率点c的三维矢量关联度
C9、计算每路频谱信号φpv(ω)在频带i内的频谱信号的三维矢量序列的平均自相关度
C10、将C5步得到的二维矢量序列的平均自相关度与C9步得到的三维矢量序列的平均自相关度相减,即得到频谱信号φpv(ω)在频带i的应变信号的近似熵特征值即
本发明采用的三向应变片和信号调理仪可选用各种现有的应变片和调理仪,如应变片可选用中航工业电测仪器股份有限公司的BE120-2CA(11)三向应变片,应变调理仪可选用北京东方振动和噪声技术研究所INV1861A便携式应变调理仪,数据采集设备选用北京东方振动和噪声技术研究所的INV3020C信号采集分析仪。
Claims (2)
1.一种基于应变信号近似熵的风洞试验天平的评估方法,其步骤如下:
A、试验前,在飞行器模型(1)的顶端内壁安装模型三向应变片(1a),在模型框架(2)上安装框架三向应变片(2a),以及在天平(3)上安装天平三向应变片(3a),模型三向应变片(1a)、框架三向应变片(2a)和天平三向应变片(3a)的信号输出端均与信号调理仪相连;信号调理仪还通过数据采集设备与计算机相连;
B、模拟实际工况,对飞行器模型进行风洞试验;试验的同时,模型三向应变片(1a)、框架三向应变片(2a)和天平三向应变片(3a)将各自测出的应变信号fpv(t)以采样频率L通过信号调理仪送入计算机;其中p=1,2,3,分别表示模型三向应变片(1a)、框架三向应变片(2a)和天平三向应变片(3a)上采集的应变信号,v=X,Y,Z,分别表示三向应变片在X、Y、Z方向上采集到的应变信号,t=1,2,3,…,N,代表应变信号fpv(t)的采样点的序列数,N为总采样点数;
C、计算机通过傅里叶变换将每路应变信号fpv(t)变换成频谱信号φpv(ω);其中ω=1,2,…,L/2为频率;并将每路频谱信号φpv(ω)在0~300Hz频率段均分成六个频带,提取频谱信号φpv(ω)在频带i的应变信号的近似熵特征值i=1,2,…,6,i为频带的序号;
D、将由模型三向应变片(1a)得到的X、Y、Z方向上的应变信号在频带i的近似熵特征值与由框架三向应变片(2a)得到的X、Y、Z方向上在频带i的应变信号的近似熵特征值分别按下式计算得到X、Y、Z方向上在频带i的模型应变与框架应变的相对差值
式中|·|表示求绝对值;
将由模型三向应变片(1a)得到的X、Y、Z方向上的应变信号在频带i的近似熵特征值与由天平三向应变片(3a)得到的X、Y、Z方向上在频带i的应变信号的近似熵特征值分别按下式计算得到X、Y、Z方向上在频带i的模型应变与天平应变的相对差值
E、若X、Y、Z方向上在频带i的模型应变与框架应变的相对差值中的最大值小于5%,且X、Y、Z方向上在频带i的模型应变与天平应变的相对差值中的最大值小于10%,则判定所有的气动力和力矩从试验模型通过模型框架传递到天平的误差在允许范围内,风洞试验的天平测试数据可信;否则,判定气动力和力矩从试验模型通过模型框架传递到天平的误差超出允许范围,风洞试验的天平测试数据不可信。
2.根据权利要求1所述的一种基于应变信号近似熵的风洞试验天平的评估方法,其特征在于:所述的步骤C中提取频谱信号φpv(ω)在频带i的应变信号的近似熵特征值的具体做法是:
C1、提取每路频谱信号φpv(ω)在频带i内的频谱信号其中m∈M,M为频谱信号φi pv(m)的频率点的集合,M={1,2,…M0},M0为频谱信号的最大频率点;并计算出频带i内的频谱信号的标准差
C2、令j∈M且j≠M0,得到频率点j的二维矢量
C3、计算频率点j的二维矢量与其余频率点b的二维矢量之间的距离djb,2,
式中,b∈M且b≠j,b≠M0;
C4、统计频率点j的二维矢量与其余频率点b的二维矢量之间的距离djb,2小于0.2倍频谱信号的标准差的个数,随后算出该个数与距离djb,2的总个数M0-2的比值,即频率点j与其余频率点b的二维矢量关联度
C5、计算每路频谱信号φpv(ω)在频带i内的频谱信号的二维矢量序列的平均自相关度
C6、令k∈M且k≠M0,k≠M0-1得到频率点k的三维矢量
C7、计算频率点k的三维矢量与其余频率点c的三维矢量之间的距离dkc,3,
其中,c∈M且c≠k,c≠M0,c≠M0-1;
C8、统计频率点k的三维矢量与其余频率点c的三维矢量之间的距离dkc,3小于0.2倍频谱信号的标准差的个数,随后算出该个数与距离dkc,3的总个数M0-3的比值,即频率点k的频率点k与其余频率点c的三维矢量关联度
C9、计算每路频谱信号φpv(ω)在频带i内的频谱信号的三维矢量序列的平均自相关度
C10、将C5步得到的二维矢量序列的平均自相关度与C9步得到的三维矢量序列的平均自相关度相减,即得到频谱信号φpv(ω)在频带i的应变信号的近似熵特征值即
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410164546.0A CN103940578B (zh) | 2014-04-22 | 2014-04-22 | 一种基于应变信号近似熵的风洞试验天平的评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410164546.0A CN103940578B (zh) | 2014-04-22 | 2014-04-22 | 一种基于应变信号近似熵的风洞试验天平的评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103940578A CN103940578A (zh) | 2014-07-23 |
CN103940578B true CN103940578B (zh) | 2017-02-01 |
Family
ID=51188342
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410164546.0A Expired - Fee Related CN103940578B (zh) | 2014-04-22 | 2014-04-22 | 一种基于应变信号近似熵的风洞试验天平的评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103940578B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1260481A (zh) * | 2000-01-19 | 2000-07-19 | 山西省气象局 | 自动检测和自动控制风洞风速的方法及装置 |
JP3787594B2 (ja) * | 2003-03-31 | 2006-06-21 | 川崎重工業株式会社 | 超音速風洞設備のノズル形状調整方法 |
CN101968399A (zh) * | 2010-10-18 | 2011-02-09 | 西南交通大学 | 一种车辆、桥梁模型风洞试验系统 |
CN102650565A (zh) * | 2012-04-24 | 2012-08-29 | 中国空气动力研究与发展中心高速空气动力研究所 | 风洞模拟实验中涡轮动力模拟器短舱唇口及其设计方法 |
RU2477460C1 (ru) * | 2011-07-14 | 2013-03-10 | Российская Федерация, от имени которой выступает Министерство промышленности и торговли Российской Федерации (Минпромторг России) | Способ определения коэффициентов аэродинамических сил и моментов при установившемся вращении модели летательного аппарата и устройство для его осуществления |
-
2014
- 2014-04-22 CN CN201410164546.0A patent/CN103940578B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1260481A (zh) * | 2000-01-19 | 2000-07-19 | 山西省气象局 | 自动检测和自动控制风洞风速的方法及装置 |
JP3787594B2 (ja) * | 2003-03-31 | 2006-06-21 | 川崎重工業株式会社 | 超音速風洞設備のノズル形状調整方法 |
CN101968399A (zh) * | 2010-10-18 | 2011-02-09 | 西南交通大学 | 一种车辆、桥梁模型风洞试验系统 |
RU2477460C1 (ru) * | 2011-07-14 | 2013-03-10 | Российская Федерация, от имени которой выступает Министерство промышленности и торговли Российской Федерации (Минпромторг России) | Способ определения коэффициентов аэродинамических сил и моментов при установившемся вращении модели летательного аппарата и устройство для его осуществления |
CN102650565A (zh) * | 2012-04-24 | 2012-08-29 | 中国空气动力研究与发展中心高速空气动力研究所 | 风洞模拟实验中涡轮动力模拟器短舱唇口及其设计方法 |
Non-Patent Citations (1)
Title |
---|
风洞虚拟飞行试验模拟方法研究;李浩 等;《实验流体力学》;20111231;第25卷(第6期);第72-76页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103940578A (zh) | 2014-07-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102539107B (zh) | 一种实现风洞试验信号精确同步的方法 | |
CN103940577B (zh) | 一种基于加速度信号能量的风洞试验天平的评估方法 | |
CN103174122B (zh) | 用于测试土体静止侧压力系数的侧向应力孔压探头 | |
CN103940575B (zh) | 一种基于应变信号能量的风洞试验天平的评估方法 | |
CN107607247B (zh) | 一种炸药爆炸冲量和风动压联合测试方法 | |
CN104343709B (zh) | 一种检测风机故障的装置及方法 | |
CN107092747A (zh) | 树木风载荷测定方法 | |
CN104792450A (zh) | 一种设备对安装基座激励力的间接测量方法 | |
RU2561829C2 (ru) | Способ и устройство для определения аэродинамических характеристик летательного аппарата | |
CN105758602A (zh) | 一种桁架梁桥断面抖振力同步测量方法 | |
CN106610488A (zh) | 一种雷达探测威力的测量方法 | |
CN107544454B (zh) | 嵌入式大气数据传感系统信息融合地面验证系统 | |
CN104089791A (zh) | 基于振动的改进损伤定位和损伤程度识别方法 | |
CN104567771B (zh) | 无约束油气管线横向运动过程中竖向位移测定装置 | |
CN106768818B (zh) | 一种激波风洞中混合气体来流运行参数获得方法 | |
CN104111157B (zh) | 一种检查风洞测量设备初读数的方法 | |
CN105092132A (zh) | 一种救生伞动态开伞力多方向的测试方法及装置 | |
CN108152130A (zh) | 三边法风电叶片静力测试方法 | |
CN105865688A (zh) | 螺杆检测用张拉仪及其应用 | |
CN207850594U (zh) | 钢弦式锚杆测力计检定装置 | |
CN103940576B (zh) | 一种基于加速度信号近似熵的风洞试验天平的评估方法 | |
CN106324097A (zh) | 一种基于声波传感测定材料碰撞恢复系数的方法 | |
CN103940578B (zh) | 一种基于应变信号近似熵的风洞试验天平的评估方法 | |
CN109470888A (zh) | 基于深度学习的高g值加速度计的标定系统及标定方法 | |
CN104792656A (zh) | 一种实时测量液氢中仲氢含量的方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170201 Termination date: 20190422 |
|
CF01 | Termination of patent right due to non-payment of annual fee |