CN105374068B - 自由曲面网格结构的网格划分方法 - Google Patents

自由曲面网格结构的网格划分方法 Download PDF

Info

Publication number
CN105374068B
CN105374068B CN201510777920.9A CN201510777920A CN105374068B CN 105374068 B CN105374068 B CN 105374068B CN 201510777920 A CN201510777920 A CN 201510777920A CN 105374068 B CN105374068 B CN 105374068B
Authority
CN
China
Prior art keywords
grid
dimensional surfaces
expansion
curved surface
free form
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.)
Expired - Fee Related
Application number
CN201510777920.9A
Other languages
English (en)
Other versions
CN105374068A (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201510777920.9A priority Critical patent/CN105374068B/zh
Publication of CN105374068A publication Critical patent/CN105374068A/zh
Application granted granted Critical
Publication of CN105374068B publication Critical patent/CN105374068B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/30—Polynomial surface description

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Analysis (AREA)
  • Algebra (AREA)
  • Computer Graphics (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • Theoretical Computer Science (AREA)
  • Image Generation (AREA)

Abstract

本发明公开了一种自由曲面网格结构的网格划分方法。它是将自由曲面的参数域u、v向等分，在曲面上形成矩形点阵，将曲面的展开等效为点的映射展开，以参数域对角线所对应的空间曲面对角线为展开基线，按照展开前后面积近似不变的基本原则，对自由曲面由中心逐步向外展开；其次，展开的点阵拟合形成二维曲面，二维曲面从中心点等角度分成6部分，对每部分采用线推进方法逐层生成三角形网格，实现二维曲面的网格生成；最后，二维曲面的网格根据拓扑不变性，映射回自由曲面形成空间曲面网格。本发明能适应大多数曲面，生成的网格大小均匀、线条流畅，既满足了建筑美学的审美要求，又达到了传力连续的效果，能实现程序化，可快速实现网格的自动生成。

Description

自由曲面网格结构的网格划分方法
技术领域
本发明涉及空间结构网格划分方法，尤其涉及一种自由曲面网格结构的网格划分方法。
背景技术
随着计算机辅助技术的进步、建筑造型艺术的提高和人们审美水平的改变，新型自由的空间曲面结构越来越多地出现在大众的视野中，该类曲面结构具有较好的视觉表现力，但因为完全自由，所以不能用解析式精确表达。传统的杆件布置方案并不可行，因此，为了对新型自由曲面进行合理的杆件布置，首先需要对它进行网格划分，目前，自由曲面的网格生成尚未有成熟完美的方法，如何生成大小均匀、线条流畅的曲面网格是空间结构研究领域的一个热点和难点。
网格划分起源于上世纪60年代，最初应用于有限元领域，在这个过程中，波前法、Delaunay法、映射法等网格划分方法相继被提出并发展，但很多时候，现有的网格划分方法所生成的网格是非结构化网格，并不能满足建筑美观的要求，即网格大小均匀、杆件线条流畅等美学指标。近年来，有不少学者在这一领域进行了研究探索，并取得了一定的成果，提出了基于主应力轨迹线的波前推进法生成曲面网格；基于映射思想提出了自定义单元法来生成网格；采用等参线分割法和改进的Delaunay点云网格生成法进行曲面网格划分；基于曲面展开和Delaunay优化来生成网格。
发明内容
本发明的目的是克服现有技术的不足，提供一种自由曲面网格结构的网格划分方法。
一种自由曲面网格结构的网格划分方法是将自由曲面的参数域u、v向等分，在曲面上形成矩形点阵，将曲面的展开等效为点的映射展开，以参数域对角线所对应的空间曲面对角线为展开基线，按照展开前后面积近似不变的基本原则，对自由曲面由中心逐步向外展开；其次，展开的点阵拟合形成二维曲面，二维曲面从中心点等角度分成6部分，对每部分采用线推进方法逐层生成三角形网格，实现二维曲面的网格生成；最后，二维曲面的网格根据拓扑不变性，映射回自由曲面形成空间曲面网格。
所述的二维曲面从中心点等角度分成6部分，对每部分采用线推进方法逐层生成三角形网格，实现二维曲面的网格生成为：在展开的二维曲面上，取参数域中u=1/2，v=1/2为它的中心点，经过中心点的三条直线将二维曲面分成六部分，每部分的顶角均为60º，将每部分单独利用线推进方法进行二维曲面上三角形网格，根据外边界的情况不同，将曲面外边界分成八种情形。
所述的二维曲面的网格根据拓扑不变性，映射回自由曲面形成空间曲面网格为：将二维曲面上的三角形网格，作为网格结构的杆件、节点的信息，反向求解出每个节点在参数域中的u、v值，将此u、v值代入空间自由曲面非均匀有理B样条数学表达式，获得相对应的空间坐标，根据映射不改变网格及节点的拓扑关系，得到自由曲面的空间网格结果。
本发明与背景技术相比具有的有益效果是：
1) 本发明能适应大多数曲面，生成的网格大小均匀、线条流畅，既满足了建筑美学的审美要求，又达到了传力连续的效果。
2) 本发明能够实现程序化，可快速实现网格的自动生成。
附图说明
图1是参数域网格示意图；
图2是空间曲面离散网格示意图；
图3是二维曲面六等角分示意图；
图4外边界的八种情形示意图；
图5是算例1曲面示意图；
图6是算例1曲面的展开结果；
图7是算例1展开后的二维平面网格；
图8是算例1网格划分侧视图；
图9是算例1网格划分俯视图；
图10是算例2曲面示意图；
图11是算例2曲面的展开结果；
图12是算例2展开后的二维平面网格；
图13是算例2网格划分三视图。
具体实施方式
一种自由曲面网格结构的网格划分方法是将自由曲面的参数域u、v向等分，在曲面上形成矩形点阵，将曲面的展开等效为点的映射展开，以参数域对角线所对应的空间曲面对角线为展开基线，按照展开前后面积近似不变的基本原则，对自由曲面由中心逐步向外展开；其次，展开的点阵拟合形成二维曲面，二维曲面从中心点等角度分成6部分，对每部分采用线推进方法逐层生成三角形网格，实现二维曲面的网格生成；最后，二维曲面的网格根据拓扑不变性，映射回自由曲面形成空间曲面网格。
所述的二维曲面从中心点等角度分成6部分，对每部分采用线推进方法逐层生成三角形网格，实现二维曲面的网格生成为：在展开的二维曲面上，取参数域中u=1/2，v=1/2为它的中心点，经过中心点的三条直线将二维曲面分成六部分，每部分的顶角均为60º，将每部分单独利用线推进方法进行二维曲面上三角形网格，根据外边界的情况不同，将曲面外边界分成八种情形。
所述的二维曲面的网格根据拓扑不变性，映射回自由曲面形成空间曲面网格为：将二维曲面上的三角形网格，作为网格结构的杆件、节点的信息，反向求解出每个节点在参数域中的u、v值，将此u、v值代入空间自由曲面非均匀有理B样条数学表达式，获得相对应的空间坐标，根据映射不改变网格及节点的拓扑关系，得到自由曲面的空间网格结果。
算例1，曲面模型图5，高度约30m，上下圆弧半径分别约30m、10m。
首先，对自由曲面近似展开处理，展开方法为：将自由曲面参数域的两个方向分别等分，等分点映射回空间曲面以在曲面上形成矩形点阵，则曲面的展开等效为矩形点阵的映射展开，先以参数域对角线所对应的空间曲面对角线为展开基线，再按照展开前后面积近似不变的基本原则，对自由曲面由中心逐层向外展开。此种展开方式能很好地反映曲面的形状，走向及面积分布，图6为曲面的展开结果。
其次，展开的点阵通过拟合形成有边界的二维平面，将二维平面从中心点等角度分成6部分，考虑不同的边界情形，依次对每部分采用线推进方法逐层生成三角形网格，最终实现完整平面的网格生成。此种平面网格划分方式更具操作性、效果也更流畅，阳光谷曲面展开后的网格划分为图7。
最后，二维平面的网格根据拓扑不变性，映射回自由曲面形成空间曲面网格。阳光谷曲面映射后的空间网格为图8，图9，网格杆件长度约2m。
算例2，曲面两个方向长度分别约60m、40m，网格杆件长度约2m，按算例1同样方法得到图10~图13结果。

Claims (1)

1.一种自由曲面网格结构的网格划分方法，其特征在于，将自由曲面的参数域u、v向等分，在曲面上形成矩形点阵，将曲面的展开等效为点的映射展开，以参数域对角线所对应的空间曲面对角线为展开基线，按照展开前后面积近似不变的基本原则，对自由曲面由中心逐步向外展开；其次，展开的点阵拟合形成二维曲面，二维曲面从中心点等角度分成6部分，对每部分采用线推进方法逐层生成三角形网格，实现二维曲面的网格生成；最后，二维曲面的网格根据拓扑不变性，映射回自由曲面形成空间曲面网格；所述的二维曲面从中心点等角度分成6部分，对每部分采用线推进方法逐层生成三角形网格，实现二维曲面的网格生成为：在展开的二维曲面上，取参数域中u=1/2，v=1/2为它的中心点，经过中心点的三条直线将二维曲面分成六部分，每部分的顶角均为60 º，将每部分单独利用线推进方法进行二维曲面上三角形网格生成，根据外边界的情况不同，将曲面外边界分成八种情形；所述的二维曲面的网格根据拓扑不变性，映射回自由曲面形成空间曲面网格为：将二维曲面上的三角形网格，作为网格结构的杆件、节点的信息，反向求解出每个节点在参数域中的u、v值，将此u、v值代入空间自由曲面非均匀有理B样条数学表达式，获得相对应的空间坐标，根据映射不改变网格及节点的拓扑关系，得到自由曲面的空间网格结果。
CN201510777920.9A 2015-11-13 2015-11-13 自由曲面网格结构的网格划分方法 Expired - Fee Related CN105374068B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510777920.9A CN105374068B (zh) 2015-11-13 2015-11-13 自由曲面网格结构的网格划分方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510777920.9A CN105374068B (zh) 2015-11-13 2015-11-13 自由曲面网格结构的网格划分方法

Publications (2)

Publication Number Publication Date
CN105374068A CN105374068A (zh) 2016-03-02
CN105374068B true CN105374068B (zh) 2018-05-01

Family

ID=55376236

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510777920.9A Expired - Fee Related CN105374068B (zh) 2015-11-13 2015-11-13 自由曲面网格结构的网格划分方法

Country Status (1)

Country Link
CN (1) CN105374068B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111070721A (zh) * 2019-12-12 2020-04-28 航天特种材料及工艺技术研究所 一种柔性化曲面泡沫结构及其铺覆方法

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107274489A (zh) * 2017-05-09 2017-10-20 浙江大学 一种基于曲面展开的引导线法自由曲面网格划分方法
CN111474898B (zh) * 2020-04-17 2022-08-09 西北工业大学 一种自由曲面加工工艺参数优化方法
CN111880305B (zh) * 2020-08-14 2022-03-08 中国科学院国家天文台南京天文光学技术研究所 基于地图投影的大型光学望远镜非球面主镜划分的方法
CN112560190A (zh) * 2020-12-28 2021-03-26 中国航空工业集团公司沈阳飞机设计研究所 一种异形曲面功能结构降维设计方法
CN117253012B (zh) * 2023-09-18 2024-03-19 东南大学 一种还原平面建筑自由曲面网格结构至三维空间的方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101383047A (zh) * 2007-09-03 2009-03-11 鸿富锦精密工业（深圳）有限公司 曲面网格化方法
CN104331932A (zh) * 2014-10-21 2015-02-04 浙江大学 一种基于自定义单元法的自由曲面建筑网格划分方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101383047A (zh) * 2007-09-03 2009-03-11 鸿富锦精密工业（深圳）有限公司 曲面网格化方法
CN104331932A (zh) * 2014-10-21 2015-02-04 浙江大学 一种基于自定义单元法的自由曲面建筑网格划分方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Topological mesh for shell structures;Rafael Jurado-Pina et al;《Applied Mathematical Modelling》;20091231;第948-958页 *
自由曲面网格的划分与优化方法;危大结 等;《建筑结构》;20131031;第43卷(第19期);第50-52页第2-4节、图4、图6 *
自由曲面网格结构的网格划分技术研究;马腾;《中国优秀硕士学位论文全文数据库 信息科技辑》;20151015(第10期);第23-26页、图3.4 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111070721A (zh) * 2019-12-12 2020-04-28 航天特种材料及工艺技术研究所 一种柔性化曲面泡沫结构及其铺覆方法
CN111070721B (zh) * 2019-12-12 2022-03-15 航天特种材料及工艺技术研究所 一种柔性化曲面泡沫结构及其铺覆方法

Also Published As

Publication number Publication date
CN105374068A (zh) 2016-03-02

Similar Documents

Publication Publication Date Title
CN105374068B (zh) 自由曲面网格结构的网格划分方法
CN102254343B (zh) 一种基于凸包与obb的三维网格模型骨架提取方法
CN101639946A (zh) 植物叶片三维模型几何描述和曲面重建方法及系统
CN107886569B (zh) 一种基于离散李导数的测度可控的曲面参数化方法及系统
CN104715508A (zh) 一种网格建模方法及装置
CN105302927A (zh) 一种任意曲面建立空间网格结构的方法
CN108038259A (zh) 基于曲率生成气动部件外形的方法
CN103412985A (zh) 一种气冷叶片尾缘劈缝参数化设计方法
CN101877145A (zh) 基于质量因子的复杂平表面网格预剖分方法
CN101546351B (zh) 一种叶轮的变复杂度形状优化的几何参数化建模方法
CN112199777A (zh) 一种适用于模化仿生前缘流场特征的方法
CN116451540A (zh) 一种建筑表面施工曲面分格方法
CN116484467A (zh) 基于grasshopper的穹顶结构网格划分方法
CN103413019A (zh) 一种不规则轮廓双曲率外形面整体壁板离散方法
CN105046727A (zh) 一种产品效果图的生成方法
CN104851125B (zh) 一种植物叶子三维模型建模方法及系统
CN113158307A (zh) 一种自由曲面非奇异建筑网格划分方法
CN104537170A (zh) 用于泵站三维流场计算的网格划分方法及装置
CN101477707A (zh) 一种插值给定若干闭合曲线的曲面造型方法
Xiong et al. Automated structured all-quadrilateral and hexahedral meshing of tubular surfaces
CN109670197B (zh) 一种基于线框模型的工艺补充的创建和调整方法
CN107274489A (zh) 一种基于曲面展开的引导线法自由曲面网格划分方法
Wang et al. An extension of advancing front technique on new target surface after virtual topology operations
Wang et al. Automated hex meshing for turbomachinery secondary air system
CN107766673B (zh) 一种半径可控的参数化三维前缘钝化设计方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180501

Termination date: 20191113

CF01 Termination of patent right due to non-payment of annual fee