CN107498390A - 一种数控切割机床切割精度的预测方法 - Google Patents
一种数控切割机床切割精度的预测方法 Download PDFInfo
- Publication number
- CN107498390A CN107498390A CN201710993224.0A CN201710993224A CN107498390A CN 107498390 A CN107498390 A CN 107498390A CN 201710993224 A CN201710993224 A CN 201710993224A CN 107498390 A CN107498390 A CN 107498390A
- Authority
- CN
- China
- Prior art keywords
- numerical control
- cutting
- control machine
- machine tool
- workpiece
- 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
Links
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
- B23Q—DETAILS, COMPONENTS, OR ACCESSORIES FOR MACHINE TOOLS, e.g. ARRANGEMENTS FOR COPYING OR CONTROLLING; MACHINE TOOLS IN GENERAL CHARACTERISED BY THE CONSTRUCTION OF PARTICULAR DETAILS OR COMPONENTS; COMBINATIONS OR ASSOCIATIONS OF METAL-WORKING MACHINES, NOT DIRECTED TO A PARTICULAR RESULT
- B23Q17/00—Arrangements for observing, indicating or measuring on machine tools
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
- B23Q—DETAILS, COMPONENTS, OR ACCESSORIES FOR MACHINE TOOLS, e.g. ARRANGEMENTS FOR COPYING OR CONTROLLING; MACHINE TOOLS IN GENERAL CHARACTERISED BY THE CONSTRUCTION OF PARTICULAR DETAILS OR COMPONENTS; COMBINATIONS OR ASSOCIATIONS OF METAL-WORKING MACHINES, NOT DIRECTED TO A PARTICULAR RESULT
- B23Q17/00—Arrangements for observing, indicating or measuring on machine tools
- B23Q17/20—Arrangements for observing, indicating or measuring on machine tools for indicating or measuring workpiece characteristics, e.g. contour, dimension, hardness
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Numerical Control (AREA)
- Automatic Control Of Machine Tools (AREA)
Abstract
本发明提供一种数控切割机床切割精度的预测方法,涉及数控机床技术领域。一种数控切割机床切割精度的预测方法,首先获得数控机床自身的一些参数,判断切割速度与待切割的板材厚度之间的关系,实现在加工前预测加工工件的精度,然后计算数控机床供电电压的幅值、频率、谐波含量以及数控机床的切割速度和待加工工件的长度对切割精度的影响因子,最后通过各影响因子计算出数控机床的切割精度。本发明提供的数控切割机床切割精度的预测方法,从供电质量上对工件的加工精度进行预测,解决了人工检测效率低,误差大的缺点,避免了加工过程中对原材料的大量浪费,节约了成本,提高了检测精度。
Description
技术领域
本发明涉及数控机床技术领域,尤其涉及一种数控切割机床切割精度的预测方法。
背景技术
工厂大量的人工质量检测降低了生产效率,且由于测量工具的限制,对工件的测量也存在很大的误差,影响工件的使用。当前数控机床的制造能力和发展水平体现在数控机床加工出工件的精度上,数控机床的加工度越高,说明数控机床的整体性能越好。但是影响数控机床加工质量的因素多种多样,目前预测数控机床的精度都是从其本身的制造精度、部件装配的匹配度,部件运动过程中的受力、受热而发生的变形,零部件之间的相互摩擦等。在数控机床的总体误差中,几何误差和热误差占到了一半以上,而且在环境温度比较恒定,温度影响较小的情况下,几何误差则为数控机床的主要误差。几何误差是由机床制造、装配及磨损等因素造成的,变化范围不大,相对稳定,易于对其进行改进。对数控机床误差的研究有助于改善加工质量,提高产品精度,对进一步提高机械加工的质量和水平意义重大。同时,针对数控机床加工误差的因素及其形成机理,形成了误差元素的检测方法,及构建了数字化传感网络的检测体系。并在此基础上,基于多体系统理论建立了数控机床误差综合数学模型,而且为了使模型得更好的应用,对模型中的误差元素进行了数值拟合和误差辨识。但以上误差分析和检测方法并没有从数控机床的供电的质量来预测数控机床的误差。
发明内容
针对现有技术的缺陷,本发明提供一种数控切割机床切割精度的预测方法,实现在工件加工过程中实时检测数控机床切割工件精度。
一种数控切割机床切割精度的预测方法,包括以下步骤:
步骤1:检测数控机床供电电压的幅值U、供电电压的频率f和供电电压的谐波含量Hi,以及数控机床的切割速度v和待加工工件的长度s;
步骤2:判断数控机床的切割速度v和待切割铁板的厚度之间是否满足关系式v≤20-ln l2,其中,l为待切割铁板的厚度;如果满足条件,则数控机床为正常工作状态,执行步骤3,对数控机床的切割精度进行检验,否则,数控机床所加工的工件不符合要求,进行错误报警;
步骤3:计算影响数控机床切割精度的各影响因子;
步骤3.1:计算数控机床供电电压的幅值U对切割精度的影响因子AU,计算公式如下所示:
步骤3.2:计算数控机床供电电压的频率f对切割精度的影响因子Af,计算公式如下所示:
其中,p为数控机床电动机的极对数;
步骤3.3:计算数控机床供电电压的谐波含量Hi对切割精度的影响因子计算公式如下所示:
其中,数控机床供电电压的谐波含量Hi的计算公式如下所示:
其中,k′为奇数,k′=3,5,…,∞为谐波次数,k=1为基波,R为整个数控机床电路的等效电阻,C为整个数控机床电路的等效电容,L为整个数控机床电路的等效电感,ω表示供电电压的角频率,j为虚数单位;
步骤3.4:计算数控机床加工工件的过程中待加工工件的长度s对切割精度的影响因子As,计算公式如下所示:
As=(s-s1)×102
其中,s1为数控机床加工完成后的工件实际长度;
步骤3.5:计算数控机床的切割速度v对切割精度的影响因子Av,计算公式如下所示:
步骤4:根据各影响因子,计算得到数控机床的切割精度;
数控机床的切割精度的计算公式如下所示:
其中,A为数控机床的切割精度,A值表示所加工工件的实际值与所加工工件的理论值差值的千分值,A值越小数控机床的切割精度越高;当A<20时,数控机床所加工工件的精度满足要求,否则,数控机床所加工工件的精度不满足要求,该加工工件为不合格产品。
由上述技术方案可知,本发明的有益效果在于:一种数控切割机床切割精度的预测方法,首先获得数控机床自身的一些参数,判断切割速度与待切割的板材厚度之间的关系能否满足精度要求,实现在加工前预测加工工件的精度,避免了加工过程中对原材料的大量浪费。然后将数控机床的精度计算从机械角度转换到供电质量的角度,计算数控机床供电电压的幅值、频率及谐波含量对切割精度的影响因子,同时计算出数控机床的切割速度和待加工工件的长度对切割精度的影响,提高了预测切割精度的准确性,解决了人工检测的效率低,误差大的缺陷,节约了加工成本。
附图说明
图1为本发明实施例提供的一种数控切割机床切割精度的预测方法的流程图。
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。
本实施例以某数控机床为例,使用本发明发明对该数控机床的切割精度进行预测。
一种数控切割机床切割精度的预测方法,如图1所示,包括以下步骤:
步骤1:检测数控机床供电电压的幅值U、供电电压的频率f和供电电压的谐波含量Hi,以及数控机床的切割速度v和待加工工件的长度s;
本实施例中,测得该数控机床供电电压的幅值U=320V,频率f=50.2Hz,该数控机床的切割速度为v=15m/min,待加工工件的长度s=200mm。
步骤2:判断数控机床的切割速度v和待切割铁板的厚度之间是否满足关系式v≤20-ln l2,其中,l为待切割铁板的厚度;如果满足条件,则数控机床为正常工作状态,执行步骤3,对数控机床的切割精度进行检验,否则,数控机床所加工的工件不符合要求,进行错误报警;
本实施例中,待切割板材厚度l=2cm,数控机床的切割速度15<20-ln 22,满足条件,该数控机床为正常工作状态,执行步骤3,对该数控机床的切割精度进行检验。
步骤3:计算影响数控机床切割精度的各影响因子;
步骤3.1:计算数控机床供电电压的幅值U对切割精度的影响因子AU,计算公式如下所示:
本实施例中,计算得到该数控机床供电电压的幅值U对切割精度的影响因子AU=1.89。
步骤3.2:计算数控机床供电电压的频率f对切割精度的影响因子Af,计算公式如下所示:
其中,p为数控机床电动机的极对数;
本实施例中,该数控机床电动机的极对数p=1,计算得到该数控机床供电电压的频率f对切割精度的影响因子Af=4.79。
步骤3.3:计算数控机床供电电压的谐波含量对切割精度的影响因子计算公式如下所示:
其中,数控机床供电电压的谐波含量Hi的计算公式如下所示:
其中,k′为奇数,k′=3,5,…,∞为谐波次数,k=1为基波,R为整个数控机床电路的等效电阻,C为整个数控机床电路的等效电容,L为整个数控机床电路的等效电感,ω表示供电电压的角频率,j为虚数单位;
本实施例中,该数控机床整个电路的等效电阻R=2kΩ,等效电容C=10μF,等效电感L=5H。计算得到该数控机床供电电压的谐波含量Hi=0.04,供电电压的谐波含量对切割精度的影响因子
步骤3.4:计算数控机床加工工件的过程中待加工工件的长度s对切割精度的影响因子As,计算公式如下所示:
As=(s-s1)×102
其中,s1为数控机床加工完成后的工件实际长度;
本实施例中,待加工工件的需求尺寸s=200mm,该数控机床加工完成后的工件实际测的尺寸s1=199.87mm,计算得到待加工工件的长度对切割精度的影响因子As=13。
步骤3.5:计算数控机床的切割速度v对切割精度的影响因子Av,计算公式如下所示:
本实施例中,待切割铁板的厚度l=2cm,数控机床的切割速度v=15m/min,计算得到数控机床的切割速度对切割精度的影响因子Av=13。
步骤4:根据各影响因子,计算得到数控机床的切割精度;
数控机床的切割精度的计算公式如下所示:
其中,A为数控机床的切割精度,A值表示所加工工件的实际值与所加工工件的理论值差值的千分值,A值越小数控机床的切割精度越高;当A<20时,数控机床所加工工件的精度满足要求,否则,数控机床所加工工件的精度不满足要求,该加工工件为不合格产品。
本实施例中,计算得到该数控机床的切割精度A=9.87<20,该数控机床所加工的工件精度满足要求。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明权利要求所限定的范围。
Claims (6)
1.一种数控切割机床切割精度的预测方法,包括以下步骤:
步骤1:检测数控机床供电电压的幅值U、供电电压的频率f和供电电压的谐波含量Hi,以及数控机床的切割速度v和待加工工件的长度s;
步骤2:判断数控机床的切割速度v和待切割铁板的厚度之间是否满足关系式v≤20-lnl2,其中,l为待切割铁板的厚度;如果满足条件,则数控机床为正常工作状态,执行步骤3,否则,数控机床所加工的工件不符合要求,进行错误报警;
步骤3:计算影响数控机床切割精度的各影响因子,包括数控机床供电电压的幅值U对切割精度的影响因子AU、供电电压的频率f对切割精度的影响因子Af、供电电压的谐波含量Hi对切割精度的影响因子数控机床加工工件的过程中待加工工件的长度s对切割精度的影响因子As、切割速度v对切割精度的影响因子Av;
步骤4:根据各影响因子,计算得到数控机床的切割精度;
数控机床的切割精度的计算公式如下所示:
<mrow>
<mi>A</mi>
<mo>=</mo>
<mfrac>
<msqrt>
<mrow>
<msubsup>
<mi>A</mi>
<mi>U</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>A</mi>
<mi>f</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>A</mi>
<mi>s</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>A</mi>
<mi>v</mi>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
<mn>4</mn>
</mfrac>
<mo>+</mo>
<msqrt>
<msub>
<mi>A</mi>
<mi>H</mi>
</msub>
</msqrt>
</mrow>
其中,A为数控机床的切割精度,A值表示所加工工件的实际值与所加工工件的理论值差值的千分值,A值越小数控机床的切割精度越高;当A<20时,数控机床所加工工件的精度满足要求,否则,数控机床所加工工件的精度不满足要求,该加工工件为不合格产品。
2.根据权利要求1所述的一种数控切割机床切割精度的预测方法,其特征在于:步骤3所述数控机床供电电压的幅值U对切割精度的影响因子AU的计算公式如下所示:
<mrow>
<msub>
<mi>A</mi>
<mi>U</mi>
</msub>
<mo>=</mo>
<mn>3</mn>
<mo>&times;</mo>
<mo>|</mo>
<mn>220</mn>
<mo>-</mo>
<mfrac>
<mi>U</mi>
<msqrt>
<mn>2</mn>
</msqrt>
</mfrac>
<mo>|</mo>
<mo>&times;</mo>
<mfrac>
<mn>1</mn>
<mn>10</mn>
</mfrac>
<mo>.</mo>
</mrow>
3.根据权利要求1所述的一种数控切割机床切割精度的预测方法,其特征在于:步骤3所述数控机床供电电压的频率f对切割精度的影响因子Af的计算公式如下所示:
<mrow>
<msub>
<mi>A</mi>
<mi>f</mi>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<mi>lg</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mn>60</mn>
<mi>f</mi>
</mrow>
<mi>p</mi>
</mfrac>
<mo>&times;</mo>
<mn>1000</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>6</mn>
<mo>&rsqb;</mo>
<mo>&times;</mo>
<mn>10</mn>
</mrow>
其中,p为数控机床电动机的极对数。
4.根据权利要求1所述的一种数控切割机床切割精度的预测方法,其特征在于:步骤3所述数控机床供电电压的谐波含量Hi对切割精度的影响因子的计算公式如下所示:
<mrow>
<msub>
<mi>A</mi>
<msub>
<mi>H</mi>
<mi>i</mi>
</msub>
</msub>
<mo>=</mo>
<mi>a</mi>
<mi>r</mi>
<mi>c</mi>
<mi>t</mi>
<mi>a</mi>
<mi>n</mi>
<mfrac>
<mrow>
<mo>|</mo>
<msub>
<mi>H</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<mo>&times;</mo>
<mn>100</mn>
</mrow>
<mn>5</mn>
</mfrac>
</mrow>
其中,数控机床供电电压的谐波含量Hi的计算公式如下所示:
<mrow>
<msub>
<mi>H</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&Sigma;</mi>
<mfrac>
<mi>U</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mi>R</mi>
</mfrac>
<mo>+</mo>
<msup>
<mi>jk</mi>
<mo>&prime;</mo>
</msup>
<mi>&omega;</mi>
<mi>C</mi>
<mo>)</mo>
<mo>&times;</mo>
<msup>
<mi>Ljk</mi>
<mo>&prime;</mo>
</msup>
<mi>&omega;</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</mfrac>
</mrow>
<mfrac>
<mi>U</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mi>R</mi>
</mfrac>
<mo>+</mo>
<mi>j</mi>
<mi>k</mi>
<mi>&omega;</mi>
<mi>C</mi>
<mo>)</mo>
<mo>&times;</mo>
<mi>L</mi>
<mi>j</mi>
<mi>k</mi>
<mi>&omega;</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</mfrac>
</mfrac>
</mrow>
其中,k′为奇数,k′=3,5,…,∞为谐波次数,k=1为基波,R为整个数控机床电路的等效电阻,C为整个数控机床电路的等效电容,L为整个数控机床电路的等效电感,ω表示供电电压的角频率,j为虚数单位。
5.根据权利要求1所述的一种数控切割机床切割精度的预测方法,其特征在于:步骤3所述待加工工件的长度s对切割精度的影响因子As的计算公式如下所示:
As=(s-s1)×102
其中,s1为数控机床加工完成后的工件实际测的尺寸。
6.根据权利要求1所述的一种数控切割机床切割精度的预测方法,其特征在于:步骤3所述数控机床的切割速度v对切割精度的影响因子Av的计算公式如下所示:
<mrow>
<msub>
<mi>A</mi>
<mi>v</mi>
</msub>
<mo>=</mo>
<mroot>
<mrow>
<msup>
<mi>l</mi>
<mn>2</mn>
</msup>
<mi>v</mi>
</mrow>
<mn>3</mn>
</mroot>
<mo>+</mo>
<mi>lg</mi>
<mi> </mi>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
<mo>.</mo>
</mrow>
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710993224.0A CN107498390B (zh) | 2017-10-23 | 2017-10-23 | 一种数控切割机床切割精度的预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710993224.0A CN107498390B (zh) | 2017-10-23 | 2017-10-23 | 一种数控切割机床切割精度的预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107498390A true CN107498390A (zh) | 2017-12-22 |
CN107498390B CN107498390B (zh) | 2019-06-18 |
Family
ID=60701932
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710993224.0A Expired - Fee Related CN107498390B (zh) | 2017-10-23 | 2017-10-23 | 一种数控切割机床切割精度的预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107498390B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108663225A (zh) * | 2018-05-31 | 2018-10-16 | 沪东中华造船(集团)有限公司 | 一种检验数控切割机切割精度的方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1117577A (zh) * | 1994-01-28 | 1996-02-28 | 三菱电机株式会社 | 绝对位置检测装置及其误差补偿方法 |
JPH08155700A (ja) * | 1994-12-05 | 1996-06-18 | Cosmo Net:Kk | プレス装置の工具類の使用回数管理装置 |
US20100154161A1 (en) * | 2006-01-20 | 2010-06-24 | Kumio Takahashi | Dust collector |
CN102501137A (zh) * | 2011-11-03 | 2012-06-20 | 西安交通大学 | 一种主轴径向回转精度在线监测装置 |
CN104482850A (zh) * | 2014-12-22 | 2015-04-01 | 深圳市五株科技股份有限公司 | Pcb钻孔精度检测方法 |
-
2017
- 2017-10-23 CN CN201710993224.0A patent/CN107498390B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1117577A (zh) * | 1994-01-28 | 1996-02-28 | 三菱电机株式会社 | 绝对位置检测装置及其误差补偿方法 |
JPH08155700A (ja) * | 1994-12-05 | 1996-06-18 | Cosmo Net:Kk | プレス装置の工具類の使用回数管理装置 |
US20100154161A1 (en) * | 2006-01-20 | 2010-06-24 | Kumio Takahashi | Dust collector |
CN102501137A (zh) * | 2011-11-03 | 2012-06-20 | 西安交通大学 | 一种主轴径向回转精度在线监测装置 |
CN104482850A (zh) * | 2014-12-22 | 2015-04-01 | 深圳市五株科技股份有限公司 | Pcb钻孔精度检测方法 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108663225A (zh) * | 2018-05-31 | 2018-10-16 | 沪东中华造船(集团)有限公司 | 一种检验数控切割机切割精度的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107498390B (zh) | 2019-06-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102854841B (zh) | 一种曲面零件的形位误差原位补偿加工方法 | |
Wu et al. | Thermal error compensation method for machine center | |
CN106808320B (zh) | 一种考虑刀具齿间角、螺旋角及偏心的铣削力预测方法 | |
Wan et al. | A novel cutting force modelling method for cylindrical end mill | |
Zhang et al. | Prediction of cutting forces in milling of circular corner profiles | |
CN103676782B (zh) | 数控铣床加工过程中能量效率在线检测方法 | |
Du et al. | Peripheral milling force induced error compensation using analytical force model and APDL deformation calculation | |
Camposeco-Negrete et al. | Sustainable machining as a mean of reducing the environmental impacts related to the energy consumption of the machine tool: a case study of AISI 1045 steel machining | |
CN105269402A (zh) | 一种基于铣削加工的钛合金材料表面粗糙度预测方法 | |
Zhang et al. | A new method for cutting force prediction in peripheral milling of complex curved surface | |
CN108490873A (zh) | 一种基于最近邻空间点的数控机床主轴热伸长预测方法 | |
CN107498390A (zh) | 一种数控切割机床切割精度的预测方法 | |
Blaser et al. | LONG-TERM THERMAL COMPENSATION OF 5-AXIS MACHINE TOOLS DUE TO THERMAL ADAPTIVE LEARNING CONTROL. | |
CN101780488B (zh) | 内螺纹冷挤压加工质量在线预测方法 | |
Wang et al. | Energy consumption model for milling processes considering auxiliary load loss and its applications | |
CN217096927U (zh) | 一种数控机床刀补检测系统 | |
Fan et al. | Error prediction and clustering compensation on shaft machining | |
Horejs et al. | Real-time compensation of machine tool thermal error including cutting process | |
Zhu et al. | A model-based monitoring and fault diagnosis methodology for free-form surface machining process | |
Bolar et al. | Experimental investigation and optimization of wall deflection and material removal rate in milling thin-wall parts | |
Zaides et al. | Automated complex for stabilized straightening of low-stiff cylindrical parts | |
Pavel et al. | Multi-constraint optimization for grinding nickel-based alloys | |
CN108733878B (zh) | 一种高猛刚材料高速加工的切削力预测模型的建立方法 | |
Hong et al. | Research on straightening process model based on iteration and self-learning | |
CN112434447B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190618 Termination date: 20191023 |
|
CF01 | Termination of patent right due to non-payment of annual fee |