CN109885946B - 一种确定复合热源的能量分配的方法及焊接模拟方法 - Google Patents

一种确定复合热源的能量分配的方法及焊接模拟方法 Download PDF

Info

Publication number
CN109885946B
CN109885946B CN201910142616.5A CN201910142616A CN109885946B CN 109885946 B CN109885946 B CN 109885946B CN 201910142616 A CN201910142616 A CN 201910142616A CN 109885946 B CN109885946 B CN 109885946B
Authority
CN
China
Prior art keywords
heat source
welding
graph
sub
model
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
Application number
CN201910142616.5A
Other languages
English (en)
Other versions
CN109885946A (zh
Inventor
蒋磊
范国成
万虹
陈根余
陈焱
高云峰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Han s Laser Technology Industry Group Co Ltd
Hans Laser Smart Equipment Group Co Ltd
Original Assignee
Han s Laser Technology Industry Group Co Ltd
Hans Laser Smart Equipment Group Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Han s Laser Technology Industry Group Co Ltd, Hans Laser Smart Equipment Group Co Ltd filed Critical Han s Laser Technology Industry Group Co Ltd
Priority to CN201910142616.5A priority Critical patent/CN109885946B/zh
Publication of CN109885946A publication Critical patent/CN109885946A/zh
Application granted granted Critical
Publication of CN109885946B publication Critical patent/CN109885946B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Arc Welding In General (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

本发明涉及焊接加工领域,具体涉及一种确定复合热源的能量分配的方法及焊接模拟方法。确定复合热源的能量分配的方法包括:选择焊接位置处的一截面,获取截面上的熔池图形;将所述熔池图形按各类几何形状划分成不同的子图形;计算所述子图形的能量分配系数;计算所述子图形的能量输入值。由于此能量输入值接近实际的能量输入量,在进行热源校核时,通过将此能量输入值作为初始值导入至热源模型中,再对比经由数值计算得到的熔池图形与实际工件的熔池图形,只需微调能量输入值等相关热源参数即可实现基本吻合,而无需繁琐地进行多次调整,极大地缩短了复合热源的校核周期,进而缩短了整个焊接仿真模拟的耗时。

Description

一种确定复合热源的能量分配的方法及焊接模拟方法
技术领域
本发明涉及焊接加工领域,具体涉及一种确定复合热源的能量分配的方法及焊接模拟方法。
背景技术
随着计算仿真技术的发展,焊接仿真技术也在快速发展,越来越多的工程师利用仿真技术来研究焊接过程和现象,通过优化焊缝结构和焊接工艺,以实现有效控制焊接变形和焊后开裂等问题。目前焊接仿真技术还不是很成熟,存在如复合热源的校核周期长、整体焊接仿真模拟时间长等问题。
发明内容
本发明要解决的技术问题在于,针对现有技术的上述缺陷,提供一种确定复合热源的能量分配的方法,克服现有仿真技术的复合热源校核周期长的缺陷。
本发明解决其技术问题所采用的技术方案是,提供一种确定复合热源的能量分配的方法,包括:步骤A、选择焊接位置处的一截面,获取截面上的熔池图形;步骤B、将所述熔池图形按各类几何形状划分成不同的子图形;步骤C、计算所述子图形的能量分配系数;步骤D、计算所述子图形的能量输入值,所述子图形的能量输入值为总输入能量值与所述子图形的能量分配系数的乘积。
本发明的更进一步优选方案是,所述步骤A中,所述截面为沿熔池深度方向选取的一平面。
本发明的更进一步优选方案是,所述子图形的能量分配系数为所述子图形的面积除以所述熔池图形的总面积的商。
本发明的更进一步优选方案是,所述步骤B具体包括:步骤1、建立焊接工件的数值模型;步骤2、将热源加载顺序、热源加载位置以及由确定复合热源的能量分配的方法得出的能量输入值导入至热源模型;步骤3、将焊接工况参数导入热源模型;步骤4、热源模型计算并分析结果。
本发明的更进一步优选方案是,所述步骤2与步骤3之间还包括:步骤20、对复合热源进行校核。
本发明的更进一步优选方案是,所述热源加载顺序为热源在每条焊缝上的移动路线或各焊缝焊接的先后顺序。
本发明的更进一步优选方案是,所述热源加载位置与时间、速度之间的关系为:
Figure BDA0001979015830000021
其中,n为组数;Sxn、Syn、Szn分别为第n组热源的起始位置在X、Y和Z轴方向的坐标;Vx、Vy、Vz分别为第n组热源在X、Y、Z轴方向对应的移动速度;t为热源加载的时间。
本发明的更进一步优选方案是,所述步骤1具体包括以下步骤:建立焊接工件的实体模型;选择网格类型以及调整网格质量,在实体模型上划分网格;将网格进行分组并对分组后的组分进行命名;生成焊接工件的数值模型。
本发明的更进一步优选方案是,所述步骤1还包括以下步骤:通过省略焊接坡口、余高、圆角、倒角实现简化焊接工件的实体模型。
本发明的更进一步优选方案是,所述步骤3具体包括以下步骤:输入热交换的形式;输入待焊工件的固定方式及位置;输入接触形式;输入待焊工件所受的外力作用。
本发明提供的一种确定复合热源的能量分配的方法,先获取焊接位置处的一个截面上的熔池图形,然后将该熔池图形划分成不同形状的子图形,并计算出各子图形对应的能量输入值。由于此能量输入值接近实际的能量输入量,在进行热源校核时,通过将此能量输入值作为初始值导入至热源模型中,再对比经由数值计算得到的熔池图形与实际工件的熔池图形,只需微调能量输入值等相关热源参数即可实现基本吻合,而无需繁琐地进行多次调整,极大地缩短了复合热源的校核周期,进而缩短了整个焊接仿真模拟的耗时。
附图说明
下面将结合附图及实施例对本发明作进一步说明,附图中:
图1是本发明实施例的确定复合热源的能量分配的方法流程示意图;
图2是本发明实施例的焊接模拟方法的流程示意图;
图3是本实施例的焊接模拟方法的步骤1的流程示意图;
图4是本实施例的焊接模拟方法的步骤3的流程示意图;
图5是本发明实施例的焊接工件的立体模型图;
图6是本发明实施例的焊接工件的实体模型示意图;
图7是本发明实施例的横截面上的熔池示意图;
图8是本发明实施例的温度云图;
图9是本发明实施例的应力云图;
图10是本发明实施例的变形云图。
具体实施方式
现结合附图,对本发明的较佳实施例作详细说明。
在本实施例中,焊接模拟软件选用Sysweld。Sysweld是完全实现了机械、热传导和金属冶金的耦合计算的计算机模拟开发系统软件。焊接模拟软件中涉及的数值模型为有限元模型。热源在有限元计算里是以一个函数形式来表示。而热源模型是对作用于焊接工件上的、在一定时间和位置上的热输入分布特点的一种数学表达。
如图1所示,本实施例的确定复合热源的能量分配的方法,包括以下步骤:
步骤A、选择焊接位置处的一截面,获取截面上的熔池图形;
步骤B、将所述熔池图形按各类几何形状划分成不同的子图形;
步骤C、计算所述子图形的能量分配系数μ;所述能量分配系数μ为各子图形的面积Mn除以熔池图形的总面积M的商,即μ=Mn/M;
步骤D、计算所述子图形的能量输入值Qn;所述能量输入值Qn为总输入能量值Q与所述子图形的能量分配系数μ的乘积,即
Qn=Q*μ=Q×(Mn/M);
其中,总输入能量值为焊接头所提供的能量。总输入能量值Q=P/V,P为焊接功率,V为焊接速度。
本发明实施例的确定复合热源的能量分配的方法,先获取焊接位置处的一个截面上的熔池图形,然后将该熔池图形划分成不同形状的子图形,并计算出各子图形对应的能量输入值。由于此能量输入值接近实际的能量输入量,在进行热源校核时,通过将此能量输入值作为初始值导入至热源模型中,再对比经由数值计算得到的熔池图形与实际工件的熔池图形,只需微调能量输入值等相关热源参数即可实现基本吻合,而无需繁琐地进行多次调整,极大地缩短了复合热源的校核周期,进而缩短了整个焊接仿真模拟的耗时,提高了生产效率。
更佳地,如图2所示,本实施例还提供一种焊接模拟方法,包括以下步骤:
步骤1、建立焊接工件的数值模型,并加载至热源模型;
步骤2、将热源加载顺序、热源加载位置以及由上述的确定复合热源的能量分配的方法得出的能量输入值Qn导入至热源模型;
步骤20、对复合热源进行校核
步骤3、将焊接工况参数导入热源模型;
步骤4、热源模型计算并分析结果。
通过建立焊接工件的数值模型并加载至热源模型,再将热源加载顺序、热源加载位置和熔池上各子图形对应的能量输入值等工艺参数导入热源模型中,并将此热源模型用以实现对数值的模拟计算。经热源校核后,将实际生产中的焊接工况导入至热源模型中,焊接模拟软件进行焊接仿真分析,得出模拟结果,用以指导实际的焊接生产,尤其适用于大尺寸、焊缝长、具备复合热源且焊接工艺复杂的生产加工,达到提高焊接质量和生产效率的效果。
如图5所示,本实施例以大尺寸的U字形焊接工件为例。第一母材1通过第一焊缝4与第二母材2连接,第二母材2通过第二焊缝5与第三母材3连接。由于U字形焊接工件的结构较为简单,可采用软件Visual-environment来进行建模及进行网格划分等。焊接头(图未示)的焊接速度V为1.2m/min,焊接功率P为2000w。
具体的,如图3所示,所述步骤1具体包括:
步骤11、建立焊接工件的实体模型,并加载至热源模型;
所述步骤11中,根据待焊工件的形状和尺寸建立出三维实体模型,如图6所示。
步骤12、简化焊接工件的实体模型;
所述步骤12中,将待焊工件的实体模型的焊接坡口、余高、圆角、倒角等作简化,以提高在软件中的响应速度和运算速度。
步骤13、选择网格类型以及调整网格质量,并在实体模型上划分网格;
所述步骤13中,网格类型选择线性六面体网格类型和8节点;网格质量的调整包括将网格密度设为2、焊缝宽度方向的单元数设为4个、熔池长度方向的单元数设为4个、焊缝与焊接母材之间过渡部分与焊缝的网格密度之比为2。
步骤14、将网格进行分组并对分组后的组分进行命名;
所述步骤14中,将网格分组为焊缝组、焊接母材组、热源加载组、热交换组及约束组等组分。各组分命名(图未示)如下:焊缝组命名为J01_WIRE、J02_WIRE;母材组命名为COMPONENT_01、COMPONENT_02、COMPONENT_03;两组热源加载组分别命名:第一组热源加载组可命名为J01_PATH、J01_REF、J01_SNO、J01_SEL、J01_ENO;第二组热源加载组可命名为J02_PATH、J02_REF、J02_SNO、J02_SEL、J02_ENO;热交换组命名为SKIN_01;约束组的命名可为CLAMP_01、CLAMP_02和CLAMP_03。
步骤15、生成焊接工件的数值模型;
所述步骤15中,根据上述的实体模型、网格类型、网格质量、划分网格、网格分组、组分命名,生成Sysweld软件能识别的数值模型,即焊接工件被转化成了可以在求解过程中能够被识别的网格。将此数值模型保存并加载至Sysweld软件中的热源模型中。
如图4所示,所述步骤2具体包括:
步骤21、确定热源加载顺序;
所述步骤21中,所述热源加载顺序为热源在每条焊缝上的移动路线或各焊缝焊接的先后顺序。移动路线如图6中的虚线箭头所示。
步骤22、建立各组热源加载位置S与时间t、速度V之间的关系;
所述步骤12中,任意一组热源加载位置加载位置Sn(v,t)与时间t、速度V(vx(t),vy(t),vz(t))之间的关系如下:
Figure BDA0001979015830000071
其中,n为组数;Sxn、Syn、Szn分别为第n组热源的起始位置在X、Y和Z轴方向的坐标;Vx、Vy、Vz分别为第n组热源在X、Y、Z轴方向对应的移动速度;t为热源加载的时间。
在本实施例中,参照图5中的坐标轴,第一组热源6的加载位置S1(v,t)与时间t、速度V(vx(t),vy(t),vz(t))之间的关系为:
Figure BDA0001979015830000072
第二组热源7的加载位置S2(v,t)与时间t、速度V(vx(t),vy(t),vz(t))之间的关系为:
Figure BDA0001979015830000073
步骤23、确定复合热源的能量分配;
所述步骤23包括以下步骤:
步骤A、选择焊接位置处的一截面,获取截面上的熔池图形;
所述步骤A中,在焊接位置处,沿熔池的深度方向选取一个截面,获取截面上的熔池图形的形状和尺寸,如图7所示。
步骤B、将所述熔池图形按各类几何形状划分成不同的子图形;
所述步骤B中,根据熔池图形的整体形状将熔池图形划分成不同几何形状类型的子图形,并计算出各个划分好的子图形的面积Mn。继续参考图7,在本实施例中,熔池图形被划分成一个类长方形子图形和一个类半圆形子图形,类长方形子图形的面积M1=3.75mm2,类半圆形子图形的面积M2=1mm2
步骤C、计算所述子图形的能量分配系数μ;能量分配系数μ为各子图形的面积Mn除以熔池图形的总面积M的商,即μ=Mn/M;
步骤D、计算所述子图形的能量输入值Qn;所述能量输入值Qn为总输入能量值Q与所述子图形的能量分配系数μ的乘积,即
Qn=Q*μ=Q×(Mn/M);其中,Q=P/V,P为焊接功率,V为焊接速度。
代入数值后可得:
Q=P/V=2000/(1.2*1000/60)=100J/mm
Q1=Q*M1/(M1+M2)=100*3.75/(3.75+1)=78.9J/mm
Q2=Q*M2/(M1+M2)=100*1/(3.75+1)=21.1J/mm
步骤24、将热源加载顺序、热源加载位置、能量输入值以及其他相关热源参数输入到热源模型中。
步骤20、对复合热源进行校核;
所述步骤20中,将经由焊接模拟软件Sysweld模拟计算得到的熔池与焊接工件上的实际熔池进行对比。若校得的每一层的热源熔池形状与实际焊接所形成的的熔池形状不吻合,则微调相关热源参数;若基本吻合,则用该组热源模型可以很好地对焊接工件进行数值模拟,保存此校正好的热源文件。
如图4所示,所述步骤3具体包括:
步骤31、输入热交换的形式;
所述步骤31中,热交换形式包括空冷、水冷和油冷。本实施例的热交换形式为空冷。
步骤32、输入待焊工件的固定方式及位置;
所述步骤32中,固定方式为自由约束,固定位置无。
步骤33、输入接触形式;
所述步骤33中,在本实施例中无接触。
步骤34、输入待焊工件所受的外力作用;
所述步骤34中,外力作用包括压力载荷、集中力作用和扭矩等。在本实施例中,无外力作用。
将实际加工生产时的上述的热交换的形式、固定方式及位置、接触形式、外力作用等焊接工况参数输入热源模型中。
一同参考图8-10,所述步骤5具体包括:
载入数值模型、校正好的热源模型及输入焊接工况后,焊接模拟软件Sysweld开始进行焊接模拟。热源模型计算并对结果进行焊接温度场、变形和应力场分析。最终分析结果显示为温度云图(图8)、应力云图(图9)和变形云图(图10)。
通过其分析的结果,参考温度云图、应力云图和变形云图,以优化焊接工艺,提高生产效率和品质。
应当理解的是,以上实施例仅用以说明本发明的技术方案,而非对其限制,对本领域技术人员来说,可以对上述实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而所有这些修改和替换,都应属于本发明所附权利要求的保护范围。

Claims (10)

1.一种确定复合热源的能量分配的方法,其特征在于,包括:
步骤A、选择焊接位置处的一截面,获取截面上的熔池图形;
步骤B、将所述熔池图形按各类几何形状划分成不同的子图形;
步骤C、计算所述子图形的能量分配系数;
步骤D、计算所述子图形的能量输入值,所述子图形的能量输入值为总输入能量值与所述子图形的能量分配系数的乘积。
2.根据权利要求1所述的确定复合热源的能量分配的方法,其特征在于,所述步骤A中,所述截面为沿熔池深度方向选取的一平面。
3.根据权利要求1所述的确定复合热源的能量分配的方法,其特征在于,所述子图形的能量分配系数为所述子图形的面积除以所述熔池图形的总面积的商。
4.一种焊接模拟方法,其特征在于,包括:
步骤1、建立焊接工件的数值模型,并加载至热源模型;
步骤2、将热源加载顺序、热源加载位置以及由权利要求1所述的确定复合热源的能量分配的方法得出的能量输入值导入至热源模型;
步骤3、将焊接工况参数导入热源模型;
步骤4、热源模型计算并分析结果。
5.根据权利要求4所述的焊接模拟方法,其特征在于,所述步骤2与步骤3之间还包括:
步骤20、对复合热源进行校核。
6.根据权利要求4所述的焊接模拟方法,其特征在于,所述热源加载顺序为热源在每条焊缝上的移动路线或各焊缝焊接的先后顺序。
7.根据权利要求4所述的焊接模拟方法,其特征在于,所述热源加载位置与时间、速度之间的关系为:
Figure FDA0003887917700000021
其中,n为组数;Sxn、Syn、Szn分别为第n组热源的起始位置在X、Y和Z轴方向的坐标;Vx、Vy、Vz分别为第n组热源在X、Y、Z轴方向对应的移动速度;t为热源加载的时间。
8.根据权利要求4所述的焊接模拟方法,其特征在于,所述步骤1具体包括以下步骤:
建立焊接工件的实体模型;
选择网格类型以及调整网格质量,在实体模型上划分网格;
将网格进行分组并对分组后的组分进行命名;
生成焊接工件的数值模型。
9.根据权利要求8所述的焊接模拟方法,其特征在于,所述步骤1还包括以下步骤:
通过省略焊接坡口、余高、圆角、倒角实现简化焊接工件的实体模型。
10.根据权利要求4至9任一项所述的焊接模拟方法,其特征在于,所述步骤3具体包括以下步骤:
输入热交换的形式;
输入待焊工件的固定方式及位置;
输入接触形式;
输入待焊工件所受的外力作用。
CN201910142616.5A 2019-02-26 2019-02-26 一种确定复合热源的能量分配的方法及焊接模拟方法 Active CN109885946B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910142616.5A CN109885946B (zh) 2019-02-26 2019-02-26 一种确定复合热源的能量分配的方法及焊接模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910142616.5A CN109885946B (zh) 2019-02-26 2019-02-26 一种确定复合热源的能量分配的方法及焊接模拟方法

Publications (2)

Publication Number Publication Date
CN109885946A CN109885946A (zh) 2019-06-14
CN109885946B true CN109885946B (zh) 2022-12-30

Family

ID=66929546

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910142616.5A Active CN109885946B (zh) 2019-02-26 2019-02-26 一种确定复合热源的能量分配的方法及焊接模拟方法

Country Status (1)

Country Link
CN (1) CN109885946B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110889213B (zh) * 2019-11-18 2023-05-02 重庆理工大学 一种基于切片分割加载体的热源焊接模拟方法
CN114799639B (zh) * 2022-05-16 2023-03-14 广州东焊智能装备有限公司 一种基于工业机器人的曲线焊接系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103049623A (zh) * 2013-01-18 2013-04-17 哈尔滨工业大学 一种激光焊接热源模型的建立方法
CN103706938A (zh) * 2013-12-11 2014-04-09 中国航空工业集团公司北京航空制造工程研究所 一种利用多热源控制残余应力的电子束焊接的方法及系统
CN106670666A (zh) * 2016-12-14 2017-05-17 南京航空航天大学 基于精准能量分配的激光‑电弧复合加工的能量分配系数模型的构建方法
CN108672893A (zh) * 2018-05-21 2018-10-19 天津工业大学 一种非对称角焊缝熔透形态及熔深的控制方法
CN109299554A (zh) * 2018-09-30 2019-02-01 西南石油大学 一种激光-电弧复合焊接热源模型的建立方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103049623A (zh) * 2013-01-18 2013-04-17 哈尔滨工业大学 一种激光焊接热源模型的建立方法
CN103706938A (zh) * 2013-12-11 2014-04-09 中国航空工业集团公司北京航空制造工程研究所 一种利用多热源控制残余应力的电子束焊接的方法及系统
CN106670666A (zh) * 2016-12-14 2017-05-17 南京航空航天大学 基于精准能量分配的激光‑电弧复合加工的能量分配系数模型的构建方法
CN108672893A (zh) * 2018-05-21 2018-10-19 天津工业大学 一种非对称角焊缝熔透形态及熔深的控制方法
CN109299554A (zh) * 2018-09-30 2019-02-01 西南石油大学 一种激光-电弧复合焊接热源模型的建立方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于形状先验活动轮廓的焊接图像熔池分割方法;罗祥等;《焊接学报》;20171125(第11期);全文 *

Also Published As

Publication number Publication date
CN109885946A (zh) 2019-06-14

Similar Documents

Publication Publication Date Title
Zhang et al. CAD-based design and pre-processing tools for additive manufacturing
Montevecchi et al. Finite element modelling of wire-arc-additive-manufacturing process
CN108022236B (zh) 基于焊缝图像提取拟合的焊接应力变形预测方法及系统
CN104504186B (zh) 一种基于3d打印数据处理软件平台的实现方式
KR101996933B1 (ko) 적층 생성된 구조물에서 응력들 및 형상 편차들을 결정하기 위한 방법, 컴퓨터-판독가능 데이터 캐리어, 컴퓨터 프로그램, 및 시뮬레이터
CN109284524B (zh) 一种创建高精度增材制造有限元模型的方法
CN108399307A (zh) 一种激光3d打印有限元模拟方法
CN108646689B (zh) 一种基于焊接的虚拟生产制造方法及系统
CN112149330B (zh) 一种风电塔筒油封平台焊接残余应力预测、焊接工艺优化方法
US6813533B1 (en) Method for simulation of laser material deposition
CN109885946B (zh) 一种确定复合热源的能量分配的方法及焊接模拟方法
CN112512729B (zh) 用于确定针对增材制造方法的构造规范的方法
CN109933007A (zh) 一种激光-mig复合焊接热源模型的建立方法
CN112883623A (zh) 一种Fe-36%Ni合金PMIG摆动焊热-机械耦合有限元建模方法
CN112380650A (zh) 工作装置结构件设计方法
CN110153582A (zh) 焊接方案生成方法、装置以及焊接系统
Xu et al. Fabrication strategy and macroscopic defect control of large-size component based on double-wire arc additive manufacturing
WO2017141070A1 (en) Method and system for rapid prototyping via virtual simulation
CN102073769B (zh) 一种快速预测铸件热节的方法
CN111950182A (zh) 一种提高大型复杂构件有限元计算效率的模型简化方法
CN112883518A (zh) 一种预测tig增材与轧制复合制造件残余应力及变形的方法
Zhou et al. Simulation and optimization framework for additive manufacturing processes
CN113283139A (zh) 基于ansys的slm单道多层扫描成形数值模拟方法
Potočnik et al. Genetic Algorithm-Based Framework for Optimization of Laser Beam Path in Additive Manufacturing
CN113849923B (zh) 一种基于abaqus的钢结构整体焊接变形的数值模拟方法

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