CN102779200A - Analytical method for contact performance of junction surface containing microcosmic surface shape - Google Patents

Abstract

The invention discloses an analytical method for a contact performance of a junction surface containing a microcosmic surface shape. According to the analytical method, the gridding dividing quality and the solving efficiency are increased and the forecasting accuracy and reliability of the contact performance of rough surfaces are improved. The analytical method comprises the following steps: measuring a practical surface through a laser confocal microscopy or obtaining the rough surfaces by utilizing a three-dimensional shape digital simulation method; extracting a height matrix of the shape of each 3D (Three-Dimensional) rough surface; utilizing a point cloud processing method to discretize each 3D rough surface into a height matrix file which is easily extracted by finite element software; importing the height matrix files into the finite element software, thereby generating a microcosmic shape of each connecting surface of the junction surface; establishing body models considering the microcosmic surface shapes by utilizing translation coordinate modifying functions of a key point of the finite element software; performing hexahedron gridding division on each body model according to a gridding control method; establishing finite element contact models in contact with the two 3D rough surfaces; gradually applying displacement and force loading boundary conditions; and analyzing the contact feature of the junction surface.

Description

A kind of faying face contact performance analytical approach that comprises the microcosmic surface pattern

Technical field

The present invention relates to a kind of faying face contact performance simulating analysis; The faying face performance simulation analytical approach of particularly a kind of consideration two 3D rough surfaces contact can be used for analyzing contact performancies such as the surperficial friction of micro-rough, wearing and tearing, contact stiffness, thermal contact resistance.

Background technology

There is a large amount of faying faces in the mechanical system; The existence of these faying faces has very big influence to the mechanical system overall permanence; And the contact of faying face actual be the contact of its rough surface, rough surface morphology all has significant effects to faying face friction, wearing and tearing and juxtaposition metamorphose thereof.

For a long time, Chinese scholars has been carried out a large amount of research to the contact of rough surface on theoretical aspect.Aspect theoretical analytic model; Made up multiple classical contact model: GW model, CEB model, MB model and KE model etc.; But all there is certain assumed condition in the foundation of these theoretical models; Its distortion is separate when all having the contact of identical summit radius-of-curvature, micro-bulge like all micro-bulges, and ignoring interaction between the micro-bulge, not having contact between large deformation, the micro-bulge all is that the peak is just contacting the peak etc.Just because of these hypothesis, make its actual micro-bulge contact condition of can't fitting between faying face, this has greatly influenced the precision of utilizing analytic model prediction faying face performance.

In order to address the above problem; In recent years; There is the scholar to propose to utilize the method for microcosmic contact between finite element numerical model simulation faying face; But it mainly concentrates on single to the contact performance analysis of micro-bulge and the contact performance analysis aspect of smooth surface-rough surface, the two 3D rough surface contact performance simulation analysis of actual contact state of more fitting of still being unrealized.The present invention starts with from a cloud treatment technology, phantom type generation technique and three aspects of grid dividing technology; Overcome the prior art bottleneck; Improve the grid dividing quality, made up two 3D rough surface contact models, thereby improved the accuracy of faying face performance prediction to a certain extent.

Summary of the invention

Fundamental purpose of the present invention provides a kind of faying face contact performance analytical approach that comprises the microcosmic surface pattern; Overcome two rough surface finite element contact models and made up the processing of mid point cloud, the generation of phantom type and grid dividing technical bottleneck; Improve the grid dividing quality and found the solution efficient, improved rough surface contact performance prediction accuracy and reliability.

Technical scheme of the present invention is:

A kind of faying face contact performance analytical approach that comprises the microcosmic surface pattern comprises following step:

1) measure real surface through laser confocal microscope or utilize the three-dimensional appearance analogy method to obtain Gauss or non-Gauss's rough surface, the height matrix that extracts the 3D rough surface morphology be Z (m, n), wherein m, n represent that respectively x, y direction get number a little;

2) be object with the microcosmic surface pattern that obtains in the step 1), choose suitable spacing, utilize some cloud disposal route that the 3D rough surface is separated into the txt file that finite element software extracts easily;

3) utilize step 2) txt file of the three-dimensional point that generates, import in the finite element software and generate point cloud model, utilize the dual loop command of finite element, connect adjacent 4 points of X/Y direction, generate microcosmic dough sheet crowd, each connects the microscopic appearance on surface to obtain faying face;

4) two groups of rough surfaces that generate with step 3) are basis, through planar stretch, and utilize translation and the coordinate modify feature of finite element software key point, make up the phantom type of one group of consideration connection surface microscopic surface topography;

5) the phantom type that generates with step 4) is the basis, through the method for gauze lattice control, the phantom type is carried out hexahedral mesh divide;

6) grid model that generates with step 5) is the basis; The contact of definition faying face is right; The fixing displacement of connector lower surface z direction down through last connector upper surface displacement constraint imposed load boundary condition, makes up the finite element contact model of two 3D rough surfaces contact; And carry out structural static and find the solution, observation connects the contact load and the contact pressure regularity of distribution of faying face.

The described three-dimensional appearance analogy method of utilizing obtains Gauss's rough surface, specifically comprises following steps:

1) utilize white noise generator randn (m, n) produce white noise sequence η (m, n), and calculate its Fourier transform A (m, n):

$A(k+1,l+1)=\underset{r-0}{\overset{m-1}{\mathrm{\Σ}}}\underset{s-0}{\overset{n-1}{\mathrm{\Σ}}}\mathrm{\η}(r+1,s+1)e^{-2\mathrm{\πi}(\mathrm{kr}/m+\mathrm{ls}/n)}\begin{array}{c}k=\mathrm{0,1,2}...m-1\\ l=\mathrm{0,1,2}...n-1\end{array}$

2) utilize given autocorrelation function generate correlation matrix R (m, n); Suppose autocorrelation function be f (x, y), β wherein _x, β _yBe respectively the auto-correlation length of X, Y direction;

$f(x,y)=\exp (-2.3\sqrt{{\left(\frac{x}{{\mathrm{\β}}_x}\right)}^2+{\left(\frac{y}{{\mathrm{\β}}_y}\right)}^2})$

Specify X, the discrete spacing of Y direction to be respectively Δ _x=1 μ m, Δ _y=1 μ m, at 0≤x≤m/2,0≤y≤n/2 disperses to autocorrelation function in the zone; (m/2+1 n/2+1), expands to-m/2≤x≤m/2 the discrete matrix R of acquisition autocorrelation function then; The autocorrelation matrix function discrete matrix R of correspondence in-n/2≤y≤n/2 zone (m, n):

R(k+1，l+1)＝f(k，l) k＝0，1，2...0.5m，l＝0，1，2...0.5n

R(0.5m+2+k，l+1)＝f(0.5m-k，l+1) k＝0，1，2...0.5m-2，l＝0，1，2...0.5n

R(k+1，0.5n+l+2)＝f(k+1，0.5n-l) k＝0，1，2...0.5m，l＝0，1，2...0.5n-2

R(0.5m+k+2，0.5n+l+2)＝f(0.5m-k，0.5n-l) k＝0，1，2...0.5m-2，l＝0，1，2...0.5n-2

3) (m n) carries out Fourier transform, obtains the power spectrum density P of simulating surface to autocorrelation function matrix R;

$P(1+k,1+l)=\left|\underset{r=0}{\overset{m-1}{\mathrm{\Σ}}}\underset{s=0}{\overset{n-1}{\mathrm{\Σ}}}R(r+1,s+1)e^{-2\mathrm{\πi}(\mathrm{kr}/m+\mathrm{ls}/n)}\right|\begin{array}{c}k=\mathrm{0,1,2}...m-1\\ l=\mathrm{0,1,2}...n-1\end{array}$

4) power spectrum density by white noise is a constant C, supposes C=1, calculate white noise to the transfer function H on simulation Gauss surface (m, n):

$H(m,n)=\sqrt{P}$

5) utilize frequency domain dot product (A.*H) negate Fourier transform method to obtain the elemental height matrix sequence z on Gauss surface:

$z(k+1,l+1)=\frac{1}{\mathrm{mn}}\underset{r=0}{\overset{m-1}{\mathrm{\&Sigma;}}}\underset{s=0}{\overset{n-1}{\mathrm{\&Sigma;}}}(A.*H){\mathrm{<figure-callout\; id="2"\; label="e"\; filenames="HDA0000120186190000013.png,HDA0000120186190000014.png"\; state="\{\{state\}\}">e</figure-callout>}}^{2\mathrm{\&pi;i}(\mathrm{kr}/m+\mathrm{ls}/n)}\begin{array}{c}k=\mathrm{0,1,2}...m-1\\ l=\mathrm{0,1,2}...n-1\end{array}$

6), ask white noise power spectrum density constant C=σ/std2 (z), thereby simulation generates the Gauss's surface elevation matrix z:z=z*C. that meets requirements such as given standard deviation according to the actual standard deviation std2 (z) of assigned altitute standard deviation and simulating surface z

The described three-dimensional appearance analogy method of utilizing obtains non-Gauss's rough surface, specifically comprises following steps:

1) directly utilizes the height sequence z that digitized simulation generated of Gauss's rough surface, adopt Pearson or the non-Gaussian transformation of Johnson system to carry out non-Gaussian transformation, generate non-gaussian sequence Z;

2) if the measure of skewness S of the non-gaussian sequence that is generated _kWith kurtosis K _uDo not satisfy accuracy requirement, then adopt new white noise sequence, repeat the simulation and the non-Gauss conversion thereof of Gauss's rough surface height sequence,, accomplish the digitized simulation of non-Gauss's rough surface up to satisfying accuracy requirement.

Described utilization point cloud disposal route is separated into the height matrix file that finite element software extracts easily with the 3D rough surface, specifically comprises following steps:

1) for actual measurement or simulating surface, can generate the dot file that comprises 3 row in the horizontal direction according to certain rules, wherein three data of each row are represented a point coordinate, and the 1st, 2,3 row are represented x, y, the z coordinate of corresponding point respectively;

2) be the some establishment order of finite element software identification through Ultraedit software with the form modifying of putting, thereby make things convenient for finite element software to read the txt dot file of this type, directly generate point cloud model.

Described translation and the coordinate modify feature that utilizes the finite element software key point makes up one group of phantom type of considering the microcosmic surface appearance, specifically comprises following steps:

1) respectively two groups of microcosmic surfaces is stretched along Z positive dirction and opposite direction, generate two groups of little hexahedron crowds with certain altitude;

2) revise the point coordinate that each organizes little hexahedron crowd, make connect the surface among each group hexahedron crowd vertex of surface is had identical Z coordinate figure, form smooth surface;

3) Z that regulates a certain connector guarantees that to the position two connection surfaces have partial points that initial contact has just taken place.

Described method through the control of gauze lattice is carried out the hexahedral mesh division to the phantom type, specifically comprises following steps:

1) control of grid cell quantity: in the phantom type of above-mentioned structure, select all lines of thickness direction, control the sizing grid or the grid number of plies of lines short transverse, pass through the size of controlling level direction grid simultaneously, accurately the quantity of control mesh;

2) in same little hexahedron; If the lines of four short transverses; Highly differ greatly; Sizing grid to its short transverse is suitably adjusted, and makes the little corresponding grid number of lines of big grid number of height value and height value differ even number, thereby guarantees to generate smoothly high-quality hexahedral mesh;

3) specify grid cell type and material properties, generate hexahedron scanning grid.

The present invention proposes the faying face performance simulation analytical approach of a kind of consideration two 3D rough surfaces contact; Start with from a cloud treatment technology, phantom type generation technique and three aspects of grid dividing technology; Overcome the prior art bottleneck; Improved the grid dividing quality and found the solution efficient, made up two 3D rough surface contact models, thereby solved the difficult convergence of finding the solution in the problem, accuracy and the reliability of having improved the faying face performance prediction to a certain extent.With this, can effectively instruct design, processing and the assembling of faying face, make it satisfy the faying face performance requirement, thereby ensure the mechanized equipment overall performance.

Description of drawings

Fig. 1 isotropy Gauss surface and height profile (β thereof _x=β _y=3, μ=0, σ=1.6 count 256 * 256)

Fig. 2 anisotropic gaussian surface and height profile (β thereof _x=2, β _y=100, μ=0, σ=1.6 count 256 * 256)

Fig. 3 isotropy non-Gauss surface and height profile (β thereof _x=β _y=5, μ=0, σ=1.6, S _k=-1, K _u=4.5, count 128 * 128)

Fig. 4 anisotropy non-Gauss surface and height profile (β thereof _x=2, β _y=100, μ=0, σ=1.6, S _Kz=0.5, K _Uz=3.5, count 256 * 256)

Fig. 5 constitutes the some cloud (128 * 128 * 2) of faying face

Fig. 6 connects the dough sheet crowd on surface

The body that Fig. 7 connects the surface up and down stretches, connector crowd's apex coordinate is revised and the Z axle translation of connector up and down

Fig. 8 considers the hexahedral mesh illustraton of model of two rough surfaces contact

Fig. 9 load constantly increases corresponding last connector Z to variation of stress

Figure 10 load constantly increases corresponding last connector equivalent stress Changing Pattern

Figure 11 load constantly increases corresponding faying face contact stress and contact area Changing Pattern

Embodiment

Digitized simulation with the non-Gauss's rough surface of 3D gaussian sum that satisfies statistical natures such as certain autocorrelation function, standard deviation, measure of skewness and kurtosis is an example, and microcosmic surface pattern acquisition methods of the present invention is described; With one group of Gauss's surface topography contact is example, adopts the ANSYS finite element software, and faying face contact performance model construction of the present invention and analytical approach are described.

1. the microcosmic surface pattern obtains:

Utilize the three-dimensional appearance analogy method to obtain Gauss or non-Gauss's rough surface, the height matrix that extracts the 3D rough surface morphology be Z (m, n), wherein m, n represent that respectively x, y direction get number a little, concrete steps are following:

(1) utilize the three-dimensional appearance analogy method to obtain Gauss's rough surface, step is following:

1) utilizes white noise generator randn (m; N) produce white noise sequence η (m; And obtain its Fourier transform n),

$n):A(k+1,l+1)=\underset{r-0}{\overset{m-1}{\mathrm{\&Sigma;}}}\underset{s-0}{\overset{n-1}{\mathrm{\&Sigma;}}}\mathrm{\&eta;}(r+1,s+1)e^{-2\mathrm{\&pi;i}(\mathrm{kr}/m+\mathrm{ls}/n)}\begin{array}{c}k=\mathrm{0,1,2}...m-1\\ l=\mathrm{0,1,2}...n-1\end{array}$

2) utilize given autocorrelation function generate correlation matrix R (m, n); Suppose that autocorrelation function is f, wherein β _x, β _yBe respectively the auto-correlation length of X, Y direction;

$f(x,y)=\exp (-2.3\sqrt{{\left(\frac{x}{{\mathrm{\&beta;}}_x}\right)}^2+{\left(\frac{y}{{\mathrm{\&beta;}}_y}\right)}^2})$

Specify X, the discrete spacing of Y direction to be respectively Δ _x=1 μ m, Δ _y=1 μ m, at 0≤x≤m/2,0≤y≤n/2 disperses to autocorrelation function in the zone; (m/2+1 n/2+1), expands to-m/2≤x≤m/2 the discrete matrix R of acquisition autocorrelation function then; The autocorrelation matrix function discrete matrix R of correspondence in-n/2≤y≤n/2 zone (m, n):

R(k+1，l+1)＝f(k，l) k＝0，1，2...0.5m，l＝0，1，2...0.5n

R(0.5m+2+k，l+1)＝f(0.5m-k，l+1) k＝0，1，2...0.5m-2，l＝0，1，2...0.5n

R(k+1，0.5n+l+2)＝f(k+1，0.5n-l) k＝0，1，2...0.5m，l＝0，1，2...0.5n-2

R(0.5m+k+2，0.5n+l+2)＝f(0.5m-k，0.5n-l) k＝0，1，2...0.5m-2，l＝0，1，2...0.5n-2

3) (m n) carries out Fourier transform, obtains the power spectrum density P of simulating surface to autocorrelation function matrix R;

$P(1+k,1+l)=\left|\underset{r=0}{\overset{m-1}{\mathrm{\&Sigma;}}}\underset{s=0}{\overset{n-1}{\mathrm{\&Sigma;}}}R(r+1,s+1)e^{-2\mathrm{\&pi;i}(\mathrm{kr}/m+\mathrm{ls}/n)}\right|\begin{array}{c}k=\mathrm{0,1,2}...m-1\\ l=\mathrm{0,1,2}...n-1\end{array}$

4) power spectrum density by white noise is a constant C, supposes C=1, calculate white noise to the transfer function H on simulation Gauss surface (m, n):

$H(m,n)=\sqrt{P}$

5) utilize frequency domain dot product (A.*H) negate Fourier transform method to obtain the elemental height matrix sequence z on Gauss surface:

$z(k+1,l+1)=\frac{1}{\mathrm{mn}}\underset{r=0}{\overset{m-1}{\mathrm{\&Sigma;}}}\underset{s=0}{\overset{n-1}{\mathrm{\&Sigma;}}}(A.*H){\mathrm{<figure-callout\; id="2"\; label="e"\; filenames="HDA0000120186190000013.png,HDA0000120186190000014.png"\; state="\{\{state\}\}">e</figure-callout>}}^{2\mathrm{\&pi;i}(\mathrm{kr}/m+\mathrm{ls}/n)}\begin{array}{c}k=\mathrm{0,1,2}...m-1\\ l=\mathrm{0,1,2}...n-1\end{array}$

6), ask white noise power spectrum density constant C=σ/std2 (z), thereby simulation generates the Gauss's surface elevation matrix z:z=z*C that meets requirements such as given standard deviation according to the actual standard deviation std2 (z) of assigned altitute standard deviation and simulating surface z.

Utilize above-mentioned given autocorrelation function; The autocorrelation function of x, y direction is all got 3 μ m; The x of simulating surface, y direction get a little be 256 * 256, the standard deviation on surface is 1.6 μ m, isotropy Gauss's surface topography and height profile thereof that simulation obtains are as shown in Figure 1.The autocorrelation function of x, y direction is got 2 μ m, 100 μ m respectively, and Gauss's surface topography that simulation obtains shows anisotropy, has certain texture, and is as shown in Figure 2.

(2) utilize the three-dimensional appearance analogy method to obtain non-Gauss's rough surface, step is following:

1) directly utilizes the height sequence z that digitized simulation generated of Gauss's rough surface, adopt Pearson or the non-Gaussian transformation of Johnson system to carry out non-Gaussian transformation, generate non-gaussian sequence Z;

2) if the measure of skewness S of the non-gaussian sequence that is generated _kWith kurtosis K _uDo not satisfy accuracy requirement, then adopt new white noise sequence, repeat the simulation and the non-Gauss conversion thereof of Gauss's rough surface height sequence,, accomplish the digitized simulation of non-Gauss's rough surface up to satisfying accuracy requirement.

Utilize above-mentioned given autocorrelation function; The autocorrelation function of x, y direction is all got 5 μ m; The x of simulating surface, y direction get a little be 128 * 128, average, standard deviation, measure of skewness and the kurtosis on surface be respectively 0,1.6 μ m ,-1 and 4.5, non-Gauss's surface topography of isotropy and the height profile thereof of simulating acquisition are as shown in Figure 3.If the autocorrelation function of x, y direction is got 2 μ m, 100 μ m respectively; The x of simulating surface, y direction get a little be 256 * 256, standard deviation, measure of skewness and the kurtosis on surface be respectively 1.6 μ m, 0.5 and 3.5; The non-Gauss's surface topography of anisotropy that simulation obtains demonstrates certain texture, and is as shown in Figure 2.

2. microscopic appearance disperses and information stores

Utilize the acquisition methods of step 1 microscopic appearance; Obtain two 128 * 128 Gauss's rough surfaces (the upper and lower surfaces standard deviation is respectively 0.8 μ m and 1.6 μ m) of counting; X, y direction spacing are chosen 1 μ m; Utilize some cloud disposal route that the 3D rough surface is separated into the txt file that finite element software extracts easily, concrete steps are following:

1) for actual measurement or simulating surface, can generate the dot file that comprises 3 row in the horizontal direction according to certain rules, wherein three data of each row are represented a point coordinate, and the 1st, 2,3 row are represented x, y, the z coordinate of corresponding point respectively;

2) be that the point of finite element software identification is created order " k,, X, Y, Z " (point that k is equivalent among the ANSYS is created order through Ultraedit software with the form modifying of point; Omitted the code name of point between ",, " two commas; X, Y, Z represent coordinate figure a little respectively, separate through comma), directly generate point cloud model thereby make things convenient for ANSYS software to read such txt dot file.

3. the surface model of microcosmic surface pattern makes up

The txt file of the three-dimensional point that step 2 is generated; Import in the ANSYS software and generate point cloud model (as shown in Figure 5), utilize ANSYS dual loop command (* do ...; * end); Connect adjacent 4 points of X/Y direction, generate microcosmic dough sheet crowd, each connects the microscopic appearance (as shown in Figure 6) on surface to obtain faying face.

4. the phantom type of microcosmic surface pattern makes up

Two groups of rough surfaces that generate with step 3 are the basis; Pass through planar stretch; Translation and coordinate modify feature in conjunction with the ANSYS key point; Make up one group of phantom type that consider to connect the surface microscopic surface topography, connect the stretching on surface up and down, modification and the connector crowd's of connector crowd summit Z coordinate Z axle translation is as shown in Figure 7 up and down, concrete steps are following:

1) respectively the positive and negative direction of upper and lower two groups of microcosmic surfaces along the Z axle stretched, generation has certain thickness two groups of little hexahedron crowds of (thickness is about 12 μ m here greater than the largest contours height of two rough surfaces);

2) revise the apex coordinate on upper and lower two groups of little hexahedron crowds' upper and lower surface respectively, make the summit on the upper and lower surface among upper and lower two groups of hexahedron crowds have identical Z coordinate figure respectively, form smooth surface;

3) Z that regulates a last connector guarantees that to the position two connection surfaces have partial points that initial contact just takes place.

5. phantom type grid dividing

The phantom type that generates with step 4 is the basis, through the method for gauze lattice control, the phantom type is carried out hexahedral mesh divide;

1) control of grid cell quantity: all lines of in the phantom type of above-mentioned structure, selecting thickness direction; The sizing grid or the grid number of plies of control lines short transverse; Pass through the size of controlling level direction grid simultaneously; The accurate quantity of control mesh is 127 * 127 * 4 * 2 like the number of grid of Fig. 8;

2) in same little hexahedron, because the lines of four short transverses highly are more or less the same, at 4 layers of grid cell of short transverse unified Definition;

3) specifying the grid cell type is Solid45, and material is 45 steel up and down, generates hexahedron scanning grid, as shown in Figure 8.

6. the faying face contact performance is analyzed

The grid model that generates with step 5 is the basis, and the contact of definition faying face is right, the fixedly displacement of lower end phantom type lower surface z direction; Through upper end phantom type upper surface displacement constraint imposed load boundary condition (0.5,0.8,11; 1.4 ..., 2.3; 2.6 μ m), make up the finite element contact model of two 3D rough surfaces contact, and distortion, the Z of micro-bulge in the contact of two rough surfaces analyzed to the contact pressure and the contact area of stress and equivalent stress and faying face.Last connector Z corresponding when increasing for load wherein shown in Figure 9 is to variation of stress; Last connector equivalent stress Changing Pattern corresponding when increasing shown in Figure 10 for load, the faying face contact pressure of correspondence and the Changing Pattern of contact area when increasing shown in Figure 11 for load.

Claims (6)

1. a faying face contact performance analytical approach that comprises the microcosmic surface pattern is characterized in that, comprises following step:

1) measure real surface through laser confocal microscope or utilize the three-dimensional appearance analogy method to obtain Gauss or non-Gauss's rough surface, the height matrix that extracts the 3D rough surface morphology be Z (m, n), wherein m, n represent that respectively x, y direction get number a little;

2) be object with the microcosmic surface pattern that obtains in the step 1), choose suitable spacing, utilize some cloud disposal route that the 3D rough surface is separated into the txt file that finite element software extracts easily;

3) utilize step 2) txt file of the three-dimensional point that generates, import in the finite element software and generate point cloud model, utilize the dual loop command of finite element, connect adjacent 4 points of X/Y direction, generate microcosmic dough sheet crowd, each connects the microscopic appearance on surface to obtain faying face;

4) two groups of rough surfaces that generate with step 3) are basis, through planar stretch, and utilize translation and the coordinate modify feature of finite element software key point, make up the phantom type of one group of consideration connection surface microscopic surface topography;

5) the phantom type that generates with step 4) is the basis, through the method for gauze lattice control, the phantom type is carried out hexahedral mesh divide;

6) grid model that generates with step 5) is the basis; The contact of definition faying face is right; The fixing displacement of connector lower surface z direction down through to last connector upper surface displacement constraint imposed load boundary condition, makes up the finite element contact model of two 3D rough surfaces contact; And carry out structural static and find the solution, observation connects the contact load and the contact pressure regularity of distribution of faying face.

2. a kind of faying face contact performance analytical approach that comprises the microcosmic surface pattern according to claim 1, the described three-dimensional appearance analogy method of utilizing obtains Gauss's rough surface, specifically comprises following steps:

1) utilize white noise generator randn (m, n) produce white noise sequence η (m, n), and obtain its Fourier transform A (m, n);

2) utilize given autocorrelation function generate correlation matrix R (m, n); Suppose autocorrelation function be f (x, y), β wherein _x, β _yBe respectively the auto-correlation length of X, Y direction;

Specify X, the discrete spacing of Y direction to be respectively Δ _x=1 μ m, Δ _y=1 μ m, at 0≤x≤m/2,0≤y≤n/2 disperses to autocorrelation function in the zone; (m/2+1 n/2+1), expands to-m/2≤x≤m/2 the discrete matrix R of acquisition autocorrelation function then; The autocorrelation matrix function discrete matrix R of correspondence in-n/2≤y≤n/2 zone (m, n):

R(k+1,l+1)=f(k,l) k=0,1,2...0.5m，l=0,1,2...0.5n

R(0.5m+2+k,l+1)=f(0.5m-k,l+1) k=0,1,2...0.5m-2，l=0,1,2...0.5n

R(k+1,0.5n+l+2)=f(k+1,0.5n-l) k=0,1,2...0.5m，l=0,1,2...0.5n-2

R(0.5m+k+2,0.5n+l+2)=f(0.5m-k,0.5n-l) k=0,1,2...0.5m-2，l=0,1,2...0.5n-2

3) (m n) carries out Fourier transform, obtains the power spectrum density P of simulating surface to autocorrelation function matrix R;

4) power spectrum density by white noise is a constant C, supposes C=1, calculate white noise to the transfer function H on simulation Gauss surface (m, n):

5) utilize frequency domain dot product (A.*H) negate Fourier transform method to obtain the elemental height matrix sequence z on Gauss surface:

6), ask white noise power spectrum density constant C=σ/std2 (z), thereby simulation generates the Gauss's surface elevation matrix z:z=z*C. that meets requirements such as given standard deviation according to the actual standard deviation std2 (z) of assigned altitute standard deviation and simulating surface z

3. a kind of faying face contact performance analytical approach that comprises the microcosmic surface pattern according to claim 1, the described three-dimensional appearance analogy method of utilizing obtains non-Gauss's rough surface, specifically comprises following steps:

1) directly utilizes the height sequence z that digitized simulation generated of Gauss's rough surface, adopt Pearson or the non-Gaussian transformation of Johnson system to carry out non-Gaussian transformation, generate non-gaussian sequence Z;

2) if the measure of skewness S of the non-gaussian sequence that is generated _kWith kurtosis K _uDo not satisfy accuracy requirement, then adopt new white noise sequence, repeat the simulation and the non-Gauss conversion thereof of Gauss's rough surface height sequence,, accomplish the digitized simulation of non-Gauss's rough surface up to satisfying accuracy requirement.

4. a kind of faying face contact performance analytical approach that comprises the microcosmic surface pattern according to claim 1, described utilization point cloud disposal route is separated into the height matrix file that finite element software extracts easily with the 3D rough surface, specifically comprises following steps:

1) for actual measurement or simulating surface, can generate the dot file that comprises 3 row in the horizontal direction according to certain rules, wherein three data of each row are represented a point coordinate, and the 1st, 2,3 row are represented x, y, the z coordinate of corresponding point respectively;

2) be the some establishment order of finite element software identification through Ultraedit software with the form modifying of putting, thereby make things convenient for finite element software to read the txt dot file of this type, directly generate point cloud model.

5. a kind of faying face contact performance analytical approach that comprises the microcosmic surface pattern according to claim 1; Described translation and the coordinate modify feature that utilizes the finite element software key point; Make up one group of phantom type of considering the microcosmic surface appearance, specifically comprise following steps:

1) respectively two groups of microcosmic surfaces is stretched along Z positive dirction and opposite direction, generate two groups of little hexahedron crowds with certain altitude;

2) revise the point coordinate that each organizes little hexahedron crowd, make connect the surface among each group hexahedron crowd vertex of surface is had identical Z coordinate figure, form smooth surface;

3) Z that regulates a certain connector guarantees that to the position two connection surfaces have partial points that initial contact has just taken place.

6. a kind of faying face contact performance analytical approach that comprises the microcosmic surface pattern according to claim 1, described method through the control of gauze lattice is carried out the hexahedral mesh division to the phantom type, specifically comprises following steps:

1) control of grid cell quantity: in the phantom type of above-mentioned structure, select all lines of thickness direction, control the sizing grid or the grid number of plies of lines short transverse, pass through the size of controlling level direction grid simultaneously, accurately the quantity of control mesh;

2) in same little hexahedron; If the lines of four short transverses; Highly differ greatly; Sizing grid to its short transverse is suitably adjusted, and makes the little corresponding grid number of lines of big grid number of height value and height value differ even number, thereby guarantees to generate smoothly high-quality hexahedral mesh;

3) specify grid cell type and material properties, generate hexahedron scanning grid.

Images