CN102609560B - 一种3d任意粗糙表面的数字化模拟方法 - Google Patents
一种3d任意粗糙表面的数字化模拟方法 Download PDFInfo
- Publication number
- CN102609560B CN102609560B CN201110421316.4A CN201110421316A CN102609560B CN 102609560 B CN102609560 B CN 102609560B CN 201110421316 A CN201110421316 A CN 201110421316A CN 102609560 B CN102609560 B CN 102609560B
- Authority
- CN
- China
- Prior art keywords
- sequence
- gaussian
- rough surface
- sigma
- high degree
- 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
Links
Landscapes
- Complex Calculations (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种3D任意粗糙表面的数字化模拟方法,提高了模拟效率与统计参数的模拟精度。通过白噪声序列的反傅里叶变换获得随机初始相角序列,利用初始相角序列获得新的白噪声及其傅里叶变换;通过指定的自相关函数,进行离散与傅里叶变换等处理,获得高斯粗糙表面高度序列的功率谱密度与传递函数;利用频域点乘并求反傅里叶变换的方法完成高斯表面高度序列的模拟;在此基础上,通过指定偏斜度与峰度等高度分布统计特征参数,利用Pearson与Johnson非高斯转换***相结合,生成非高斯粗糙表面;若偏斜度与峰度模拟精度不合格,则更新相角序列与白噪声的傅里叶变换,重新进行高斯滤波与非高斯转换,直到满足给定精度要求。
Description
技术领域
本发明涉及一种3D微观任意粗糙表面的数字模拟方法,更具体地说,本发明,根据连接表面点的高度分布前四阶矩统计特征与水平方向自相关函数统计规律,可以高效准确地构建具有相同统计特征的粗糙表面。
背景技术
高速、精密加工技术的发展及产品可靠度要求的不断提高,推动着粗糙表面的加工制造、检测、数字化模拟与性能预测等领域的快速发展。
近年来,国内外在粗糙表面的数字模拟领域的研究,主要集中在针对粗糙高度符合高斯分布的粗糙表面方面,在非高斯粗糙表面的数字化模拟方面,已能成功将高斯分布的数据转化为非高斯数据,出现了基于Johnson、Pearson等非高斯转换***。但在非高斯粗糙表面的数字化模拟当中,经非高斯转换后的非高斯序列还需要开展滤波,给模拟序列的偏斜度、峰度等统计特征参数引入了新的误差,同时也改变了偏斜度-峰度的对应的模拟区域。从而导致有些统计特征的粗糙表面不能通过Johnson、Pearson等转换***转换与滤波生成相应的非高斯序列,即使能转换,也会在滤波阶段引入新的误差,影响滤波的非高斯表面数字化模拟的精度。
考虑到滤波(滤波系数为h(k,l))的影响,Johnson、Pearson等非高斯转换***中生成的非高斯序列粗糙表面的偏斜度Skη、峰度Kuη与滤波后的粗糙表面的偏斜度SkZ、峰度KuZ存在如下关系:
即使目标偏斜度与峰度满足非高斯转换必要条件:但在考虑滤波系数与自相关函数的相互影响后,Skη、Kuη有可能不再满足非高斯转换的必要条件:从而导致不能通过Johnson或Pearson等非高斯转换***生成非高斯序列;即使能转换,也会在滤波阶段引入新的误差,影响滤波的非高斯表面数字化模拟的精度。
现实中,不同的加工方法会产生具有不同统计特征的粗糙表面,粗糙表面的数字化模拟模型作为结合面性能分析的载体,其模拟的准确性直接影响相应结合面的接触性能预测。本研究利用一定的相角序列生成白噪声,通过白噪声的高斯滤波与高斯表面高度序列的非高斯转换***相结合,形成满足高度分布统计特征与自相关函数的非高斯高度序列,高斯或非高斯粗糙表面的数字化模拟方法为结合面性能研究奠定了基础。
发明内容
本发明的主要目的是提供一种高效高精度实现任意粗糙表面的数字化模拟方法,减少了粗糙表面数字化模拟中的滤波误差环节,在满足自相关函数、高度分布前四阶矩等统计特征的前提下,提高了模拟粗糙表面高度分布统计特征的精度与数字化模拟的效率。
本发明通过以下技术方案实现,主要步骤包括:
一种3D任意粗糙表面的数字化模拟方法,其特征在于,包含下述步骤:
1)当任意粗糙表面的偏斜度值约等于3且峰度值约等于0时,其表面形貌符合高斯分布,否则,其表面形貌符合非高斯分布;
2)如果步骤1)中粗糙表面形貌为高斯分布,利用三维高斯粗糙表面的模拟方法获得粗糙表面高度序列z,完成粗糙表面的数字化模拟;
3)如果步骤1)中粗糙表面形貌为非高斯分布,利用三维非高斯粗糙表面的模拟方法获得粗糙表面高度序列z,完成粗糙表面的数字化模拟。
所述的三维高斯粗糙表面的模拟方法,具体包含以下步骤:
1)通过给定的自相关函数生成离散的自相关函数矩阵R(m,n),其中m、n分别表示x、y方向取点的数目;
2)对离散的自相关函数矩阵R(m,n)进行快速傅里叶变换,获得高斯粗糙表面的功率谱密度P(m,n);
3)假设白噪声的功率谱密度常数C等于1,利用获得的功率谱密度,计算传递函数H(m,n);
4)利用白噪声生成器,获得白噪声序列η(m,n),利用反傅里叶变换生成相角序列Φ;
5)以步骤4)获得的相角序列Φ,计算相角相关序列exp(iΦ),并对其进行反傅里叶变换,更新白噪声序列η(m,n)及其傅里叶变换A(m,n);
6)以步骤3)获得的传递函数H(m,n),步骤5)获得的白噪声傅里叶变换序列A(m,n),求序列H与序列A的点积运算,并对点积运算结果进行反傅里叶变换,获得高斯表面的高度序列z0(m,n);
7)以步骤6)生成的高斯高度序列z0,对其进行高度方向平移与缩放,获得满足均值与标准偏差的高斯高度序列z,完成高斯粗糙表面的数字化模拟。
所述的三维非高斯粗糙表面的模拟方法,具体包含以下步骤:
1)以获得的高斯高度序列z,利用高斯序列的标准化方法,在高度方向进行平移与缩放,获得标准化正态分布序列Z0;
2)以步骤1)中获得的标准化正态分布序列Z0与高度序列的前四阶矩统计特征作为输入,利用Johnson或Pearson转换***进行非高斯变换,获得非高斯序列Z1;
3)利用步骤2)生成的非高斯序列Z1,对其进行平移与缩放,获得满足均值与标准偏差高度统计特征非高斯序列Z,同时,计算其偏斜度与峰度高度统计特征,判断其精度是否满足要求;如满足精度要求,则完成非高斯粗糙表面的数字化模拟,否则,继续执行后面步骤;
4)以步骤3)获得的非高斯序列Z1,利用反傅里叶变换生成相角序列Φ,更新白噪声序列η(m,n)及其傅里叶变换A(m,n);
5)以步骤5)获得序列A与权利要求1所获得的传递函数H(m,n),通过求序列A与H的点积运算,并对点积运算结果进行反傅里叶变换,获得新的高斯表面的高度序列z0(m,n);
6)以步骤5)生成的高斯高度序列z0,利用高斯序列的标准化方法,获得满足均值与标准偏差分别为0与1的新标准化正态分布序列Z0;
7)重复步骤2)~3),直到非高斯序列Z满足偏斜度与峰度的精度要求为止。
所述对z0高度方向进行平移与缩放的方法,具体包含以下步骤:
1)对高度序列z0进行平移,使其均值等于给定μ;
z0=z0-mean2(z0)+μ
2)对经平移后的高度序列z0在高度方向进行缩放,获得标准偏差为σ的高斯高度序列z。
z=σ·z0/std2(z0)
所述高斯序列的标准化方法,具体包含以下步骤:
1)对高度序列z或z0进行平移,使其均值等于0;
z0=z0-mean2(z0)或z=z-mean2(z)
2)对经平移后的高度序列z或z0在高度方向进行缩放,获得标准偏差1的标准高斯序列。
Z0=σ·z0/std2(z0)或Z0=σ·z/std2(z)
所述的利用Johnson或Pearson转换***进行非高斯变换,具体步骤如下:
1)优先采用Johnson非高斯转换***;
对数正态拟合方式SL:Z1=ξ+λexp[(Z0-γ)δ] (ξ<η)
有界***拟合方式SB:z=ξ+λe(η-γ)δ/[1+e(η-γ)/δ] (ξ<η<ξ+λ)
无界***拟合方式SU:z=ξ+λsinh[(η-γ)/δ]
γ、δ为形状参数,λ、ξ分别为比例参数与位置参数,γ、δ、λ与ξ均由均值、标准偏差、偏斜度与峰度等高度统计特征输入参数确定。
2)如果有界***SB拟合方式无法收敛,则采用Pearson非高斯转换***,以提高非高斯转换***的精度:
利用概率密度函数P满足以下微分方程,根据均值、方差、偏斜度与峰度等高度统计参数获得概率密度函数与高度分布函数。
式中a、c0、c1与c2是由均值、标准偏差、偏斜度与峰度等高度统计特征输入参数确定。
所述的利用反傅里叶变换生成相角序列Φ,更新白噪声序列η(m,n)及其傅里叶变换A(m,n),具体包含以下步骤:
1)以权利要求3所获得的非高斯序列Z1,通过反傅里叶变换,生成新的白噪声相角序列Φ(m,n);
2)以步骤1)获得的相角序列Φ(m,n),计算相角相关序列exp(iΦ),并对其进行反傅里叶变换,更新白噪声序列η(m,n)及其傅里叶变换A(m,n);
3)更新高斯粗糙表面数字化模拟中的白噪声序列η(m,n)的傅里叶变换序列A(m,n)。
本发明提出一种任意3D粗糙表面的数字化模拟方法,从高斯滤波、高精度非高斯变换、变相角循环迭代技术三个方面入手,克服了现有技术瓶颈,提高了3D微观粗糙表面的模拟精度和模拟效率,在一定程度上保证了整个偏斜度-峰度平面对应高度统计特征的任意非高斯表面的数字化模拟的可行性、准确性与高效性。可为结合面的接触性能研究提供高精度统计特征的形貌模型。
附图说明
图1各向同性高斯表面及其高度分布(βx=βy=5,μ=0,σ=1.6,点数128×128)
图2各向异性高斯表面及其高度分布(βx=5,βy=50,μ=0,σ=1.6,点数512×512)
图3各向同性非高斯粗糙表面及其高度分布(βx=βy=5,μ=0,σ=1.6,Skz=-1,Kuz=4,点数128×128)
图4各向同性粗糙表面的x、y方向平均自相关函数与理论自相关函数的对比:(a)x方向自相关函数(βx=5),(b)y方向自相关函数(βx=5)
图5各向异性非高斯粗糙表面及其高度分布(βx=5,βy=50,μ=0,σ=1.6,Skz=1,Kuz=5点数512×512)
图6各向异性粗糙表面的x、y方向平均自相关函数与理论自相关函数的对比:(a)x方向自相关函数(βx=5),(b)y方向自相关函数(βy=50)
具体实施方式
以满足一定的自相关函数、均值、标准偏差、偏斜度与峰度等统计特征的3D高斯或非高斯粗糙表面的数字化模拟为例,对本发明的任意3D粗糙表面的数字化模拟方法进行说明。
1.高斯粗糙表面的数字化模拟:
1)利用给定自相关函数生成相关函数离散矩阵R(m,n)
a)假定自相关函数为f(x,y)(σ为高斯序列的标准偏差,βx、βy分别为X、Y方向的自相关长度);
b)指定X、Y方向的离散间距均为1μm,并在-m/2≤x≤m/2与-n/2≤y≤n/2区域内对自相关函数进行离散,获得自相关函数的离散矩阵R(m+1,n+1);
c)所得R(m+1,n+1)序列中的第m/2+1,m/2+2行重复,删除其中一行;第n/2+1与n/2+2列重复,删除其中一列;获得与模拟表面高度序列具有相同行列数的自相关函数矩阵R(m,n)。
2)根据自相关函数求高斯表面的功率谱密度与传递函数;
a)对自相关函数离散矩阵进行傅里叶变换,并求解高斯粗糙表面的功率谱密度P;
b)根据白噪声的功率谱密度为常数C,假定C=1,则传递函数为;
3)白噪声序列的生成与傅里叶变换。
a)利用白噪声生成器randn(m,n)产生白噪声序列η(m,n);
b)求白噪声序列η的傅里叶变换A(m,n):
4)利用频域点乘(A.*H)取反傅里叶变换的方法获得高斯表面的高度序列;
5)根据给定高度均值μ、标准偏差σ与获得序列z0的实际均值mean2(z0)(mean2(z0)求序列均值的函数)与标准偏差std2(z)(std2为求序列标准偏差函数)的关系,对高斯高度序列z0进行高度方向平移与缩放,消除白噪声功率谱密度长度C不等于1对高度序列统计特征的影响,获得满足均值与标准偏差分布为μ与σ的高斯高度序列z,完成高斯粗糙表面的数字化模拟。
z0=z0-mean2(z0)+μ z=σ·z0/std2(z0)
利用上述给定的自相关函数,x、y方向的自相关函数均取5μm,模拟表面的x,y方向取点为128×128、表面的均值为0,标准偏差为1.6μm,模拟获得的各向同性高斯表面形貌及其高度分布如图1所示,其前四阶矩高度统计特征误差如表1所示。x、y方向的自相关函数分别取5μm、50μm,模拟获得的高斯表面形貌表现出各向异性,具有一定的纹理。模拟获得的各向同性高斯表面形貌及其高度分布如图2所示,其前四阶矩高度统计特征误差如表2所示。
表1 各向同性高斯粗糙表面(βx=βy=5,μ=0,σ=1.6,点数128×128)的高度统计前四阶矩仿真精度
表2 各向异性高斯粗糙表面(βx=5,βy=50,μ=0,σ=1.6,点数512×512)的高度统计前四阶矩仿真精度计算
2.非高斯粗糙表面的数字化模拟:
6)对上述模拟的高斯高度序列z平移与缩放,获得均值与标准偏差分别为0与1的标准化正态分布序列Z0;
z=z-mean2(z) Z0=z/std2(z)
7)采用Johnson或Pearson非高斯转换***进行非高斯转换,生成非高斯序列Z1;
a)优先采用Johnson非高斯转换***;
对数正态拟合方式SL:Z1=ξ+λexp[(Z0-γ)δ] (ξ<η)
有界***拟合方式SB:
无界***拟合方式SU:Z1=ξ+λsinh[(Z0-γ)/δ]
γ、δ为形状参数,λ、ξ分别为比例参数与位置参数,γ、δ、λ与ξ均由均值、标准偏差、偏斜度与峰度等高度统计特征输入参数确定。
b)如果有界***SB拟合方式无法收敛,则采用Pearson非高斯转换***,以提高非高斯转换***的精度:
利用概率密度函数P满足以下微分方程,根据均值、方差、偏斜度与峰度等高度统计参数获得概率密度函数与高度分布函数。
式中a、c0、c1与c2是由均值、标准偏差、偏斜度与峰度等高度统计特征输入参数确定。
8)非高斯高度序列偏斜度与峰度精度控制与循环迭代策略:
a)对非高斯高度序列Z1进行平移与缩放,获得均值与标准偏差分别为μ与σ的非高斯高度序列Z;
Z1=Z1-mean2(Z1)+μ Z=σ·Z1/std2(Z1)
b)计算非高斯序列Z的偏斜度Skz与峰度Kuz,并判断其精度能否满足要求,若满足精度要求,则完成非高斯粗糙表面的模拟,否则继续执行后面的步骤;
c)利用非高斯序列Z1的快速反傅里叶变换,更新的相角序列Φ(m,n);
d)利用相角序列,计算相角相关序列exp(iΦ),对其进行快速反傅里叶变换更新高斯粗糙表面数字化模拟中的白噪声序列η,则更新的白噪声序列的傅里叶变换A约等于mn·exp(iΦ);
e)重复步骤4)~8),直到非高斯序列Z满足偏斜度与峰度的精度要求。
利用模拟的各向同性高斯粗糙表面高度序列z,其x、y方向的自相关函数均取5μm,x、y方向取点为128×128、表面的均值为0,标准偏差为1.6μm,通过非高斯模拟获得偏斜度Skz与峰度Kuz分别为-1与4的各向同性非高斯序列,对应的非高斯表面形貌及其高度分布如图3所示。产生非高斯粗糙表面的迭代循环次数为2,模拟时间小于0.5秒。在高度统计特征方面:其前四阶矩高度统计特征误差如表3所示,其误差小于1.00%;在自相关函数统计特征方面:如图4所示,产生的高斯粗糙粗糙表面与非高斯粗糙表面对应的x方向与y方向平均自相关函数(已经标准化,最大值为1)接近重合,且与相应理论自相关函数基本吻合,在下部稍有偏离,可以通过增加模拟的点数来提高其吻合程度。
利用模拟的各向异性高斯粗糙表面高度序列z,其x、y方向的自相关函数分别为5μm与50μm,x、y方向取点为512×512、表面的均值为0,标准偏差为1.6μm,通过非高斯模拟获得偏斜度Skz与峰度Kuz分别为1与5的各向异性非高斯序列,对应的非高斯表面形貌及其高度分布如图5所示。产生非高斯粗糙表面的迭代循环次数为3,模拟时间小于0.5秒。在高度统计特征方面:其前四阶矩高度统计特征误差如表4所示,其误差小于1.00%;在自相关函数统计特征方面:如图6所示,由于x方向模拟点数512(反映模拟长度,原来点数为128)与x方向自相关长度5μm的比值明显增大,模拟产生的高斯粗糙粗糙表面与非高斯粗糙表面的x方向平均自相关函数(已经标准化,最大值为1)重合,几乎与相应理论自相关函数曲线重合;由于y方向模拟点数512(反映模拟长度,原来点数为128)与y方向自相关长度50μm的比值稍有减小,导致在自相关函数值较小的情况下,模拟产生的高斯粗糙粗糙表面与非高斯粗糙表面的y方向平均自相关函数曲线与相应的理论自相关曲线有一定的偏离,这是由于模拟点数反映模拟长度,某方向的模拟长度与相应自相关长度比值的减小,会导致该方向的相似波数目变少,导致该方向的统计特征偏离。如要提高自相关函数的统计精度,可以通过增加点数与减小自相关长度。
表3 非高斯表面(βx=βy=5,μ=0,σ=1.6,Skz=-1,Kuz=4,点数128×128)的高度统计前四阶矩仿真精度计算
表4 非高斯表面(βx=5,βy=50,μ=0,σ=1.6,Skz=1,Kuz=5点数512×512)的高度统计前四阶矩仿真精度计算
Claims (6)
1.一种3D任意粗糙表面的数字化模拟方法,其特征在于,包含下述步骤:
1)当任意粗糙表面的偏斜度值等于0且峰度值约等于3时,其表面形貌符合高斯分布,否则,其表面形貌符合非高斯分布;
2)如果步骤1)中粗糙表面形貌为高斯分布,利用三维高斯粗糙表面的模拟方法获得粗糙表面高度序列z,完成粗糙表面的数字化模拟;
3)如果步骤1)中粗糙表面形貌为非高斯分布,利用三维非高斯粗糙表面的模拟方法获得粗糙表面高度序列Z,完成粗糙表面的数字化模拟;
其中,步骤2)中三维高斯粗糙表面的模拟方法,具体包含以下步骤:
(1)通过给定的自相关函数生成离散的自相关函数矩阵R(m,n),其中m、n分别表示x、y方向取点的数目;
(2)对离散的自相关函数矩阵R(m,n)进行快速傅里叶变换,获得高斯粗糙表面的功率谱密度P(m,n);
(3)假设白噪声的功率谱密度常数C等于1,利用获得的功率谱密度,计算传递函数H(m,n);
(4)利用白噪声生成器,获得白噪声序列η(m,n),利用反傅里叶变换生成相角序列Φ;
(5)以步骤(4)获得的相角序列Φ,计算相角相关序列exp(iΦ),并对其进行反傅里叶变换,更新白噪声序列η(m,n)及其傅里叶变换A(m,n);
(6)以步骤(3)获得的传递函数H(m,n),步骤5)获得的白噪声傅里叶变换序列A(m,n),求序列H与序列A的点积运算,并对点积运算结果进行反傅里叶变换,获得高斯表面的高度序列z0(m,n);
(7)以步骤(6)生成的高斯高度序列z0,对其进行高度方向平移与缩放,获得满足均值与标准偏差的高斯高度序列z,完成高斯粗糙表面的数字化模拟。
2.根据权利要求1所述的一种3D任意粗糙表面的数字化模拟方法,所述对z0高度方向进行平移与缩放的方法,具体包含以下步骤:
1)对高度序列z0进行平移,使其均值等于给定μ;
z0=z0-mean2(z0)+μ
其中mean2(z0)为求序列z0均值的函数,μ为高度均值;
2)对经平移后的高度序列z0在高度方向进行缩放,获得标准偏差为σ的高斯高度序
列z,
z=σ·z0/std2(z0)
其中std2(z0)为求序列z0的标准偏差函数。
3.根据权利要求1所述的一种3D任意粗糙表面的数字化模拟方法,所述的三维非高斯粗糙表面的模拟方法,具体包含以下步骤:
1)以三维高斯粗糙表面的模拟方法中所获得的高斯高度序列z,利用高斯序列的标准化方法,获得标准化正态分布序列Z0;
2)以步骤1)中获得的标准化正态分布序列Z0与高度序列的前四阶矩统计特征作为输入,利用Johnson或Pearson转换***进行非高斯变换,获得非高斯序列Z1;
3)利用步骤2)生成的非高斯序列Z1,对其进行平移与缩放,获得满足均值与标准偏差高度统计特征非高斯序列Z,同时,计算其偏斜度与峰度高度统计特征,判断其精度是否满足要求;如满足精度要求,则完成非高斯粗糙表面的数字化模拟,否则,继续执行后面步骤;
4)以步骤3)获得的非高斯序列Z1,利用反傅里叶变换生成相角序列Φ,更新白噪声序列η(m,n)及其傅里叶变换A(m,n);
5)以步骤5)获得序列A与权利要求1所获得的传递函数H(m,n),通过求序列A与H的点积运算,并对点积运算结果进行反傅里叶变换,获得新的高斯表面的高度序列z0(m,n);
6)以步骤5)生成的高斯高度序列z0,利用高斯序列的标准化方法,获得满足均值与标准偏差分别为0与1的新标准化正态分布序列Z0;
7)重复步骤2)~3),直到非高斯序列Z满足偏斜度与峰度的精度要求为止。
4.根据权利要求3所述的一种3D任意粗糙表面的数字化模拟方法,所述高斯序列的标准化方法,具体包含以下步骤:
1)对高度序列z0进行平移,使其均值等于0;
z0=z0-mean2(z0)
2)对经平移后的高度序列z0在高度方向进行缩放,获得标准偏差1的标准高斯序列,Z0=z0/std2(z0)。
5.根据权利要求3所述的一种3D任意粗糙表面的数字化模拟方法,所述的利用Johnson或Pearson转换***进行非高斯变换,具体步骤如下:
1)优先采用Johnson非高斯转换***;
对数正态拟合方式SL:Z1=ξ+λexp[(Z0-γ)/δ] (ξ<η)
有界***拟合方式SB:z=ξ+λe(η-γ)/δ/[1+e(η-γ)/δ] (ξ<η<ξ+λ)
无界***拟合方式SU:z=ξ+λsinh[(η-γ)/δ]
γ、δ为形状参数,λ、ξ分别为比例参数与位置参数,γ、δ、λ与ξ均由均值、标准偏差、偏斜度与峰度等高度统计特征输入参数确定;
2)如果有界***SB拟合方式无法收敛,则采用Pearson非高斯转换***,以提高非高斯转换***的精度:
利用概率密度函数P满足以下微分方程,根据均值、方差、偏斜度与峰度等高度统计参数获得概率密度函数与高度分布函数:
式中a、c0、c1与c2是由均值、标准偏差、偏斜度与峰度等高度统计特征输入参数确定。
6.根据权利要求3所述的一种3D任意粗糙表面的数字化模拟方法,所述的利用反傅里叶变换生成相角序列Φ,更新白噪声序列η(m,n)及其傅里叶变换A(m,n),具体包含以下步骤:
1)以三维非高斯粗糙表面的模拟方法中所获得的非高斯序列Z1,通过反傅里叶变换,生成新的白噪声相角序列Φ(m,n);
2)以步骤1)获得的相角序列Φ(m,n),计算相角相关序列exp(iΦ),并对其进行反傅里叶变换,更新白噪声序列η(m,n)及其傅里叶变换A(m,n);
3)更新高斯粗糙表面数字化模拟中的白噪声序列η(m,n)的傅里叶变换序列A(m,n);
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110421316.4A CN102609560B (zh) | 2011-12-14 | 2011-12-14 | 一种3d任意粗糙表面的数字化模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110421316.4A CN102609560B (zh) | 2011-12-14 | 2011-12-14 | 一种3d任意粗糙表面的数字化模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102609560A CN102609560A (zh) | 2012-07-25 |
CN102609560B true CN102609560B (zh) | 2015-07-01 |
Family
ID=46526928
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110421316.4A Expired - Fee Related CN102609560B (zh) | 2011-12-14 | 2011-12-14 | 一种3d任意粗糙表面的数字化模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102609560B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3372954A1 (en) * | 2017-03-09 | 2018-09-12 | Kabushiki Kaisha Toyota Chuo Kenkyusho | Arithmetic unit, method and program for computing a contact condition between two surfaces |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105956300B (zh) * | 2016-05-10 | 2019-05-07 | 清华大学 | 一种分层粗糙表面的数字化重构方法 |
CN106021661B (zh) * | 2016-05-10 | 2019-03-22 | 清华大学 | 一种表面分析、仿真与重构*** |
CN106021660B (zh) * | 2016-05-10 | 2018-12-28 | 清华大学 | 一种分层粗糙表面的分析方法 |
JP6597716B2 (ja) * | 2017-03-09 | 2019-10-30 | 株式会社豊田中央研究所 | 演算装置 |
CN110622046B (zh) * | 2017-05-15 | 2022-07-08 | 日本电气硝子株式会社 | 透明物品以及透明物品的制造方法 |
CN107063162A (zh) * | 2017-06-02 | 2017-08-18 | 浙江水利水电学院 | 一种弹性接触表面形貌形成方法 |
CN109359333B (zh) * | 2018-09-12 | 2021-09-24 | 大连理工大学 | 一种包含多尺度形貌特征的体模型构建方法 |
CN112797891B (zh) * | 2020-12-28 | 2021-11-05 | 大连理工大学 | 基于传递函数的白光扫描干涉测量法高频形貌补偿方法 |
CN114611356B (zh) * | 2022-03-11 | 2024-07-09 | 西安交通大学 | 用于激光超声仿真分析的增材制造零件粗糙表面建模方法 |
-
2011
- 2011-12-14 CN CN201110421316.4A patent/CN102609560B/zh not_active Expired - Fee Related
Non-Patent Citations (5)
Title |
---|
Numerical generation of arbitrarily oriented non-Gaussian;Vasilios Bakolas;《Wear》;20030331;第254卷(第5-6期);全文 * |
Simulation of rough surfaces with FFT;Jiunn-Jong Wu;《Tribology International》;20000131;第33卷(第1期);全文 * |
利用FFT实现对非高斯随机粗糙表面的模拟;宋俊杰 等;《西安工业大学学报》;20080229;第28卷(第1期);全文 * |
粗糙表面计算机模拟;陈辉 等;《润滑与密封》;20061031(第10期);第1、2节 * |
非高斯随机粗糙表面的数字模拟;田爱玲 等;《***仿真学报》;20090531;第21卷(第10期);全文 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3372954A1 (en) * | 2017-03-09 | 2018-09-12 | Kabushiki Kaisha Toyota Chuo Kenkyusho | Arithmetic unit, method and program for computing a contact condition between two surfaces |
Also Published As
Publication number | Publication date |
---|---|
CN102609560A (zh) | 2012-07-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102609560B (zh) | 一种3d任意粗糙表面的数字化模拟方法 | |
CN101413824B (zh) | 一种基于随机传声器阵列的运动物体声场测量方法 | |
CN104730518B (zh) | 一种基于高斯拟合的雷达多普勒谱估计海面流场的方法 | |
CN103954953B (zh) | 一种基于数据驱动的机载激光雷达盲源误差补偿方法 | |
CN101438184B (zh) | 一种跟踪移动电子设备的状态的方法 | |
CN101604019B (zh) | 一种海洋环境与声场不确实性的表征和传递的快速计算方法 | |
CN101980044B (zh) | 未知测量噪声分布下的多目标跟踪方法 | |
CN104569911B (zh) | 一种obu定位方法、rsu及etc*** | |
CN103759726B (zh) | 一类循环平稳泊松信号快速模拟方法及其硬件*** | |
CN104035083A (zh) | 一种基于量测转换的雷达目标跟踪方法 | |
CN103472450A (zh) | 基于压缩感知的非均匀空间构形分布式sar动目标三维成像方法 | |
CN104052557A (zh) | 一种Nakagami复衰落信道建模方法 | |
CN102253385A (zh) | 基于合成孔径雷达图像和内波模型的海洋内波预测方法 | |
CN106583509A (zh) | 不规则零件的折弯工艺 | |
CN109038675A (zh) | 基于风电波动多尺度分解的建模方法 | |
CN103699810A (zh) | 一种粗糙面微波段双向反射分布函数的建模方法 | |
CN103575298A (zh) | 基于自调节的ukf失准角初始对准方法 | |
CN101854216A (zh) | 一种基于修正电离层信道模型进行信道对等性研究的方法 | |
CN103852648A (zh) | 获取空间电磁强度数据的方法 | |
Jeong et al. | A shipyard simulation system using the process-centric simulation modeling methodology: Case study of the simulation model for the shipyard master plan validation | |
CN106342186B (zh) | 地磁导航定位匹配误差确定方法 | |
CN109001670B (zh) | 一种联合时差和角度的分布式无源定位方法及装置 | |
CN101598788B (zh) | 合成孔径声纳信号的快速仿真方法 | |
CN110532659A (zh) | 一种考虑接触表面形貌的有限元节点建模方法 | |
CN103294522A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150701 Termination date: 20201214 |
|
CF01 | Termination of patent right due to non-payment of annual fee |