CN110991009A - 一种基于上覆荷载作用下埋地管线下方土体流失对管线受力形变的确定方法 - Google Patents

一种基于上覆荷载作用下埋地管线下方土体流失对管线受力形变的确定方法 Download PDF

Info

Publication number
CN110991009A
CN110991009A CN201911092349.1A CN201911092349A CN110991009A CN 110991009 A CN110991009 A CN 110991009A CN 201911092349 A CN201911092349 A CN 201911092349A CN 110991009 A CN110991009 A CN 110991009A
Authority
CN
China
Prior art keywords
pipeline
soil
stress
pipe
deformation
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
Application number
CN201911092349.1A
Other languages
English (en)
Other versions
CN110991009B (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.)
Ningbo University
Original Assignee
Ningbo University
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 Ningbo University filed Critical Ningbo University
Priority to CN201911092349.1A priority Critical patent/CN110991009B/zh
Publication of CN110991009A publication Critical patent/CN110991009A/zh
Application granted granted Critical
Publication of CN110991009B publication Critical patent/CN110991009B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

一种基于上覆荷载作用下埋地管线下方土体流失对管线受力形变的确定方法,根据埋地管线从受力平衡到失衡的受力特征,将管线和土体抽离出来构建管‑土耦合模型进行分析;提出了基于角点法和Winkler弹性地基梁理论,给出了管线受力变形方程,计算出了管线的位移、剪力和弯矩。保证了计算结果的准确性,克服了传统方法因考虑管线受力因素过于单一,导致计算结果不能正确反映实际工程中管线的受力特征的局限性。

Description

一种基于上覆荷载作用下埋地管线下方土体流失对管线受力 形变的确定方法
技术领域
本发明属于工程技术领域,涉及一种基于上覆荷载作用条件下埋地管线下方土体流失对管线受力与形变的定量计算方法,特别是本发明基于角点法和Winkler弹性地基梁理论,将管与管周土体脱离区段划分为管土分离区和管土半分离区,确定了管线的位移、剪力和弯矩的计算方法,进而提出了解决了实际工程中在上覆荷载和土体流失的共同作用下管线的受力和变形问题的技术方案,为埋地管线的保护措施和安全性分析提供了参考依据。
背景技术
埋地管线作为天然气、石油、自来水等资源的运输工具,具有运输快速、便捷、不占地面空间等优点,它是城市供水、供气、热力供给、排污等地下生命线系统的基本组成部分。管线因施工方式不恰当、材料老化、养护不到位、管壁腐蚀等问题,容易引起管线的渗漏。其次,管线周围的土体发生不均匀沉降或者邻近位置的管线破裂以及地下水的过度开采,都会导致管线周围土体流失、力学特征发生改变。另外,由于降水、人类的活动等原因也会引起管线的上覆荷载发生变化,改变管线原有的受力平衡。这些因素容易引起管线发生破坏,给人民的生产、生活带来不便,也造成了一定的经济损失。目前,国内外诸多学者采用不同的方法对埋地管线的力学特征展开研究,一些代表性研究如下:
T.E.Vorster等考虑了隧道的施工、土的非线性问题等因素对既有管线的影响,提出了一种计算连续(或刚性连接)管道最大弯矩的方法;
Wang Chenqi等采用Winkler弹性地基梁理论,提出了管-土分离模型并分析了非均质滑坡驱动力下管线的位移、弯矩等受力特征,得出了随着分离段长度的增加,管道的挠度和应力变化的规律;
龚晓南等利用Winkler弹性地基短梁理论,在有限差分法的基础上建立了地面超载作用下对埋地管道的位移影响分析模型,探讨了地面超载下邻近埋地管道位移的规律。
刘金梅等基于Drucker-Prager弹塑性准则,对目前常用的大管径直埋管道进行了沉降变形模拟;
史江伟等探讨了不均匀的土体位移对管线变形影响的规律,用ABAQUS有限元软件对管线间的相互作用进行辅助计算,为隧道开挖工程提供了计算管线弯曲变形的方法;
吴张中等分析了采空塌陷作用下油气管道的受力特征,结合急倾斜煤层开采实例,利用FLAC3D进行数值模拟,验证了管道在受采空塌陷作用下经过的各个阶段及受力特征。
李新亮等基于Boussinesq解和Winkler弹性地基梁模型,运用光纤布拉格光栅(Fiber Bragg Grating,FBG)传感器技术监测管道的位移,研究了交通荷载和土体各项性能指标对埋地管道力学性状的影响规律;
朱治齐等制作管线试验模型研究了工程荷载对地埋管线纵向响应的影响机理,通过检测得到了地埋管线的位移和弯矩的变化规律。
上述研究中考虑管线受力因素过于单一,只研究了上覆荷载对管线的影响或者管线下方土体流失对管线的影响,目前没有提出能够综合考虑了上覆荷载和土体流失的因素对管线受力变形影响的方法。因此,迫切需要提出一种能够正确反映管线在地表荷载和土体荷载综合作用下的受力与变形的计算方法。
发明内容
为了克服现有的技术不能准确的计算地表荷载和土体荷载综合作用下管线的受力与变形,本发明提出了基于角点法和Winkler弹性地基梁理论,给出了管线受力变形方程,计算出了管线的位移、剪力和弯矩。保证了计算结果的准确性,避免了传统方法考虑管线受力因素过于单一,导致计算结果不能正确反映实际工程中管线的受力特征。
本发明解决其技术问题所采用的技术方案是:
一种基于上覆荷载作用下埋地管线下方土体流失对管线受力形变的确定方法,包括以下步骤:
1)根据埋地管线从受力平衡到失衡的受力特征,将管线和土体抽离出来构建管-土耦合模型进行分析;
2)运用角点法求得地表荷载作用下管线的附加应力σz’
3)地基土假定为各向同性、均质线性变形体,不考虑土体和边界效应的影响,运用自重应力公式计算土体作用于管线的有效自重应力σc
4)考虑到管线直径比较小,管顶和管底的受力基本相同,梁的长与管线相同,宽为管线的直径,作用在埋地管线上的上覆附加荷载为:q(x)=(σz’c)D;
5)埋地管线埋置于土体中的情况比较复杂,受到的影响因素比较多,将埋地管线的上覆附加荷载假定为q(x)=ax2+bx+c,当管线的附加荷载系数a,b为0时,可求得q(x)=c,能够反映步骤4)中管线的受力情况;
6)根据埋地管线在土体中承载环境的差别,将管线形变段划分为管土分离区的变形和管土半分离区的变形;
7)通过管土分离区的受力与变形方程计算出管土分离区管线的位移y1、剪力Q1和弯矩M1
8)对于管土半分离区,将管线视为放置于半无限弹性地基上的半无限梁,运用Winkler地基梁理论推导出管线变形方程计算管土半分离区的位移y2、剪力Q2和弯矩M2
9)在算例的基础上,分析管线的位移、剪力以及弯矩的变化情况。
进一步,所述步骤2)中,管线的附加应力方程形式表达为:
Figure BDA0002267215820000031
其中,σz’为平行于z’坐标轴的应力,单位:kPa;p为作用于地表的竖向均布荷载,单位:kPa;l为地表荷载的长度;b为地表荷载的宽度;z’为地表荷载距离埋地管线的深度。
再进一步,所述步骤3)中,管线的有效自重应力方程形式表达为:
Figure BDA0002267215820000032
其中,σc为天然地面下任意深度z’处的竖向有效自重应力,单位:kPa;n为深度z’范围内的土层总数;hi为第i层土的厚度,单位:m;γi为第i层土的天然重度,单位:kN/m3
更进一步,所述步骤5)中,管土分离区的受力与变形方程形式表达为:
Figure BDA0002267215820000033
Figure BDA0002267215820000034
Figure BDA0002267215820000035
Figure BDA0002267215820000041
其中,y1为管土分离区管线的位移;M1为管土分离区管线的弯矩;Q1为管土分离区管线的剪力;E为管线的弹性模量;I为埋地管线的截面惯性矩;a,b,c为管线附加荷载系数;β为特征系数;yP为管线在P4点受到的位移,单位:m;MP为管线在P4点受到的弯矩,单位;kN·m;QP为管线在P4点受到的剪力,单位:kN。
更进一步,所述步骤6)中,管土分离区的受力与变形方程形式表达为:
Figure BDA0002267215820000042
Figure BDA0002267215820000043
Figure BDA0002267215820000044
其中,y2为管土半分离区管线的位移;M2为管土半分离区管线的弯矩;Q2为管土半分离区管线的剪力;k为地基弹性阻力系数。
本发明的有益效果为:保证了计算结果的准确性,避免了传统方法考虑管线受力因素过于单一,导致计算结果不能正确反映实际工程中管线的受力特征。
附图说明
图1为管-土耦合模型示意图,其中,1是地表载荷,2是埋地管道,3是土体流失区,4是地基土;
图2为管土分离段受力模型示意图;
图3为管土半分离段受力模型示意图;
图4为管线竖向位移图;
图5为管线剪力图;
图6为管线弯矩图。
具体实施方式
下面结合附图对本发明作进一步描述。
参照图1~图6,一种基于上覆荷载作用下埋地管线下方土体流失对管线受力形变的确定方法,包括以下步骤:
1)根据埋地管线从受力平衡到失衡的受力特征,将管线和土体抽离出来构建管-土耦合模型进行分析;
2)运用角点法求得地表荷载作用下管线的附加应力σz’
3)地基土假定为各向同性、均质线性变形体,不考虑土体和边界效应的影响,运用自重应力公式计算土体作用于管线的有效自重应力σc
4)考虑到管线直径比较小,管顶和管底的受力基本相同,梁的长与管线相同,宽为管线的直径,作用在埋地管线上的上覆附加荷载为:q(x)=(σz’c)D;
5)埋地管线埋置于土体中的情况比较复杂,受到的影响因素比较多,将埋地管线的上覆附加荷载假定为q(x)=ax2+bx+c,当管线的附加荷载系数a,b为0时,可求得q(x)=c,能够反映步骤4)中管线的受力情况;
6)根据埋地管线在土体中承载环境的差别,将管线形变段划分为管土分离区的变形和管土半分离区的变形;
7)通过管土分离区的受力与变形方程计算出管土分离区管线的位移y1、剪力Q1和弯矩M1
8)对于管土半分离区,将管线视为放置于半无限弹性地基上的半无限梁,运用Winkler地基梁理论推导出管线变形方程计算管土半分离区的位移y2、剪力Q2和弯矩M2
9)在算例的基础上,分析管线的位移、剪力以及弯矩的变化情况。
实例:一种基于上覆荷载作用下埋地管线下方土体流失对管线受力形变的确定方法,包括以下步骤:
1)通过管-土耦合实例,土体的长宽高分别取5m、1m、2.4m,管线长度为5m,管径D=0.12m,管壁厚度为0.002m,管线距离土体表面h=0.5m,管土分离区长度l1=0.373m。以软土为例,对土体参数进行取值,密度ρ=1400kg/m3,弹性模量E=5MPa,泊松比υ=0.3,抗弯刚度EI=1032N·m2。管线密度ρ=800kg/m3,弹性模量E=800MPa,泊松比υ=0.45。地表荷载位于管线正上方,长度l和宽度b均为0.4m,荷载大小为7.620KPa的均布荷载,计算模型如图1所示;
2)运用角点法求得地表荷载作用下管线的附加应力σz’=1.06KPa;
3)地基土假定为各向同性、均质线性变形体,不考虑土体和边界效应的影响,运用自重应力公式计算土体作用于管线的有效自重应力σc=7KPa;
4)考虑到管线直径比较小,管顶和管底的受力基本相同,梁的长与管线相同,宽为管线的直径,作用在埋地管线上的上覆附加荷载为:q(x)=(σz’c)D;
5)埋地管线埋置于土体中的情况比较复杂,受到的影响因素比较多,将埋地管线的上覆附加荷载假定为q(x)=ax2+bx+c,当管线的附加荷载系数a,b为0时,可求得q(x)=c,能够反映步骤4)中管线的受力情况;
6)根据埋地管线在土体中承载环境的差别,将管线形变段划分为管土分离区的变形和管土半分离区的变形;
7)通过管土分离区的受力与变形方程计算出管土分离区管线的位移y1、剪力Q1和弯矩M1为:
Figure BDA0002267215820000061
Figure BDA0002267215820000062
Q1=0.157-cx1
8)对于管土半分离区,将管线视为放置于半无限弹性地基上的半无限梁,运用Winkler地基梁理论推导出管线变形方程计算管土半分离区的位移y2、剪力Q2和弯矩M2为:
Figure BDA0002267215820000063
Figure BDA0002267215820000064
Figure BDA0002267215820000065
9)在算例的基础上,分析管线的位移、剪力以及弯矩的变化情况,如图3~图6所示。

Claims (5)

1.一种基于上覆荷载作用下埋地管线下方土体流失对管线受力形变的确定方法,其特征在于,所述方法包括以下步骤:
1)根据埋地管线从受力平衡到失衡的受力特征,将管线和土体抽离出来构建管-土耦合模型进行分析;
2)运用角点法求得地表荷载作用下管线的附加应力σz’
3)地基土假定为各向同性、均质线性变形体,不考虑土体和边界效应的影响,运用自重应力公式计算土体作用于管线的有效自重应力σc
4)考虑到管线直径比较小,管顶和管底的受力基本相同,梁的长与管线相同,宽为管线的直径,作用在埋地管线上的上覆附加荷载为:q(x)=(σz’c)D;
5)埋地管线埋置于土体中的情况比较复杂,受到的影响因素比较多,将埋地管线的上覆附加荷载假定为q(x)=ax2+bx+c,当管线的附加荷载系数a,b为0时,可求得q(x)=c,能够反映步骤4)中管线的受力情况;
6)根据埋地管线在土体中承载环境的差别,将管线形变段划分为管土分离区的变形和管土半分离区的变形;
7)通过管土分离区的受力与变形方程计算出管土分离区管线的位移y1、剪力Q1和弯矩M1
8)对于管土半分离区,将管线视为放置于半无限弹性地基上的半无限梁,运用Winkler地基梁理论推导出管线变形方程计算管土半分离区的位移y2、剪力Q2和弯矩M2
9)在算例的基础上,分析管线的位移、剪力以及弯矩的变化情况。
2.如权利要求1所述的基于上覆荷载作用下埋地管线下方土体流失对管线受力形变的确定方法,其特征在于,所述步骤2)中,管线的附加应力方程形式表达为:
Figure FDA0002267215810000011
其中,σz’为平行于z’坐标轴的应力,单位:kPa;p为作用于地表的竖向均布荷载,单位:kPa;l为地表荷载的长度;b为地表荷载的宽度;z’为地表荷载距离埋地管线的深度。
3.如权利要求1或2所述的基于上覆荷载作用下埋地管线下方土体流失对管线受力形变的确定方法,其特征在于,所述步骤3)中,管线的有效自重应力方程形式表达为:
Figure FDA0002267215810000012
其中,σc为天然地面下任意深度z’处的竖向有效自重应力,单位:kPa;n为深度z’范围内的土层总数;hi为第i层土的厚度,单位:m;γi为第i层土的天然重度,单位:kN/m3
4.如权利要求1或2所述的基于上覆荷载作用下埋地管线下方土体流失对管线受力形变的确定方法,其特征在于,所述步骤5)中,管土分离区的受力与变形方程形式表达为:
Figure FDA0002267215810000013
Figure FDA0002267215810000021
Figure FDA0002267215810000022
Figure FDA0002267215810000023
其中,y1为管土分离区管线的位移;M1为管土分离区管线的弯矩;Q1为管土分离区管线的剪力;E为管线的弹性模量;I为埋地管线的截面惯性矩;a,b,c为管线附加荷载系数;β为特征系数;yP为管线在P4点受到的位移,单位:m;MP为管线在P4点受到的弯矩,单位;kN·m;QP为管线在P4点受到的剪力,单位:kN。
5.如权利要求1或2所述的基于上覆荷载作用下埋地管线下方土体流失对管线受力形变的确定方法,其特征在于,所述步骤6)中,管土分离区的受力与变形方程形式表达为:
Figure FDA0002267215810000024
Figure FDA0002267215810000025
Figure FDA0002267215810000026
其中,y2为管土半分离区管线的位移;M2为管土半分离区管线的弯矩;Q2为管土半分离区管线的剪力;k为地基弹性阻力系数。
CN201911092349.1A 2019-11-11 2019-11-11 一种埋地管线下方土体流失对管线受力形变的确定方法 Active CN110991009B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911092349.1A CN110991009B (zh) 2019-11-11 2019-11-11 一种埋地管线下方土体流失对管线受力形变的确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911092349.1A CN110991009B (zh) 2019-11-11 2019-11-11 一种埋地管线下方土体流失对管线受力形变的确定方法

Publications (2)

Publication Number Publication Date
CN110991009A true CN110991009A (zh) 2020-04-10
CN110991009B CN110991009B (zh) 2023-05-23

Family

ID=70083683

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911092349.1A Active CN110991009B (zh) 2019-11-11 2019-11-11 一种埋地管线下方土体流失对管线受力形变的确定方法

Country Status (1)

Country Link
CN (1) CN110991009B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112597681A (zh) * 2020-12-22 2021-04-02 上海勘察设计研究院(集团)有限公司 一种大面积堆载引起的水平位移计算方法及计算系统
CN115435743A (zh) * 2022-08-24 2022-12-06 中山大学 一种管道形变监测方法及装置

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105184102A (zh) * 2015-10-10 2015-12-23 天津大学 一种管土动态耦合作用分析方法
CN106250618A (zh) * 2016-07-28 2016-12-21 浙江大学城市学院 地面出入式盾构施工引起地下管线弯矩及应变的计算方法
JP2017053725A (ja) * 2015-09-09 2017-03-16 株式会社クボタ 断層横断埋設管路の挙動推定方法及び断層横断埋設管路の挙動推定装置
CN108362583A (zh) * 2017-12-28 2018-08-03 武汉科技大学 一种模拟爆破荷载对临近埋地油气管道影响的试验装置及方法
CN109543286A (zh) * 2018-11-20 2019-03-29 中南大学 隧道开挖及爆破荷载作用下邻近埋地管道的力学分析模型及振速控制标准确定方法
CN110008613A (zh) * 2019-04-12 2019-07-12 交通运输部天津水运工程科学研究所 一种网格式软土地基承载特性计算方法
CN110135125A (zh) * 2019-01-30 2019-08-16 天津大学前沿技术研究院 一种水平力下钻井船插桩对邻近群桩影响的数值分析方法
CN110147587A (zh) * 2019-04-29 2019-08-20 浙江大学城市学院 一种隧道内集中荷载引起既有隧道位移变形的计算方法
CN110263448A (zh) * 2019-06-24 2019-09-20 广西路桥工程集团有限公司 一种双隧道开挖对邻近地理管线影响的预测方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017053725A (ja) * 2015-09-09 2017-03-16 株式会社クボタ 断層横断埋設管路の挙動推定方法及び断層横断埋設管路の挙動推定装置
CN105184102A (zh) * 2015-10-10 2015-12-23 天津大学 一种管土动态耦合作用分析方法
CN106250618A (zh) * 2016-07-28 2016-12-21 浙江大学城市学院 地面出入式盾构施工引起地下管线弯矩及应变的计算方法
CN108362583A (zh) * 2017-12-28 2018-08-03 武汉科技大学 一种模拟爆破荷载对临近埋地油气管道影响的试验装置及方法
CN109543286A (zh) * 2018-11-20 2019-03-29 中南大学 隧道开挖及爆破荷载作用下邻近埋地管道的力学分析模型及振速控制标准确定方法
CN110135125A (zh) * 2019-01-30 2019-08-16 天津大学前沿技术研究院 一种水平力下钻井船插桩对邻近群桩影响的数值分析方法
CN110008613A (zh) * 2019-04-12 2019-07-12 交通运输部天津水运工程科学研究所 一种网格式软土地基承载特性计算方法
CN110147587A (zh) * 2019-04-29 2019-08-20 浙江大学城市学院 一种隧道内集中荷载引起既有隧道位移变形的计算方法
CN110263448A (zh) * 2019-06-24 2019-09-20 广西路桥工程集团有限公司 一种双隧道开挖对邻近地理管线影响的预测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
雍睿: "基于概率统计的结构面粗糙度各向异性和尺寸效应分析" *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112597681A (zh) * 2020-12-22 2021-04-02 上海勘察设计研究院(集团)有限公司 一种大面积堆载引起的水平位移计算方法及计算系统
CN112597681B (zh) * 2020-12-22 2022-11-08 上海勘察设计研究院(集团)有限公司 一种大面积堆载引起的水平位移计算方法及计算系统
CN115435743A (zh) * 2022-08-24 2022-12-06 中山大学 一种管道形变监测方法及装置
CN115435743B (zh) * 2022-08-24 2024-03-08 中山大学 一种管道形变监测方法及装置

Also Published As

Publication number Publication date
CN110991009B (zh) 2023-05-23

Similar Documents

Publication Publication Date Title
Livneh et al. Axial testing and numerical modeling of square shaft helical piles under compressive and tensile loading
CN111058855B (zh) 一种盾构下穿的结构物的形变控制方法及评估系统
Rankin Ground movements resulting from urban tunnelling: predictions and effects
CN103243747B (zh) 预应力混凝土管桩桩基的纠偏方法
CN101251498B (zh) 一种基于电磁辐射原理的围岩松动圈测试及评价方法
CN102221332A (zh) 松散地层的光纤光栅多点传感装置及监测系统与监测方法
CN104727828A (zh) 基于变形控制的隧道超前支护管理方法
CN105220718A (zh) 一种建筑施工对邻近既有隧道的影响预测定方法
CN110046470A (zh) 一种盾构掘进引起工后地表沉降量的确定方法
CN108842821A (zh) 一种钻爆法修建海底隧道合理埋深的计算方法
CN103437318A (zh) 悬浮式海堤的工后沉降计算预测方法
CN110991009A (zh) 一种基于上覆荷载作用下埋地管线下方土体流失对管线受力形变的确定方法
CN110245426B (zh) 一种用于管廊结构顶管施工的有限元精细化模拟方法
Miśkiewicz et al. Monitoring system of the road embankment
Kunecki Field Test and Three-Dimensional Numerical Analysis of Soil–Steel Tunnel during Backfilling
Mao et al. Water disaster susceptible areas in loess multi-arch tunnel construction under the lateral recharge condition
CN106682267B (zh) 一种潮汐荷载下沉管隧道管节环向应变的计算方法
CN106096162A (zh) 一种确定盾构土仓压力的方法及其数学模型和构建方法
Tappenden Predicting the axial capacity of screw piles installed in Western Canadian soils
Chaiyaput et al. Evaluation of the road failure using resistivity and screw driving sounding testing techniques: A case study in Ang Thong province, Thailand
Parker et al. Field performance and analysis of 3-m-diameter induced trench culvert under a 19.4-m soil cover
Pathak Analysis of static lateral load test of battered pile group at I-10 twin span bridge
CN112307606A (zh) 浅埋公路隧道二次衬砌施做时间的确定方法
Liedberg Load reduction on a rigid pipe: pilot study of a soft cushion installation
Sun et al. Numerical simulation of a deep excavation near a shield tunnel

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