CN101947693B - A Tailored Laser Welded Blank Process Optimization System and Method Based on Performance Prediction - Google Patents
A Tailored Laser Welded Blank Process Optimization System and Method Based on Performance Prediction Download PDFInfo
- Publication number
- CN101947693B CN101947693B CN 201010274947 CN201010274947A CN101947693B CN 101947693 B CN101947693 B CN 101947693B CN 201010274947 CN201010274947 CN 201010274947 CN 201010274947 A CN201010274947 A CN 201010274947A CN 101947693 B CN101947693 B CN 101947693B
- Authority
- CN
- China
- Prior art keywords
- module
- prediction
- tailor
- welded
- pls
- 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
- 238000000034 method Methods 0.000 title claims abstract description 89
- 238000005457 optimization Methods 0.000 title claims abstract description 35
- 238000003466 welding Methods 0.000 claims abstract description 53
- 238000007781 pre-processing Methods 0.000 claims abstract description 14
- 238000012805 post-processing Methods 0.000 claims abstract description 9
- 239000002184 metal Substances 0.000 claims description 10
- 229910000831 Steel Inorganic materials 0.000 claims description 7
- 239000000203 mixture Substances 0.000 claims description 7
- 239000010959 steel Substances 0.000 claims description 7
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000004088 simulation Methods 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 2
- 238000012360 testing method Methods 0.000 claims description 2
- 229910000679 solder Inorganic materials 0.000 claims 6
- 239000010953 base metal Substances 0.000 abstract description 12
- 238000005516 engineering process Methods 0.000 abstract description 4
- 238000013386 optimize process Methods 0.000 description 6
- CMRJPMODSSEAPL-FYQPLNBISA-N Trenbolone acetate Chemical compound C1CC2=CC(=O)CCC2=C2[C@@H]1[C@@H]1CC[C@H](OC(=O)C)[C@@]1(C)C=C2 CMRJPMODSSEAPL-FYQPLNBISA-N 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 5
- 230000001186 cumulative effect Effects 0.000 description 5
- 238000007689 inspection Methods 0.000 description 5
- 238000000513 principal component analysis Methods 0.000 description 5
- 230000008676 import Effects 0.000 description 4
- 239000000463 material Substances 0.000 description 4
- 229910001335 Galvanized steel Inorganic materials 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 239000008397 galvanized steel Substances 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000012827 research and development Methods 0.000 description 2
- 229910000922 High-strength low-alloy steel Inorganic materials 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 238000002790 cross-validation Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000010238 partial least squares regression Methods 0.000 description 1
- 230000036314 physical performance Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000002436 steel type Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000002054 transplantation Methods 0.000 description 1
Images
Landscapes
- Laser Beam Processing (AREA)
Abstract
Description
技术领域 technical field
本发明涉及激光焊接技术领域,具体是一种针对拼焊板产品的性能要求,制定优化激光焊接工艺的系统和方法。The invention relates to the technical field of laser welding, in particular to a system and method for formulating and optimizing a laser welding process according to the performance requirements of tailor welded blank products.
背景技术 Background technique
近年来,随着国民经济的增长,经济实用的汽车以其质量轻、耗油少、安全性高而备受关注。由于激光拼焊板的应用减低了车身质量及生产成本,因而被广泛应用于汽车生产企业。但是,长期以来关于各种高性能激光拼焊板的研发工作均局限于实验观察和一般性理论探讨,特别是涉及到各种激光拼焊板结构的工艺设计和工艺优化时,很少能给出定量的工艺模拟及其预测结果,由此大大增加了其研发费用及周期。若能在激光焊接领域引入先进的计算机模拟和预测技术,则可迅速而准确地预测激光拼焊板的力学性能并及时调整和优化激光焊接工艺,对高性能激光拼焊板的研发具有十分重要的工程意义。然而,此类研究的相关成果在国内外均鲜见报道。In recent years, with the growth of the national economy, economical and practical cars have attracted much attention because of their light weight, low fuel consumption and high safety. Since the application of tailor-welded blanks reduces the quality of the body and production costs, it is widely used in automobile manufacturers. However, for a long time, the research and development work on various high-performance tailor-welded blanks has been limited to experimental observation and general theoretical discussion, especially when it comes to the process design and process optimization of various tailor-welded blank structures, it is rarely given Quantitative process simulation and its prediction results are produced, which greatly increases its research and development costs and cycle. If advanced computer simulation and prediction technology can be introduced in the field of laser welding, the mechanical properties of tailor-welded blanks can be predicted quickly and accurately, and the laser welding process can be adjusted and optimized in time, which is very important for the development of high-performance tailor-welded blanks engineering significance. However, the relevant results of such research are rarely reported at home and abroad.
最新检索工作表明,日本Kobe Steel公司在2007年向美国专利局申请了专利“Weld metalexcellent in toughness and SR cracking resistance”(专利号:US07597841),该专利的内容主要为通过对焊接金属成份的调整,以获得优良的焊缝力学性能。此外,新日本制铁公司的JosidaKhirosi在2009年向欧洲专利局申请了专利Device to forecast rupture of part subjected to pointwelding method to this end computer software and machine-readable data carrier(专利号:RU2006013994820050412),该专利的主要内容为发明了一种基于终端计算机处理的预测点焊件破裂的装置。因为点焊时,随着焊点数目的增加,电极头部产生的塑性变形所导致的电极头直径增大等原因会产生脱焊,需要进行人工干预。但关于激光拼焊板力学性能预测和工艺优化方法的专利成果未见报道。The latest retrieval work shows that Kobe Steel Corporation of Japan applied for the patent "Weld metalexcellent in toughness and SR cracking resistance" (Patent No.: US07597841) to the US Patent Office in 2007. The content of this patent is mainly through the adjustment of the welding metal composition, In order to obtain excellent weld mechanical properties. In addition, Josida Khirosi of Nippon Steel Corporation applied to the European Patent Office in 2009 for a patent Device to forecast disruption of part subjected to pointwelding method to this end computer software and machine-readable data carrier (patent number: RU2006013994820050412). The main content is to invent a device for predicting the rupture of spot weldments based on terminal computer processing. Because during spot welding, as the number of welding points increases, the plastic deformation of the electrode head causes the diameter of the electrode head to increase and other reasons will cause desoldering, which requires manual intervention. However, there is no report on the patent results of the mechanical property prediction and process optimization methods of laser tailor welded blanks.
发明内容 Contents of the invention
本发明的目的是以常用激光拼焊钢板为研究对象,具体涉及冷轧深冲板系列St12及其镀锌板;高强度镀锌钢DOGAL800DP/超级拉延钢BUSD、高强度低合金钢板系列B240/390DP等,提供一种针对拼焊板产品的性能要求,制定优化激光焊接工艺的系统和方法。The purpose of the present invention is to take the commonly used laser tailor welded steel plate as the research object, specifically related to the cold-rolled deep-drawing plate series St12 and its galvanized plate; high-strength galvanized steel DOGAL800DP/super-drawn steel BUSD, high-strength low-alloy steel plate series B240 /390DP, etc., provide a system and method for formulating and optimizing the laser welding process according to the performance requirements of tailor-welded blank products.
实现本发明一个目的的技术方案是:一种基于性能预测的激光拼焊板工艺优化系统,该系统由SQL数据库、前处理模块、力学性能预测模块、工艺优化模块、后处理模块组成;The technical solution to achieve one of the objectives of the present invention is: a tailor-welded blank process optimization system based on performance prediction, the system is composed of SQL database, pre-processing module, mechanical performance prediction module, process optimization module, and post-processing module;
所述SQL数据库包含了拼焊板母材的基本信息;The SQL database contains the basic information of the tailor-welded blank base metal;
所述前处理模块,用于从数据库中读取激光焊接过程所需的焊件及母材基本信息和工艺参数,为后续流程提供初始条件;The pre-processing module is used to read the basic information and process parameters of the weldment and base metal required for the laser welding process from the database, and provide initial conditions for the subsequent process;
所述力学性能预测模块,通过对不同工艺状态的模拟得出目标相对误差≤5%的性能预测值以及相应的最优PLS预测模型组;The mechanical performance prediction module obtains a performance prediction value with a target relative error ≤ 5% and a corresponding optimal PLS prediction model group by simulating different process states;
所述工艺优化模块,以系统数据库中已内置的优化性能范围为依据,该范围即为拼焊板屈强比的极小值~极大值,由用户在输入相关焊接母材信息后,调出系统内置的该拼焊板的所有屈强比值,并对其进行选择;对于所选择的某一具体的屈强比值,则可由最优的PLS预测模型组反求得到使拼焊板力学性能最优的工艺方案;The process optimization module is based on the built-in optimization performance range in the system database, which is the minimum value to the maximum value of the yield ratio of tailor-welded blanks, and is adjusted by the user after inputting relevant welding base material information. All the yield strength ratios of the tailor-welded blanks built in the system are obtained and selected; for a specific yield ratio selected, the mechanical properties of the tailor-welded blanks can be obtained by inverse calculation by the optimal PLS prediction model group Optimal process plan;
所述后处理模块,完成结果的显示输出,采用表格或图表方式对工艺优化模块的结果进行输出,同时输出普通文本分析报告。The post-processing module completes the display and output of the results, outputs the results of the process optimization module in the form of tables or graphs, and outputs ordinary text analysis reports at the same time.
具体如下:details as follows:
SQL数据库功能:该数据库包含了拼焊板母材的基本信息,基本信息包括材料牌号、化学成分、母材厚度、焊件尺寸等物理性能参数、焊接工艺参数及其它系统所需数据。SQL database function: This database contains the basic information of the tailor-welded blank base metal, including material grade, chemical composition, base metal thickness, weldment size and other physical performance parameters, welding process parameters and other data required by the system.
前处理模块具体流程为:首先从SQL数据库中读取焊接工艺参数,包括:焊接的母材牌号、成分、母材厚度、焊件尺寸、激光功率、焊接速度、光斑直径、线能量、离焦量、焦距、吸热系数等焊接工艺参数。如读取无误,则转入力学性能预测模块;如有错误,则可返回重新读取工艺参数。The specific process of the pre-processing module is: first read the welding process parameters from the SQL database, including: welding base metal grade, composition, base metal thickness, weldment size, laser power, welding speed, spot diameter, line energy, defocus Welding process parameters such as weight, focal length, heat absorption coefficient, etc. If the reading is correct, it will be transferred to the mechanical performance prediction module; if there is an error, it can return to re-read the process parameters.
力学性能预测模块,利用钢成分、拼焊板母材厚度、焊接工艺参数,建立偏最小二乘预测(Partial Least-Squares Regression,简称“PLS”)公式,并对模型精度的检验与控制,以期得到精度较高的PLS力学性能预测公式。The mechanical performance prediction module uses steel composition, tailor-welded blank base metal thickness, and welding process parameters to establish a partial least squares prediction (Partial Least-Squares Regression, referred to as "PLS") formula, and to test and control the accuracy of the model, in order to The prediction formula of PLS mechanical properties with high precision is obtained.
工艺优化预测模块,由用户在输入相关焊接母材信息后,即可调出系统内置的该拼焊板的所有屈强比值,并对其进行选择。对于所选择的某一具体的屈强比值,则可由最优的PLS预测模型组反求得到使拼焊板力学性能最优的工艺方案。当用户选择项为缺省时,工艺优化模块则将拼焊板屈强比极值范围的中限值默认为优化性能值,并由此得出相应的最优工艺方案因素水平组合。In the process optimization prediction module, after the user enters the relevant welding base metal information, all the yield strength ratios of the tailored welded blank built in the system can be called out and selected. For a selected specific yield strength ratio, the process scheme for optimizing the mechanical properties of tailor-welded blanks can be obtained from the optimal PLS prediction model group. When the user selects the default option, the process optimization module defaults the middle limit of the extreme value range of the yield strength ratio of the tailored welded blank as the optimal performance value, and thus obtains the corresponding optimal process plan factor level combination.
后处理模块,用于显示输出计算结果,采用表格、图表等方式对工艺优化模块的结果进行输出,还包括普通文本分析报告。The post-processing module is used to display and output calculation results, output the results of the process optimization module by means of tables, charts, etc., and also includes ordinary text analysis reports.
实现本发明的另一发明目的的技术方案是:一种基于性能预测的激光拼焊板工艺优化方法,该方法包括以下步骤:The technical solution to achieve another object of the invention is: a method for optimizing the laser tailor welded blank process based on performance prediction, the method includes the following steps:
(1)从数据库中读取激光焊接过程所需的母材信息和工艺参数等,为后续流程提供初始条件的步骤;(1) Read the base metal information and process parameters required for the laser welding process from the database, and provide initial conditions for the subsequent process;
(2)结合待测焊件的基本信息,基本信息包括成分、厚度和工艺参数,建立PLS力学性能预测模型组,并对PLS公式进行精度的检验与控制,得到精度最高的PLS力学性能预测模型组的步骤;(2) Combining the basic information of the weldment to be tested, including composition, thickness and process parameters, a PLS mechanical performance prediction model group is established, and the accuracy of the PLS formula is checked and controlled to obtain the most accurate PLS mechanical performance prediction model group steps;
(3)以系统数据库中已内置的优化性能范围为依据,该范围即为拼焊板屈强比的极小值~极大值,由用户在输入相关焊接母材信息后,调出系统内置的该拼焊板的所有屈强比值,并对其进行选择,并根据所选择的某一具体的屈强比值,求得到使拼焊板力学性能最优的工艺方案。(3) Based on the optimized performance range already built in the system database, this range is the minimum value to maximum value of the yield ratio of tailor-welded blanks. All the yield strength ratios of the tailored welded blank are selected and selected, and the process plan for optimizing the mechanical properties of the tailor welded blank is obtained according to a specific selected yield strength ratio.
(4)输出优化工艺参数。(4) Output optimized process parameters.
所述步骤(1)进一步包括,当焊接工艺参数读取无误,则转入下一步骤;如有错误,则可返回重新读取工艺参数;所述的焊接工艺参数是本领域技术人员所熟知的,包括:焊接板的钢种及其成分、待焊板材厚度规格、激光功率、焊接速度和光斑直径等。The step (1) further includes, when the welding process parameters are read correctly, then go to the next step; if there is an error, then return to re-read the process parameters; the welding process parameters are well known to those skilled in the art It includes: the steel type and its composition of the welded plate, the thickness specification of the plate to be welded, the laser power, the welding speed and the spot diameter, etc.
所述步骤(2)中具体包括:采用交叉有效性确定主成分的个数,以获得PRESSh最小化的预测方程,利用目标相对误差5%的精度检验与控制,获得精度最高的PLS力学性能预测模型组以及逆映射反求得到的最优激光焊接方案。The step (2) specifically includes: using the cross-validity to determine the number of principal components to obtain the prediction equation for the minimization of PRESSh, and using the precision inspection and control of the target relative error of 5% to obtain the highest precision PLS mechanical performance prediction The model group and the optimal laser welding scheme obtained by inverse mapping.
所述步骤(3)中进一步包括,对于所选择的某一具体的屈强比值,由最优的PLS预测模型组反推求得到使拼焊板力学性能最优的工艺方案;当用户选择项为缺省时,工艺优化模块则将拼焊板屈强比极值范围的中限值默认为优化性能值,并由此得出相应的最优工艺方案。The step (3) further includes, for the selected specific yield strength ratio, the optimal PLS prediction model group is reversely obtained to obtain the optimal process plan for the mechanical properties of the tailored welded blank; when the user selection item is By default, the process optimization module defaults the middle limit of the extreme value range of the yield strength ratio of tailor-welded blanks as the optimal performance value, and thus obtains the corresponding optimal process plan.
所述步骤(4)中输出的方式可以是采用表格、曲线或智能报告方式对各预测结果进行显示输出。The way of outputting in the step (4) may be displaying and outputting each prediction result in a form, a curve or an intelligent report.
本发明优点是:Advantage of the present invention is:
1、能根据用户提出的拼焊板目标力学性能,提供相应的最优激光焊接工艺方案。1. According to the target mechanical properties of tailor-welded blanks proposed by users, the corresponding optimal laser welding process scheme can be provided.
2、能根据最优激光焊接工艺方案,进行相应的工艺模拟并预测拼焊板力学性能。利用本发明系统,可使拼焊板的屈服强度、抗拉强度和延伸率的预测精度均达到95%以上,从而保证了拼焊板材质控制的精度。2. According to the optimal laser welding process plan, the corresponding process simulation can be carried out and the mechanical properties of tailor-welded blanks can be predicted. By using the system of the invention, the prediction accuracy of the yield strength, tensile strength and elongation of the tailored welded blank can all reach more than 95%, thereby ensuring the precision of material control of the tailored welded blank.
3、本发明提供的数据库具有大量的拼焊板母材基本信息及激光焊接生产工艺参数,界面友好、输入、输出均与生产过程保持一致,易于操作。系统各模块工作在标准配置的计算机上运行,实现了计算与结果输出的分离,便于程序的调试、升级、维护和移植。3. The database provided by the present invention has a large amount of basic information of tailor-welded blank base metal and laser welding production process parameters, with friendly interface, input and output consistent with the production process, and easy to operate. Each module of the system runs on a computer with standard configuration, which realizes the separation of calculation and result output, and facilitates program debugging, upgrading, maintenance and transplantation.
4、本发明具有优异的普适性,可以推广应用于各种新型、超级强韧化材料的激光焊接应用领域,其优化、预测结果有助于技术人员改进现有生产的工艺,为提高各种新型、超级强韧化材料拼焊板的最终性能提供可靠依据。4. The present invention has excellent universality, and can be popularized and applied to the laser welding application field of various new and super toughened materials. It provides a reliable basis for the final performance of a new type of super-toughened tailor-welded blank.
附图说明 Description of drawings
图1为本发明的系统整体框图;Fig. 1 is a system overall block diagram of the present invention;
图2为力学性能预测模块建立流程框图;Figure 2 is a block diagram of the establishment of the mechanical performance prediction module;
图3为工艺优化模块建立流程框图。Figure 3 is a block diagram of the establishment process of the process optimization module.
具体实施方式 Detailed ways
下面结合附图对本发明作进一步详细说明。The present invention will be described in further detail below in conjunction with the accompanying drawings.
如图1所示,本发明以常用激光拼焊板为研究对象,建立了基于性能预测的激光拼焊板工艺优化的方法和系统。整体模型由数据库、前处理模块、力学性能预测模块、工艺优化模块、后处理模块组成。其中,前处理模块、后处理模块为辅模块,以实现力学性能预测模块的建立和工艺优化模块的优化结果的显示输出。As shown in Fig. 1, the present invention takes commonly used tailor-welded blanks as the research object, and establishes a method and system for process optimization of tailor-welded blanks based on performance prediction. The overall model consists of a database, a pre-processing module, a mechanical performance prediction module, a process optimization module, and a post-processing module. Among them, the pre-processing module and the post-processing module are auxiliary modules to realize the establishment of the mechanical property prediction module and the display output of the optimization results of the process optimization module.
在建立如图2所示的拼焊板力学性能预测模块,为便于控制激光拼焊板力学性能预测精度,必须结合待测焊件的厚度、工艺参数建立PLS力学性能预测模型组,并对PLS式进行精度的检验与控制,得到精度最高的PLS力学性能预测模型组。When establishing the mechanical performance prediction module of tailor welded blanks as shown in Figure 2, in order to facilitate the control of the prediction accuracy of the mechanical performance of tailor welded blanks, it is necessary to establish a PLS mechanical performance prediction model group in combination with the thickness of the weldment to be measured and the process parameters, and PLS The accuracy is tested and controlled by the formula, and the PLS mechanical performance prediction model group with the highest accuracy is obtained.
为使上述拼焊板的PLS力学性能预测模块输出值与实际值的拟合度最好,应使得输出值的误差平方和PRESSh最小化;同时为消除各输入变量的多重相关性影响,需采用交叉有效性确定主成分的个数,以获得PRESSh最小化的预测方程,计算工具采用MATLAB软件。In order to make the best fitting degree between the output value of the PLS mechanical property prediction module of tailor-welded blanks and the actual value, the error square sum PRESS h of the output value should be minimized; at the same time, in order to eliminate the influence of multiple correlations of each input variable, it is necessary to The cross validity is used to determine the number of principal components to obtain the prediction equation for the minimization of PRESS h , and the calculation tool uses MATLAB software.
另外,为使所建模块具备与典型多元线性相关分析相类似的解释输入、输出变量空间的解释能力并使其达到较高的水平,应使得主成分th对输入变量X和输出变量Y的累积解释能力——RdXtt和RdYtt最大化。In addition, in order to make the built module have the ability to explain the input and output variable space similar to the typical multiple linear correlation analysis and make it reach a higher level, the principal component t h should be made to have the input variable X and the output variable Y. Cumulative explanatory power - RdXtt and RdYtt maximized.
对所建模块预测精度的控制,是利用模块的预测值与真实值的目标相对误差加以控制(目标相对误差≤5%),如果符合该精度要求,则输出结果;若不符合时,则返回前处理模块重复图2的建模工作,直至达到目标相对误差为止。The control of the prediction accuracy of the built module is to control the target relative error between the predicted value of the module and the real value (target relative error ≤ 5%). If the accuracy requirement is met, the result will be output; if not, it will return The pre-processing module repeats the modeling work in Figure 2 until the target relative error is reached.
在进行如图3所示的工艺优化模块是以系统数据库中已内置的优化性能范围为依据,该范围即为拼焊板屈强比的极小值~极大值,由用户在输入相关焊接母材信息后,即可调出系统内置的该拼焊板的所有屈强比值,并对其进行选择。对于所选择的某一具体的屈强比值(用户目标屈强比值),则可由最优的PLS预测模型组反求得到相应的优化工艺因素水平组合,意即得到使拼焊板力学性能最优的工艺方案。当用户选择项为缺省时,系统则将拼焊板屈强比极值范围的中限值默认为优化目标性能值,并由此得出相应的最优工艺方案因素水平组合。这种默认值是基于工程应用中的拼焊板均需后续的成形工艺而设,其后续成形工艺要求拼焊板具有较高的成形性能,即屈强比不能过高。若屈强比过高,则势必使成形性能显著降低,继而易产生各种成形缺陷。The process optimization module shown in Figure 3 is based on the optimized performance range built in the system database. This range is the minimum value to maximum value of the yield ratio of tailor-welded blanks. The user inputs the relevant welding After the parent metal information is obtained, all the yield strength ratios of the tailor-welded blank built in the system can be recalled and selected. For a specific selected yield strength ratio (user target yield strength ratio), the optimal PLS prediction model group can be inversely calculated to obtain the corresponding optimized process factor level combination, which means that the optimal mechanical properties of tailor welded blanks can be obtained process plan. When the user selects the default option, the system defaults the middle limit of the extreme value range of the yield strength ratio of tailor-welded blanks as the optimal target performance value, and thus obtains the corresponding optimal process plan factor level combination. This default value is based on the fact that tailor-welded blanks in engineering applications require subsequent forming processes, and the subsequent forming processes require that tailor-welded blanks have high formability, that is, the yield strength ratio should not be too high. If the yield ratio is too high, the formability will be significantly reduced, and various forming defects will easily occur.
下面结合附图,通过4个实施例进一步对本发明的实施过程进行逐步说明。Below in conjunction with the accompanying drawings, the implementation process of the present invention will be further described step by step through four embodiments.
实施例1Example 1
以1.5mm St12板/镀锌板拼焊板为例,焊接工艺采用单面焊双面成形,焊接工艺参数为:功率P=1525~1850W,焊接速度1.6~2.0m/min,光斑直径Φ0.3mm~1mm,吸收率为0.7,焊接用透镜的焦距为127mm。Taking 1.5mm St12 sheet/galvanized sheet tailor-welded blank as an example, the welding process adopts single-sided welding and double-sided forming. The welding process parameters are: power P=1525~1850W, welding speed 1.6~2.0m/min, spot diameter Φ0. 3mm to 1mm, the absorption rate is 0.7, and the focal length of the welding lens is 127mm.
应用MATLAB软件对已有的焊接工艺参数与力学性能的映射关系 j=0,1,…k,i=1,2)进行PLS计算,其算法步骤如下。Using MATLAB software to map the existing welding process parameters and mechanical properties j=0, 1,...k, i=1, 2) to perform PLS calculation, the algorithm steps are as follows.
步骤1:焊件试样(≥9)相关工艺数据,如表1所示。Step 1: Relevant process data of weldment samples (≥9), as shown in Table 1.
表1 焊件试样相关工艺数据Table 1 Process data related to weldment samples
步骤2:考察输入变量之间的多重相关性问题,由表2可见,输入变量之间的多重相关性较明显。Step 2: Investigate the issue of multiple correlations between input variables. It can be seen from Table 2 that the multiple correlations between input variables are obvious.
表2 焊缝区输入变量、输出变量间相关系数矩阵Table 2 Correlation coefficient matrix between input variables and output variables in weld area
步骤3:经交叉有效性原则处理后,得到屈服强度的最优主成分个数,如见表3所示,当h=3时,PRESShmin=0.324679,h=1,2,3;同理得到抗拉强度的最优主成分个数,如见表4所示,当h=3时,PRESShmin=0.376854;同理得到延伸率的最优主成分个数,如见表5所示,当h=3时,PRESShmin=0.489121;故本发明最终取三个主成分建立屈服强度、抗拉强度和延伸率预测模块,分别如下。Step 3: After being processed by the principle of cross validity, the optimal number of principal components of the yield strength is obtained, as shown in Table 3, when h=3, PRESS hmin =0.324679, h=1, 2, 3; the same reason Obtain the optimal number of principal components of tensile strength, as shown in Table 4, when h=3, PRESS hmin =0.376854; similarly, obtain the optimal number of principal components of elongation, as shown in Table 5, When h=3, PRESS hmin =0.489121; therefore, the present invention finally takes three main components to establish yield strength, tensile strength and elongation prediction modules, respectively as follows.
屈服强度y1为Yield strength y1 is
y1=-0.0289x1+36.9053x2-20.9374x3+175.3353y 1 =-0.0289x 1 +36.9053x 2 -20.9374x 3 +175.3353
抗拉强度y2为:The tensile strength y2 is:
y2=-0.0121x1+17.6536x2-38.8001x3+358.1600y 2 =-0.0121x 1 +17.6536x 2 -38.8001x 3 +358.1600
延伸率y3为:The elongation y3 is:
y3=0.0028x1-3.4113x2-6.6175x3+28.6674y 3 =0.0028x 1 -3.4113x 2 -6.6175x 3 +28.6674
以上两式中的x1一激光功率(W),x2-焊接速度(m/min),x3-光斑直径(mm)。In the above two formulas, x 1 - laser power (W), x 2 - welding speed (m/min), x 3 - spot diameter (mm).
步骤4:PLS方程的主成分分析Step 4: Principal Component Analysis of the PLS Equation
对以上的y1、y2、y3分别进行主成分分析后,得到表3、表4、表5。After performing principal component analysis on the above y 1 , y 2 , and y 3 respectively, Table 3, Table 4, and Table 5 are obtained.
表3中,符号RdXt表示主成分th对输入变量X的解释能力;符号RdXtt表示主成分th对输入变量X的累计解释能力;符号RdYt表示主成分th对输出变量Y的解释能力;其中,符号RdYtt表示主成分th对输出变量Y的累计解释能力;PRESSh为Y的预测误差平方和,是交叉有效性方法判断主成分选取个数的依据。In Table 3, the symbol RdXt represents the explanatory ability of the principal component t h to the input variable X; the symbol RdXtt represents the cumulative explanatory ability of the principal component t h to the input variable X; the symbol RdYt represents the explanatory ability of the principal component t h to the output variable Y; Among them, the symbol RdYtt represents the cumulative explanatory ability of the principal component t h to the output variable Y; PRESS h is the sum of squared prediction errors of Y, which is the basis for judging the number of principal components selected by the cross-validation method.
表3表明,应选取三个主成分,此时PRESShmin=0.324679,一共解释了原输入变量系统中99.9384%的变异信息,解释了输出变量系统中97.6633%的变异信息,对输入变量和输出变量的累计解释能力均很高。Table 3 shows that three principal components should be selected. At this time, PRESS hmin = 0.324679, which explains 99.9384% of the variation information in the original input variable system and 97.6633% of the variation information in the output variable system. The cumulative explanatory power is high.
表3 屈服强度的主成分分析数据Table 3 Principal component analysis data of yield strength
表4表明,应选取三个主成分,此时PRESShmin=0.376854,一共解释了原输入变量系统中93.4899%的变异信息,解释了输出变量系统中96.4832%的变异信息,对输入变量和输出变量的累计解释能力均很高。Table 4 shows that three principal components should be selected. At this time, PRESS hmin = 0.376854, which explains 93.4899% of the variation information in the original input variable system and 96.4832% of the variation information in the output variable system. The cumulative explanatory power is high.
表4 抗拉强度的主成分分析数据Table 4 Principal component analysis data of tensile strength
同理,表5表明,应选取三个主成分,此时PRESShmin=0.489121,一共解释了原输入变量系统中94.3812%的变异信息,解释了输出变量系统中97.2396%的变异信息,对输入变量和输出变量的累计解释能力均很高。Similarly, Table 5 shows that three principal components should be selected. At this time, PRESS hmin = 0.489121, which explains 94.3812% of the variation information in the original input variable system and 97.2396% of the variation information in the output variable system. The accumulative explanatory power of both the output variable and the output variable is high.
表5 延伸率的主成分分析数据Table 5 Principal component analysis data of elongation
步骤5:预测精度检验与控制Step 5: Prediction accuracy inspection and control
将前处理模块的工艺数据导入力学性能模块,进行力学性能预测,预测结果如表6所示。由表6可见,屈服强度相对误差为0.1813~4.6043%,抗拉强度相对误差为0.0274~3.7452%,延伸率预测相对误差为0.0658~4.9447%。Import the process data of the pre-processing module into the mechanical performance module to predict the mechanical performance. The prediction results are shown in Table 6. It can be seen from Table 6 that the relative error of yield strength is 0.1813-4.6043%, the relative error of tensile strength is 0.0274-3.7452%, and the relative error of elongation prediction is 0.0658-4.9447%.
表6 抗拉强度、延伸率预测值与实际值Table 6 Predicted and actual values of tensile strength and elongation
步骤6:用户输入期望性能Step 6: User enters desired performance
该拼焊板屈强比的极值范围为:0.48~0.55,用户确定屈强比值为0.50。The extreme value range of the yield ratio of the tailored welded blank is 0.48 to 0.55, and the user determines the yield ratio to be 0.50.
步骤7:输出优化工艺参数Step 7: Output optimized process parameters
x1=1735W,x2-1.6m/min,x3=1mm,Rel=163,Rm=327MPa,A=21%x 1 =1735W, x 2 -1.6m/min, x 3 =1mm, Rel=163, Rm=327MPa, A=21%
实施例2Example 2
以0.8mm St12板/1.5mm镀锌板拼焊板为例,焊接工艺与实例1相同,其建模和优化流程及步骤与实例1相同Taking 0.8mm St12 plate/1.5mm galvanized sheet tailor-welded blank as an example, the welding process is the same as that of Example 1, and its modeling and optimization process and steps are the same as those of Example 1
步骤1:焊件试样(≥9)相关工艺数据,如表7所示。Step 1: Relevant process data of weldment samples (≥9), as shown in Table 7.
表7 焊件试样相关工艺数据Table 7 Process data related to weldment samples
经过实施例2所述的步骤2~步骤4之后,得到了力学性能的PLS预测模型组,如下所示。After steps 2 to 4 described in Example 2, a PLS prediction model group of mechanical properties is obtained, as shown below.
屈服强度y1为:The yield strength y1 is:
y1=-0.0471x1+39.1667x2-21.4414x3+202.0455y 1 =-0.0471x 1 +39.1667x 2 -21.4414x 3 +202.0455
抗拉强度y2为:The tensile strength y2 is:
y2=-0.0472x1+38.3333x2-23.8288x3+359.7018y 2 =-0.0472x 1 +38.3333x 2 -23.8288x 3 +359.7018
延伸率y3为:The elongation y3 is:
y3=0.0103x1-6.6667x2+4.2793x3+16.6625y 3 =0.0103x 1 -6.6667x 2 +4.2793x 3 +16.6625
步骤5:预测精度检验与控制Step 5: Prediction accuracy inspection and control
将前处理模块的工艺数据导入力学性能预测模块,进行力学性能预测,预测结果如表8所示。由表8可见,屈服强度相对误差为0.6948~4.8694%,抗拉强度相对误差为0.1224~4.5038%,延伸率预测相对误差为1.3132~4.4867%。Import the process data of the pre-processing module into the mechanical property prediction module to predict the mechanical properties. The prediction results are shown in Table 8. It can be seen from Table 8 that the relative error of yield strength is 0.6948-4.8694%, the relative error of tensile strength is 0.1224-4.5038%, and the relative error of elongation prediction is 1.3132-4.4867%.
表8 抗拉强度、延伸率预测值与实际值Table 8 Predicted and actual values of tensile strength and elongation
步骤6:用户输入期望性能Step 6: User enters desired performance
该拼焊板屈强比的极值范围为:0.50~0.58,用户确定屈强比值为0.54。The extreme value range of the yield strength ratio of the tailored welded blank is: 0.50~0.58, and the user determines the yield strength ratio to be 0.54.
步骤7:输出优化工艺参数Step 7: Output optimized process parameters
x1=1735W,x2-1.6m/min,x3=1mm,Rel=163,Rm=327MPa,A=21%x 1 =1735W, x 2 -1.6m/min, x 3 =1mm, Rel=163, Rm=327MPa, A=21%
实施例3Example 3
以1.5mm高强度镀锌钢DOGAL800DP/超级拉延钢BUSD拼焊板为例,焊接工艺采用单面焊双面成形,焊接工艺参数为:激光功率为900~1400W,焊接速度为1~2m/min,光斑直径为0.3~1.0mm,吸收率为0.7,焊接用透镜的焦距为127mm。其建模和优化流程及步骤与实施例1相同。Taking 1.5mm high-strength galvanized steel DOGAL800DP/super-drawn steel BUSD tailor-welded blank as an example, the welding process adopts single-sided welding and double-sided forming. The welding process parameters are: laser power 900-1400W, welding speed 1-2m/ min, the spot diameter is 0.3-1.0mm, the absorption rate is 0.7, and the focal length of the welding lens is 127mm. Its modeling and optimization process and steps are the same as those in Example 1.
步骤1:随机调用一组(≥9)相关工艺数据,如表9所示。Step 1: Call a group (≥9) of relevant process data randomly, as shown in Table 9.
表9 焊件试样相关工艺数据Table 9 Process data related to weldment samples
经过实施例1所述的步骤2~步骤4之后,得到了力学性能的PLS预测模型组,如下所示。After steps 2 to 4 described in Example 1, a PLS prediction model group of mechanical properties is obtained, as shown below.
拼焊板屈服强度y1为:Yield strength y1 of tailor welded blank is:
y1=-0.0187x1+9.0000x2-13.8739x3+250.9757y 1 =-0.0187x 1 +9.0000x 2 -13.8739x 3 +250.9757
拼焊板抗拉强度y2为:The tensile strength y2 of tailor welded blank is:
y2=-0.0193x1+10.3333x2-13.8739x3+339.7423y 2 =-0.0193x 1 +10.3333x 2 -13.8739x 3 +339.7423
拼焊板延伸率y3为:Tailored welded blank elongation y 3 is:
y3=0.0047x1-2.6667x2+3.7838x3+29.9036y 3 =0.0047x 1 -2.6667x 2 +3.7838x 3 +29.9036
以上两式中的x1-激光功率(W),x2-焊接速度(m/min),x3-光斑直径(mm)。In the above two formulas, x 1 - laser power (W), x 2 - welding speed (m/min), x 3 - spot diameter (mm).
步骤5:预测精度检验与控制Step 5: Prediction accuracy inspection and control
将前处理模块的工艺数据导入力学性能模块,进行力学性能预测,预测结果如表10所示。由表10可见,抗拉强度相对误差为0.0549~3.3265%,抗拉强度相对误差为0.0900~3.8579%,延伸率预测相对误差为0.5213~3.9148%。Import the process data of the pre-processing module into the mechanical performance module to predict the mechanical performance. The prediction results are shown in Table 10. It can be seen from Table 10 that the relative error of tensile strength is 0.0549-3.3265%, the relative error of tensile strength is 0.0900-3.8579%, and the relative error of elongation prediction is 0.5213-3.9148%.
表10 抗拉强度、延伸率预测值与实际值Table 10 Predicted and actual values of tensile strength and elongation
步骤6:用户输入期望性能Step 6: User enters desired performance
该拼焊板屈强比的极值范围为:0.68~0.75,用户确定屈强比值为缺省。The extreme value range of the yield ratio of the tailored welded blank is: 0.68~0.75, and the default value of the yield ratio is determined by the user.
步骤7:输出优化工艺参数Step 7: Output optimized process parameters
x1=1400W,x2-1.0m/min,x3=0.7mm,Rel=224,Rm=313MPa,A=37%x 1 =1400W, x 2 -1.0m/min, x 3 =0.7mm, Rel=224, Rm=313MPa, A=37%
实施例4Example 4
以1.5mm高强度镀锌钢DOGAL800DP/2.0mm超级拉延钢BUSD拼焊板拼焊板为例,焊接工艺与实例3相同,建模和优化流程及步骤与实例1相同。Taking 1.5mm high-strength galvanized steel DOGAL800DP/2.0mm super-drawn steel BUSD tailor-welded blank as an example, the welding process is the same as Example 3, and the modeling and optimization process and steps are the same as Example 1.
步骤1:焊件试样(≥9)相关工艺数据,如表11所示。Step 1: Relevant process data of weldment samples (≥9), as shown in Table 11.
表11 焊件试样相关工艺数据Table 11 Process data related to weldment samples
经过实施例2所述的步骤2~步骤4之后,得到了力学性能的PLS预测模型组,如下所示。After steps 2 to 4 described in Example 2, a PLS prediction model group of mechanical properties is obtained, as shown below.
抗拉强度y1为:The tensile strength y1 is:
y1=-0.0271x1+13.2568x2-19.4452x3+267.8816y 1 =-0.0271x 1 +13.2568x 2 -19.4452x 3 +267.8816
抗拉强度y2为:The tensile strength y2 is:
y2=-0.0132x1+7.8176x2-22.8733x3+347.6128y 2 =-0.0132x 1 +7.8176x 2 -22.8733x 3 +347.6128
延伸率y3为:The elongation y3 is:
y3=0.0033x1-2.1053x2+5.8724x3+22.4861y 3 =0.0033x 1 -2.1053x 2 +5.8724x 3 +22.4861
步骤5:预测精度检验与控制Step 5: Prediction accuracy inspection and control
将前处理模块的工艺数据导入力学性能模块,进行力学性能预测,预测结果如表12所示。由表12可见,屈服强度相对误差为0.3931~3.9920%,抗拉强度相对误差为0.3931~3.9920%,延伸率预测相对误差为1.6324~3.1255%。Import the process data of the pre-processing module into the mechanical performance module to predict the mechanical performance. The prediction results are shown in Table 12. It can be seen from Table 12 that the relative error of yield strength is 0.3931-3.9920%, the relative error of tensile strength is 0.3931-3.9920%, and the relative error of elongation prediction is 1.6324-3.1255%.
表12 抗拉强度、延伸率预测值与实际值Table 12 Predicted and actual values of tensile strength and elongation
步骤6:用户输入期望性能Step 6: User enters desired performance
该拼焊板屈强比的极值范围为:0.70~0.78,用户确定屈强比值为缺省。The extreme value range of the yield ratio of the tailored welded blank is: 0.70~0.78, and the default value of the yield ratio is determined by the user.
步骤7:输出优化工艺参数Step 7: Output optimized process parameters
x1=1016W,x2-1.3m/min,x3=1mm,Rel=238,Rm=321MPa,A=29%。x 1 =1016W, x 2 -1.3m/min, x 3 =1mm, Rel=238, Rm=321MPa, A=29%.
Claims (6)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010274947 CN101947693B (en) | 2010-09-07 | 2010-09-07 | A Tailored Laser Welded Blank Process Optimization System and Method Based on Performance Prediction |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010274947 CN101947693B (en) | 2010-09-07 | 2010-09-07 | A Tailored Laser Welded Blank Process Optimization System and Method Based on Performance Prediction |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101947693A CN101947693A (en) | 2011-01-19 |
CN101947693B true CN101947693B (en) | 2013-10-23 |
Family
ID=43451388
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010274947 Expired - Fee Related CN101947693B (en) | 2010-09-07 | 2010-09-07 | A Tailored Laser Welded Blank Process Optimization System and Method Based on Performance Prediction |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101947693B (en) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102937564B (en) * | 2012-10-18 | 2015-02-04 | 江苏大学 | Prediction method for corrosion resistance of laser welding plate, and system thereof |
CN104899416B (en) * | 2015-04-23 | 2018-08-10 | 江苏大学 | A kind of prediction technique and its system of the anti-lamellar tearing performance of Steels for High Rise Buildings |
CN106274069A (en) * | 2015-05-21 | 2017-01-04 | 上海超铂信息系统技术有限公司 | A kind of mark processing technique method of testing and system |
CN106312318A (en) * | 2015-06-23 | 2017-01-11 | 姚伟 | Laser welding process for St12 steel plate and welding quality control system for laser welding process |
CN105136589A (en) * | 2015-08-14 | 2015-12-09 | 南京工程学院 | Method for testing forming limit of tailor welded blank during metal plate incremental forming |
CN109581201B (en) * | 2018-11-30 | 2020-11-20 | 北京卫星制造厂有限公司 | Method for judging cold solder joint based on regularization apparent heat absorption coefficient |
CN112149303A (en) * | 2020-09-24 | 2020-12-29 | 扬州浩辰电力设计有限公司 | Carbon fiber wire structure optimization system and method based on tension prediction |
CN112990650B (en) * | 2021-01-12 | 2024-12-17 | 北京钢研新材科技有限公司 | Raw material quality filtering method and system |
CN113203641B (en) * | 2021-03-25 | 2023-01-17 | 武汉钢铁有限公司 | Processing parameter determination method, sample processing method, device and equipment |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6281465B1 (en) * | 1996-10-02 | 2001-08-28 | Siemens Aktiengesellschaft | Device for predictive diagnosis of the current quality of technical production from a technical installation, particularly of the current quality of welding points of a spot welding robot |
CN101119430A (en) * | 2006-07-31 | 2008-02-06 | 富士胶片株式会社 | Template generating apparatus, image layout apparatus, modified template generating apparatus, and programs therefor |
CN101559512A (en) * | 2009-05-21 | 2009-10-21 | 山东大学 | Welding track detection and control method of plate butt weld based on laser ranging |
CN101722371A (en) * | 2009-12-10 | 2010-06-09 | 江苏大学 | Laser welding weld joint/ heat affected zone shape and crystalline grain size prediction method |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5133490B2 (en) * | 2000-09-21 | 2013-01-30 | マサチューセッツ インスティテュート オブ テクノロジー | Spot welding apparatus and method for detecting welding conditions in real time |
EP1642366B1 (en) * | 2003-06-03 | 2019-05-22 | Esab AB | Laser-weld process control system and method |
-
2010
- 2010-09-07 CN CN 201010274947 patent/CN101947693B/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6281465B1 (en) * | 1996-10-02 | 2001-08-28 | Siemens Aktiengesellschaft | Device for predictive diagnosis of the current quality of technical production from a technical installation, particularly of the current quality of welding points of a spot welding robot |
CN101119430A (en) * | 2006-07-31 | 2008-02-06 | 富士胶片株式会社 | Template generating apparatus, image layout apparatus, modified template generating apparatus, and programs therefor |
CN101559512A (en) * | 2009-05-21 | 2009-10-21 | 山东大学 | Welding track detection and control method of plate butt weld based on laser ranging |
CN101722371A (en) * | 2009-12-10 | 2010-06-09 | 江苏大学 | Laser welding weld joint/ heat affected zone shape and crystalline grain size prediction method |
Also Published As
Publication number | Publication date |
---|---|
CN101947693A (en) | 2011-01-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101947693B (en) | A Tailored Laser Welded Blank Process Optimization System and Method Based on Performance Prediction | |
Shehabeldeen et al. | A novel method for predicting tensile strength of friction stir welded AA6061 aluminium alloy joints based on hybrid random vector functional link and henry gas solubility optimization | |
CN101722371B (en) | Laser welding weld joint/ heat affected zone shape and crystalline grain size prediction method, and device for relizing the method | |
CN101354731B (en) | Method for establishing and using difference thick laser assembly solder plate forming limit picture | |
CN111199072B (en) | All-aluminum body riveting system and its implementation method based on online process library self-learning | |
CN110614416B (en) | Method for optimizing welding parameters of dual-phase high-strength steel | |
CN103273207A (en) | Press machine body welding deformation eliminating method based on residual stress quantitative analysis | |
CN102218642B (en) | Manufacturing method of thick wall steel pipe with small radius-thickness ratio used for building structure | |
CN103914593A (en) | Virtual prediction method for mechanical behaviors of laminated composites | |
CN110362861B (en) | Multi-objective optimization method for mold structural parameters in consideration of energy efficiency | |
Equbal et al. | Deformation behavior of micro-alloyed steel by using thermo mechanical simulator and finite element method | |
Zhou et al. | Process optimization of aluminum/steel resistance spot welding based on dynamic resistance analysis | |
CN102937564B (en) | Prediction method for corrosion resistance of laser welding plate, and system thereof | |
CN116206712A (en) | A Design Method for Improving the Elongation of Variable-thickness Metal Thin-walled Structural Parts | |
Li et al. | Improving weldability of press hardened steel through combining stepped current pulse and magnetically assisted resistance spot welding process | |
CN104899416B (en) | A kind of prediction technique and its system of the anti-lamellar tearing performance of Steels for High Rise Buildings | |
Xu et al. | Simultaneous reduction in laser welding energy consumption and strength improvement of aluminum alloy joint by addition of trace carbon nanotubes | |
CN103084811A (en) | High strength and light weight vehicle drive shaft assembly manufacturing method | |
CN105603272A (en) | Aluminum alloy material for automobile alternator and thermal treatment method of aluminum alloy material | |
Ghulman et al. | An application of a genetic algorithm based on Particle Swarm Optimization to a multiple responses problem arising in the Tube Hydroforming Process | |
CN106312318A (en) | Laser welding process for St12 steel plate and welding quality control system for laser welding process | |
CN116890166A (en) | A multi-objective response evaluation method for laser welding systems oriented to carbon emissions and welding quality | |
CN117947365A (en) | A method for optimizing the heat treatment process of CMT arc cladding 9Cr1Mo heat-resistant coating on Q235 steel surface | |
CN103480959A (en) | Spot welding process of thin-walled tube and thin plate | |
JP7578213B1 (en) | Method for predicting welded part properties, method for generating prediction model, and prediction model |
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: 20131023 Termination date: 20150907 |
|
EXPY | Termination of patent right or utility model |