CN105631065B - 一种基于背景网格的动网格方法 - Google Patents
一种基于背景网格的动网格方法 Download PDFInfo
- Publication number
- CN105631065B CN105631065B CN201410601962.2A CN201410601962A CN105631065B CN 105631065 B CN105631065 B CN 105631065B CN 201410601962 A CN201410601962 A CN 201410601962A CN 105631065 B CN105631065 B CN 105631065B
- Authority
- CN
- China
- Prior art keywords
- mesh
- grid
- point
- bcmv
- coordinate information
- 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
Links
Landscapes
- Complex Calculations (AREA)
- Image Generation (AREA)
Abstract
本发明属于针对结构网格、结构/非结构混合网格以及点云的动网格技术领域,具体涉及一种基于背景网格的动网格技术。本发明采用径向基函数法对稀疏的背景网格进行变形,然后将背景网格的变形量插值到原始体网格的网格点上,从而达到提高动网格效率的目的。本发明的背景网格仅用于将变形量回插到原始体网格上,所以并不要求背景网格与原始体网格的边界吻合,以至于可以采用简单的可自动化生成的直角网格作为背景网格,大大地减少了人工干预,且变形后的网格质量仍然较高。
Description
技术领域
本发明属于针对结构网格、结构/非结构混合网格以及点云的动网格技术领域,具体涉及一种基于背景网格的动网格技术。
背景技术
动网格技术被广泛运用于气动外形优化设计、非定常流动仿真以及气动弹性力学仿真中。高效的动网格技术是解决此类问题的关键技术之一,而且也是实现工程应该中最重要的技术“瓶颈”。动网格技术主要有三种类型:刚性运动网格技术、重叠动网格技术、网格变形技术。本发明属于网格变形类的动网格技术。
目前,网格变形类的动网格技术一般分为结构网格变形技术和结构/非结构混合网格变形技术。常用的结构网格变形技术有:超限插值网格变形技术、径向基函数网格变形技术等。常用的结构/非结构混合网格变形技术有:弹簧拉伸网格变形技术、有限元网格变形技术、径向基函数网格变形技术等。超限插值网格变形技术属于代数方法,效率极高,但在网格拓扑结构较为复杂或者变形幅度较大时适用性较差,易出现网格交叉等问题。弹簧拉伸网格变形技术和有限元网格变形技术均需求解由全体网格节点形成的大型稀疏矩阵,此类方法在网格尺度较大、网格各向同性以及网格单元数较小时得到了广泛的应用,但应用于求解Navier-Stokes方程的各向异性的网格,则效率十分低下,且易出现网格交叉等情况。径向基函数网格变形技术是一类不依赖网格拓扑的网格变形方法,能够应用于任意网格类型的网格变形技术。其网格变形能力极强,当网格规模较大时,计算量十分巨大,目前已经提出了贪心算法来精简数据,大大的提高了计算的效率。
采用径向基函数网格变形技术时,如果体网格数达到百万甚至千万级时,将边界网格的变形通过径向基函数插值到体网格上时,计算量将十分庞大,难以在工程中应用。
因此,亟需研制一种基于背景网格的动网格技术,以有效的解决此问题。
发明内容
本发明要解决的技术问题是提出的基于背景网格的动网格技术,旨在大幅减小动网格技术中的计算量,提高网格生成的速度,提出一种能够满足工程需求的高效的动网格技术,突破气动外形优化设计、非定常流动仿真以及气动弹性力学仿真中的动网格技术“瓶颈”。
为了实现这一目的,本发明采取的技术方案是:
一种基于背景网格的动网格技术,包括如下步骤:
步骤1、采集体网格点坐标信息V(x,y,z)以及边界网格点坐标信息BC(x,y,z);
步骤2、根据步骤1所采集的体网格点坐标信息,创建一个完全覆盖体网格区域的长方体,并对此长方体进行网格剖分,完成背景网格Vbg(x,y,z)的生成;
步骤3、从步骤1所采集的边界网格点坐标信息中采集参与网格变形的网格点BCmv(x,y,z)old,并根据工程实际需求,给出变形后边界网格点的新坐标信息BCmv(x,y,z)new,然后计算边界网格点位置的改变量:
ΔBCmv(x,y,z)=BCmv(x,y,z)new-BCmv(x,y,z)old
步骤4、剔除步骤3中采集的BCmv(x,y,z)old中的重复点,并将步骤3中采集的BCmv(x,y,z)new以及ΔBCmv(x,y,z)中相应的网格点的信息剔除;
步骤5、根据步骤4中剔除重复点后的BCmv(x,y,z)old以及ΔBCmv(x,y,z),采用基于贪心算法的径向基函数法计算插值系数,从BCmv(x,y,z)old中获取经过贪心算法精简后的插值基函数支撑点集合U(P1,P2,…,Pn)以及相应的权系数ωi(i=1,2,…,n),其中n为集合U中的元素个数;
步骤6、根据步骤5所获得的支撑点集合U以及权系数ωi计算步骤2生成的背景网格Vbg(x,y,z)的变形量ΔVbg(x,y,z),并生成变形后的背景网格;
步骤7、根据步骤6获得的ΔVbg(x,y,z)插值出步骤1中体网格V(x,y,z)的变形量ΔV(x,y,z),本步骤可分为以下几步:
步骤7.1、在体网格V(x,y,z)选取任意网格点Pi(x,y,z),从步骤2的背景网格Vbg(x,y,z)中搜索出包含网格点Pi(x,y,z)的网格单元Ci;
网格单元Ci由网格点Pi1,Pi2,Pi3,Pi4,Pi5,Pi6,Pi7,Pi8八个节点连接而成,根据步骤6的ΔVbg(x,y,z)可得网格点Pi1~Pi8的变形量分别为(Δxi1,Δyi1,Δzi1),(Δxi2,Δyi2,Δzi2),(Δxi3,Δyi3,Δzi3),(Δxi4,Δyi4,Δzi4),(Δxi5,Δyi5,Δzi5),(Δxi6,Δyi6,Δzi6),(Δxi7,Δyi7,Δzi7),(Δxi8,Δyi8,Δzi8);
步骤7.2、根据点Pi的坐标信息,计算Pi在网格单元Ci中的相对位置Pi(ξ,η,ζ);
步骤7.3、根据步骤7.2获得的Pi(ξ,η,ζ)以及步骤7.1中获得的Pi1~Pi8的变形量,采用三线性插值计算出Pi(x,y,z)的变形量ΔPi(x,y,z)=(Δxi,Δyi,Δzi);具体的计算公式如下:
DX1=Δxi1
DX2=-Δxi1+Δxi2
DX3=-Δxi1+Δxi4
DX4=-Δxi1+Δxi5
DX5=Δxi1-Δxi2+Δxi3-Δxi4
DX6=Δxi1-Δxi2-Δxi5+Δxi6
DX7=Δxi1-Δxi4-Δxi5+Δxi8
DX8=-Δxi1+Δxi2-Δxi3+Δxi4+Δxi5-Δxi6+Δxi7-Δxi8
对DY1~DY8,DZ1~DZ8,都采用上述方法求得;
Δxi=DX1+DX2·ξ+DX3·η+DX4·ζ+DX5·ξη+DX6·ξζ+DX7·ηζ++DX8·ξηζ
Δyi=DY1+DY2·ξ+DY3·η+DY4·ζ+DY5·ξη+DX6·ξζ+DY7·ηζ++DY8·ξηζ
Δzi=DZ1+DZ2·ξ+DZ3·η+DZ4·ζ+DZ5·ξη+DZ6·ξζ+DZ7·ηζ++DZ8·ξηζ。
进一步的,如上所述的一种基于背景网格的动网格技术,步骤1中,采集体网格点坐标信息V(x,y,z)以及边界网格点坐标信息BC(x,y,z)时,对结构网格,直接采集体网格点坐标信息V(x,y,z)以及边界网格点坐标信息BC(x,y,z)时;对非结构网格,不采集网格单元连接信息。
本发明技术方案基于背景网格的动网格技术是指将动网格技术应用于网格较为稀疏的背景网格上,实现高效的网格变形,然后将背景网格的变形量插值到原始的体网格上,这样可以大大的提高网格变形的速度。
附图说明
图1是Z=0的截面的原始网格;
图2是Z=0的截面的变形后的网格;
图3是最小长方体;
图4是长方体网格剖分;
图5是变形后的背景网格A;
图6是变形后的背景网格B;
图7是Pi与Ci的位置关系。
具体实施方式
下面结合附图和具体实施例对本发明技术方案进行详细说明。
本发明采用径向基函数法对稀疏的背景网格进行变形,然后将背景网格的变形量插值到原始体网格的网格点上,从而达到提高动网格效率的目的。本发明的背景网格仅用于将变形量回插到原始体网格上,所以并不要求背景网格与原始体网格的边界吻合,以至于可以采用简单的可自动化生成的直角网格作为背景网格,大大地减少了人工干预,且变形后的网格质量仍然较高。
本发明一种基于背景网格的动网格技术,包括如下步骤:
步骤1、采集体网格点坐标信息V(x,y,z)以及边界网格点坐标信息BC(x,y,z);对非结构网格,不采集网格单元连接信息。
步骤2、根据步骤1所采集的体网格点坐标信息,创建一个完全覆盖体网格区域的长方体(如图3所示),并对此长方体进行网格剖分(如图4所示),完成背景网格Vbg(x,y,z)的生成;
步骤3、从步骤1所采集的边界网格点坐标信息中采集参与网格变形的网格点BCmv(x,y,z)old,并根据工程实际需求,给出变形后边界网格点的新坐标信息BCmv(x,y,z)new,然后计算边界网格点位置的改变量:
ΔBCmv(x,y,z)=BCmv(x,y,z)new-BCmv(x,y,z)old
步骤4、剔除步骤3中采集的BCmv(x,y,z)old中的重复点,并将步骤3中采集的BCmv(x,y,z)new以及ΔBCmv(x,y,z)中相应的网格点的信息剔除;
步骤5、根据步骤4中剔除重复点后的BCmv(x,y,z)old以及ΔBCmv(x,y,z),采用基于贪心算法的径向基函数法计算插值系数,从BCmv(x,y,z)old中获取经过贪心算法精简后的插值基函数支撑点集合U(P1,P2,…,Pn)以及相应的权系数ωi(i=1,2,…,n),其中n为集合U中的元素个数;
步骤6、根据步骤5所获得的支撑点集合U以及权系数ωi计算步骤2生成的背景网格Vbg(x,y,z)的变形量ΔVbg(x,y,z),并生成变形后的背景网格;如图5、6所示;
步骤7、根据步骤6获得的ΔVbg(x,y,z)插值出步骤1中体网格V(x,y,z)的变形量ΔV(x,y,z),本步骤可分为以下几步:
步骤7.1、在体网格V(x,y,z)选取任意网格点Pi(x,y,z),从步骤2的背景网格Vbg(x,y,z)中搜索出包含网格点Pi(x,y,z)的网格单元Ci;其位置关系如图7所示;
网格单元Ci由网格点Pi1,Pi2,Pi3,Pi4,Pi5,Pi6,Pi7,Pi8八个节点连接而成,根据步骤6的ΔVbg(x,y,z)可得网格点Pi1~Pi8的变形量分别为(Δxi1,Δyi1,Δzi1),(Δxi2,Δyi2,Δzi2),(Δxi3,Δyi3,Δzi3),(Δxi4,Δyi4,Δzi4),(Δxi5,Δyi5,Δzi5),(Δxi6,Δyi6,Δzi6),(Δxi7,Δyi7,Δzi7),(Δxi8,Δyi8,Δzi8);
步骤7.2、根据点Pi的坐标信息,计算Pi在网格单元Ci中的相对位置Pi(ξ,η,ζ);
步骤7.3、根据步骤7.2获得的Pi(ξ,η,ζ)以及步骤7.1中获得的Pi1~Pi8的变形量,采用三线性插值计算出Pi(x,y,z)的变形量ΔPi(x,y,z)=(Δxi,Δyi,Δzi);具体的计算公式如下:
DX1=Δxi1
DX2=-Δxi1+Δxi2
DX3=-Δxi1+Δxi4
DX4=-Δxi1+Δxi5
DX5=Δxi1-Δxi2+Δxi3-Δxi4
DX6=Δxi1-Δxi2-Δxi5+Δxi6
DX7=Δxi1-Δxi4-Δxi5+Δxi8
DX8=-Δxi1+Δxi2-Δxi3+Δxi4+Δxi5-Δxi6+Δxi7-Δxi8
对DY1~DY8,DZ1~DZ8,都采用上述方法求得;
Δxi=DX1+DX2·ξ+DX3·η+DX4·ζ+DX5·ξη+DX6·ξζ+DX7·ηζ++DX8·ξηζ
Δyi=DY1+DY2·ξ+DY3·η+DY4·ζ+DY5·ξη+DX6·ξζ+DY7·ηζ++DY8·ξηζ
Δzi=DZ1+DZ2·ξ+DZ3·η+DZ4·ζ+DZ5·ξη+DZ6·ξζ+DZ7·ηζ++DZ8·ξηζ。
图1为某体网格在Z=0截面的网格,采用本发明进行变形后,Z=0的截面的网格为图2所示。对比可见,变形后的网格质量较高。本发明与其它动网格技术相比的有益效果:
(1)本发明技术方案采用径向基函数法仅计算背景网格点的变形量,然后将背景网格的变形量插值到体网格上,因为背景网格非常稀疏,网格点较少,可以大大的提高动网格的效率,尤其对于求解Navier-Stokes方程所使用的较大规模的网格,效率提高非常明显。
(2)本发明技术方案不依赖于网格拓扑,能够适用于能任意网格类型的网格变形(结构网格、结构/非结构混合网格、点云等)。
(3)本发明技术方案所采用的背景网格为完全覆盖体网格区域的简单的直角网格,背景网格点不需要与边界吻合。
(4)本发明技术方案所采用的背景网格可以自动生成,省去了为了进行网格变形而专门生成的一套较为稀疏的背景网格。
(5)本发明技术方案所采用的背景网格为规整的直角网格,体网格点与背景网格之间的映射关系非常容易获得,大大缩短了映射搜索的时间。
Claims (2)
1.一种基于背景网格的动网格方法,其特征在于:包括如下步骤:
步骤1、采集体网格点坐标信息V(x,y,z)以及边界网格点坐标信息BC(x,y,z);
步骤2、根据步骤1所采集的体网格点坐标信息,创建一个完全覆盖体网格区域的长方体,并对此长方体进行网格剖分,完成背景网格Vbg(x,y,z)的生成;
步骤3、从步骤1所采集的边界网格点坐标信息中采集参与网格变形的网格点BCmv(x,y,z)old,并根据工程实际需求,给出变形后边界网格点的新坐标信息BCmv(x,y,z)new,然后计算边界网格点位置的改变量:
ΔBCmv(x,y,z)=BCmv(x,y,z)new-BCmv(x,y,z)old
步骤4、剔除步骤3中采集的BCmv(x,y,z)old中的重复点,并将步骤3中采集的BCmv(x,y,z)new以及ΔBCmv(x,y,z)中与BCmv(x,y,z)old中的重复的网格点的信息剔除;
步骤5、根据步骤4中剔除重复点后的BCmv(x,y,z)old以及ΔBCmv(x,y,z),采用基于贪心算法的径向基函数法计算插值系数,从BCmv(x,y,z)old中获取经过贪心算法精简后的插值基函数支撑点集合U(P1,P2,…,Pn)以及相应的权系数ωi(i=1,2,…,n),其中n为集合U中的元素个数;
步骤6、根据步骤5所获得的支撑点集合U以及权系数ωi计算步骤2生成的背景网格Vbg(x,y,z)的变形量ΔVbg(x,y,z),并生成变形后的背景网格;
步骤7、根据步骤6获得的ΔVbg(x,y,z)插值出步骤1中体网格V(x,y,z)的变形量ΔV(x,y,z),本步骤可分为以下几步:
步骤7.1、在体网格V(x,y,z)选取任意网格点Pi(x,y,z),从步骤2的背景网格Vbg(x,y,z)中搜索出包含网格点Pi(x,y,z)的网格单元Ci;
网格单元Ci由网格点Pi1,Pi2,Pi3,Pi4,Pi5,Pi6,Pi7,Pi8八个节点连接而成,根据步骤6的ΔVbg(x,y,z)可得网格点Pi1~Pi8的变形量分别为(Δxi1,Δyi1,Δzi1),(Δxi2,Δyi2,Δzi2),(Δxi3,Δyi3,Δzi3),(Δxi4,Δyi4,Δzi4),(Δxi5,Δyi5,Δzi5),(Δxi6,Δyi6,Δzi6),(Δxi7,Δyi7,Δzi7),(Δxi8,Δyi8,Δzi8);
步骤7.2、根据点Pi的坐标信息,计算Pi在网格单元Ci中的相对位置Qi(ξ,η,ζ);
步骤7.3、根据步骤7.2获得的Qi(ξ,η,ζ)以及步骤7.1中获得的Pi1~Pi8的变形量,采用三线性插值计算出Pi(x,y,z)的变形量ΔPi(x,y,z)=(Δxi,Δyi,Δzi);具体的计算公式如下:
DX1=Δxi1
DX2=-Δxi1+Δxi2
DX3=-Δxi1+Δxi4
DX4=-Δxi1+Δxi5
DX5=Δxi1-Δxi2+Δxi3-Δxi4
DX6=Δxi1-Δxi2-Δxi5+Δxi6
DX7=Δxi1-Δxi4-Δxi5+Δxi8
DX8=-Δxi1+Δxi2-Δxi3+Δxi4+Δxi5-Δxi6+Δxi7-Δxi8
对DY1~DY8,DZ1~DZ8,都采用上述方法求得;
Δxi=DX1+DX2·ξ+DX3·η+DX4·ζ+DX5·ξη+DX6·ξζ+DX7·ηζ++DX8·ξηζ
Δyi=DY1+DY2·ξ+DY3·η+DY4·ζ+DY5·ξη+DX6·ξζ+DY7·ηζ++DY8·ξηζ
Δzi=DZ1+DZ2·ξ+DZ3·η+DZ4·ζ+DZ5·ξη+DZ6·ξζ+DZ7·ηζ++DZ8·ξηζ。
2.如权利要求1所述的一种基于背景网格的动网格方法,其特征在于:步骤1中,采集体网格点坐标信息V(x,y,z)以及边界网格点坐标信息BC(x,y,z)时,对结构网格,直接采集体网格点坐标信息V(x,y,z)以及边界网格点坐标信息BC(x,y,z)时;对非结构网格,不采集网格单元连接信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410601962.2A CN105631065B (zh) | 2014-10-31 | 2014-10-31 | 一种基于背景网格的动网格方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410601962.2A CN105631065B (zh) | 2014-10-31 | 2014-10-31 | 一种基于背景网格的动网格方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105631065A CN105631065A (zh) | 2016-06-01 |
CN105631065B true CN105631065B (zh) | 2019-05-10 |
Family
ID=56045997
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410601962.2A Active CN105631065B (zh) | 2014-10-31 | 2014-10-31 | 一种基于背景网格的动网格方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105631065B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108180918B (zh) * | 2017-11-29 | 2021-04-30 | 西南科技大学 | 一种点云测地路径正向跟踪生成方法及装置 |
CN112632869B (zh) * | 2020-12-23 | 2022-09-27 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种基于网格框架的非结构附面层网格生成方法 |
CN113111612B (zh) * | 2021-06-15 | 2021-08-10 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种基于自适应空间剖分的离散点云重复点快速查找方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005216131A (ja) * | 2004-01-30 | 2005-08-11 | Digital Fashion Ltd | 化粧シミュレーション装置、化粧シミュレーション方法、化粧シミュレーションプログラム |
CN1828672A (zh) * | 2006-03-28 | 2006-09-06 | 中山大学 | 一种基于六面体插值体细分的几何数据细分方法 |
CN102682152A (zh) * | 2011-03-01 | 2012-09-19 | 利弗莫尔软件技术公司 | 使用扩展的无网格有限元法对结构性能的数值仿真 |
-
2014
- 2014-10-31 CN CN201410601962.2A patent/CN105631065B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005216131A (ja) * | 2004-01-30 | 2005-08-11 | Digital Fashion Ltd | 化粧シミュレーション装置、化粧シミュレーション方法、化粧シミュレーションプログラム |
CN1828672A (zh) * | 2006-03-28 | 2006-09-06 | 中山大学 | 一种基于六面体插值体细分的几何数据细分方法 |
CN102682152A (zh) * | 2011-03-01 | 2012-09-19 | 利弗莫尔软件技术公司 | 使用扩展的无网格有限元法对结构性能的数值仿真 |
Non-Patent Citations (2)
Title |
---|
Efficient mesh motion using radial basis functions with data;T.C.S. Rendall et al;《Journal of Computational Physics》;20090518;第6231-6249页 |
基于非结构背景网格的动网格方法研究;张加乐等;《航空计算技术》;20120331;第42卷(第2期);第95-99页 |
Also Published As
Publication number | Publication date |
---|---|
CN105631065A (zh) | 2016-06-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105631065B (zh) | 一种基于背景网格的动网格方法 | |
CN102819647A (zh) | 一种非均质材料随机微观结构有限元建模方法 | |
CN103353907A (zh) | 一种法兰连接螺栓极限强度校核的计算方法 | |
CN103310058B (zh) | 基于三维模型的尺寸链校核方法 | |
CN103065059B (zh) | 一种基于变量代换的辐射型配电网潮流计算方法 | |
CN104751006B (zh) | 一种计及变量相关性的概率潮流计算方法 | |
CN105631068A (zh) | 一种非结构网格cfd计算的网格边界条件处理方法 | |
CN102890703A (zh) | 一种网络异质多维标度方法 | |
CN103530812B (zh) | 一种基于局部敏感哈希的电网状态相似度量化分析方法 | |
CN114676522B (zh) | 融合gan和迁移学习的气动形状优化设计方法及系统及设备 | |
CN104408773B (zh) | 一种结构化网格非匹配界面插值的方法 | |
CN105932673B (zh) | 一种含大规模新能源的大型互联电网的外网等值方法 | |
CN104361625A (zh) | 一种基于射线原理的带边界保留的云数据精简算法 | |
CN105224726B (zh) | 结构网格动网格技术用于非结构网格流场求解器的方法 | |
CN113673186A (zh) | 一种基于stl文件的笛卡尔网格快速生成方法 | |
CN105678844A (zh) | 一种基于地物散乱点逐点增加的轮廓构建算法 | |
CN104535040A (zh) | 用于叶片的有限元单元划分方法和叶片的检测方法 | |
CN105550385A (zh) | 一种含分布式电源配电网的小步长暂态仿真方法及系统 | |
CN103942368A (zh) | 一种激光切割机床的结构设计方法 | |
CN106712029B (zh) | 小阻抗支路pq端点变雅可比矩阵的牛顿法潮流计算方法 | |
CN109241546A (zh) | 一种风机塔筒连接法兰疲劳强度的校核方法 | |
CN114943167A (zh) | 一种结构网格壁面距离的计算方法、系统、介质和设备 | |
CN109049715B (zh) | 一种基于多相材料拓扑优化的3d打印模型数据处理方法 | |
CN106156443A (zh) | 一种在冷却塔上布设加劲环和子午肋的方法及装置 | |
CN106156537A (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 |