CN102609560B - Digitalized simulation method of any three-dimensional (3D) rough surface - Google Patents
Digitalized simulation method of any three-dimensional (3D) rough surface 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
The invention discloses a digitalized simulation method of any three-dimensional (3D) rough surface, which improves simulation efficiency and simulation accuracy of statistical parameters. Random initial phase angle sequences are obtained through inverse Fourier transform of white noise sequences, and novel white noise and Fourier transform are obtained by utilizing the initial phase angle sequences. Processing including dispersion and Fourier transform and the like is performed through appointed autocorrelation functions to obtain power spectral density and transfer functions of Gaussian rough surface height sequences. The simulation of the Gaussian rough surface height sequences is finished by utilizing frequency domain dot product and a method for obtaining the inverse Fourier transform. On the basis, non-Gaussian rough surfaces are generated through height distribution statistical characteristic parameters including appointed skewness degree, kurtosis and the like and by utilizing combined Pearson and Johnson non-Gaussian transition systems. If the simulation accuracy of the skewness degree and the kurtosis is below standard, the phase angle sequences and the Fourier transform of the white noise are updated to re-perform Gaussian filtering and non-Gaussian transition till the given accuracy requirement is met.
Description
Technical field
The present invention relates to the digital simulation method of any rough surface of a kind of 3D microcosmic, more particularly, the present invention, according to Fourth-order moment statistical nature and horizontal direction autocorrelation function statistical law before the height distribution of connection table millet cake, the rough surface with identical statistical nature can be built efficiently and accurately.
Background technology
At a high speed, improving constantly of requiring of the development of precision processing technology and production reliability, promote the fast development in the fields such as the processing and manufacturing of rough surface, detection, digitized simulation and performance prediction.
In recent years, research both at home and abroad in the digital simulation field of rough surface, mainly concentrate on the rough surface aspect meeting Gaussian distribution for roughness height, in the digitized simulation of non-gaussian rough surface, can be successfully non-gaussian data by the data transformations of Gaussian distribution, occur based on non-gaussian converting systems such as Johnson, Pearson.But in the middle of the digitized simulation of non-gaussian rough surface, non-gaussian sequence after non-gaussian conversion also needs to carry out filtering, introduce new error to statistical nature parameters such as the measure of skewness of simulated series, kurtosis, have also been changed the simulated domain of the correspondence of measure of skewness-kurtosis simultaneously.Thus cause the rough surface of some statistical nature not generate corresponding non-gaussian sequence by converting system conversions such as Johnson, Pearson to filtering, i.e. enable transition, also can introduce new error at filtering stage, affect the precision of the non-gaussian surface digitizing simulation of filtering.
Consider the impact of filtering (filter factor is h (k, l)), the measure of skewness S of the non-gaussian sequence rough surface generated in the non-gaussian converting systems such as Johnson, Pearson
k η, kurtosis K
u ηwith the measure of skewness S of filtered rough surface
kZ, kurtosis K
uZthere is following relation:
Even if target deflection degree and kurtosis meet non-gaussian change necessary condition:
but after considering the influencing each other of filter factor and autocorrelation function, S
k η, K
u ηlikely no longer meet the necessary condition of non-gaussian conversion:
thus cause not generating non-gaussian sequence by non-gaussian converting systems such as Johnson or Pearson; I.e. enable transition, also can introduce new error at filtering stage, affects the precision of the non-gaussian surface digitizing simulation of filtering.
In reality, different job operations can produce the rough surface with different statistical nature, the carrier that the digitized simulation model of rough surface is analyzed as interface properties, and accuracy of its simulation directly affects the contact performance prediction of corresponding faying face.The sequence of phase angles generation white noise that this research and utilization is certain, combined with the non-gaussian converting system of Gaussian surface high degree of sequence by the gaussian filtering of white noise, form the non-gaussian high degree of sequence meeting height distribution statistics feature and autocorrelation function, the digitized simulation method of Gauss or non-gaussian rough surface is that interface properties research is laid a good foundation.
Summary of the invention
Fundamental purpose of the present invention is to provide a kind of digitized simulation method that high-efficiency high-accuracy realizes any rough surface, decrease the filtering error link in rough surface digitized simulation, before meeting autocorrelation function, highly distribution the statistical nature such as Fourth-order moment prerequisite under, improve and simulate the precision of rough surface height distribution statistics feature and the efficiency of digitized simulation.
The present invention is achieved through the following technical solutions, and key step comprises:
A digitized simulation method for any rough surface of 3D, is characterized in that, comprises following step:
1) when the deflection angle value of any rough surface approximates 3 and kurtosis value approximates 0, its surface topography meets Gaussian distribution, otherwise its surface topography meets non-gaussian distribution;
2) if rough surface morphology is Gaussian distribution in step 1), utilize the analogy method of three-dimensional Gaussian rough surface to obtain rough surface high degree of sequence z, complete the digitized simulation of rough surface;
3) if rough surface morphology is non-gaussian distribution in step 1), utilize the analogy method of three-dimensional non-gaussian rough surface to obtain rough surface high degree of sequence z, complete the digitized simulation of rough surface.
The analogy method of described three-dimensional Gaussian rough surface, specifically comprises following steps:
1) generate discrete autocorrelation function matrix R (m, n) by given autocorrelation function, wherein m, n represent that number is a little got in x, y direction respectively;
2) Fast Fourier Transform (FFT) is carried out to discrete autocorrelation function matrix R (m, n), obtain the power spectrum density P (m, n) of Gauss's rough surface;
3) suppose that the power spectrum density constant C of white noise equals 1, utilize the power spectrum density obtained, calculation of transfer function H (m, n);
4) utilize white noise maker, obtain white noise sequence η (m, n), utilize inversefouriertransform to generate sequence of phase angles Φ;
5) with the sequence of phase angles Φ that step 4) obtains, calculate phase angle correlated series exp (i Φ), and inversefouriertransform is carried out to it, upgrade white noise sequence η (m, n) and Fourier transform A (m, n) thereof;
6) with the transfer function H (m, n) that step 3) obtains, the white noise Fourier transform sequence A (m that step 5) obtains, n), ask the dot-product operation of sequence H and sequence A, and inversefouriertransform is carried out to dot-product operation result, obtain the high degree of sequence z of Gaussian surface
0(m, n);
7) with Gauss's high degree of sequence z that step 6) generates
0, short transverse translation and convergent-divergent are carried out to it, obtain the Gauss's high degree of sequence z meeting average and standard deviation, complete the digitized simulation of Gauss's rough surface.
The analogy method of described three-dimensional non-gaussian rough surface, specifically comprises following steps:
1) with the Gauss's high degree of sequence z obtained, utilize the standardized method of gaussian sequence, carry out translation and convergent-divergent in short transverse, obtain standardized normal distribution series Z
0;
2) with the standardized normal distribution series Z obtained in step 1)
0with Fourth-order moment statistical nature before high degree of sequence as input, utilize Johnson or Pearson converting system to carry out non-gaussian conversion, obtain non-gaussian sequence Z
1;
3) step 2 is utilized) the non-gaussian sequence Z that generates
1, translation and convergent-divergent are carried out to it, obtain and meet average and standard deviation height statistical nature non-gaussian sequence Z, meanwhile, calculate its measure of skewness and kurtosis height statistical nature, judge whether its precision meets the demands; As met accuracy requirement, then complete the digitized simulation of non-gaussian rough surface, otherwise, continue to perform later step;
4) with the non-gaussian sequence Z that step 3) obtains
1, utilize inversefouriertransform to generate sequence of phase angles Φ, upgrade white noise sequence η (m, n) and Fourier transform A (m, n) thereof;
5) with the transfer function H (m, n) that step 5) acquisition sequence A and claim 1 obtain, by asking the dot-product operation of sequence A and H, and inversefouriertransform being carried out to dot-product operation result, obtaining the high degree of sequence z of new Gaussian surface
0(m, n);
6) with Gauss's high degree of sequence z that step 5) generates
0, utilize the standardized method of gaussian sequence, obtain and meet the new standardized normal distribution sequence Z that average and standard deviation are respectively 0 and 1
0;
7) step 2 is repeated) ~ 3), until non-gaussian sequence Z meets the accuracy requirement of measure of skewness and kurtosis.
Described to z
0short transverse carries out the method for translation and convergent-divergent, specifically comprises following steps:
1) to high degree of sequence z
0carry out translation, make its average equal given μ;
z
0=z
0-mean2(z
0)+μ
2) to the high degree of sequence z after translation
0carry out convergent-divergent in short transverse, obtain Gauss's high degree of sequence z that standard deviation is σ.
z=σ·z
0/std2(z
0)
The standardized method of described gaussian sequence, specifically comprises following steps:
1) to high degree of sequence z or z
0carry out translation, make its average equal 0;
Z
0=z
0-mean2 (z
0) or z=z-mean2 (z)
2) to high degree of sequence z or z after translation
0carry out convergent-divergent in short transverse, obtain the standard gaussian sequence of standard deviation 1.
Z
0=σ z
0/ std2 (z
0) or Z
0=σ z/std2 (z)
Described Johnson or the Pearson converting system that utilizes carries out non-gaussian conversion, and concrete steps are as follows:
1) preferential employing Johnson non-gaussian converting system;
Lognormality fit approach S
l: Z
1=ξ+λ exp [(Z
0-γ) δ] (ξ < η)
Bounded system fit approach S
b: z=ξ+λ e
(η-γ) δ/ [1+e
(η-γ)/δ] (ξ < η < ξ+λ)
Boundless system fit approach S
u: z=ξ+λ sinh [(η-γ)/δ]
γ, δ are form parameter, and λ, ξ are respectively scale parameter and location parameter, and γ, δ, λ and ξ determine by average, standard deviation, measure of skewness and kurtosis equal altitudes statistical nature input parameter.
2) if bounded system S
bfit approach cannot be restrained, then adopt Pearson non-gaussian converting system, to improve the precision of non-gaussian converting system:
Utilize probability density function P to meet the following differential equation, obtain probability density function and height distribution function according to average, variance, measure of skewness and kurtosis equal altitudes statistical parameter.
A, c in formula
0, c
1with c
2determined by average, standard deviation, measure of skewness and kurtosis equal altitudes statistical nature input parameter.
The described inversefouriertransform that utilizes generates sequence of phase angles Φ, upgrades white noise sequence η (m, n) and Fourier transform A (m, n) thereof, specifically comprises following steps:
1) with the non-gaussian sequence Z that claim 3 obtains
1, by inversefouriertransform, generate new white noise sequence of phase angles Φ (m, n);
2) with the sequence of phase angles Φ (m, n) that step 1) obtains, phase angle correlated series exp (i Φ) is calculated, and inversefouriertransform is carried out to it, upgrade white noise sequence η (m, n) and Fourier transform A (m, n) thereof;
3) the Fourier transform sequence A (m, n) of the white noise sequence η (m, n) in Gauss's rough surface digitized simulation is upgraded.
The present invention proposes a kind of digitized simulation method of any 3D rough surface, start with from gaussian filtering, the conversion of high precision non-gaussian, covert loop iteration technology three aspects, angle, overcome prior art bottleneck, improve simulation precision and the simulation precision of 3D micro-roughened surface, ensure that the feasibility of the digitized simulation on any non-gaussian surface of whole measure of skewness-kurtosis plane respective heights statistical nature, accuracy and high efficiency to a certain extent.The contact behaviours that can be faying face provides the pattern model of high precision statistical nature.
Accompanying drawing explanation
Fig. 1 isotropy Gaussian surface and the (β that highly distributes thereof
x=β
y=5, μ=0, σ=1.6, count 128 × 128)
Fig. 2 anisotropic Gaussian surface and the (β that highly distributes thereof
x=5, β
y=50, μ=0, σ=1.6, count 512 × 512)
Fig. 3 isotropy non-gaussian rough surface and the (β that highly distributes thereof
x=β
y=5, μ=0, σ=1.6, S
kz=-1, K
uz=4, count 128 × 128)
X, y orientation average autocorrelation function of Fig. 4 isotropy rough surface and the contrast of theoretical autocorrelation function: (a) x direction autocorrelation function (β
x=5), (b) y direction autocorrelation function (β
x=5)
Fig. 5 anisotropy non-gaussian rough surface and the (β that highly distributes thereof
x=5, β
y=50, μ=0, σ=1.6, Skz=1, K
uz=5 count 512 × 512)
X, y orientation average autocorrelation function of Fig. 6 anisotropic rough surface and the contrast of theoretical autocorrelation function: (a) x direction autocorrelation function (β
x=5), (b) y direction autocorrelation function (β
y=50)
Embodiment
To meet the 3D Gauss of statistical nature or the digitized simulations of non-gaussian rough surface such as certain autocorrelation function, average, standard deviation, measure of skewness and kurtosis, the digitized simulation method of any 3D rough surface of the present invention is described.
1.
the digitized simulation of Gauss's rough surface:
1) given autocorrelation function is utilized to generate related function discrete matrix R (m, n)
A) suppose that autocorrelation function is that (σ is the standard deviation of gaussian sequence to f (x, y), β
x, β
ybe respectively the auto-correlation length of X, Y-direction);
B) the discrete spacing of appointment X, Y-direction is 1 μm, and carries out discrete to autocorrelation function in-m/2≤x≤m/2 with-n/2≤y≤n/2 region, obtains the discrete matrix R (m+1, n+1) of autocorrelation function;
C) the capable repetition of the m/2+1 in gained R (m+1, n+1) sequence, m/2+2, deletes wherein a line; N-th/2+1 and n/2+2 column weight is multiple, deletes wherein row; Obtain autocorrelation function matrix R (m, n) with simulation surface elevation sequence with identical ranks number.
2) power spectrum density and the transport function of Gaussian surface is asked according to autocorrelation function;
A) Fourier transform is carried out to autocorrelation function discrete matrix, and solve the power spectrum density P of Gauss's rough surface;
B) be constant C according to the power spectrum density of white noise, assuming that C=1, then transport function is;
3) generation of white noise sequence and Fourier transform.
A) white noise maker randn (m, n) is utilized to produce white noise sequence η (m, n);
B) the Fourier transform A (m, n) of white noise sequence η is asked:
4) method of frequency domain dot product (A.*H) negate Fourier transform is utilized to obtain the high degree of sequence of Gaussian surface;
5) according to assigned altitute average μ, standard deviation and acquisition sequence z
0actual average mean2 (z
0) (mean2 (z
0) ask the function of serial mean) and with the relation of standard deviation std2 (z) (std2 is for asking sequence criteria departure function), to Gauss's high degree of sequence z
0carry out short transverse translation and convergent-divergent, eliminate Power Spectrum of White Noise density length C and be not equal to 1 impact on high degree of sequence statistical nature, obtain and meet Gauss's high degree of sequence z that average and distribution of standard deviation are μ and σ, complete the digitized simulation of Gauss's rough surface.
z
0=z
0-mean2(z
0)+μ z=σ·z
0/std2(z
0)
Utilize above-mentioned given autocorrelation function, the autocorrelation function in x, y direction all gets 5 μm, the x on simulation surface, it is a little 128 × 128 that y gets in direction, the average on surface is 0, standard deviation is 1.6 μm, as shown in Figure 1, its front Fourth-order moment height statistical nature error is as shown in table 1 for the isotropy Gaussian surface pattern that simulation obtains and highly distribution.The autocorrelation function in x, y direction gets 5 μm, 50 μm respectively, and the Gaussian surface pattern that simulation obtains shows anisotropy, has certain texture.As shown in Figure 2, its front Fourth-order moment height statistical nature error is as shown in table 2 for the isotropy Gaussian surface pattern that simulation obtains and highly distribution.
Table 1 isotropy Gauss rough surface (β
x=β
y=5, μ=0, σ=1.6, count 128 × 128) height statistics before Fourth-order moment simulation accuracy
Table 2 anisotropic Gaussian rough surface (β
x=5, β
y=50, μ=0, σ=1.6, count 512 × 512) height statistics before Fourth-order moment simulation accuracy calculate
2. the digitized simulation of non-gaussian rough surface:
6) to the z translation of Gauss's high degree of sequence and the convergent-divergent of above-mentioned simulation, acquisition average and standard deviation are respectively the standardized normal distribution series Z of 0 and 1
0;
z=z-mean2(z) Z
0=z/std2(z)
7) adopt Johnson or Pearson non-gaussian converting system to carry out non-gaussian conversion, generate non-gaussian sequence Z
1;
A) preferential employing Johnson non-gaussian converting system;
Lognormality fit approach S
l: Z
1=ξ+λ exp [(Z
0-γ) δ] (ξ < η)
Bounded system fit approach S
b:
Boundless system fit approach S
u: Z
1=ξ+λ sinh [(Z
0-γ)/δ]
γ, δ are form parameter, and λ, ξ are respectively scale parameter and location parameter, and γ, δ, λ and ξ determine by average, standard deviation, measure of skewness and kurtosis equal altitudes statistical nature input parameter.
If b) bounded system S
bfit approach cannot be restrained, then adopt Pearson non-gaussian converting system, to improve the precision of non-gaussian converting system:
Utilize probability density function P to meet the following differential equation, obtain probability density function and height distribution function according to average, variance, measure of skewness and kurtosis equal altitudes statistical parameter.
A, c in formula
0, c
1with c
2determined by average, standard deviation, measure of skewness and kurtosis equal altitudes statistical nature input parameter.
8) non-gaussian high degree of sequence measure of skewness and kurtosis precision controlling and loop iteration strategy:
A) to non-gaussian high degree of sequence Z
1carry out translation and convergent-divergent, acquisition average and standard deviation are respectively the non-gaussian high degree of sequence Z of μ and σ;
Z
1=Z
1-mean2(Z
1)+μ Z=σ·Z
1/std2(Z
1)
B) the measure of skewness S of non-gaussian sequence Z is calculated
kzwith kurtosis K
uz, and judge that can its precision meet the demands, if meet accuracy requirement, then complete the simulation of non-gaussian rough surface, otherwise continue to perform step below;
C) non-gaussian sequence Z is utilized
1quick inversefouriertransform, the sequence of phase angles Φ (m, n) of renewal;
D) sequence of phase angles is utilized, calculate phase angle correlated series exp (i Φ), carry out the white noise sequence η in quick inversefouriertransform renewal Gauss rough surface digitized simulation to it, then the Fourier transform A of the white noise sequence upgraded approximates mnexp (i Φ);
E) step 4 is repeated) ~ 8), until non-gaussian sequence Z meets the accuracy requirement of measure of skewness and kurtosis.
Utilize the isotropy Gauss rough surface high degree of sequence z of simulation, the autocorrelation function in its x, y direction all gets 5 μm, and it is a little 128 × 128 that x, y get in direction, the average on surface is 0, and standard deviation is 1.6 μm, obtains measure of skewness S by non-gaussian simulation
kzwith kurtosis K
uzbe respectively the isotropy non-gaussian sequence of-1 and 4, corresponding non-gaussian surface topography and highly distribution are as shown in Figure 3.The iterative loop number of times producing non-gaussian rough surface is 2, and simulated time is less than 0.5 second.In height statistical nature: its front Fourth-order moment height statistical nature error is as shown in table 3, and its error is less than 1.00%; In autocorrelation function statistical nature: as shown in Figure 4, the x direction that the coarse rough surface of Gauss produced is corresponding with non-gaussian rough surface and (the standardization of y orientation average autocorrelation function, maximal value is 1) close to overlapping, and it is substantially identical with corresponding theory autocorrelation function, being slightly offset in bottom, its degree of agreement can being improved by increasing counting of simulation.
Utilize the anisotropic Gaussian rough surface high degree of sequence z of simulation, the autocorrelation function in its x, y direction is respectively 5 μm and 50 μm, it is a little 512 × 512 that x, y get in direction, the average on surface is 0, and standard deviation is 1.6 μm, obtains measure of skewness S by non-gaussian simulation
kzwith kurtosis K
uzbe respectively the anisotropy non-gaussian sequence of 1 and 5, corresponding non-gaussian surface topography and highly distribution are as shown in Figure 5.The iterative loop number of times producing non-gaussian rough surface is 3, and simulated time is less than 0.5 second.In height statistical nature: its front Fourth-order moment height statistical nature error is as shown in table 4, and its error is less than 1.00%, in autocorrelation function statistical nature: as shown in Figure 6, because 512(reflection simulation length of counting is simulated in x direction, originally counted is 128) obviously increase with the ratio of x direction auto-correlation length 5 μm, x orientation average autocorrelation function (the standardization of the coarse rough surface of Gauss that simulation produces and non-gaussian rough surface, maximal value is 1) overlap, almost with corresponding theory autocorrelation function curve co-insides, because 512(reflection simulation length of counting is simulated in y direction, originally counted is 128) slightly reduce with the ratio of y direction auto-correlation length 50 μm, cause when auto-correlation function value is less, the coarse rough surface of Gauss of simulation generation has certain departing to the y orientation average autocorrelation function curve of non-gaussian rough surface and corresponding theoretical autocorrelator trace, this counts due to simulation to reflect simulation length, the simulation length in certain direction and the reduction of corresponding auto-correlation lenth ratio, the similar wave number order in this direction can be caused to tail off, the statistical nature in this direction is caused to depart from.As the statistical precision of autocorrelation function will be improved, can be counted by increase and reduce auto-correlation length.
Table 3 non-gaussian surface (β
x=β
y=5, μ=0, σ=1.6, S
kz=-1, K
uz=4, count 128 × 128) height statistics before Fourth-order moment simulation accuracy calculate
Table 4 non-gaussian surface (β
x=5, β
y=50, μ=0, σ=1.6, S
kz=1, K
uz=5 count 512 × 512) height statistics before Fourth-order moment simulation accuracy calculate
Claims (6)
1. a digitized simulation method for any rough surface of 3D, is characterized in that, comprises following step:
1) when the deflection angle value of any rough surface equals 0 and kurtosis value approximates 3, its surface topography meets Gaussian distribution, otherwise its surface topography meets non-gaussian distribution;
2) if step 1) in rough surface morphology be Gaussian distribution, utilize the analogy method of three-dimensional Gaussian rough surface to obtain rough surface high degree of sequence z, complete the digitized simulation of rough surface;
3) if step 1) in rough surface morphology be non-gaussian distribution, utilize the analogy method of three-dimensional non-gaussian rough surface to obtain rough surface high degree of sequence Z, complete the digitized simulation of rough surface;
Wherein, step 2) in the analogy method of three-dimensional Gaussian rough surface, specifically comprise following steps:
(1) generate discrete autocorrelation function matrix R (m, n) by given autocorrelation function, wherein m, n represent that number is a little got in x, y direction respectively;
(2) Fast Fourier Transform (FFT) is carried out to discrete autocorrelation function matrix R (m, n), obtain the power spectrum density P (m, n) of Gauss's rough surface;
(3) suppose that the power spectrum density constant C of white noise equals 1, utilize the power spectrum density obtained, calculation of transfer function H (m, n);
(4) utilize white noise maker, obtain white noise sequence η (m, n), utilize inversefouriertransform to generate sequence of phase angles Φ;
(5) with the sequence of phase angles Φ that step (4) obtains, calculate phase angle correlated series exp (i Φ), and inversefouriertransform is carried out to it, upgrade white noise sequence η (m, and Fourier transform A (m, n) n);
(6) with the transfer function H (m that step (3) obtains, n), step 5) the white noise Fourier transform sequence A (m that obtains, n), ask the dot-product operation of sequence H and sequence A, and inversefouriertransform is carried out to dot-product operation result, obtain the high degree of sequence z of Gaussian surface
0(m, n);
(7) with Gauss's high degree of sequence z that step (6) generates
0, short transverse translation and convergent-divergent are carried out to it, obtain the Gauss's high degree of sequence z meeting average and standard deviation, complete the digitized simulation of Gauss's rough surface.
2. the digitized simulation method of any rough surface of a kind of 3D according to claim 1, described to z
0short transverse carries out the method for translation and convergent-divergent, specifically comprises following steps:
1) to high degree of sequence z
0carry out translation, make its average equal given μ;
z
0=z
0-mean2(z
0)+μ
Wherein mean2 (z
0) for asking sequence z
0the function of average, μ is height average;
2) to the high degree of sequence z after translation
0carry out convergent-divergent in short transverse, obtain Gauss's height sequence that standard deviation is σ
Row z,
z=σ·z
0/std2(z
0)
Wherein std2 (z
0) for asking sequence z
0standard deviation function.
3. the digitized simulation method of any rough surface of a kind of 3D according to claim 1, the analogy method of described three-dimensional non-gaussian rough surface, specifically comprises following steps:
1) with the Gauss's high degree of sequence z obtained in the analogy method of three-dimensional Gaussian rough surface, utilize the standardized method of gaussian sequence, obtain standardized normal distribution series Z
0;
2) with step 1) the middle standardized normal distribution series Z obtained
0with Fourth-order moment statistical nature before high degree of sequence as input, utilize Johnson or Pearson converting system to carry out non-gaussian conversion, obtain non-gaussian sequence Z
1;
3) step 2 is utilized) the non-gaussian sequence Z that generates
1, translation and convergent-divergent are carried out to it, obtain and meet average and standard deviation height statistical nature non-gaussian sequence Z, meanwhile, calculate its measure of skewness and kurtosis height statistical nature, judge whether its precision meets the demands; As met accuracy requirement, then complete the digitized simulation of non-gaussian rough surface, otherwise, continue to perform later step;
4) with step 3) the non-gaussian sequence Z that obtains
1, utilize inversefouriertransform to generate sequence of phase angles Φ, upgrade white noise sequence η (m, n) and Fourier transform A (m, n) thereof;
5) with step 5) obtain the transfer function H (m, n) that sequence A and claim 1 obtain, by asking the dot-product operation of sequence A and H, and inversefouriertransform being carried out to dot-product operation result, obtaining the high degree of sequence z of new Gaussian surface
0(m, n);
6) with step 5) Gauss's high degree of sequence z of generating
0, utilize the standardized method of gaussian sequence, obtain and meet the new standardized normal distribution sequence Z that average and standard deviation are respectively 0 and 1
0;
7) step 2 is repeated) ~ 3), until non-gaussian sequence Z meets the accuracy requirement of measure of skewness and kurtosis.
4. the digitized simulation method of any rough surface of a kind of 3D according to claim 3, the standardized method of described gaussian sequence, specifically comprises following steps:
1) to high degree of sequence z
0carry out translation, make its average equal 0;
z
0=z
0-mean2(z
0)
2) to the high degree of sequence z after translation
0carry out convergent-divergent in short transverse, obtain the standard gaussian sequence of standard deviation 1, Z
0=z
0/ std2 (z
0).
5. the digitized simulation method of any rough surface of a kind of 3D according to claim 3, described Johnson or the Pearson converting system that utilizes carries out non-gaussian conversion, and concrete steps are as follows:
1) preferential employing Johnson non-gaussian converting system;
Lognormality fit approach S
l: Z
1=ξ+λ exp [(Z
0-γ)/δ] (ξ < η)
Bounded system fit approach S
b: z=ξ+λ e
(η-γ)/δ/ [1+e
(η-γ)/δ] (ξ < η < ξ+λ)
Boundless system fit approach S
u: z=ξ+λ sinh [(η-γ)/δ]
γ, δ are form parameter, and λ, ξ are respectively scale parameter and location parameter, and γ, δ, λ and ξ determine by average, standard deviation, measure of skewness and kurtosis equal altitudes statistical nature input parameter;
2) if bounded system S
bfit approach cannot be restrained, then adopt Pearson non-gaussian converting system, to improve the precision of non-gaussian converting system:
Utilize probability density function P to meet the following differential equation, obtain probability density function and height distribution function according to average, variance, measure of skewness and kurtosis equal altitudes statistical parameter:
A, c in formula
0, c
1with c
2determined by average, standard deviation, measure of skewness and kurtosis equal altitudes statistical nature input parameter.
6. the digitized simulation method of any rough surface of a kind of 3D according to claim 3, the described inversefouriertransform that utilizes generates sequence of phase angles Φ, upgrades white noise sequence η (m, n) and Fourier transform A (m thereof, n), specifically following steps are comprised:
1) with the non-gaussian sequence Z obtained in the analogy method of three-dimensional non-gaussian rough surface
1, by inversefouriertransform, generate new white noise sequence of phase angles Φ (m, n);
2) with step 1) the sequence of phase angles Φ (m, n) that obtains, calculate phase angle correlated series exp (i Φ), and inversefouriertransform is carried out to it, upgrade white noise sequence η (m, n) and Fourier transform A (m, n) thereof;
3) the Fourier transform sequence A (m, n) of the white noise sequence η (m, n) in Gauss's rough surface digitized simulation is upgraded;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110421316.4A CN102609560B (en) | 2011-12-14 | 2011-12-14 | Digitalized simulation method of any three-dimensional (3D) rough surface |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110421316.4A CN102609560B (en) | 2011-12-14 | 2011-12-14 | Digitalized simulation method of any three-dimensional (3D) rough surface |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102609560A CN102609560A (en) | 2012-07-25 |
CN102609560B true CN102609560B (en) | 2015-07-01 |
Family
ID=46526928
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110421316.4A Expired - Fee Related CN102609560B (en) | 2011-12-14 | 2011-12-14 | Digitalized simulation method of any three-dimensional (3D) rough surface |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102609560B (en) |
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 |
---|---|---|---|---|
CN106021661B (en) * | 2016-05-10 | 2019-03-22 | 清华大学 | A kind of surface analysis, emulation and reconfiguration system |
CN106021660B (en) * | 2016-05-10 | 2018-12-28 | 清华大学 | A kind of analysis method being layered rough surface |
CN105956300B (en) * | 2016-05-10 | 2019-05-07 | 清华大学 | A kind of digital reconstruction method being layered rough surface |
JP6597716B2 (en) * | 2017-03-09 | 2019-10-30 | 株式会社豊田中央研究所 | Arithmetic unit |
WO2018212146A1 (en) * | 2017-05-15 | 2018-11-22 | 日本電気硝子株式会社 | Transparent product and method for producing transparent product |
CN107063162A (en) * | 2017-06-02 | 2017-08-18 | 浙江水利水电学院 | A kind of resilient contact surfaces morphogenesis method |
CN109359333B (en) * | 2018-09-12 | 2021-09-24 | 大连理工大学 | Body model construction method containing multi-scale morphology features |
CN112797891B (en) * | 2020-12-28 | 2021-11-05 | 大连理工大学 | High-frequency morphology compensation method of white light scanning interferometry based on transfer function |
CN114611356B (en) * | 2022-03-11 | 2024-07-09 | 西安交通大学 | Modeling method for rough surface of additive manufactured part for laser ultrasonic simulation analysis |
-
2011
- 2011-12-14 CN CN201110421316.4A patent/CN102609560B/en 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 (en) | 2012-07-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102609560B (en) | Digitalized simulation method of any three-dimensional (3D) rough surface | |
CN101413824B (en) | Moving body acoustic field measuring method based on random microphone array | |
CN104730518B (en) | A kind of method in the RADOP Power estimation sea flow field based on Gauss curve fitting | |
CN103954953B (en) | The blind source error compensation method of a kind of airborne laser radar based on data-driven | |
CN101604019B (en) | Method for quickly calculating characterization and transfer of uncertainties of marine environment and sound fields | |
CN101980044B (en) | Method for tracking multiple targets under unknown measurement noise distribution | |
CN104569911B (en) | OBU positioning method, RSU and ETC system | |
CN103759726B (en) | One class cyclo-stationary Poisson signal rapid simulation method and hardware system thereof | |
CN104035083A (en) | Radar target tracking method based on measurement conversion | |
CN103472450A (en) | Non-uniform space configuration distributed SAR moving target three-dimensional imaging method based on compressed sensing | |
CN104077440A (en) | Junction surface contact area and rigidity confirming method based on surface fitting | |
CN104052557A (en) | Method for modeling Nakagami repeated fading channel | |
CN105353408A (en) | Wigner higher-order spectrum seismic signal spectral decomposition method based on matching pursuit | |
CN104516771A (en) | Efficient non-stationary random process simulation method | |
CN102073754B (en) | Comprehensive electromechanical analysis method of reflector antenna based on error factor | |
CN102253385A (en) | Ocean internal wave forecast method based on synthetic aperture radar image and internal wave model | |
CN104318593A (en) | Simulation method and system of radar sea clusters | |
CN106583509A (en) | Bending process for irregular part | |
CN103699810A (en) | Modeling method of rough surface microwave band bidirectional reflectance distribution function | |
CN101854216A (en) | Channel parity researching method based on ionosphere correction layer channel model | |
CN103852648A (en) | Method for obtaining space electromagnetic intensity data | |
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 (en) | Earth-magnetism navigation position matching error is determined method | |
CN109001670B (en) | Distributed passive positioning method and device combining time difference and angle | |
CN101598788B (en) | Rapid simulation method of synthetic aperture sonar signal |
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 |