CN103049623B - 一种激光焊接热源模型的建立方法 - Google Patents
一种激光焊接热源模型的建立方法 Download PDFInfo
- Publication number
- CN103049623B CN103049623B CN201310020198.5A CN201310020198A CN103049623B CN 103049623 B CN103049623 B CN 103049623B CN 201310020198 A CN201310020198 A CN 201310020198A CN 103049623 B CN103049623 B CN 103049623B
- Authority
- CN
- China
- Prior art keywords
- heat source
- laser
- model
- ellipsoid
- heat
- 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
Landscapes
- Laser Beam Processing (AREA)
Abstract
一种激光焊接热源模型的建立方法,本发明涉及焊接热源模型的建立方法。本发明是要解决大型复杂构件的激光焊接模拟中存在的激光热源实现难度大,计算效率低的问题,而提供的一种激光焊接热源模型的建立方法。第一步:建立三维有限元网格模型;第二步:建立热源表面的高斯热源模型;第三步:将焊接能量沿着激光焊接熔池深度方向进行拓展;第四步:基于有限元计算软件求解控制方程,进行热-机耦合计算,即完成了激光焊接热源模型的建立。本发明应用于焊接领域。
Description
技术领域
本发明涉及焊接热源模型的建立方法。
背景技术
激光焊接作为一种高效率、高质量、高精度和小变形的焊接方法得到日益广泛的应用和不断的开发,已经在机械制造、航空航天、汽车工业以及微电子行业等领域得到飞速的发展。尤其是在航空航天领域,激光焊接因其能够实现多种类型材料的连接,具有其他熔焊工业无法比拟的优越性等优点在航空工业的应用中做出了开创性贡献。采用激光焊接的焊接结构从降低飞行器制造成本、缩短生产周期、减轻重量、提高构件的机械性能,提高运载能力、增强稳定性方面,都有良好的效果。目前,激光焊接技术已应用于波音、空中客车等主要航空企业的飞机制造中。空中客车公司在A318飞机机身壁板的桁条上用激光焊接代替铆接工艺。该工艺不仅可以节约成本,而且可以提高抗腐蚀能力,减重并减少裂纹增长。德国飞机蒙皮与筋条的连接采用激光焊接技术,用于A330/340飞机壁板的制造,其结果与铆接相比,机身的重量减轻18%左右,制造成本降低了近25%。航空领域应用一般结构件比较复杂,且材料造价比较高,制造工艺不成熟,如果采用物理实验或者是依靠经验的方法需要做很多次试验,试验周期很长,这将会耗费大量的不必要成本,严重消耗大量物力人力资源,严重的影响到构件的投入生产使用的时间。所以焊接模拟和仿真具有十分重要的意义。采用数值模拟的方法来代替前期不必要的物理实验,可以节约大量的成本和物力人力资源,缩短研究周期。将生产车间的反复性试验交由计算机处理,可减少实际试验的反复性和盲目性,并且焊接在模拟过程中调整方便,便于控制。
激光焊接模拟首要解决的问题是激光焊接热源模型的问题。目前国内外文献中所提到的比较成熟的激光热源模型多为组合热源模型,解析公式较为复杂,商用有限元计算软件中只有高斯面热源和双椭球热源模型,没有与之对应的激光热源模型。所以要实现激光焊接模拟,必须人工编写含有激光热源相关代码的子程序,然后通过有限元软件调用该程序进行计算。这样便带来了诸多弊端:其一是热源子程序的编写复杂多变,通用性较差,热源移动路径难于实现,对编程者水平的要求较高;其二有限元软件调用子程序会耗用许多额外时间,大大降低了计算效率,若有限元模型网格数量较大,无疑会进一步加大计算量,甚至会达到计算机无法承受的程度;此外激光热源的校核比较复杂,重复性调试操作量较大。利用本发明所提供的热源模型无需人工编程,利用有限元计算软件自带热源模型模块便可实现,降低了热源实现难度,提高了计算效率,可完成大型复杂构件的激光焊接过程。另外该热源构建形式简单,热源形状调节方便,可控性强,与实验结果相似程度高。
公开日期2004年4月30日的不锈钢焊接温度场的三维数值模拟,但该文献仅仅公开了利用网格划分的高斯热源模型,并未公开有限元计算软件自带激光热源模型模块。
发明内容
本发明是要解决大型复杂构件的激光焊接模拟中存在的激光热源实现难度大,计算效率低的问题,而提供的一种激光焊接热源模型的建立方法。
本发明的一种激光焊接热源模型的建立方法按以下步骤实现:
第一步:建立三维有限元网格模型:
在三维造型软件中建立工件的几何模型,利用网格划分软件或者有限元计算软件对工件进行网格划分;
第二步:建立热源表面的高斯热源模型:
式中q(x,y,t)为面热源在工件表面上的热流密度分布,R为面热源有效加热半径,Q为焊接功率,参数η1为面热源能量分配系数,t为时间;
第三步:将焊接能量沿着激光焊接熔池深度方向进行拓展:
复制双椭球热源模型,将焊接能量分布到激光焊接熔池深度方向,在表面的高斯热源的下方叠加一个或者多个热源尺寸较小的双椭球热源可实现深宽比较大的激光熔池形貌;其中,所述双椭球热源模型将焊接熔池的前半部分作为一个1/4椭球,后半部分作为另一个1/4椭球,面热源下方第i个双椭球热源的公式为:
前半部分椭球内热源分布函数:
后半部分椭球内热源分布函数:
式中q(x,y,z,t)为双椭球热源在工件内部的热流密度分布,ff和fr分别为总的输入功率在熔池前、后两部分的分配指数,通常ff+fr=2,ηi为第i个热源的能量分配系数,参数ai1,ai2,bi,ci为激光焊接热源内部双椭球热源参数,其中(n-1)*i≤ci≤n*i;
对于整个激光热源来说
Q=ηUI
h为激光焊接热源总的有效作用深度,η为激光热源功率有效系数;
第四步:基于有限元计算软件求解控制方程,进行热-机耦合计算:
将所建立的激光热源模型加载到网格模型上,设置热源边界条件及控制热源位置,提交运算并导出计算结果,建立三维具有内热源和瞬态温度场的固体导热微分方程:
式中,T为物体的瞬态温度,t为过程进行的时间,k为材料的导热系数,ρ为材料的密度,cp为材料的定压比热,qv为材料的内热源强度,x,y,z分别为坐标系中的三向坐标,建立热-弹塑性控制方程:
dσ=[Dep]dε=([De]-[Dp])dε(5)
式中,[De]为弹性矩阵,[Dp]为塑性矩阵,[Dep]为弹塑性矩阵,即完成了激光焊接热源模型的建立。
有益效果:本发明提供的一种激光焊接热源模型的建立方法,针对激光焊接特别是深熔型激光焊接过程,建立了有限元网格模型,添加合理的初始条件和边界条件,利用有限元计算软件进行求解,计算结果与实验结果吻合良好。通过试验证明,整个计算过程中加载本发明的热源模型的计算总时间为764秒,而加载根据国内外文献常用热源模型所编制的子程序的计算总时间为1532秒,所以本发明所提供的热源模型运用到激光焊接的模拟过程较常规的计算方法可提高效率50%以上,而且当计算的构件尺寸越大,焊缝总长度越长时,本发明所提供的热源模型的高效性效果越明显。
附图说明
图1为高斯热源模型;
图2为双椭球热源模型;
图3为激光焊接件网格划分示意图;
图4为激光焊接宏观温度场示意图;
图5为激光焊缝横截面剖面图;
图6为根据目前国内外通用激光热源模型编写相关子程序,并加载到网格模型上进行计算得到的焊缝熔池剖面图;
图7为利用本发明所提供的热源模型进行计算所得到的焊缝熔池剖面图。
具体实施方式
具体实施方式一:本实施方式的一种激光焊接热源模型的建立方法按以下步骤实现:
第一步:建立三维有限元网格模型:
在三维造型软件中建立工件的几何模型,利用网格划分软件或者有限元计算软件对工件进行网格划分;
第二步:建立热源表面的高斯热源模型:
式中q(x,y,t)为面热源在工件表面上的热流密度分布,R为面热源有效加热半径,Q为焊接功率,参数η1为面热源能量分配系数,t为时间;
第三步:将焊接能量沿着激光焊接熔池深度方向进行拓展:
复制双椭球热源模型,将焊接能量分布到激光焊接熔池深度方向,在表面的高斯热源的下方叠加一个或者多个热源尺寸较小的双椭球热源可实现深宽比较大的激光熔池形貌;其中,所述双椭球热源模型将焊接熔池的前半部分作为一个1/4椭球,后半部分作为另一个1/4椭球,面热源下方第i个双椭球热源的公式为:
前半部分椭球内热源分布函数:
后半部分椭球内热源分布函数:
式中q(x,y,z,t)为双椭球热源在工件内部的热流密度分布,ff和fr分别为总的输入功率在熔池前、后两部分的分配指数,通常ff+fr=2,ηi为第i个热源的能量分配系数,参数ai1,ai2,bi,ci为激光焊接热源内部双椭球热源参数,其中(n-1)*i≤ci≤n*i;
对于整个激光热源来说
Q=ηUI
h为激光焊接热源总的有效作用深度,η为激光热源功率有效系数;
第四步:基于有限元计算软件求解控制方程,进行热-机耦合计算:
将所建立的激光热源模型加载到网格模型上,设置热源边界条件及控制热源位置,提交运算并导出计算结果,建立三维具有内热源和瞬态温度场的固体导热微分方程:
式中,T为物体的瞬态温度,t为过程进行的时间,k为材料的导热系数,ρ为材料的密度,cp为材料的定压比热,qv为材料的内热源强度,x,y,z分别为坐标系中的三向坐标,建立热-弹塑性控制方程:
dσ=[Dep]dε=([De]-[Dp])dε(5)
式中,[De]为弹性矩阵,[Dp]为塑性矩阵,[Dep]为弹塑性矩阵,即完成了激光焊接热源模型的建立。
本实施方式中激光焊接熔池的横截面程一个钉子形状,熔池表面有一个“钉子头”,该特征是由于激光热源在焊接熔池表面形成一个半球状的光致等离子球体,该球体对焊缝熔池表面进行加热,导致熔池表面熔化尺寸较熔池内部要大。因而激光热源模型表面的高斯热源的热源参数较大,从而形成焊接熔池的“钉子头”;
激光焊接熔池在“钉子头”下方的熔化区域尺寸较小,熔深随着焊接功率的增大而增大。
本实施方式的商用有限元计算软件为Marc,购买自MSC。
有益效果:本实施方式提供的一种激光焊接热源模型的建立方法,针对激光焊接特别是深熔型激光焊接过程,建立了有限元网格模型,添加合理的初始条件和边界条件,利用有限元计算软件进行求解,计算结果与实验结果吻合良好。通过试验证明,整个计算过程中加载本发明的热源模型的计算总时间为764秒,而加载根据国内外文献常用热源模型所编制的子程序的计算总时间为1532秒,所以本实施方式所提供的热源模型运用到激光焊接的模拟过程较常规的计算方法可提高效率50%以上,而且当计算的构件尺寸越大,焊缝总长度越长时,本实施方式所提供的热源模型的高效性效果越明显。
通过以下试验验证本发明有益效果:
下面以某材料的平板与肋板组成的T型接头,在接头的背面即平板一侧的激光束焊接过程为例进行模拟计算并分析:
T型接头肋板厚度都为1.4mm,底板厚度为0.8mm,模型底板的大小为5×16mm,肋板的高度为3mm,长16mm,根据实际情况,在T型背面平板一侧施加激光束进行焊接;
第一步:基于MSC.Marc软件的网格画分模块对几何模型划分网格,网格划分统一采用六面体单元,靠近焊缝的区域采用单元较小的加密型网格,离焊缝越远,网格尺寸越大,共73600个单元,81116个节点;T型接头的材料为GH3128镍基合金;
焊接所用单束激光的有效功率为1.2kW,激光束移动速度为2m/min;
工件底面与工作台之间存在传导换热,其他表面的热边界条件为热对流和热辐射;
第二步:加载本发明所提供的热源模型,在平板表面添加高斯面热源,平板和肋板上添加四次双椭球热源,每个热源的作用深度ci为4mm,合理分配每个热源的焊接热输入,最后提交运算,提取结果;
第三步:为对比本发明的热源模型的计算精度和计算效率,我们将上述运算过程采用国内外常用的激光热源模型进行加载并再次运算:
为了验证模型计算结果的准确性,将计算结果与实际焊缝进行比较,通过比较可以看出,本发明所提供的热源模型同样得到了与实际激光焊缝相似度很高的熔池形貌,但是本发明所提供的热源模型构建形式简洁,仅通过高斯面热源和双椭球热源组合而成,与其他激光热源模型相比简单明了,此外该模型加载方式便捷,利用有限元计算软件内部的模块便可实现加载,避免了复杂热源的人工编程,也节省了软件调用子程序文件的时间。
上述两种情况的计算过程均在同一计算机上进行计算,整个计算过程中加载本发明的热源模型的计算总时间为764秒,而加载根据国内外文献常用热源模型所编制的子程序的计算总时间为1532秒,所以本发明所提供的热源模型运用到激光焊接的模拟过程较常规的计算方法可提高效率50%以上,而且当计算的构件尺寸越大,焊缝总长度越长时,本发明所提供的热源模型的高效性效果越明显。
Claims (1)
1.一种激光焊接热源模型的建立方法,其特征在于激光焊接热源模型的建立方法按以下步骤实现:
第一步:建立三维有限元网格模型:
在三维造型软件中建立工件的几何模型,利用网格划分软件或者有限元计算软件对工件进行网格划分;
第二步:建立热源表面的高斯热源模型:
式中q(x,y,t)为面热源在工件表面上的热流密度分布,R为面热源有效加热半径,Q为焊接功率,参数η1为面热源能量分配系数,t为时间;
第三步:将焊接能量沿着激光焊接熔池深度方向进行拓展:
复制双椭球热源模型,将焊接能量分布到激光焊接熔池深度方向,在表面的高斯热源的下方叠加一个或者多个热源尺寸较小的双椭球热源可实现深宽比较大的激光熔池形貌;其中,所述双椭球热源模型将焊接熔池的前半部分作为一个1/4椭球,后半部分作为另一个1/4椭球,面热源下方第i个双椭球热源的公式为:
前半部分椭球内热源分布函数:
后半部分椭球内热源分布函数:
式中q(x,y,z,t)为双椭球热源在工件内部的热流密度分布,ff和fr分别为总的输入功率在熔池前、后两部分的分配指数,ff+fr=2,ηi为第i个热源的能量分配系数,参数ai1,ai2,bi,ci为激光焊接热源内部双椭球热源参数,其中(n-1)*i≤ci≤n*i;
对于整个激光热源来说
Q=ηUI
h为激光焊接热源总的有效作用深度,η为激光热源功率有效系数;
第四步:基于有限元计算软件求解控制方程,进行热-机耦合计算:
将所建立的激光热源模型加载到网格模型上,设置热源边界条件及控制热源位置,提交运算并导出计算结果,建立三维具有内热源和瞬态温度场的固体导热微分方程:
式中,T为物体的瞬态温度,t为过程进行的时间,k为材料的导热系数,ρ为材料的密度,cp为材料的定压比热,qv为材料的内热源强度,x,y,z分别为坐标系中的三向坐标,建立热-弹塑性控制方程:
dσ=[Dep]dε=([De]-[Dp])dε(5)
式中,dε为应变分量,dσ为应力增量,[De]为弹性矩阵,[Dp]为塑性矩阵,[Dep]为弹塑性矩阵,即完成了激光焊接热源模型的建立。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310020198.5A CN103049623B (zh) | 2013-01-18 | 2013-01-18 | 一种激光焊接热源模型的建立方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310020198.5A CN103049623B (zh) | 2013-01-18 | 2013-01-18 | 一种激光焊接热源模型的建立方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103049623A CN103049623A (zh) | 2013-04-17 |
CN103049623B true CN103049623B (zh) | 2016-02-17 |
Family
ID=48062260
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310020198.5A Active CN103049623B (zh) | 2013-01-18 | 2013-01-18 | 一种激光焊接热源模型的建立方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103049623B (zh) |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104809291A (zh) * | 2015-04-27 | 2015-07-29 | 江苏金通灵流体机械科技股份有限公司 | 一种基于ansys的双相不锈钢与异种钢焊接变形预测方法 |
CN105138772B (zh) * | 2015-08-24 | 2018-12-07 | 北京星航机电装备有限公司 | 变截面构件电子束焊的有限元模拟方法 |
CN105718690A (zh) * | 2016-01-26 | 2016-06-29 | 南京航空航天大学 | 基于时间和空间主动跟踪的激光3d打印熔池凝固行为数值模拟方法 |
CN106529047B (zh) * | 2016-11-15 | 2019-07-16 | 河海大学常州校区 | 一种串列双丝埋弧焊数值模拟热源模型的建模方法 |
CN106529051B (zh) * | 2016-11-15 | 2019-02-19 | 河海大学常州校区 | 一种单丝埋弧焊数值模拟热源模型参数的确定方法 |
CN106709176A (zh) * | 2016-11-29 | 2017-05-24 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种激光熔化沉积成形熔池动力学数值模拟技术 |
CN106670666A (zh) * | 2016-12-14 | 2017-05-17 | 南京航空航天大学 | 基于精准能量分配的激光‑电弧复合加工的能量分配系数模型的构建方法 |
CN107330181A (zh) * | 2017-06-28 | 2017-11-07 | 上海交通大学 | 预测激光焊接变形量的实现方法 |
CN107695535B (zh) * | 2017-10-11 | 2019-09-10 | 广东正业科技股份有限公司 | 一种硬脆性材料激光切割辅助方法、模型创建方法及系统 |
CN108062432A (zh) * | 2017-11-17 | 2018-05-22 | 西安铂力特增材技术股份有限公司 | 一种激光选区熔化过程的数值模拟方法 |
CN108038266A (zh) * | 2017-11-17 | 2018-05-15 | 西安铂力特增材技术股份有限公司 | 一种选择性激光修复过程的数值模拟方法 |
CN108304641B (zh) * | 2018-01-29 | 2021-06-25 | 北京航空航天大学 | 一种高能容干式摩擦元件耐热性能设计方法 |
CN109299554A (zh) * | 2018-09-30 | 2019-02-01 | 西南石油大学 | 一种激光-电弧复合焊接热源模型的建立方法 |
CN109933007A (zh) * | 2019-01-30 | 2019-06-25 | 中国科学院金属研究所 | 一种激光-mig复合焊接热源模型的建立方法 |
CN109885946B (zh) * | 2019-02-26 | 2022-12-30 | 大族激光科技产业集团股份有限公司 | 一种确定复合热源的能量分配的方法及焊接模拟方法 |
CN110598358A (zh) * | 2019-09-26 | 2019-12-20 | 华中科技大学 | 一种增材制造应力变形仿真方法、装置、设备及存储介质 |
CN112214863B (zh) * | 2020-07-31 | 2024-01-02 | 中国科学院金属研究所 | 一种奥氏体不锈钢焊接热影响区粗晶区特征热循环曲线的获取方法 |
CN112016228B (zh) * | 2020-08-14 | 2022-08-19 | 河海大学常州校区 | 一种基于水冷压缩系数的水下焊接热源模型的建模方法 |
CN111950183A (zh) * | 2020-08-20 | 2020-11-17 | 南京中车浦镇城轨车辆有限责任公司 | 一种适合于复杂大构件的热源加载方法及系统 |
CN111950182A (zh) * | 2020-08-20 | 2020-11-17 | 南京中车浦镇城轨车辆有限责任公司 | 一种提高大型复杂构件有限元计算效率的模型简化方法 |
CN113268898A (zh) * | 2021-03-25 | 2021-08-17 | 南京航空航天大学 | 一种tc4钛合金h型结构双激光束双侧同步焊接的变形模拟方法 |
CN115630512A (zh) * | 2022-10-28 | 2023-01-20 | 江苏徐工工程机械研究院有限公司 | 一种摆动激光焊接等效热源建模方法及摆动激光焊接仿真方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101324469A (zh) * | 2008-08-06 | 2008-12-17 | 中国航空工业第一集团公司北京航空制造工程研究所 | 一种激光焊接温度场三维测量方法 |
CN101559511A (zh) * | 2009-05-22 | 2009-10-21 | 清华大学 | 一种以温度为控制变量的焊接数值模拟计算方法 |
-
2013
- 2013-01-18 CN CN201310020198.5A patent/CN103049623B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101324469A (zh) * | 2008-08-06 | 2008-12-17 | 中国航空工业第一集团公司北京航空制造工程研究所 | 一种激光焊接温度场三维测量方法 |
CN101559511A (zh) * | 2009-05-22 | 2009-10-21 | 清华大学 | 一种以温度为控制变量的焊接数值模拟计算方法 |
Non-Patent Citations (2)
Title |
---|
5A06铝合金间断焊温度场的数值模拟;曾志等;《天津大学学报》;20080731;第41卷(第7期);第849-853页 * |
不锈钢焊接温度场的三维数值模拟;董志波等;《焊接学报》;20040430;第25卷(第2期);第9-15页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103049623A (zh) | 2013-04-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103049623B (zh) | 一种激光焊接热源模型的建立方法 | |
CN105138772B (zh) | 变截面构件电子束焊的有限元模拟方法 | |
CN104809291A (zh) | 一种基于ansys的双相不锈钢与异种钢焊接变形预测方法 | |
CN103246774B (zh) | 数值模拟p92钢管焊接热影响区宽度的方法 | |
CN109513924B (zh) | 一种激光选区熔化过程中表面粗糙度控制方法 | |
CN104899345A (zh) | 用于确定复杂曲面形状工件激光喷丸成形工艺参数的方法 | |
CN107194057A (zh) | 壁板铆接扭翘变形预测方法 | |
Arora et al. | Prediction of temperature distribution and displacement of carbon steel plates by FEM | |
CN106055766A (zh) | 一种随机损伤对点阵夹层板动力学特性影响的分析方法 | |
Fossati et al. | Numerical modelling of crack growth profiles in integral skin-stringer panels | |
Piekarska et al. | Analytical and numerical model of laser welding phenomena with the initial preheating | |
CN113065272A (zh) | 一种焊接变形数值模拟快速实现方法 | |
CN110895634B (zh) | 一种面向精确变形控制的2.5mm厚度铝锂合金T型接头焊接结构集成模拟方法 | |
Liu et al. | Optimization design of drawbead in drawing tools of autobody cover panel | |
Miah et al. | Friction stir welding on temperature field for aluminum alloy based on combined HSM | |
CN109885946A (zh) | 一种确定复合热源的能量分配的方法及焊接模拟方法 | |
Papadakis et al. | Numerical investigation of the influence of preliminary manufacturing processes on the crash behaviour of automotive body assemblies | |
Lyu et al. | Strategies for solving high-fidelity aerodynamic shape optimization problems | |
Shen et al. | Investigation of molten pool and heat transfer mechanisms after moving annular laser incidence metal by lattice Boltzmann method | |
CN111914383A (zh) | 一种基于壳-实体耦合单元的t型结构激光焊接模拟方法 | |
Merle et al. | Gradient-based optimization of a flexible long-range transport aircraft using a high-dimensional CAD-ROM parameterization | |
Wang et al. | Optimization design of towbarless aircraft tractor frame based on ANSYS Workbench | |
Mohan et al. | Induction welding simulations of curvature-based thermoplastic composites | |
Shen et al. | Mesoscopic simulation of molten pool heat transfer and fluid flow with moving annular laser via lattice Boltzmann method | |
Tong et al. | Research on the Application of Structural Topology Optimisation in the High-Precision Design of a Press Machine Frame |
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 |