CN109668960B - 一种点声源房屋群衰减计算方法 - Google Patents

一种点声源房屋群衰减计算方法 Download PDF

Info

Publication number
CN109668960B
CN109668960B CN201811408549.9A CN201811408549A CN109668960B CN 109668960 B CN109668960 B CN 109668960B CN 201811408549 A CN201811408549 A CN 201811408549A CN 109668960 B CN109668960 B CN 109668960B
Authority
CN
China
Prior art keywords
point
building
sound source
attenuation
sound
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
CN201811408549.9A
Other languages
English (en)
Other versions
CN109668960A (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.)
Sun Yat Sen University
Original Assignee
Sun Yat Sen 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 Sun Yat Sen University filed Critical Sun Yat Sen University
Priority to CN201811408549.9A priority Critical patent/CN109668960B/zh
Publication of CN109668960A publication Critical patent/CN109668960A/zh
Application granted granted Critical
Publication of CN109668960B publication Critical patent/CN109668960B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04Analysing solids
    • G01N29/11Analysing solids by measuring attenuation of acoustic waves
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/15Correlation function computation including computation of convolution operations
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/01Indexing codes associated with the measuring variable
    • G01N2291/015Attenuation, scattering

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Biochemistry (AREA)
  • Acoustics & Sound (AREA)
  • Immunology (AREA)
  • Computing Systems (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Algebra (AREA)
  • Pathology (AREA)
  • Chemical & Material Sciences (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Analytical Chemistry (AREA)
  • Building Environments (AREA)

Abstract

本发明提供一种点声源房屋群衰减计算方法,该方法从考虑点声源和接收点高度出发,利用声传播路线受遮挡段与完整长度的比值,进行房屋群衰减的计算;该方法能够应用于三维空间下的点声源房屋群衰减计算,有效的解决了工程应用中点声源衰减未考虑高度、难以计算、处理复杂和适用性差等问题。

Description

一种点声源房屋群衰减计算方法
技术领域
本发明涉及点声源声音排放和立体几何技术领域,更具体地,涉及一种点声源房屋群衰减计算方法。
背景技术
随着城市化进程的加快,噪声污染越来越严重,人们对于城市声环境质量也愈发的关注。城市内部存在着众多声源,如高音喇叭、道路交通流、工厂厂房等,这些声源产生的声音在城市内部传播,会遇到许多结构不一的障碍物,如何计算不同类型声源在城市内部的传播衰减量成为了学者们亟待解决的问题。实际生活中,城市房屋群的结构、表面材料、布局复杂多样,很难严格按照数学模型对声音传播的衰减量进行精确的计算。因此,在实际工程应用中,通常依据声音传播路线与房屋群的重叠部分,来进行衰减量的简化计算。
在GBT 17247.2-1998《声学户外声传播的衰减第2部分:一般计算方法》,对点声源房屋群衰减计算做出了规定,但标准中的点声源和噪声接收点均局限于平面,没有考虑到声源点和接收点的高度对衰减量带来的影响,无法计算前排建筑物对后排建筑物的遮挡效应,且标准中提到的沿声传播路线上的建筑物密度也难以计算,适用性较差。其他的点声源房屋群衰减计算模型也存在未考虑高度因素、难以应用等问题。对于在现场或用模型测量的方式,其测量点位有限,需要耗费大量的成本,且得到的结果难以用于噪声地图绘制或其他噪声模型的构建
发明内容
本发明提供一种三维空间下的考虑声源点和接收点高度、适用性强的点声源房屋群衰减计算方法。
为了达到上述技术效果,本发明的技术方案如下:
一种点声源房屋群衰减计算方法,包括以下步骤:
S1:获取房屋群在水平面上的投影轮廓面,将点声源和接收点的声传播路线投影到水平面上,筛选出与声传播路线相交的所有建筑物轮廓面;
S2:求出点声源与接收点的高度差,利用高度差与声传播路线长度的比值求出声传播路线与水平面的夹角;
S3:求出点声源与所有相交的轮廓面中所有边的最短距离;
S4:利用夹角与最短距离,求得声传播路线与建筑物最外侧表面的交点高度,若该高度高于建筑物的高度,则说明声传播路线没有受建筑物衰减,即无需计算衰减量;若该高度小于最外侧建筑物的高度,说明声传播路线受建筑物遮挡,需要计算衰减;
S5:利用相似三角形原理计算声传播路线在建筑物内部的传播距离与声传播路线总长度的比值,将比值乘以声传播路线总长度再乘以折减系数,得到点声源传播经过建筑物的衰减量。
进一步地,所述步骤S1中,与声传播路线相交是指建筑物轮廓面与声传播路线在水平面投影上相交。
进一步地,所述步骤S2中,点声源与接收点的高度差,无论为正值或负值,直接用于求声传播路线与水平面的夹角,具体为:
Figure GDA0002447586630000021
其中,θ表示为声传播路线与水平面的夹角,hp、hps分别表示接收点与点声源的高度,l表示声传播路线的长度。
进一步地,所述步骤S3中,点声源与所有相交轮廓面中所有边的最短距离,表示为求出点声源与所有边的距离集合,并从中找出最小值,具体为:
lmin=min{dis(ps,edge)}
其中,dis表示求解点与线段距离的函数,ps表示点声源,edge表示相交轮廓面中的所有边集合,lmin表示最短距离。
进一步地,所述步骤S4中,最外侧建筑物和建筑物最外侧表面分别指的是距离点声源最近边所属的建筑物及所在侧面,利用夹角和最短距离求声传播路线与最外侧表面的交点高度,具体为:
hc=lmin·tanθ+hps
其中,lmin表示最短距离,θ表示为声传播路线与水平面的夹角,hc、hps分别表示交点高度和点声源高度。
进一步地,所述步骤S4中,声传播路线受建筑物遮挡的情况分为两类情况,一类是传播路线由建筑物侧面穿入,由顶面穿出,另一类是传播路线由建筑物侧面穿入,由侧面穿出,这两类情况分别对应衰减类型为部分受房屋遮挡以及完全受房屋遮挡。
进一步地,所述步骤S5中,对于部分受房屋遮挡的情况而言,点声源传播经过建筑物的衰减量,具体为:
Figure GDA0002447586630000031
其中,Atten表示衰减量,α表示折减系数,hb表示最外侧建筑物高度,hc表示交点高度,θ表示声传播路线与水平面的夹角,lxy表示声传播路线在水平面上投影的长度,l表示声传播路线的长度。
进一步地,所述步骤S5中,对于完全受房屋遮挡的情况而言,点声源传播经过建筑物的衰减量,具体为:
Figure GDA0002447586630000032
式中,Atten表示衰减量,α表示折减系数,lcross表示声传播路线与建筑物轮廓面在水平面上相交线段的长度,lxy表示声传播路线在水平面上投影的长度,l表示声传播路线的长度。
进一步地,所述步骤S5中,折减系数在计算城市内部点声源受房屋群衰减时取值为0.1。
与现有技术相比,本发明技术方案的有益效果是:
本发明方法从考虑点声源和接收点高度出发,利用声传播路线受遮挡段与完整长度的比值,进行房屋群衰减的计算;该方法能够应用于三维空间下的点声源房屋群衰减计算,有效的解决了工程应用中点声源衰减未考虑高度、难以计算、处理复杂和适用性差等问题。
附图说明
图1是本发明流程示意图;
图2是本发明中点声源声传播场景示意图;
图3是本发明中点声源声传播场景俯视图;
图4是声传播未受建筑物遮挡情况示意图;
图5是声传播部分受建筑物遮挡情况示意图;
图6是声传播完全受建筑物遮挡情况示意图。
具体实施方式
附图仅用于示例性说明,不能理解为对本专利的限制;
为了更好说明本实施例,附图某些部件会有省略、放大或缩小,并不代表实际产品的尺寸;
对于本领域技术人员来说,附图中某些公知结构及其说明可能省略是可以理解的。
下面结合附图和实施例对本发明的技术方案做进一步的说明。
实施例1
如图1所示,一种点声源房屋群衰减计算方法,通过以下步骤实现:
S1:获取房屋群在水平面上的投影轮廓面,将点声源和接收点的声传播路线投影到水平面上,筛选出与声传播路线相交的所有建筑物轮廓面;
与声传播路线相交的所有建筑物轮廓面是指在水平面投影上二者相交;
S2:求出点声源与接收点的高度差,利用高度差与声传播路线长度的比值求出声传播路线与水平面的夹角;
S3:求出点声源与所有相交的轮廓面中所有边的最短距离;
点声源与所有相交轮廓面中所有边的最短距离,可表示为求出点声源与所有边的距离集合,并从中找出最小值;
S4:利用夹角与最短距离,可求得声传播路线与建筑物最外侧表面的交点高度,若该高度高于建筑物的高度,则说明声传播路线没有受建筑物衰减,即无需计算衰减量;若该高度小于最外侧建筑物的高度,说明声传播路线受建筑物遮挡,需要计算衰减;
最外侧建筑物和建筑物最外侧表面分别指的是距离点声源最近边所属的建筑物及所在侧面;
声传播路线受建筑物遮挡的情况可分为两类情况,一类是传播路线由建筑物侧面穿入,由顶面穿出,另一类是传播路线由建筑物侧面穿入,由侧面穿出,这两类情况分别对应衰减类型为部分受房屋遮挡以及完全受房屋遮挡;
S5:利用相似三角形原理计算声传播路线在建筑物内部的传播距离与声传播路线总长度的比值,将比值乘以声传播路线总长度再乘以折减系数,得到点声源传播经过建筑物的衰减量。
实例1
本发明适用于三维空间下点声源房屋群衰减量的计算。声音衰减量可由声传播路线在建筑物内部的传播距离与声传播路线总长度的比值乘以折减系数得到。
首先,对当前计算场景中(如图2所示)所有建筑物进行遍历,获取所有建筑物在水平面上的投影轮廓面,计算声传播路线与建筑物轮廓面在水平面上是否相交,将相交的建筑物筛选出来,如图3所示;
求出点声源与筛选得到的相交建筑物水平面轮廓中的每一条边的垂直距离,找到最短垂直距离lmin对应的边,该边为离点声源最近的边,该边所在侧面为离点声源最近的面,lmin可按下式进行计算:
lmin=min{dis(ps,edge)}
式中,dis表示求解点与线段距离的函数,ps表示点声源,edge表示相交轮廓面中的所有边集合。
然后,转换观察视角,所选视角应为侧面且垂直于上一步求得的最短距离,如图4所示。利用点声源和接收点的高度差,求取声传播路线与水平面的夹角θ,具体可按下式进行计算:
Figure GDA0002447586630000051
式中,hp、hps分别表示接收点与点声源的高度,l表示声传播路线的长度。
然后,利用最短距离lmin和水平夹角θ计算声传播路线与建筑物交点的高度hc,具体如下式进行计算:
hc=lmin·tanθ+hps
式中,lmin表示最短距离,θ表示为声传播路线与水平面的夹角,hps表示点声源高度。
然后,若交点高度hc大于最外侧建筑物的高度hb,这说明声传播路线并没有穿过建筑物,即无需计算衰减量,该情况如图4所示。若交点高度hc小于最外侧建筑物的高度hb,则说明声传播路线受到了建筑物遮挡,具体可分为两类情况,一类是传播路线由建筑物侧面穿入,由顶面穿出,另一类是传播路线由建筑物侧面穿入,由侧面穿出,这两类情况分别对应衰减类型为部分受房屋遮挡以及完全受房屋遮挡;
对于部分受房屋遮挡的情况,如图5所示。利用相似三角形原理计算声传播路线在建筑物内部的传播距离与声传播路线总长度的比值,将比值乘以声传播路线总长度再乘以折减系数,得到点声源传播经过建筑物的衰减量Atten,具体图下式计算:
Figure GDA0002447586630000061
式中,α表示折减系数,在本例中取值0.1,hb表示最外侧建筑物高度,hc表示交点高度,θ表示声传播路线与水平面的夹角,lxy表示声传播路线在水平面上投影的长度,l表示声传播路线的长度。
对于完全受房屋遮挡的情况,如图6所示,其衰减量具体计算如下式所示:
Figure GDA0002447586630000062
式中,α表示折减系数,在本例中取值0.1,lcross表示声传播路线与建筑物轮廓面在水平面上相交线段的长度,lxy表示声传播路线在水平面上投影的长度,l表示声传播路线的长度。
相同或相似的标号对应相同或相似的部件;
附图中描述位置关系的用于仅用于示例性说明,不能理解为对本专利的限制;
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明权利要求的保护范围之内。

Claims (7)

1.一种点声源房屋群衰减计算方法,其特征在于,包括以下步骤:
S1:获取房屋群在水平面上的投影轮廓面,将点声源和接收点的声传播路线投影到水平面上,筛选出与声传播路线相交的所有建筑物轮廓面;
S2:求出点声源与接收点的高度差,利用高度差与声传播路线长度的比值求出声传播路线与水平面的夹角;
S3:求出点声源与所有相交的轮廓面中所有边的最短距离;
S4:利用夹角与最短距离,求得声传播路线与建筑物最外侧表面的交点高度,若该高度高于建筑物的高度,则说明声传播路线没有受建筑物衰减,即无需计算衰减量;若该高度小于最外侧建筑物的高度,说明声传播路线受建筑物遮挡,需要计算衰减;
S5:利用相似三角形原理计算声传播路线在建筑物内部的传播距离与声传播路线在水平面上投影的长度的比值,将比值乘以声传播路线总长度再乘以折减系数,得到点声源传播经过建筑物的衰减量;
所述步骤S5中,对于部分受房屋遮挡的情况而言,点声源传播经过建筑物的衰减量,具体为:
Figure FDA0002447586620000011
其中,Atten表示衰减量,α表示折减系数,hb表示最外侧建筑物高度,hc表示交点高度,θ表示声传播路线与水平面的夹角,lxy表示声传播路线在水平面上投影的长度,l表示声传播路线的长度;
所述步骤S5中,对于完全受房屋遮挡的情况而言,点声源传播经过建筑物的衰减量,具体为:
Figure FDA0002447586620000012
式中,Atten表示衰减量,α表示折减系数,lcross表示声传播路线与建筑物轮廓面在水平面上相交线段的长度,lxy表示声传播路线在水平面上投影的长度,l表示声传播路线的长度。
2.根据权利要求1所述的点声源房屋群衰减计算方法,其特征在于,所述步骤S1中,与声传播路线相交是指建筑物轮廓面与声传播路线在水平面投影上相交。
3.根据权利要求2所述的点声源房屋群衰减计算方法,其特征在于,所述步骤S2中,点声源与接收点的高度差,无论为正值或负值,直接用于求声传播路线与水平面的夹角,具体为:
Figure FDA0002447586620000021
其中,θ表示为声传播路线与水平面的夹角,hp、hps分别表示接收点与点声源的高度,l表示声传播路线的长度。
4.根据权利要求3所述的点声源房屋群衰减计算方法,其特征在于,所述步骤S3中,点声源与所有相交轮廓面中所有边的最短距离,表示为求出点声源与所有边的距离集合,并从中找出最小值,具体为:
lmin=min{dis(ps,edge)}
其中,dis表示求解点与线段距离的函数,ps表示点声源,edge表示相交轮廓面中的所有边集合,lmin表示最短距离。
5.根据权利要求4所述的点声源房屋群衰减计算方法,其特征在于,所述步骤S4中,最外侧建筑物和建筑物最外侧表面分别指的是距离点声源最近边所属的建筑物及所在侧面,利用夹角和最短距离求声传播路线与最外侧表面的交点高度,具体为:
hc=lmin·tanθ+hps
其中,lmin表示最短距离,θ表示为声传播路线与水平面的夹角,hc、hps分别表示交点高度和点声源高度。
6.根据权利要求5所述的点声源房屋群衰减计算方法,其特征在于,所述步骤S4中,声传播路线受建筑物遮挡的情况分为两类情况,一类是传播路线由建筑物侧面穿入,由顶面穿出,另一类是传播路线由建筑物侧面穿入,由侧面穿出,这两类情况分别对应衰减类型为部分受房屋遮挡以及完全受房屋遮挡。
7.根据权利要求6所述的点声源房屋群衰减计算方法,其特征在于,所述步骤S5中,折减系数在计算城市内部点声源受房屋群衰减时取值为0.1。
CN201811408549.9A 2018-11-23 2018-11-23 一种点声源房屋群衰减计算方法 Active CN109668960B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811408549.9A CN109668960B (zh) 2018-11-23 2018-11-23 一种点声源房屋群衰减计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811408549.9A CN109668960B (zh) 2018-11-23 2018-11-23 一种点声源房屋群衰减计算方法

Publications (2)

Publication Number Publication Date
CN109668960A CN109668960A (zh) 2019-04-23
CN109668960B true CN109668960B (zh) 2020-09-29

Family

ID=66142249

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811408549.9A Active CN109668960B (zh) 2018-11-23 2018-11-23 一种点声源房屋群衰减计算方法

Country Status (1)

Country Link
CN (1) CN109668960B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106546661A (zh) * 2016-10-24 2017-03-29 河海大学常州校区 一种干涉超声合成孔径成对混凝土内部成像方法
CN106652446A (zh) * 2016-11-30 2017-05-10 中山大学 一种基于离线‑在线模式的道路交通噪声动态模拟方法
CN108842642A (zh) * 2018-06-13 2018-11-20 中山大学 一种频率等效的声屏障的道路交通噪声衰减计算方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5166991B2 (ja) * 2008-06-25 2013-03-21 パナソニック株式会社 浮遊粒子測定システム
JP2010038696A (ja) * 2008-08-04 2010-02-18 Toyama Prefecture 金属疲労損傷度の非破壊評価方法と超音波金属疲労損傷度測定装置
US8847813B2 (en) * 2010-06-15 2014-09-30 Stolar Research Corporation Unsynchronized radio imaging
CN102622518B (zh) * 2012-03-08 2014-11-26 中山大学 一种基于建筑物群密度的室外声预测方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106546661A (zh) * 2016-10-24 2017-03-29 河海大学常州校区 一种干涉超声合成孔径成对混凝土内部成像方法
CN106652446A (zh) * 2016-11-30 2017-05-10 中山大学 一种基于离线‑在线模式的道路交通噪声动态模拟方法
CN108842642A (zh) * 2018-06-13 2018-11-20 中山大学 一种频率等效的声屏障的道路交通噪声衰减计算方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
《Areawide dynamic traffic noise simulation in urban built-up area using beam tracing approach》;Haibo Wang等;《Sustainable Cities and Society》;20171231;全文 *
《Comprehensive Evaluation Of Traffic Noise Pollution Based On Population Exposure》;XIE Lin-hua等;《Science Direct》;20131231;全文 *
《Road traffic noise mapping in Guangzhou using GIS and GPS》;Ming Cai等;《Applied Acoustics》;20151231;全文 *
《The 3D Attenuation Calculation of Traffic Noise among Building Groups by Using Beam Tracing Method》;Hai-bo Wang等;《Science Direct》;20131231;全文 *
《一种道路交通噪声建筑物群三维衰减计算方法》;王海波等;《环境科学与技术》;20130831;全文 *
《单体建筑在点声源环境中对预测点噪声等级的影响》;郝崑等;《沈阳航空工业学院学报》;20041231;全文 *

Also Published As

Publication number Publication date
CN109668960A (zh) 2019-04-23

Similar Documents

Publication Publication Date Title
AU2020101587A4 (en) Coupling method of 1D and 2D hydrodynamic models based on spatial topology
US11532123B2 (en) Method for visualizing large-scale point cloud based on normal
Hornikx et al. Noise abatement schemes for shielded canyons
CN102930509B (zh) 一种机载激光点云数据的智能化滤波方法
GB2456301A (en) Surface Representation from a Point Cloud
Jang et al. Scale-model method for measuring noise reduction in residential buildings by vegetation
Farcaş et al. Road traffic noise: GIS tools for noise mapping and a case study for Skåne region
WO2022083108A1 (zh) 一种施工质量检测方法、装置、终端和存储介质
CN109668960B (zh) 一种点声源房屋群衰减计算方法
CN106096321B (zh) 结合声线跟踪法与声束跟踪法的室内外声音传播模拟方法
Van Renterghem et al. Road traffic noise reduction by vegetated low noise barriers in urban streets
CN109286794A (zh) 一种基于相对空间位置的室内告警设备关系计算方法
CN110175345B (zh) 一种三维场景下机动车尾气扩散的模拟方法
Schröder et al. Real-time hybrid simulation method including edge diffraction
CN102622518B (zh) 一种基于建筑物群密度的室外声预测方法
CN110967737B (zh) 一种构造约束的初始模型构建方法
CN111859612B (zh) 激光雷达仿真方法以及装置
Suhanek et al. Combination of boundary element method and genetic algorithm for optimization of T-shape noise barrier
CN112700824A (zh) 纤维混凝土细观模型的建立方法
CN112668768A (zh) 基于导航点和rvo模型的人群疏散仿真方法
CN115205433B (zh) 基于glTF模型和建筑物轮廓拓展的三维渲染融合方法
CN106815430B (zh) 电磁环境仿真中建筑物模型构建方法及装置
CN110415360B (zh) 一种家装设计中地面区域的图形处理方法及系统
EP3937125B1 (en) Method, apparatus for superimposing laser point clouds and high-precision map and electronic device
CN114664106B (zh) 一种交叉口车道分配方法、装置、电子设备和存储介质

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
OL01 Intention to license declared