CN108733892A - 一种关于混凝土中非凸形颗粒的三维数值构造方法 - Google Patents

一种关于混凝土中非凸形颗粒的三维数值构造方法 Download PDF

Info

Publication number
CN108733892A
CN108733892A CN201810376224.0A CN201810376224A CN108733892A CN 108733892 A CN108733892 A CN 108733892A CN 201810376224 A CN201810376224 A CN 201810376224A CN 108733892 A CN108733892 A CN 108733892A
Authority
CN
China
Prior art keywords
vertex
polygon
convex
coordinate
sphere
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
CN201810376224.0A
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.)
Hohai University HHU
Original Assignee
Hohai University HHU
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 Hohai University HHU filed Critical Hohai University HHU
Priority to CN201810376224.0A priority Critical patent/CN108733892A/zh
Publication of CN108733892A publication Critical patent/CN108733892A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation

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)
  • Processing Or Creating Images (AREA)

Abstract

本发明公开了一种关于混凝土中非凸形颗粒的三维数值构造方法。该方法具体包括如下几个步骤:步骤1:在局部坐标系下的某个坐标面内生成多边形,获取其顶点信息;步骤2:将步骤1生成的多边形扩展至三维空间并判断各个顶点的凹凸性;步骤3:确定拓展球体的半径和球心所在位置;步骤4:确定扩展球的滚动方向,并令扩展球沿非凸多边形的边滚动;步骤5:确定新的表面顶点坐标。本发明在计算机图形学的基础上,构造出三维的非凸形颗粒。本发明具有操作简单,对计算机硬件要求低等优点。

Description

一种关于混凝土中非凸形颗粒的三维数值构造方法
技术领域
本发明属于混凝土力学性能分析领域,具体涉及一种关于混凝土中非凸形颗粒的三维数值构造方法。
背景技术
混凝土作为一种颗粒增强复合材料,它的承载能力很大程度上受制于骨料,特别是它的形状和大小。长期以来,人们主要通过实验来分析骨料特征。但是,由于实验过程中出现的一系列不可控因素,经常会导致实验结果离散、可重复性差。近几年,随着计算机的迅猛发展,计算机仿真模拟骨料的几何形状已经在混凝土研究中得到广泛地应用。
关于混凝土骨料形状的相关数值构造工作,前期研究大部分是集中在球体、椭球体、凸面体等规则的几何凸体构造。实际混凝土中的骨料的形状并非完全是完美的凸体颗粒,而是复杂的非凸形颗粒组成。因此,构建非凸形颗粒的几何形貌是实施精确分析混凝土和其他颗粒材料的物理力学性能的基础和保证。相关研究人员提出了球拼装、像素扫描等方法来构造复杂的非凸形颗粒,但是它们存在着一些缺陷。对于球拼装法,很大程度上受到球体数量的限制,不能够精确的模拟出所要构造的非凸形颗粒,而像素扫描法则很大程度上受制于成像设备的分辨率。
发明内容
针对现有技术的不足,本发明的目的在于提供一种关于混凝土中非凸形颗粒的三维数值构造方法,用此方法模拟混凝土中的骨料有效克服了骨料与现实的差距,为混凝土的力学性能提供更有利的保障。
为解决现有技术问题,本发明采取的技术方案为:
一种关于混凝土中非凸形颗粒的三维数值构造方法,包括以下步骤:
步骤1,在O-XYZ全局坐标系中随机生成点P0(x0,y0,z0),利用公式(1)平移矩阵
将原点(0,0,0)平移至P0(x0,y0,z0)形成局部坐标系O1-X1Y1Z1,P0(x0,y0,z0)即为局部坐标系O1-X1Y1Z1的原点,将局部坐标系X1O1Y1面内(0,0)点作为极坐标的极心,若生成多边形的边数为N,则在区间[0,2π)内随机的取N个不同的极角值,可以通过公式(2)生成第i个顶点所对应的极角,将所取的极角按照从小到大的顺序排列,
θi=ηi×2π
(2),其中,θi是第i个顶点所对应的极角,ηi是[0,1)内的一个随机数;
在生成多边形各个顶点对应的极角的同时,随机地生成对应的极径,极径随机的分布在(A0-A1)到(A0+A1)之间,通过公式(3)可以随机地生成第i个顶点所对应的极径ri
ri=A0i×A1 (3),
其中,A0是颗粒的平均粒径;αi是[-1,1]内的一个随机数,A1是颗粒的最小粒径,ri是第i个顶点所对应的极径;
根据每个顶点所对应的极角和极径,通过公式(4)可以计算出在局部直角坐标系O1-X1Y1Z1下的第i个顶点的顶点坐标P′i(x′i,y′i,z′i),i=1,2,...,N;
xi=x0+ricos(θi)
yi=y0+risin(θi)
zi=z0
(4),
其中,ri和θi分别是第i个顶点对应的极径和极角;Pi′(x′i,y′i,z′i)是在局部直角坐标系下的第i个顶点的顶点坐标,i=1,2…N;
步骤2,将步骤1生成的多边形通过平移变换和欧拉变换扩展至三维空间中得多边形,并判断所得多边形中所有顶点的凹凸性
根据步骤1所得顶点坐标,通过平移变换和欧拉变换,将其在各自范围内随机地生成三个欧拉角α、β和γ。其中,α和γ的取值区间在[-π,π],β的取值区间在[0,π],计算其余弦和正弦值,得到如公式(5)所示的旋转矩阵,顺时针旋转的角度为正值,逆时针时旋转的角度为负值;
其中,R是旋转矩阵,通过公式(6),将在局部坐标O1-X1Y1Z1求得的第i个顶点的坐标Pi′(x′i,y′i,z′i),变换成在全局坐标系O-XYZ下的顶点坐标Pi(xi,yi,zi)
X=TRX′
(6),其中,X′是局部坐标系下的坐标,X是全局坐标系下的坐标,
通过多边形各个顶点转向来判断顶点的凹凸性,可以通过公式(7)计算出过顶点A(x1,y1,z1)、点B(x2,y2,z2)和点C(x3,y3,z3)的多边形的法向量n0,通过公式(8)来判断多边形的顶点A(x1,y1,z1)的凹凸性,其中,公式(7)中的点B(x2,y2,z2)和点C(x3,y3,z3)是与顶点A(x1,y1,z1)相邻的两个顶点,依次判断完多边形的每个顶点的凹凸性。若多边形的所有顶点为凸顶点,此时生成的多边形为无效多边形,重新生成;
其中,n0是过顶点(x1,y1,z1)的多边形法向量,V0是以(x1,y1,z1),(x0,y0,z0)为端点的多边形的棱边的方向向量,V1是以(x1,y1,z1),(x2,y2,z2)为端点的非凸多边形的棱边的方向向量;
其中,n是非多边形的法向量,n0是过顶点(x1,y1,z1)的多边形法向量。λ是判定顶点凹凸性的参数:当λ>0时,该顶点为凸顶点;当λ<0时,该顶点为凹顶点;
步骤3,确定扩展球半径和球心的位置
根据粒子的厚度确定扩展球的半径,根据步骤2得到的非凸多边形的顶点坐标,任意取一个顶点作为扩展球体的球心;
步骤4,确定扩展球的滚动方向,令扩展球沿非凸多边形的边滚动
根据步骤2得到的顶点坐标,确定非凸多边形各个边的方向,以任一顶点为扩展球体的球心,让扩展球体沿着经过该点的边滚动,滚动至球心位于邻近顶点位置,滚动的轨迹构成一个球柱体,按照相同的滚动方式,依次滚完非凸多边形的其余的边;
步骤5,确定新的表面顶点坐标
根据步骤4滚动所得的球柱体,首先通过公式(9)计算得出初始状态的法向量n(x0,y0,z0),基于该法向量通过公式(10)计算得出新的表面顶点坐标(x,y,z)
x0=y1z2+y2z3+y3z1-y1z3-y2z1-y3z2
y0=x1z3+x2z1+x3z2-x1z2-x2z3-x3z1
z0=x1y2+x2y3+x3y1-x1y3-x2y1-x3y2 (9),
其中,(x1,y1,z1)、(x2,y2,z2)、(x3,y3,z3)为非凸多边形平面内任意不共线的三点坐标,通常随机地取其三个顶点,计算出初始状态的法向量;(x0,y0,z0)是初始状态的法向量n,即非凸多边形的法向量;
其中,(X,Y,Z)为扩展球的球心坐标;r是扩展球体的半径;(x0,y0,z0)是初始状态的法向量n。
有益效果
与现有技术相比,本发明在计算机图形学基础上,构造出三维状态下的非凸形颗粒,克服了模拟出的颗粒形状和现实相差很大的问题。
本发明构造方法简单,操作简便,且对硬件要求低,成本低,为混凝土性能的计算提供了更有利的证据。
附图说明
图1为本发明的三维非凸形颗粒构造的流程图;
图2为模拟生成的在三维空间中的一个非凸多边形单元示意图;
图3为扩展球体所在非凸多边形单元某个顶点的示意图;
图4为扩展球体沿着非凸多边形单元的棱边滚动示意图;
图5为生成的新的表面顶点坐标示意图;
图6为实现的边数为8的三维非凸形颗粒的示例图。
具体实施方式
下面结合附图和具体实施例,进一步阐明本发明,应理解这些实施例仅用于说明本发明而不用于限制本发明的范围。在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
本发明的思路是在现有的计算机图形学理论的基础上,实现了三维非凸形颗粒的构造。具体而言,本发明采用以下技术方案解决上述技术问题。
对于需要构造的三维非凸形颗粒,将在三维空间中的非凸多边形作为基本单元体。扩展单元体设为半径为r的球体,且球心在基本单元体表面上。具体的构造方法图如流程图1所示。
一种关于混凝土中非凸形颗粒的三维数值构造方法,包括以下步骤:
步骤1,在O-XYZ全局坐标系中随机生成点P0(x0,y0,z0),利用公式(1)平移矩阵(1),将原点(0,0,0)平移至P0(x0,y0,z0)形成局部坐标系O1-X1Y1Z1,P0(x0,y0,z0)即为局部坐标系O1-X1Y1Z1的原点,将局部坐标系X1O1Y1面内(0,0)点作为极坐标的极心,若生成多边形的边数为N,则在区间[0,2π)内随机的取N个不同的极角值,可以通过公式(2)生成第i个顶点所对应的极角,将所取的极角按照从小到大的顺序排列,
θi=ηi×2π
(2),其中,θi是第i个顶点所对应的极角,ηi是[0,1)内的一个随机数;
在生成多边形各个顶点对应的极角的同时,随机地生成对应的极径,极径随机的分布在(A0-A1)到(A0+A1)之间,通过公式(3)可以随机地生成第i个顶点所对应的极径ri
ri=A0i×A1 (3),
其中,A0是颗粒的平均粒径;αi是[-1,1]内的一个随机数,A1是颗粒的最小粒径,ri是第i个顶点所对应的极径;
根据每个顶点所对应的极角和极径,通过公式(4)可以计算出在局部直角坐标系O1-X1Y1Z1下的第i个顶点的顶点坐标P'i(x'i,y'i,z'i),i=1,2,...,N;
其中,ri和θi分别是第i个顶点对应的极径和极角;Pi′(x′i,y′i,z′i)是在局部直角坐标系下的第i个顶点的顶点坐标,i=1,2…N;
步骤2,,将步骤1生成的多边形通过平移变换和欧拉变换扩展至三维空间中得多边形,并判断所得多边形中所有顶点的凹凸性
根据步骤1所得顶点坐标,通过平移变换和欧拉变换,将其在各自范围内随机地生成三个欧拉角α、β和γ。其中,α和γ的取值区间在[-π,π],β的取值区间在[0,π],计算其余弦和正弦值,得到如公式(5)所示的旋转矩阵。顺时针旋转的角度为正值,逆时针时旋转的角度为负值;
其中,R是旋转矩阵,通过公式(6),将在局部坐标O1-X1Y1Z1求得的第i个顶点的坐标P′i(x′i,y′i,z′i),变换成在全局坐标系O-XYZ下的顶点坐标Pi(xi,yi,zi)
X=TRX′ (6),
其中,X′是局部坐标系下的坐标,X是全局坐标系下的坐标,
通过多边形各个顶点转向来判断顶点的凹凸性,可以通过公式(7)计算出过顶点A(x1,y1,z1)、点B(x2,y2,z2)和点C(x3,y3,z3)的多边形的法向量n0,通过公式(8)来判断多边形的顶点A(x1,y1,z1)的凹凸性,其中,公式(7)中的点B(x2,y2,z2)和点C(x3,y3,z3)是与顶点A(x1,y1,z1)相邻的两个顶点,依次判断完多边形的每个顶点的凹凸性。若多边形的所有顶点为凸顶点,此时生成的多边形为无效多边形,重新生成;
其中,n0是过顶点(x1,y1,z1)的多边形法向量,V0是以(x1,y1,z1),(x0,y0,z0)为端点的多边形的棱边的方向向量,V1是以(x1,y1,z1),(x2,y2,z2)为端点的非凸多边形的棱边的方向向量;
其中,n是非多边形的法向量,n0是过顶点(x1,y1,z1)的多边形法向量。λ是判定顶点凹凸性的参数:当λ>0时,该顶点为凸顶点;当λ<0时,该顶点为凹顶点,如图2所示;
步骤3,设扩展球体的半径为r,以某个顶点为扩展球体的球心,如图3所示;
步骤4,扩展球体以该顶点位置为出发点,沿着棱边滚动至球心位于另一顶点位置,滚动的轨迹构成一个球柱体,球柱体的上、下顶点分别对应于该棱边的两个顶点,则该圆柱体表露在颗粒外围的区域即三维非凸颗粒的边,如图4所示;
步骤5,根据所得的球柱体,计算出新的表面顶点位置。具体细节如下:假设,非凸多边形的三个顶点坐标为(x1,y1,z1)、(x2,y2,z2)、(x3,y3,z3),通过公式(1)可以求出,基本单元非凸多边形的平面法向量n(x0,y0,z0)
x0=y1z2+y2z3+y3z1-y1z3-y2z1-y3z2
y0=x1z3+x2z1+x3z2-x1z2-x2z3-x3z1
z0=x1y2+x2y3+x3y1-x1y3-x2y1-x3y2 (9)
其中,(x1,y1,z1)、(x2,y2,z2)、(x3,y3,z3)是非凸多边形的三个顶点坐标;(x0,y0,z0)是基本单元的法向量,即非凸多边形的法向量。
然后通过公式(10)可以计算出现的表面顶点位置,从而构造出三维的非凸形颗粒。
其中,(X,Y,Z)为扩展球的球心坐标;r是扩展球体的半径;(x0,y0,z0)是初始状态的法向量n。
实施例2
图6展示了,基于边数N=8的非凸边形构造出三维的非凸形颗粒。将扩展球体的半径r设为1,其球心为非凸多边形的某个顶点;让扩展球沿着非凸多边形的棱边滚动;而后计算其新的表面的顶点坐标。最终得到图6最右侧的三维非凸形颗粒。

Claims (5)

1.一种关于混凝土中非凸形颗粒的三维数值构造方法,其特征在于,包括以下步骤:
步骤1,在局部坐标系的任一坐标平面内生成多边形,获取其顶点坐标信息在全局坐标系O-XYZ上随机生成一个点P0(x0,y0,z0),并将全局坐标系平移至该点,形成以P0为原点的局部坐标系O1-X1Y1Z1,随机生成欧拉角,计算旋转矩阵和平移矩阵,将局部坐标系X1O1Y1面内(0,0)点作为极坐标的极心,基于该极心,随机生成各个顶点对应的极角和极径,计算在局部坐标系下的顶点坐标;
步骤2,将步骤1生成的多边形扩展至三维空间,并判断其任一顶点的凹凸性根据步骤1计算的顶点坐标,通过平移变换和欧拉变换,将其扩展至一个三维空间中得,并判断其顶点的凹凸性,若生成的多边形中至少有1个顶点为凹性,即多边形为非凸多边形,则进行步骤3,若生成的多边形所有顶点中均为凸形,则回到步骤1重新生成多边形,直至满足条件;
步骤3,确定扩展球体半径和球心所在位置
根据步骤2中非凸多边形的顶点坐标,任意选取一个顶点作为扩展球体的球心;
步骤4,确定扩展球的滚动方向,令扩展求沿边滚动,直至扩展求滚动完所有边
根据步骤2得到的顶点坐标,确定非凸多边形各个边的方向,以任一顶点为扩展球体的球心,让扩展球体沿着经过该点的边滚动,滚动至球心位于邻近顶点位置,滚动的轨迹构成一个球柱体,按照相同的滚动方式,依次滚完非凸多边形的其余的边;
步骤5,确定新的表面顶点坐标
依据步骤4滚动得到的球柱体,计算出初始状态的法向量n(x0,y0,z0),再公式(10)计算出新的表面顶点坐标(x,y,z),其中,公式(10)如下所示
其中,(X,Y,Z)为扩展球的球心坐标;r是扩展球体的半径;(x0,y0,z0) 是初始状态的法向量n。
2.根据权利要求1所述的一种关于混凝土中非凸形颗粒的三维数值构造方法,其特征在于,步骤1中各个顶点对应的极角和极径的生成方法为:若多边形的边数为N,通过在[0,2π)内随机地生成极角θi,i=1,2,…,N,并将其按照从小达到地顺序排列,可以通过公式(2)生成第i个顶点所对应的极角,θi=ηi×2π (2),
其中,ηi是[0,1)内随机生成的一个随机数,在生成多边形各个顶点对应的极角的同时,通过公式(3)来生成第i个顶点所对应的极径ri
ri=A0i×A1 (3),
其中,A0是颗粒的平均粒径;αi是[-1,1]内的一个随机数,A1是颗粒的最小粒径,ri是第i个顶点所对应的极径。
3.根据权利要求1所述的一种关于混凝土中非凸形颗粒的三维数值构造方法,其特征在于,步骤2中判断所得多边形中所有顶点的凹凸性的方法,具体步骤如下:根据步骤1所得顶点坐标,通过平移变换和欧拉变换,将其在各自范围内随机地生成三个欧拉角α、β和γ,其中,α和γ的取值区间在[-π,π],β的取值区间在[0,π],计算其余弦和正弦值,得到如公式(5)所示的旋转矩阵,顺时针旋转的角度为正值,逆时针时旋转的角度为负值;
其中,R是旋转矩阵,通过公式(6),将在局部坐标O1-X1Y1Z1求得的第i个顶点的坐标P′i(x′i,y′i,z′i),变换成在全局坐标系O-XYZ下的顶点坐标Pi(xi,yi,zi)X=TRX′ (6),
其中,X′是局部坐标系下的坐标,X是全局坐标系下的坐标,通过多边形各个顶点转向来判断顶点的凹凸性,可以通过公式(7)计算出过顶点A(x1,y1,z1)、点B(x2,y2,z2)和点C(x3,y3,z3)的多边形的法向量n0,通过公式(8)来判断多边形的顶点A(x1,y1,z1)的凹凸性,其中,公式(7)中的点B(x2,y2,z2)和点C(x3,y3,z3)是与顶点A(x1,y1,z1)相邻的两个顶点,依次判断完多边形的每个顶点的凹凸性,若多边形的所有顶点为凸顶点,此时生成的多边形为无效多边形,重新生成;
其中,n0是过顶点(x1,y1,z1)的多边形法向量,V0是以(x1,y1,z1),(x0,y0,z0)为端点的多边形的棱边的方向向量,V1是以(x1,y1,z1),(x2,y2,z2)为端点的非凸多边形的棱边的方向向量;
其中n是非多边形的法向量,n0是过顶点(x1,y1,z1)的多边形法向量,λ是判定顶点凹凸性的参数:当λ>0时,该顶点为凸顶点;当λ<0时,该顶点为凹顶点。
4.根据权利要求3所述的一种关于混凝土中非凸形颗粒的三维数值构造方法,其特征在于,
5.根据权利要求1所述的一种关于混凝土中非凸形颗粒的三维数值构造方法,其特征在于,步骤5中法向量n(x0,y0,z0)的计算公式(9)如下所示
x0=y1z2+y2z3+y3z1-y1z3-y2z1-y3z2
y0=x1z3+x2z1+x3z2-x1z2-x2z3-x3z1
z0=x1y2+x2y3+x3y1-x1y3-x2y1-x3y2 (9),
其中,(x1,y1,z1)、(x2,y2,z2)、(x3,y3,z3)为非凸多边形平面内任意不共线的三点坐标。
CN201810376224.0A 2018-04-25 2018-04-25 一种关于混凝土中非凸形颗粒的三维数值构造方法 Pending CN108733892A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810376224.0A CN108733892A (zh) 2018-04-25 2018-04-25 一种关于混凝土中非凸形颗粒的三维数值构造方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810376224.0A CN108733892A (zh) 2018-04-25 2018-04-25 一种关于混凝土中非凸形颗粒的三维数值构造方法

Publications (1)

Publication Number Publication Date
CN108733892A true CN108733892A (zh) 2018-11-02

Family

ID=63939848

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810376224.0A Pending CN108733892A (zh) 2018-04-25 2018-04-25 一种关于混凝土中非凸形颗粒的三维数值构造方法

Country Status (1)

Country Link
CN (1) CN108733892A (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110083936A (zh) * 2019-04-26 2019-08-02 中煤科工集团重庆研究院有限公司 考虑形状参数的可破碎碎石颗粒生成方法
CN111080790A (zh) * 2019-11-28 2020-04-28 河海大学 一种基于凸包算法的凹体生成方法和装置
CN111159927A (zh) * 2019-11-24 2020-05-15 浙江大学 基于体素矩阵的三维不规则形状颗粒投放的数值建模方法
CN112632814A (zh) * 2020-12-04 2021-04-09 西安理工大学 一种基于椭球体随机骨料的再生混凝土三维模型构建方法
CN114969892A (zh) * 2022-04-20 2022-08-30 西安建筑科技大学 碎石粒料颗粒接触方式的界定量化方法、系统及存储介质
CN115861399A (zh) * 2022-09-29 2023-03-28 优脑银河(浙江)科技有限公司 基于医学影像图像的配准方法、可读存储介质及电子设备
CN116486967A (zh) * 2023-05-12 2023-07-25 哈尔滨工业大学 可控颗粒间距及面积占比的带壳颗粒随机分布的生成方法

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110083936A (zh) * 2019-04-26 2019-08-02 中煤科工集团重庆研究院有限公司 考虑形状参数的可破碎碎石颗粒生成方法
CN110083936B (zh) * 2019-04-26 2023-06-16 中煤科工集团重庆研究院有限公司 考虑形状参数的可破碎碎石颗粒生成方法
CN111159927A (zh) * 2019-11-24 2020-05-15 浙江大学 基于体素矩阵的三维不规则形状颗粒投放的数值建模方法
CN111080790A (zh) * 2019-11-28 2020-04-28 河海大学 一种基于凸包算法的凹体生成方法和装置
CN111080790B (zh) * 2019-11-28 2022-04-05 河海大学 一种基于凸包算法的凹体生成方法和装置
CN112632814A (zh) * 2020-12-04 2021-04-09 西安理工大学 一种基于椭球体随机骨料的再生混凝土三维模型构建方法
CN114969892A (zh) * 2022-04-20 2022-08-30 西安建筑科技大学 碎石粒料颗粒接触方式的界定量化方法、系统及存储介质
CN114969892B (zh) * 2022-04-20 2024-05-17 西安建筑科技大学 碎石粒料颗粒接触方式的界定量化方法、系统及存储介质
CN115861399A (zh) * 2022-09-29 2023-03-28 优脑银河(浙江)科技有限公司 基于医学影像图像的配准方法、可读存储介质及电子设备
CN115861399B (zh) * 2022-09-29 2023-11-28 优脑银河(浙江)科技有限公司 基于医学影像图像的配准方法、可读存储介质及电子设备
CN116486967A (zh) * 2023-05-12 2023-07-25 哈尔滨工业大学 可控颗粒间距及面积占比的带壳颗粒随机分布的生成方法
CN116486967B (zh) * 2023-05-12 2023-10-03 哈尔滨工业大学 可控颗粒间距及面积占比的带壳颗粒随机分布的生成方法

Similar Documents

Publication Publication Date Title
CN108733892A (zh) 一种关于混凝土中非凸形颗粒的三维数值构造方法
TWI709106B (zh) 基於深度學習網路的室內場景結構估測系統及其估測方法
Lee et al. Complete swept volume generation, Part I: Swept volume of a piecewise C1-continuous cutter at five-axis milling via Gauss map
CN106777822B (zh) 二维高含量颗粒增强复合材料三相细观模型及建立方法
CN103489214A (zh) 增强现实系统中基于虚拟模型预处理的虚实遮挡处理方法
CN102982551B (zh) 空间三条不平行直线线性求解抛物折反射摄像机内参数
CN102411794A (zh) 一种基于球谐变换的三维模型的二维投影的输出方法
CN103955960B (zh) 一种基于单幅输入图像的图像视点变换方法
CN104217435A (zh) 两个相互遮挡的球线性确定拋物折反射摄像机内参数方法
CN104732565A (zh) 一种基于监控范围可视化设计的安防布局绘图方法
CN106500626A (zh) 一种手机立体成像方法及三维成像手机
CN109583377A (zh) 一种管路模型重建的控制方法、装置及上位机
Wang et al. Hexagonal meshes with planar faces
Sadoc et al. Phyllotaxis on surfaces of constant Gaussian curvature
CN101609563A (zh) 一种三维模型形状特征二叉树的构建方法
CN101196988A (zh) 三维掌纹身份鉴别系统的手掌定位及中心区提取方法
Bae et al. NURBS surface fitting using orthogonal coordinate transform for rapid product development
CN105303570A (zh) 使用双球的相交像及圆环点的像标定拋物折反射摄像机
CN106023314A (zh) 一种基于回转轴方向映射的b样条母曲线拟合方法
CN109360247A (zh) 单个球的公共自极三角形及正交消失点标定抛物摄像机
CN109523598A (zh) 单个球的公共自极三角形及圆环点标定抛物摄像机的方法
CN102760305B (zh) 采用旋转tin网和距离与夹角双权内插法的三维复杂建模方法
CN104700454B (zh) 一种混凝土非均质各向异性集料周围界面过渡区几何拓扑结构的构造方法
CN102663822A (zh) 基于细分方法的分子表面网格生成方法
CN108734181A (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
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20181102