CN101979678A - 一种残余应力振动定位均化方法 - Google Patents
一种残余应力振动定位均化方法 Download PDFInfo
- Publication number
- CN101979678A CN101979678A CN2010105327755A CN201010532775A CN101979678A CN 101979678 A CN101979678 A CN 101979678A CN 2010105327755 A CN2010105327755 A CN 2010105327755A CN 201010532775 A CN201010532775 A CN 201010532775A CN 101979678 A CN101979678 A CN 101979678A
- Authority
- CN
- China
- Prior art keywords
- stress
- exciting
- vibration
- timeliness
- residual stress
- 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
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
- Y02P10/00—Technologies related to metal processing
- Y02P10/20—Recycling
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
一种残余应力振动定位均化方法,它包含如下步骤:(1)模拟时效件加工过程,定量获得工件内残余应力场分布;(2)依据(1)确定需要均化残余应力的具体部位;(3)对工件进行模态分析,优选振型,使模态高应力区与工件高残余应力区对应;(4)依据所选振型确定支撑点、激振点、激振频率;(5)对振动时效进行动力学分析,优选激振力;(6)对振动时效进行疲劳分析,优选激振力、激振时间;(7)利用优选的振动时效工艺参数对工件实施振动;(8)确定下一处需时效的高残余应力部位;(9)重复(3)至(8),直至应力场得到均化。该方法运行成本低,周期短、节能环保,时效效果稳定,造成的疲劳损伤小。
Description
技术领域
本发明涉及一种利用振动定位消除残余应力的方法,尤其是涉及到一种残余应力振动定位均化方法,属于机械制造中的振动时效技术领域。
背景技术
在机械制造领域,金属工件经过铸造、锻造、轧制、切削等机械加工后,在金属内部通常会因为不均匀的塑型变形而产生残余应力,该残余应力的存在会导致金属工件尺寸变化,降低工件尺寸精度,因此,在生产过程中必须增加消除金属工件内残余应力的工艺环节。常用消除残余应力方法包括自然时效、热时效和振动时效。其中,自然时效占用场地大,周期长,不适应大批量生产模式;热时效投资大、耗能高、污染环境且温度难以控制易产生二次应力;振动时效是在工件的共振频率进行振动处理,使工件获得动应力,当动应力与残余应力叠加超过材料屈服极限时,工件产生少量塑性变形,残余应力得到松弛和均化,从而使尺寸精度稳定的一种方法,其具有投资少、效果好、生产周期短、便于携带、对工件形状没有限制、节能环保等优点,近年来呈现逐渐替代前两种时效方法的趋势。
尽管振动时效已经在机械制造领域得到了广泛应用,但传统的振动消除应力技术经过几十年的发展,时效效果不稳定,判断方法不被用户信服等缺点已经暴露出来,具体体现以下几点:
(1)在振动时效实践中,由于受现有残余应力检测技术(盲孔法、磁测法、X衍射法等)发展的限制,人们无法获得时效件三维初始残余应力分布,且检测费用高、耗时长,不宜作为生产现场检测手段使用。而时效件的初始残余应力分布状况未知,则需降低和均化应力的部位无法确定,振动工艺参数(支撑点、激振点、激振频率、激振力、激振时间)的选取就无的放矢,实施振动时效操作存在极大盲目性。现有振动时效在实施前完全不考虑工件初始残余应力分布以及工件各阶振动模态的动应力分布,而仅以若干不同共振频率或谐波频率进行激振,激振力与激振时间的选取也是因人而异,缺乏科学依据,因此,难以有效地全方位或定位降低和均化工件内残余应力水平;
(2)振动消除时效件残余应力工艺参数的选择完全依靠对工件进行重复振动试验,试验周期长、成本高。工件的形状、尺寸、加工工艺的任一改变,确定时效工艺参数的试验就将重新进行,因此,严重阻碍了振动时效技术这种高效节能技术的推广应用;
(3)由于缺乏经济、实时、便捷的定量判据,不同振动工艺参数下时效效果到底如何,对工件产生多大的疲劳损伤我们无法预知,而事实上在时效实施过程中也经常出现由于激振力偏小时效效果不理想,加工变形依旧无法控制;或有时甚至出现工件被振裂的现象。
发明内容
1、目的
为了解决上述问题,本发明提供一种残余应力振动定位均化方法,以解决现有振动时效技术中操作实施盲目、获取时效工艺参数实验周期长、成本高和时效效果无法控制等问题。
2、技术方案
为实现上述发明目的,本发明采用了如下技术方案:
对获得时效工件的加工工艺过程进行模拟,得到时效件残余应力分布,并据此确定需均化残余应力的部位;
对时效件进行模态分析,优选某阶模态,使其模态高应力区与残余应力高应力区对应,依据选取的模态确定振动时效时工件的支撑点、激振点及激振频率;
对振动消除残余应力进行动力学模拟分析,建立激振力与残余应力降低之间的关系,优选激振力大小,使残余应力消除比在设定范围;
对振动消除残余应力进行疲劳寿命模拟分析,建立激振力、激振时间与疲劳寿命之间的关系,优选激振力、激振时间,使振动对工件产生的疲劳损伤小于设定的疲劳损伤容限;
在上述优选的支撑点、激振点、激振频率、激振力、激振时间等振动时效工艺参数下,对工件实施振动时效,至达到设定时间后结束操作。
所述通过对工件加工工艺过程的模拟,定量获得时效工件三维残余应力分布状况,据此确定残余应力大而需要进行均化的部位;
为获得所需模态应力分布,支撑在保持工件平稳前提下,放置于所选模态振型节点位置,激振点在所选振型波峰,激振频率为该阶固有频率;
疲劳损伤容限的具体数值人为确定;
振动后高残余应力在一定区域内得以降低,区域范围取决于激振力的大小;
确定激振时间,使产生疲劳寿命在疲劳损伤容限范围内。
综上所述,本发明一种残余应力振动定位均化方法,该方法具体步骤如下:
步骤一:建立时效件仿真模型,对铸造、锻造、轧制、热处理、切削等加工工艺过程进行模拟,仿真所得时效件具有一定的结构形状特性、应力分布规律以及残余应力大小;
步骤二:依据残余应力的分布特点,确定残余应力高而需降低和均化应力的具体部位;
步骤三:分析不同结构形状时效件的振动模态,优选振型使模态高应力区对应于时效件高残余应力区。对于结构尺寸小,激振器无法安装的时效件,采取与振动台组合振动时效,时效件的模态应力分布与组合振动的振型有关,还与时效件在振动台的放置位置有关,在激振器安装方便的前提下,将时效件刚性连接在振动台高模态应力处,使时效件获得较高动应力;
步骤四:在保持时效件平稳的前提下,所选振动模态的节点即为振动时效时工件的支撑点;所选振动模态的波峰(除强度薄弱处外)即为振动时效时工件的激振点;激振频率为该阶固有频率;
步骤五:在所选振动模态下,对振动时效进行动力学分析,比较不同大小激振力产生的不同动应力场,选择激振力大小使激振产生的动应力与残余应力叠加超过材料的屈服极限,高残余应力得以降低和均化。在实施数值模拟时用激振产生的位移来表征激振力大小,这样可以避免由于阻尼系数选取不准确带来计算动应力的偏差。优选激振位移,使残余应力均化达到设定范围;
步骤六:在所选振动模态下,对时效件进行振动疲劳寿命分析,比较不同激振力、激振时间下,振动对时效件造成的疲劳损伤,合理选择激振力、激振时间,使产生的疲劳损伤不超过设定的疲劳损伤容限;
步骤七:利用优选的振动时效工艺参数对时效件实施振动时效;
步骤八:对步骤七振动时效过程进行模拟,依据振动后残余应力分布确定下一处需振动时效的高残余应力部位;
步骤九:对步骤八确定的高残余应力部位重复步骤三至步骤八,直至时效件残余应力场得到均化。
3、优点及效果
本发明的有益效果是:
(1)本发明在实施振动时效前用仿真模拟方法获知时效工件残余应力的三维分布情况和大小,克服现有残余应力测试技术的不足,一方面可确定需要降低和均化残余应力的具体部位,只需选择对降低和均化残余应力有效的振型进行振动,避免进行不必要的振动时效,造成资源浪费且易对工件产生疲劳损伤,另一方面基于残余应力大小,为合理选取激振力提供理论依据;
(2)本发明不是用传统的扫描或频谱方法寻找工件的共振频率,而是用有限元方法获得工件共振频率,快捷而且准确,同时还可获得各阶频率下工件的振动模态;
(3)本发明用模拟现实技术分析各振动工艺参数对时效效果的影响,避免了对时效件进行重复试验,周期大大缩短,成本也大大降低;
(4)本发明在优选出的振动工艺参数下振动获得的动应力场分布与残余应力场分布一致,即高动应力区域与高残余应力区域重合,从而实现定位降低和均化残余应力;
(5)本发明建立激振力与工件疲劳寿命的关系,使振动对时效件的疲劳寿命影响控制在设定范围内;
(6)本发明的实现使时效工件得到快速有效的处理,而且可以预知振动消除残余应力的效果以及对时效件疲劳寿命的影响。
附图说明
图1a是细长梁型工件前两阶振型,其中1为一阶振型,2为二阶振型;
图1b是细长梁型工件前两阶振型模态应力示意图,其中1’为一阶模态应力分布,2’为二阶模态应力分布;
图2是本发明振动定位降低和均化残余应力方法步骤流程图
图3是利用梁一阶振型振动时效装置示意图
具体实施方式
如图1a所示1为细长梁一阶模态振型,2为细长梁二阶模态振型,如图1b所示1’为一阶模态应力分布,2’为二阶模态应力分布,不同的频率对应不同的振型,不同的振型对应不同的动应力场,重点消除不同部位的残余应力。一阶振型梁中部的弯曲应力最大,重点消除梁中部的残余应力;二阶振型弯曲应力峰值移至梁1/4处,重点消除该部位残余应力,通过选择不同振型消除不同部位的残余应力。
如图2所示,是本发明所述的利用振动定位降低和均化残余应力方法,以轧制梁型件为例介绍该方法。
本发明一种残余应力振动定位均化方法,其具体步骤流程如下:
步骤一:建立时效件仿真模型,对轧制工艺过程进行模拟,得到梁型时效件轧制残余应力数值大小δ残以及残余应力沿梁厚度方向不均匀分布;
步骤二:依据步骤一,确定需整体降低梁内部轧制残余应力;
步骤三:分析梁振动模态,得梁前几阶固有频率f1、f2、f3...fn,梁一阶振型中部弯曲应力最大,利用该振型可消除梁中部轧制残余应力;
步骤四:将橡胶垫、泡沫塑料等弹性物体放置梁件一阶振型的节点处,激振器安装在一阶振型的波峰处,激振频率为一阶固有频率f1。如图3所示装置示意图。
步骤五:在一阶振动模态下,对振动时效进行动力学分析,比较不同大小激振力产生的不同动应力场,选择激振力大小使激振产生的最低动应力δ动与残余应力δ残叠加超过材料的屈服极限,进而找出激振力大小与残余应力降低之间的关系曲线。在实施数值模拟时用激振力产生的位移来表征激振力大小,这样可以避免由于阻尼系数选取不准确带来计算动应力的偏差。优选激振位移,使残余应力均化达到设定范围;
步骤六:在一阶振动模态下,对振动时效进行疲劳寿命分析,比较不同激振力、激振时间振动对时效件疲劳寿命的影响,合理选择激振力、激振时间,使产生的疲劳损伤不超过设定的疲劳损伤容限;
步骤七:利用优选的振动时效工艺参数对时效件实施振动时效;
步骤八:对步骤七振动时效过程进行模拟,得到一阶振型振动时效后残余应力大小及分布情况;
步骤九:利用一阶振型时效件中部残余应力得以降低,对剩下高残余应力部分,利用二阶振型进行振动时效,重复步骤三至步骤七,相应一阶振动模态改为二阶振动模态,最终使整个梁残余应力场得到降低和均化。
Claims (1)
1.一种残余应力振动定位均化方法,其特征在于:该方法具体步骤如下:
步骤一:建立时效件仿真模型,对铸造、锻造、轧制、热处理、切削加工工艺过程进行模拟,仿真所得时效件具有结构形状特性、应力分布规律以及残余应力大小;
步骤二:依据残余应力的分布特点,确定残余应力高而需降低和均化应力的具体部位;
步骤三:分析不同结构形状时效件的振动模态,优选振型使模态高应力区对应于时效件高残余应力区;对于结构尺寸小,激振器无法安装的时效件,采取与振动台组合振动时效,时效件的模态应力分布与组合振动的振型有关,还与时效件在振动台的放置位置有关,在激振器安装方便的前提下,将时效件刚性连接在振动台高模态应力处,使时效件获得较高动应力;
步骤四:在保持时效件平稳的前提下,所选振动模态的节点即为振动时效时工件的支撑点;所选振动模态的波峰即为振动时效时工件的激振点;激振频率为该阶固有频率;
步骤五:在所选振动模态下,对振动时效进行动力学分析,比较不同大小激振力产生的不同动应力场,选择激振力大小使激振产生的动应力与残余应力叠加超过材料的屈服极限,高残余应力得以降低和均化;在实施数值模拟时用激振产生的位移来表征激振力大小,这样可以避免由于阻尼系数选取不准确带来计算动应力的偏差;优选激振位移,使残余应力均化达到设定范围;
步骤六:在所选振动模态下,对时效件进行振动疲劳寿命分析,比较在不同激振力、激振时间下,振动对时效件造成的疲劳损伤,合理选择激振力、激振时间,使产生的疲劳损伤不超过设定的疲劳损伤容限;
步骤七:利用优选的振动时效工艺参数对时效件实施振动时效;
步骤八:对步骤七振动时效过程进行模拟,依据振动后残余应力分布确定下一处需振动时效的高残余应力部位;
步骤九:对步骤八确定的高残余应力部位重复步骤三至步骤八,直至时效件残余应力场得到均化。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010105327755A CN101979678B (zh) | 2010-11-01 | 2010-11-01 | 一种残余应力振动定位均化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010105327755A CN101979678B (zh) | 2010-11-01 | 2010-11-01 | 一种残余应力振动定位均化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101979678A true CN101979678A (zh) | 2011-02-23 |
CN101979678B CN101979678B (zh) | 2013-03-13 |
Family
ID=43600217
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010105327755A Expired - Fee Related CN101979678B (zh) | 2010-11-01 | 2010-11-01 | 一种残余应力振动定位均化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101979678B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102146504A (zh) * | 2011-03-14 | 2011-08-10 | 中国科学院力学研究所 | 局部热时效和振动辅助局部热时效方法 |
CN103589855A (zh) * | 2013-12-03 | 2014-02-19 | 北京航空航天大学 | 一种冷振复合残余应力均化的方法 |
CN103981359A (zh) * | 2014-04-08 | 2014-08-13 | 长治清华机械厂 | 一种频谱谐波时效技术激振力的选配方法 |
CN104259360A (zh) * | 2014-07-29 | 2015-01-07 | 张家港市广大机械锻造有限公司 | Cr9SiMn轴承钢锻前加热工艺 |
CN105486598A (zh) * | 2015-11-26 | 2016-04-13 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种铝锂合金表面残余应力分析方法 |
CN106197915A (zh) * | 2016-06-24 | 2016-12-07 | 北京航空航天大学 | 一种定位振动时效控制时效件疲劳寿命的方法 |
CN109554532A (zh) * | 2018-11-30 | 2019-04-02 | 湖北武铁山桥轨道装备有限公司 | 降低与均化道岔尖轨淬火后残余应力的振动时效方法 |
CN109821722A (zh) * | 2019-03-25 | 2019-05-31 | 贵州瑞吉科技有限公司 | 一种超磁致伸缩振动装置和方法 |
CN110578049A (zh) * | 2019-10-18 | 2019-12-17 | 北京航空航天大学 | 一种基于模态控制的高低频振动-热时效复合应力均化设备 |
CN110586941A (zh) * | 2019-08-26 | 2019-12-20 | 华中科技大学 | 一种金属零件增材制造过程中变形控制系统及方法 |
CN110643912A (zh) * | 2019-10-30 | 2020-01-03 | 北京航空航天大学 | 一种热轧铝合金环件振动热处理工艺方法 |
CN113732617A (zh) * | 2021-11-08 | 2021-12-03 | 滨州学院 | 一种工件损伤修复方法、装置及电子设备 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101225466A (zh) * | 2008-01-24 | 2008-07-23 | 郝俊山 | 一种新型振动时效的方法 |
CN101567021A (zh) * | 2009-04-09 | 2009-10-28 | 中国人民解放军国防科学技术大学 | 用于振动发电的矩形悬臂梁压电振子优化设计有限元方法 |
-
2010
- 2010-11-01 CN CN2010105327755A patent/CN101979678B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101225466A (zh) * | 2008-01-24 | 2008-07-23 | 郝俊山 | 一种新型振动时效的方法 |
CN101567021A (zh) * | 2009-04-09 | 2009-10-28 | 中国人民解放军国防科学技术大学 | 用于振动发电的矩形悬臂梁压电振子优化设计有限元方法 |
Non-Patent Citations (2)
Title |
---|
张洪伟等: "航空整体结构件加工变形校正技术研究", 《兵工学报》, vol. 31, no. 8, 31 August 2010 (2010-08-31), pages 1072 - 1077 * |
赵晓慈等: "基于有限元的铝合金预拉伸板振动时效仿真分析", 《辽宁工程技术大学学报(自然科学版)》, vol. 27, no. 6, 31 December 2008 (2008-12-31), pages 924 - 926 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102146504A (zh) * | 2011-03-14 | 2011-08-10 | 中国科学院力学研究所 | 局部热时效和振动辅助局部热时效方法 |
CN102146504B (zh) * | 2011-03-14 | 2013-03-13 | 中国科学院力学研究所 | 局部热时效和振动辅助局部热时效方法 |
CN103589855A (zh) * | 2013-12-03 | 2014-02-19 | 北京航空航天大学 | 一种冷振复合残余应力均化的方法 |
CN103981359A (zh) * | 2014-04-08 | 2014-08-13 | 长治清华机械厂 | 一种频谱谐波时效技术激振力的选配方法 |
CN103981359B (zh) * | 2014-04-08 | 2016-07-06 | 长治清华机械厂 | 一种频谱谐波时效技术激振力的选配方法 |
CN104259360A (zh) * | 2014-07-29 | 2015-01-07 | 张家港市广大机械锻造有限公司 | Cr9SiMn轴承钢锻前加热工艺 |
CN105486598A (zh) * | 2015-11-26 | 2016-04-13 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种铝锂合金表面残余应力分析方法 |
CN106197915A (zh) * | 2016-06-24 | 2016-12-07 | 北京航空航天大学 | 一种定位振动时效控制时效件疲劳寿命的方法 |
CN109554532A (zh) * | 2018-11-30 | 2019-04-02 | 湖北武铁山桥轨道装备有限公司 | 降低与均化道岔尖轨淬火后残余应力的振动时效方法 |
CN109554532B (zh) * | 2018-11-30 | 2020-10-09 | 湖北武铁山桥轨道装备有限公司 | 降低与均化道岔尖轨淬火后残余应力的振动时效方法 |
CN109821722A (zh) * | 2019-03-25 | 2019-05-31 | 贵州瑞吉科技有限公司 | 一种超磁致伸缩振动装置和方法 |
CN110586941A (zh) * | 2019-08-26 | 2019-12-20 | 华中科技大学 | 一种金属零件增材制造过程中变形控制系统及方法 |
CN110578049A (zh) * | 2019-10-18 | 2019-12-17 | 北京航空航天大学 | 一种基于模态控制的高低频振动-热时效复合应力均化设备 |
CN110643912A (zh) * | 2019-10-30 | 2020-01-03 | 北京航空航天大学 | 一种热轧铝合金环件振动热处理工艺方法 |
CN113732617A (zh) * | 2021-11-08 | 2021-12-03 | 滨州学院 | 一种工件损伤修复方法、装置及电子设备 |
CN113732617B (zh) * | 2021-11-08 | 2022-03-15 | 滨州学院 | 一种工件损伤修复方法、装置及电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN101979678B (zh) | 2013-03-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101979678B (zh) | 一种残余应力振动定位均化方法 | |
CN103602801B (zh) | 一种热振复合残余应力均化的方法 | |
CN103488850B (zh) | 一种热振复合残余应力定位均化的方法 | |
CN105825029B (zh) | 一种用于优化设计高频振动能量放大装置的方法 | |
CN103106307B (zh) | 激光切割机的横梁设计方法 | |
CN103528779B (zh) | 一种电机整体机座振动特性的试验系统及试验方法 | |
CN103589855B (zh) | 一种冷振复合残余应力均化的方法 | |
CN101623690B (zh) | 一种产品可靠性筛选装置及方法 | |
CN108456772B (zh) | 一种确定超声振动时效工艺参数的方法 | |
CN102799729A (zh) | 一种快速消除异质构件残余应力有效方法 | |
CN103955576A (zh) | 一种车床卡盘动平衡配重的方法及装置 | |
CN102564870A (zh) | 裂纹扩展试验方法及装置 | |
CN103255268A (zh) | 一种优化双面激光同时冲击合金厚度的方法 | |
CN203323985U (zh) | 一种模拟涌潮作用下隧道动力响应的试验装置 | |
CN108984951B (zh) | 一种优化机床整机动态特性的方法 | |
CN112464401B (zh) | 一种金属材料焊点的精准建模方法 | |
CN106650155A (zh) | 一种大规模振动工程控制装置优化布置方法 | |
CN106197915A (zh) | 一种定位振动时效控制时效件疲劳寿命的方法 | |
CN106124401B (zh) | 锚杆粘结强度测试方法 | |
CN109182729B (zh) | 一种基于变频调速与数值模拟技术的振动时效系统 | |
CN106980713B (zh) | 一种基于底载分配的大型结构件调频方法 | |
CN105385842B (zh) | 大型薄壁零件制造过程残余应力释放工艺装置及释放方法 | |
CN108229055B (zh) | 一种质刚重合的三向等刚度气浮式隔振系统优化设计方法 | |
CN108776746B (zh) | 一种提高机床动态特性的动刚度优化方法 | |
CN109740200B (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: 20130313 Termination date: 20141101 |
|
EXPY | Termination of patent right or utility model |