CN115831296B - 一种材料摩擦系数的不确定度的计算方法及应用 - Google Patents
一种材料摩擦系数的不确定度的计算方法及应用 Download PDFInfo
- Publication number
- CN115831296B CN115831296B CN202310138721.8A CN202310138721A CN115831296B CN 115831296 B CN115831296 B CN 115831296B CN 202310138721 A CN202310138721 A CN 202310138721A CN 115831296 B CN115831296 B CN 115831296B
- Authority
- CN
- China
- Prior art keywords
- uncertainty
- friction
- force
- normal force
- standard
- 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
Links
- 239000000463 material Substances 0.000 title claims abstract description 44
- 238000004364 calculation method Methods 0.000 title claims abstract description 16
- 238000000034 method Methods 0.000 claims abstract description 62
- 238000005259 measurement Methods 0.000 claims abstract description 58
- 238000012360 testing method Methods 0.000 claims description 25
- 238000005498 polishing Methods 0.000 claims description 18
- 238000012544 monitoring process Methods 0.000 claims description 11
- 239000004065 semiconductor Substances 0.000 claims description 10
- 230000002194 synthesizing effect Effects 0.000 claims description 6
- 229910000831 Steel Inorganic materials 0.000 claims description 5
- 239000010959 steel Substances 0.000 claims description 5
- 238000000227 grinding Methods 0.000 claims description 4
- 230000003068 static effect Effects 0.000 claims description 2
- 238000012545 processing Methods 0.000 abstract description 6
- 230000008859 change Effects 0.000 description 9
- 230000008569 process Effects 0.000 description 7
- 238000002474 experimental method Methods 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 239000000758 substrate Substances 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 241000577395 Thenus Species 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000005672 electromagnetic field Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000012854 evaluation process Methods 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- GBPVMEKUJUKTBA-UHFFFAOYSA-N methyl 2,2,2-trifluoroethyl carbonate Chemical compound COC(=O)OCC(F)(F)F GBPVMEKUJUKTBA-UHFFFAOYSA-N 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000010561 standard procedure Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000013076 uncertainty analysis Methods 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明涉及一种材料摩擦系数的不确定度的计算方法及应用,其包括:获取待测材料的摩擦力F的测量值及法向力N的测量值,以计算摩擦系数μ;分别获取摩擦力F的标准不确定度u F 及法向力N的标准不确定度u N ;通过全微分形式建立摩擦系数μ与摩擦力F和摩擦力F的标准不确定度u F 、以及法向力N和法向力N的标准不确定度u N 的关系,得到摩擦系数的不确定度u μ 。本发明能实现测量结果的可溯源性;可以在电数字数据处理领域中应用。
Description
技术领域
本发明涉及一种电数字数据处理技术领域,特别是关于一种材料摩擦系数的不确定度的计算方法及应用。
背景技术
材料的摩擦系数测定在基础科学与工业生产中有着重要的意义,涉及到系统能量的转移与耗散,是进行能量传递控制的重要参数。在半导体基材加工过程中,结构表面超精密加工、加工设备的主承载面的设计与加工都需要准确的摩擦力、摩擦系数测定,特别是在基于摩擦扭矩的基材抛光终点检测中更为重要。摩擦系数的测定在不同的领域有着自己的国家测量标准,但对于测量结果的不确定度评估标准方法却没有完全建立,更多的不确定度评估是基于误差的分析,无法在分布、溯源链条给出更为准确的指导方法,导致各行业中关于摩擦系数的测量结果都有较大的浮动范围。
发明内容
针对上述问题,本发明的目的是提供一种材料摩擦系数的不确定度的计算方法及应用,其能实现测量结果的可溯源性。
为实现上述目的,一方面,本发明采取的技术方案为:一种材料摩擦系数的不确定度的计算方法,其包括:获取待测材料的摩擦力F的测量值及法向力N的测量值,以计算摩擦系数μ;分别获取摩擦力F的标准不确定度u F 及法向力N的标准不确定度u N ;通过全微分形式建立摩擦系数μ与摩擦力F和摩擦力F的标准不确定度u F 、以及法向力N和法向力N的标准不确定度u N 的关系,得到摩擦系数的不确定度u μ 。
进一步,所述获取待测材料的摩擦力F的测量值及法向力N的测量值,包括:由设置的测试系统获取;
所述测试系统包括加载盘系统、上摩擦副、法向加载装置和输出监测系统;所述加载盘系统用于盛放所述待测材料,所述上摩擦副与所述加载盘系统接触以获取所述待测材料的摩擦力F,所述法向加载装置与所述加载盘系统连接以获取施加在所述待测材料上的法向力N,所述输出监测系统用于输出摩擦力F及法向力N的测量值。
进一步,所述计算摩擦力F的标准不确定度u F ,包括:
进一步,所述摩擦力F的标准不确定度u F 为:
进一步,所述计算法向力N的标准不确定度u N ,包括:
进一步,所述法向力N的标准不确定度u N 为:
进一步,所述摩擦力F与所述法向力N的相关系数为零时,所述摩擦系数的不确定度u μ 为:
另一方面,本发明采取的技术方案为:基于上述的方法在半导体材料制备中的应用。
本发明由于采取以上技术方案,其具有以下优点:
1、本发明针对材料摩擦系数测量需求,建立相应的不确定度分析方法,进而实现测量结果的可溯源性。建立针对材料摩擦系数的溯源链条,实现测量的量值可以溯源至国际单位制的基本物理量,并最终溯源至物理常数。
2、本发明通过分解影响摩擦系数测量准确度的因素,考虑各因素分布特征,得到完整的不确定度计算方法。
附图说明
图1是本发明一实施例中材料摩擦系数的不确定度的计算方法流程图;
图2是本发明一实施例中测试系统的结构示意图;
图3是本发明一实施例中测试系统中对摩擦系数不确定度的具体影响量;
图4是本发明一实施例中不确定度分量量值溯源体系图;
附图标记:
1-上抛磨盘,2-下抛磨盘,3-待抛磨材料,4-侧向力传感器,5-连接钢带,6-法向力传感器。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例的附图,对本发明实施例的技术方案进行清楚、完整地描述。显然,所描述的实施例是本发明的一部分实施例,而不是全部的实施例。基于所描述的本发明的实施例,本领域普通技术人员所获得的所有其他实施例,都属于本发明保护的范围。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本申请的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
在本发明一实施例中,提供一种材料摩擦系数的不确定度的计算方法。在本实施例中,如图1所示,该计算方法包括以下步骤:
1)获取待测材料的摩擦力F的测量值及法向力N的测量值,以计算摩擦系数μ;
2)分别获取摩擦力F的标准不确定度u F 及法向力N的标准不确定度u N ;
3)通过全微分形式建立摩擦系数μ与摩擦力F和摩擦力F的标准不确定度u F 、以及法向力N和法向力N的标准不确定度u N 的关系,得到摩擦系数的不确定度u μ 。
上述步骤1)中,获取待测材料的摩擦力F的测量值及法向力N的测量值,具体为:由设置的测试系统获取。使用时,由测试系统对待抛磨材料的摩擦系数进行测量,并根据测量结果及测试系统中所涉及的各个不确定度分量计算出最终待抛磨材料摩擦系数的不确定度。
具体的,如图2所示,本实施例中的测试系统包括加载盘系统、上摩擦副、法向加载装置和输出监测系统;加载盘系统用于盛放待测材料,上摩擦副与加载盘系统连接以获取待测材料的摩擦力F,法向加载装置与加载盘系统接触以获取施加在待测材料上的法向力N,输出监测系统用于输出摩擦力F及法向力N的测量值。
在一个可行的实施方式中,加载盘系统包括上抛磨盘1和下抛磨盘2;待抛磨材料3放置在下抛磨盘2上并随下抛磨盘2转动。上摩擦副包括侧向力传感器4和连接钢带5,其中连接钢带5直接与上抛磨盘1接触,防止上抛磨盘1跟随待抛磨材料3发生转动,并通过侧向力传感器4对静止的或有相对运动趋势的上抛磨盘1进行力值测量,并通过输出监测系统进行输出。法向加载装置包括施加法向力的装置(在图中未示出)和法向力传感器6,对转动的下抛磨盘2进行法向力力值测量,并通过输出监测系统进行输出。
上述步骤1)中,根据摩擦系数的定义,摩擦系数μ的数值可通过式(1)计算。
在实验过程中,F、N均为实际测量得到的数值;例如,法向力N 为150 N,摩擦力F为50 N。
根据摩擦系数计算方式,其表达式与实际直接测量值存在显函数关系,可通过全微分形式获得实际测量量的标准不确定度与摩擦系数不确定度的关系,如式(2)所示。
由于摩擦力F,法向力N是通过两套相互独立的系统进行测量,因此,这两个量的相关系数r(F,N)为零。
其中u μ 为摩擦系数标准不确定度,u F 为摩擦力测量标准不确定度,u N 为法向力测量标准不确定度。根据式(1),摩擦系数相对不确定度计算式为:
上述步骤2)中,计算摩擦力F的标准不确定度u F ,包括以下步骤:
u Ftemp为温度引起的测量标准不确定度,具体表达式为:
u Fzerodrift为零点漂移引起的测量标准不确定度,具体表达式为:
u Fcal为校准引起的测量标准不确定度,属于B类计算方法,源自证书记录。例如,经校准,得到力值测量结果不确定度为0.3% (k=2),则u Fcal=0.15 %。
u Fhest为滞后引起的测量标准不确定度,具体表达式为:
其中为相对滞后变化量,通过对侧向力传感器施加载荷,记录传感器在加载和卸载过程中对同一名义力值的输出值,并计算差值,取试验得到的最大相对变化量,为同载荷点进程与回程的示值差,为当量示值。如:对于准确度等级为0.3级标准测力仪,相对滞后可取,则。
u Fst为稳定性引起的测量标准不确定度,具体表达式为:
u Fresolution为分辨率引起的测量标准不确定度,具体表达式为:
u FEMC为EMC电磁兼容引起的测量标准不确定度,具体表达式为:
1)Power voltage variation(DC参考标准:IEC 61000-4-29\IEC61000-4-1;AC参考标准:IEC 61000-2-1\IEC61000-4-1);
2)Short-time power reductions(DC参考标准:IEC 61000-4-29\IEC61000-4-1,AC参考标准:IEC 61000-4-11\IEC61000-6-1\6-2);
3)Bursts(electrical fast transients)(参考标准:IEC 61000-4-4);
4)Surge(参考标准:IEC 61000-4-5);
5)Electrostatic discharge(参考标准:IEC 61000-4-2);
6)Exposure to radiated RF electromagnetic fields(参考标准:IEC 61000-4-3)
7)Exposure to conducted currents generated by RF EM fields(参考标准:IEC 61000-4-6)
U Fcreep为蠕变引起的测量标准不确定度具体表达式为:
其中为传感器满量程加载测试时间,实际测量读数时间为,蠕变相对输出为:。如对于准确度等级为0.3级标准测力仪,传感器满量程=3分钟=180秒蠕变相对输出为=40/700,考虑实际测量读数时间为=40秒,考虑到蠕变读数的随机性,包含因子使用,根据测量过程四个测量用传感器与一个校准用传感器蠕变特性相互独立,则
u Fre为测量系统中示值重复性引起的测量标准不确定度,通过标准传感器进行短时间内的多次标定,获取对应示值,根据bessel公式,计算测量重复性,属于A类计算方法,并可以得到:
U Frp为测量系统中示值复现性引起的测量标准不确定度,通过标准传感器多次不同时间相同控制条件下多次标定,获取对应示值,利用bessel公式计算复现性,属于A类计算方法,并可以得出:
其中,摩擦力F的标准不确定度u F 为:
上述步骤2)中,计算法向力N的标准不确定度u N ,包括以下步骤:
在本实施例中,法向力N的测量以及标准不确定度u N 的计算方法与摩擦力F的测量及标准不确定度u F 中描述一致,例如:u Ftemp与u Ntemp的计算方法一致,其他分量值依此类推,得法向力N的相对不确定度为:
其中,u Norient为法向力测量引起的标准不确定度;
u Ntemp为温度引起的测量标准不确定度;
u Nzerodrift为零点漂移引起的测量标准不确定度;
u Ncal为校准引起的测量标准不确定度;
u Nhest为输出监测系统中滞后引起的测量标准不确定度;
u Nst为稳定性引起的测量标准不确定度;
u Nresolution为示值分辨率引起的测量标准不确定度;
u NEMC为EMC电磁兼容引起的测量标准不确定度;
u Ncreep为蠕变引起的测量标准不确定度;
u Nre为示值重复性引起的测量标准不确定度;
u Nrp为示值复现性引起的测量标准不确定度。
其中,法向力N的标准不确定度u N 为:
上述各实施例中,计算影响摩檫力F的各个因素的不确定度分量,以及计算影响法向力N的各个因素的不确定度分量,具体为:
如图3所示,根据测试系统中每一部分所引入的不确定度分量按照主导量、影响量的划分进行归类。根据各环节中对不确定度的大小的影响机理,对其进行分类;能够用统计的方法得到的不确定度分量采用A类计算方法,不能用统计的方法得到的不确定度分量采用B类计算方法。
其中,A类、B类计算方法参见JCGM 100:2008以及 ISO/IEC Guide 99:2007。
上述步骤3)中,在计算出摩擦力F的不确定度u F 及法向力N的不确定度u N 后,将计算值及摩擦力F、法向力N的测量带入步骤3)中摩擦系数的不确定度计算公式,最终计算得出待测材料摩擦系数的不确定度。
具体的,在摩擦力F与所述法向力N的相关系数为零时,摩擦系数的不确定度u μ 为:
在本实施例中,为保障摩擦系数的不确定度u μ 评估过程中的不确定度分量能够逐一溯源,建立如图4所示的量值溯源体系。
从体系上看,测试系统的各个环节,分别对应国家计量系统相应的计量标准:
上摩擦副所涉及的关键物理量长度、角度及力值分别溯源至国家长度标准、角度标准及力值标准;
加载盘系统所涉及的的关键物理量值溯源至国家力值标准;
法向加载装置所涉及的关键物理量溯源至国家力值标准、质量标准;
输出监测系统所涉及的关键物理量溯源至国家力值标准;
同时所有环节中环境温度的测定与分析将溯源至国家温度标准。
上述国家标准通过量传溯源至国家基准,根据国际计量系统对基本单位的物理常数重新定义,本测试系统所测试的摩擦系数的不确定度最终将溯源至基本物理常数。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (8)
1.一种半导体材料摩擦系数的不确定度的计算方法,该方法应用于测试系统中的待测材料,其特征在于,包括:
获取待测材料的摩擦力F的测量值及法向力N的测量值,以计算摩擦系数μ;
分别获取摩擦力F的标准不确定度u F 及法向力N的标准不确定度u N ;
通过全微分形式建立摩擦系数μ与摩擦力F和摩擦力F的标准不确定度u F 以及法向力N和法向力N的标准不确定度u N 的关系,得到摩擦系数的不确定度u μ ,以将不确定度分量逐一溯源;所述获取待测材料的摩擦力F的测量值及法向力N的测量值,由设置的测试系统获取;
所述测试系统包括加载盘系统、上摩擦副、法向加载装置和输出监测系统;所述加载盘系统用于盛放所述待测材料,所述上摩擦副与所述加载盘系统接触以获取所述待测材料的摩擦力F,所述法向加载装置与所述加载盘系统连接以获取施加在所述待测材料上的法向力N,所述输出监测系统用于输出摩擦力F及法向力N的测量值;
所述加载盘系统包括上抛磨盘和下抛磨盘;待抛磨材料放置在所述下抛磨盘上并随所述下抛磨盘转动;所述上摩擦副包括侧向力传感器和连接钢带,其中所述连接钢带直接与所述上抛磨盘接触,并通过所述侧向力传感器对静止的或有相对运动趋势的所述上抛磨盘进行力值测量,并通过所述输出监测系统进行输出;所述法向加载装置包括施加法向力的装置和法向力传感器,对转动的所述下抛磨盘进行法向力力值测量,并通过所述输出监测系统进行输出。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310138721.8A CN115831296B (zh) | 2023-02-21 | 2023-02-21 | 一种材料摩擦系数的不确定度的计算方法及应用 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310138721.8A CN115831296B (zh) | 2023-02-21 | 2023-02-21 | 一种材料摩擦系数的不确定度的计算方法及应用 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115831296A CN115831296A (zh) | 2023-03-21 |
CN115831296B true CN115831296B (zh) | 2023-05-05 |
Family
ID=85521915
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310138721.8A Active CN115831296B (zh) | 2023-02-21 | 2023-02-21 | 一种材料摩擦系数的不确定度的计算方法及应用 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115831296B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116469793B (zh) * | 2023-06-19 | 2023-08-18 | 北京特思迪半导体设备有限公司 | 晶圆激光改制中改制深度测量不确定度的确定方法及设备 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109187330A (zh) * | 2018-08-09 | 2019-01-11 | 西安理工大学 | 一种变工况条件下摩擦系数测量方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9862070B2 (en) * | 2011-11-16 | 2018-01-09 | Applied Materials, Inc. | Systems and methods for substrate polishing end point detection using improved friction measurement |
CN111597620A (zh) * | 2020-05-25 | 2020-08-28 | 海南热带汽车试验有限公司 | 一种道路载荷测量的不确定度评估方法 |
JP7512170B2 (ja) * | 2020-10-30 | 2024-07-08 | 公益財団法人鉄道総合技術研究所 | 摩擦係数を測定する方法 |
CN114969630B (zh) * | 2022-05-20 | 2024-07-12 | 武汉大学 | 预测法向载荷作用下粗糙分形表面摩擦系数的方法及装置 |
-
2023
- 2023-02-21 CN CN202310138721.8A patent/CN115831296B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109187330A (zh) * | 2018-08-09 | 2019-01-11 | 西安理工大学 | 一种变工况条件下摩擦系数测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115831296A (zh) | 2023-03-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115831296B (zh) | 一种材料摩擦系数的不确定度的计算方法及应用 | |
CN115344814B (zh) | 一种基于柔性气囊的下压力不确定度计算方法 | |
CN109855583B (zh) | 三坐标测量机测量不确定度的确定方法 | |
CN103513195A (zh) | 霍尔传感器测量系统及温度补偿方法 | |
CN110940265A (zh) | 一种应用于静强度试验快速加载采集系统的大应变修正测量方法 | |
CN113176294A (zh) | 一种防护热流计法导热系数测定仪的校准方法 | |
Prakosa et al. | Development of simple method for quality testing of pt100 sensors due to temperature coefficient of resistance measurement | |
Wang et al. | Research and calibration experiment of characteristic parameters of high temperature resistance strain gauges | |
Bartel et al. | Creep and creep recovery response of load cells tested according to US and international evaluation procedures | |
Sekatskii et al. | Two-point calibration of coating thickness gauges: what needs to be considered to improve measurement accuracy | |
Nowicki et al. | Temperature error of Hall-effect and magnetoresistive commercial magnetometers | |
Schiering et al. | Uncertainty evaluation in industrial pressure measurement | |
CN114689950A (zh) | 一种用于静电压检测的温湿度、时间、气压补偿方法 | |
Hucl et al. | Automatic unit for measuring refractive index of air based on Ciddor equation and its verification using direct interferometric measurement method | |
CN108534918A (zh) | 一种传感器自动隔离检测系统 | |
Yang et al. | Uncertainty Assessment of Energy Meter Calibration Devices Based on MID Certification of Energy Meters and Its Application in Evaluating Test Results | |
Zakharov et al. | Peculiarities of measurement uncertainty evaluation at calibrating a ring gauge | |
Castrup et al. | Uncertainty analysis for alternative calibration scenarios | |
Klimov et al. | C1. 3-Integrated Nanophotonics for One One-Step Traceability for Temperature Measurements | |
Wang et al. | The Identification Coding Rules and the Acceptance Criteria of Measuring Instruments | |
Parker et al. | Accurate measurements of electrical conductivity of metals in the range from 2 MS/m to 14 MS/m | |
CN115628845A (zh) | 一种能够对静态电阻应变仪进行自校准的实验方法 | |
Armstrong et al. | Water Vapor Measurements of Ultra-Pure Gases for Critical Manufacturing Processes | |
Leonov et al. | Investigation of the influence of load cell diameter on mass measurement error | |
CN116175396A (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 | ||
REG | Reference to a national code |
Ref country code: HK Ref legal event code: DE Ref document number: 40085663 Country of ref document: HK |