CN106096149A - 一种微观路表形貌的数字化生成方法 - Google Patents

一种微观路表形貌的数字化生成方法 Download PDF

Info

Publication number
CN106096149A
CN106096149A CN201610427759.7A CN201610427759A CN106096149A CN 106096149 A CN106096149 A CN 106096149A CN 201610427759 A CN201610427759 A CN 201610427759A CN 106096149 A CN106096149 A CN 106096149A
Authority
CN
China
Prior art keywords
dimensional
table pattern
microcosmic
road
point cloud
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.)
Pending
Application number
CN201610427759.7A
Other languages
English (en)
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.)
Southeast University
Original Assignee
Southeast 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 Southeast University filed Critical Southeast University
Priority to CN201610427759.7A priority Critical patent/CN106096149A/zh
Publication of CN106096149A publication Critical patent/CN106096149A/zh
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Tires In General (AREA)

Abstract

本发明公开了一种微观路表形貌的数字化生成方法，包括如下步骤：采用三维扫描仪扫描路面表面，获取路表形貌的三维点云坐标（X0,Y0,Z0）；对三维点云平面坐标（X0,Y0）进行栅格化，获得栅格格式的三维点云平面坐标（X1,Y1）；采用散点插值方法，获得与点云Z0坐标对应的插值Z1坐标；创建自由格式的多边形网格，获得三维点云坐标（X1,Y1,Z1）对应的三维网格；转换三维网格为三维多重曲面；导入有限元软件，获得用于数值模拟的微观路表形貌。本发明为基于微观路表形貌的轮胎/路面抗滑性能分析和轮胎/路面噪声性能分析提供了有利条件，为路表形貌的微观调控提供了技术可能性。

Description

一种微观路表形貌的数字化生成方法
技术领域
本发明涉及道路路面，特别是微观路表形貌的数字化生成。
背景技术
轮胎/路面抗滑性能与轮胎/路面噪声性能均是路面所应提供的两项重要功能。为进行这两项性能的深入研究，一方面依赖于对轮胎的精确建模(如轮胎组成、花纹等)，另一方面依赖于对微观路表形貌的精确表征。
目前，对路表形貌，大多仍基于一维维度上的考虑，如路表纹理波长、平均纹理深度(MTD)、平均断面深度(MPD)等；也有基于二维维度上的考虑，如研究随机路表不平度对噪声、抗滑性能的影响。
申请人研究后发现：这些参数不能完整地描述路表的微观形貌特征(导致形貌信息的丢失)，更为重要的是，这些参数无法满足抗滑和噪声研究对微观路表形貌的精度要求，也无法为路表形貌的微观调控(满足抗滑的条件下，降低路面噪声)创造良好的技术条件。
因而，迫切需要微观路表形貌的数字化生成方法，以满足抗滑和噪声研究的精度要求。
发明内容
发明目的：为了克服现有技术中存在的不足，本发明提供一种微观路表形貌的数字化生成方法，以满足抗滑和噪声研究对微观路表形貌的精度要求。
技术方案：为实现上述目的，本发明利用路表形貌的三维扫描信息，通过插值和网格化，获得用于数值模拟的数字化微观路表形貌。
本发明采用的具体的技术方案为：
一种微观路表形貌的数字化生成方法，包括顺序执行的以下步骤：
步骤1、扫描路面，获取表征路表形貌的三维点云坐标(X0,Y0,Z0)；
步骤2、对三维点云平面坐标(X0,Y0)进行栅格化，获得栅格格式的三维点云平面坐标(X1,Y1)；
步骤3、采用散点插值方法，获得与点云Z0坐标对应的插值Z1坐标；
步骤4、创建自由格式的多边形网格，获得三维点云坐标(X1,Y1,Z1)对应的三维网格；
步骤5、将三维网格转换为三维多重曲面；
步骤6、将数据导入有限元软件，获得用于数值模拟的数字化微观路表形貌。
进一步的，在本发明中，步骤2中，路表形貌格栅格式的三维点云平面坐标(X1,Y1)按照如下步骤获得：将路表形貌三维点云读入科学计算软件，根据每个点云平面坐标(X0,Y0)的上下限值和期望的间隔数(即所需的精度)，采用网格采样点生成函数，获得栅格格式的三维点云平面坐标(X1,Y1)。
进一步的，在本发明中，步骤3中插值Z1坐标按照如下步骤获得：在科学计算软件中，采用散点插值函数和散点插值方法，获得基于栅格格式点云平面坐标(X1,Y1)的、与点云Z0坐标对应的插值Z1坐标。
进一步的，在本发明中，步骤4中的三维点云坐标(X1,Y1,Z1)对应的三维网格按照如下步骤获得：在CAD软件中，采用三维网格命令，创建基于格栅格式的点云坐标(X1,Y1,Z1)的、自由格式的多边形三维网格。
进一步的，在本发明中，步骤5中的三维多重曲面按照如下步骤获得：在三维建模软件中，采用网格转换方法，将三维网格转换为三维多重曲面，保存为有限元软件支持的多种文件格式。
有益效果：微观路表形貌的数字化生成，为基于微观路表形貌的轮胎/路面抗滑性能分析和轮胎/路面噪声性能分析提供了有利条件，为路表形貌的微观调控提供了技术可能性。
附图说明
图1是本发明实施例的选定区域的路表形貌图；
图2是本发明实施例的用于数值模拟的数字化二维微观路表形貌图；
图3是本发明实施例的用于数值模拟的数字化三维微观路表形貌图。
具体实施方式
下面结合附图对本发明作更进一步的说明。
如图1所示，在已有的路面表面(混合料类型为SMA13)上，选定15cm×10cm的形貌测试区域，喷洒一层反光增强剂；采用三维扫描仪扫描该形貌测试区域，获取路表形貌的三维点云坐标(X0,Y0,Z0)，保存为*.asc文件。
将上述路表形貌的三维点云坐标文件(*.asc)读入科学计算软件MATLAB中，根据每个点云平面坐标(X0,Y0)的上下限值和期望的间隔数(本实施例中X方向间隔数为219，Y方向间隔数为134)，采用网格采样点生成函数meshgrid，获得栅格格式的三维点云平面坐标(X1,Y1)。
在科学计算软件MATLAB中，采用散点插值函数TriScatteredInterp和线性插值方法，获得基于栅格格式点云平面坐标(X1,Y1)的、与点云Z0坐标对应的插值Z1坐标，按照列向量形式，将全部三维点云坐标(X1,Y1,Z1)保存为*.txt文件。
在CAD软件AutoCAD中，采用三维网格命令3dmesh，输入M方向上的网格数量219，输入N方向上的网格数量134，再依次输入*.txt文件中的网格点的三维坐标(X1,Y1,Z1)(共29346个)，创建基于格栅格式的点云坐标(X1,Y1,Z1)的、自由格式的多边形三维网格，保存为*.dwg文件。
在三维建模软件Rhinoceros中，打开*.dwg文件，采用网格转换方法MeshToNurb，将三维网格转换为三维多重曲面，保存为*.sat文件。
在有限元软件Abaqus中，导入*.sat文件，获得如图2、图3所示的数字化二维与三维微观路表形貌图。
采用标准配置的计算机实施上述方案，所需的时间如下：
在上述实施例中，第1步骤采用扫描仪扫描试件表面，获取三维点云坐标(X0,Y0,Z0)，这个过程耗时约为5min。
第2和第3步骤，采用MATLAB软件进行点云坐标栅格化和散点插值，获得栅格点云坐标(X1,Y1,Z1)，耗时不超过1min。
第4步骤中，采用AutoCAD创建三维点云的三维网格，只需要1个命令即可完成，完成时间不超过1min。
第5步骤，采用Rhinoceros进行三维网格到三维多重曲面的转换，也只需要1个命令即可完成，该步骤约需1min。
第6步骤中，将三维多重曲面导入有限元中，约需2min。所有上述6个步骤，需要的总时间约为：5+1+1+1+2＝10min。
总体而言，本专利申请对计算机的要求不高(标准配置的计算机即可)，资源消耗不大，处理速度快。
在研究过程中，申请人发现：路表的形貌具有独特的结构，如何快速地获得仿真度较高的数字化的三维模型是一个具有挑战性的问题。为此，申请人对各种三维模型构建方法进行了认真的研究，分析各算法或方法的优缺点，设计了上述实施方案。需要说明的是，上述实施例是采用各软件中的相关模块来实现的，也可以根据需要，自行设计专用的软件来实现上述过程。采用上述软件，免去了软件开发的时间和费用，使用方便、成本低。技术人员可以根据本发明公开的思路和原理对上述实施例进行改进，通过一个软件来实现实现上述过程。
需要说明的是，本发明的设计构思不仅能够应用于微观道路的三维结构重建，也可以应用于其他三维景物场景，例如宏观道路三维重建，以解决现有三维景物重建中存在的耗时耗力、效率低下和仿真度不高的问题。
在一些微观三维场景重建的方案中，需要针对三维扫描仪的特点设计专门的重建算法，开发成本高，而且通用性不好。同样，现有的三维扫描获取数据之后，通过MATLAB软件进行三维重建，其不能反映研究对象的真实状态，仿真度不高。而且，如何对MATLAB重建后的三维数据进行深度处理，当前并没有很好的方案。而采用本发明的技术方案可以解决上述问题，提高仿真度和处理速度，从而更高更好地模拟相关场景，以更准确地研究相关问题。
本发明为基于微观路表形貌的轮胎/路面抗滑性能分析和轮胎/路面噪声性能分析提供了有利条件，为路表形貌的微观调控(通过调整混合料级配来实现)提供了技术可能性。以上详细描述了本发明的优选实施方式，但是，本发明并不限于上述实施方式中的具体细节，在本发明的技术构思范围内，可以对本发明的技术方案进行多种等同变换，这些等同变换均属于本发明的保护范围。

Claims (4)

1.一种微观路表形貌的数字化生成方法，其特征在于：包括以下步骤：
步骤1、扫描路面，获取表征路表形貌的三维点云坐标（X0,Y0,Z0）；
步骤2、对三维点云平面坐标（X0,Y0）进行栅格化，获得栅格格式的三维点云平面坐标（X1,Y1）；
步骤3、采用散点插值方法，获得与点云Z0坐标对应的插值Z1坐标；
步骤4、创建自由格式的多边形网格，获得三维点云坐标（X1,Y1,Z1）对应的三维网格；
步骤5、将三维网格转换为三维多重曲面；
步骤6、将数据导入有限元软件，获得用于数值模拟的微观路表形貌。
2.如权利要求1所述的微观路表形貌的数字化生成方法，其特征在于：步骤2中，按所需精度对三维点云平面坐标（X0,Y0）进行栅格化。
3.如权利要求1所述的微观路表形貌的数字化生成方法，其特征在于：步骤3中，所述散点插值方法包括线性散点插值方法、自然临近散点插值方法和最临近散点插值方法。
4.如权利要求1所述的微观路表形貌的数字化生成方法，其特征在于：步骤5进一步为，采用网格转换方法将转换三维网格为三维多重曲面，将其保存为有限元软件支持的文件格式。
CN201610427759.7A 2016-06-15 2016-06-15 一种微观路表形貌的数字化生成方法 Pending CN106096149A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610427759.7A CN106096149A (zh) 2016-06-15 2016-06-15 一种微观路表形貌的数字化生成方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610427759.7A CN106096149A (zh) 2016-06-15 2016-06-15 一种微观路表形貌的数字化生成方法

Publications (1)

Publication Number Publication Date
CN106096149A true CN106096149A (zh) 2016-11-09

Family

ID=57235825

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610427759.7A Pending CN106096149A (zh) 2016-06-15 2016-06-15 一种微观路表形貌的数字化生成方法

Country Status (1)

Country Link
CN (1) CN106096149A (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107705363A (zh) * 2017-10-20 2018-02-16 北京世纪高通科技有限公司 一种道路三维可视化建模方法及装置
CN107818208A (zh) * 2017-10-25 2018-03-20 东南大学 一种路表激励的获取方法
CN108763625A (zh) * 2018-04-12 2018-11-06 广东省水利水电科学研究院 一种flac3d软件前处理的方法及系统
CN109003332A (zh) * 2018-06-25 2018-12-14 重庆交通大学 沥青路面表面纹理仿真系统及其仿真方法
CN112613218A (zh) * 2020-12-26 2021-04-06 北京工业大学 一种基于abaqus的真实粗糙表面建模方法
CN115825410A (zh) * 2022-11-22 2023-03-21 重庆大学 基于路面构造预估轮胎路面作用噪声水平的方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101825444A (zh) * 2010-04-09 2010-09-08 上海辉格科技发展有限公司 基于面结构光的车载路谱测试系统
CN104282040A (zh) * 2014-09-29 2015-01-14 北京航空航天大学 一种用于三维实体模型重建的有限元前处理方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101825444A (zh) * 2010-04-09 2010-09-08 上海辉格科技发展有限公司 基于面结构光的车载路谱测试系统
CN104282040A (zh) * 2014-09-29 2015-01-14 北京航空航天大学 一种用于三维实体模型重建的有限元前处理方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
KP KUMARA: "Reconstructing Solid Model from 2D Scanned Images of Biological Organs for Finite Element Simulation", 《CENTRE FOR PRODUCT DESIGN & MANUFACTURING》 *
徐新泉: "路面形貌表征及固体颗粒对轮胎湿附着性能的影响", 《中国优秀硕士学位论文全文数据库（工程科技II辑）》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107705363A (zh) * 2017-10-20 2018-02-16 北京世纪高通科技有限公司 一种道路三维可视化建模方法及装置
CN107705363B (zh) * 2017-10-20 2021-02-23 北京世纪高通科技有限公司 一种道路三维可视化建模方法及装置
CN107818208A (zh) * 2017-10-25 2018-03-20 东南大学 一种路表激励的获取方法
CN107818208B (zh) * 2017-10-25 2021-06-01 东南大学 一种路表激励的获取方法
CN108763625A (zh) * 2018-04-12 2018-11-06 广东省水利水电科学研究院 一种flac3d软件前处理的方法及系统
CN109003332A (zh) * 2018-06-25 2018-12-14 重庆交通大学 沥青路面表面纹理仿真系统及其仿真方法
CN112613218A (zh) * 2020-12-26 2021-04-06 北京工业大学 一种基于abaqus的真实粗糙表面建模方法
CN112613218B (zh) * 2020-12-26 2024-03-15 北京工业大学 一种基于abaqus的真实粗糙表面建模方法
CN115825410A (zh) * 2022-11-22 2023-03-21 重庆大学 基于路面构造预估轮胎路面作用噪声水平的方法
CN115825410B (zh) * 2022-11-22 2023-12-22 重庆大学 基于路面构造预估轮胎路面作用噪声水平的方法

Similar Documents

Publication Publication Date Title
CN106096149A (zh) 一种微观路表形貌的数字化生成方法
WO2005010669A3 (en) Method for creating single 3d surface model from a point cloud
CN105354873A (zh) 用于多孔介质三维重构的模式密度函数模拟算法
CN105205861A (zh) 基于Sphere-Board的树木三维可视化模型实现方法
CN106127847B (zh) 一种生成实时山水画特征线条以及适当渲染的方法
CN109003322A (zh) 一种动力定位船舶海上作业的三维海浪仿真优化方法
CN113421338A (zh) 一种将二维矢量进行三维渲染的方法
JPWO2008032593A1 (ja) シボ柄印刷用のシボ柄並びにそのシボ柄の創成方法および創成プログラム、並びにそのシボ柄を印刷した住宅建材製品、自動車用内装部品、家庭電化製品および情報機器
CN112509105A (zh) 一种纹理连续的建筑物表面铺贴方法及装置
Graham et al. TravellerSim: growing settlement structures and territories with agent-based modeling
Talgorn et al. Real-time sketch-based terrain generation
Wang et al. An improving algorithm for generating real sense terrain and parameter analysis based on fractal
Steinø et al. Developing a parametric urban design tool. Some structural challenges and possible ways to overcome them
Koltsova et al. A case study of script-based techniques in urban planning
CN107767451A (zh) 一种基于gocad栅格模型的3dec复杂地形建模方法
CN106570934A (zh) 针对大场景的空间隐函数建模方法
CN102867325B (zh) 植物叶片病斑渲染方法
Carlson Multi-scale truchet patterns
Adams et al. Feature-based interactively sketched terrain
Backhaus Improved representation of topographic effects by a vertical adaptive grid in vector-ocean-model (VOM). Part I: Generation of adaptive grids
Wang et al. A study of leaf modeling technology based on morphological features
Fletcher et al. Challenges and perspectives of procedural modelling and effects
CN110750927A (zh) 基于有限元结果的连续体表面轻量化结构生成方法及系统
CN104732582B (zh) 基于gpu的三维海洋矢量场动态脉面可视化方法
Sakurai et al. Procedural Modeling of Leather Texture with Structural Elements.

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20161109

RJ01 Rejection of invention patent application after publication