CN109063328B - Method for acquiring rough surface bearing area rate curve - Google Patents
Method for acquiring rough surface bearing area rate curve Download PDFInfo
- Publication number
- CN109063328B CN109063328B CN201810860990.4A CN201810860990A CN109063328B CN 109063328 B CN109063328 B CN 109063328B CN 201810860990 A CN201810860990 A CN 201810860990A CN 109063328 B CN109063328 B CN 109063328B
- Authority
- CN
- China
- Prior art keywords
- height
- area ratio
- smc
- bearing
- bearing area
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
Abstract
The invention relates to a method for acquiring a rough surface bearing area rate curve. The method comprises the following steps: measuring the three-dimensional geometrical morphology of the rough surface and carrying out data pretreatment; determining the bearing surface area ratio of the height of the surface to be calculatedp i ，i=1,2,3, \ 8230; for each onep i Calculating the corresponding surface height by using a search algorithmc i ，i=1,2,3, \ 8230; from the calculated series (p i ,c i ) And obtaining a bearing area rate curve. The invention is based on the bearing area ratiop i Calculating surface heightc i In the core area of the obtained bearing area rate curve, under the condition of the same calculated amount, the resolution of the corresponding relation between the bearing area rate and the surface height is far higher than that of the traditional method; the invention adopts a search algorithm to calculate the surface height at the position of the given bearing area rate, and the calculation precision is higher.
Description
Technical Field
The invention belongs to the technical field of analysis and characterization of the geometric morphology of a machined surface, and particularly relates to a method for acquiring a rough surface bearing area rate curve.
Background
The surface of a mechanical part has various surface geometries that have very important effects on various physical and chemical properties of the surface, such as friction, wear, sealing, lubrication, thermal conduction, electrical conduction, light reflection and refraction, and the like. Researches show that only the geometric characteristics of the rough surface can be completely reflected by the three-dimensional geometric shape, and the prediction and control of the surface performance can be realized, so that the three-dimensional geometric shape measurement and characterization of the rough surface are more and more widely and deeply applied in scientific researches and industries along with the progress and development of scientific technologies.
In surface three-dimensional characterization parameter systems such as Birminghan14, STANDSURF and ISO 25178, a surface bearing area rate curve is defined, and the curve and characterization parameters defined on the basis of the curve can well reflect the bearing performance, the running-in characteristic, the wear resistance, the lubricating oil film forming capacity and the lubricating oil storage capacity of a rough surface.
Currently, the method for obtaining the surface bearing area ratio curve is as follows: the method comprises the steps of dividing a surface at equal intervals along the height direction, calculating corresponding surface supporting area rates at all heights, obtaining a series of (supporting area rates and surface heights), and finally drawing a surface supporting area rate curve by taking the supporting area rates as a horizontal axis and the surface heights as a vertical axis. However, after machining (especially finishing) of a rough surface of a part, the heights of most surface sampling points tend to be distributed in a narrow interval (core region). However, there may be places where individual, particularly prominent high points, or particularly deep valleys, often do not have much material effect on the surface properties. The main drawbacks with the existing methods are therefore:
disadvantage 1: the prior art method determines the bearing area ratio by the surface height. However, if a reverse solution is required, that is, the surface height is obtained by specifying the supporting area ratio, the functional relationship between the supporting area ratio and the surface height needs to be established by using an interpolation method on the basis of obtaining all data points (supporting area ratio, surface height), and finally the surface height can be obtained. Therefore, the method is not flexible enough, and the calculation accuracy of the surface height is related to the data point quantity, so that the relationship between the calculation accuracy and the calculation speed is difficult to adjust and balance according to the requirement.
And (2) disadvantage: the surface core area is an important area for representing the surface performance, is divided at equal intervals along the height direction, and the method for calculating the supporting area rate usually ensures that the number of the representing points of a supporting area rate curve in the core area is small, and the resolution of the variation relation between the representing supporting area rate and the surface height is insufficient; while the resolution of the characterization of the apparent incidental salient locations is too great, wasting computational resources.
Disclosure of Invention
The invention aims to provide a method for acquiring a rough surface bearing area rate curve, which can acquire the bearing area rate curve by dividing the bearing area rate at equal intervals or at any intervals by taking the bearing area rate as an independent variable and the surface height as a dependent variable.
In order to realize the purpose, the technical scheme of the invention is as follows: a method for acquiring a rough surface bearing area ratio curve comprises the following steps:
s1, calculating corresponding surface height when the bearing area ratio of the surface is a preset value by adopting a search algorithm;
and S2, taking the supporting area rate as an independent variable and the corresponding surface height as a dependent variable, and acquiring a supporting area rate curve.
In an embodiment of the present invention, the step S1 specifically includes:
step S11, inputting the parameter as the bearing area ratio of the surfacep i The result being the corresponding surface heightc i ；
Step S12, setting the surface heightc i Variable for left end point of search space of (2)clRepresenting that the initial value is the maximum value of the height of the surface sampling point;
step S13, setting the surface heightc i Variable for right end point of search space (2)crRepresenting that the initial value is the minimum value of the height of the surface sampling point;
step S14, setting the height of the surface to be obtainedc i By means of variablessmcIndicates and sets its initial value to (cl + cr) / 2；
Step S15, entering a circulation search, if the surface height issmcCalculated value of area ratio of bearing andp i are equal to each other, thensmcWhen the bearing area ratio isp i Height of surface of timec i Ending the calculation; otherwise, executing step S16;
step S16, if in step S15, the surface height issmcThe calculated value of the area ratio of the bearing is larger thanp i Then set the right end point of the search intervalcrIs equal tosmc(ii) a If in step S15, the surface height issmcThe calculated value of the area ratio of the bearing is less thanp i Then set the left end point of the search intervalclIs equal tosmc(ii) a Order tosmcIs updated with a value of (cl + cr) 2; step S15 is repeatedly performed.
In an embodiment of the present invention, the step S2 is specifically implemented as follows:
obtaining a curve of rough surface bearing area ratio as bearing area ratiopAs independent variable, surface heightcIs a dependent variable; area ratio of supportpEqually spaced or divided at arbitrary intervals, and the division points are marked asp i ，iIs a natural number.
Compared with the prior art, the invention has the following beneficial effects:
(1) When a surface bearing area rate curve is obtained, the invention can carry out equal interval or arbitrary interval division on the bearing area rate;
(2) When the surface height is calculated according to the given surface bearing area ratio, the calculation accuracy can be adjusted according to the requirement, and the balance between the calculation accuracy and the calculation speed is achieved.
Drawings
Fig. 1 is a general flow chart of the present invention.
FIG. 2 is a flow chart of the present invention for calculating surface height based on bearing area ratio.
FIG. 3 is a graph of sample point height data for a polished surface after data preprocessing.
FIG. 4 is a graph of surface bearing area ratios obtained using the method of the present invention.
Fig. 5 is a graph of surface bearing area ratio obtained by a conventional method.
Detailed Description
The technical scheme of the invention is specifically explained by combining the attached drawings 1-5.
The invention provides a method for acquiring a rough surface bearing area rate curve, which comprises the following steps:
s1, calculating corresponding surface height when the bearing area ratio of the surface is a preset value by adopting a search algorithm;
and S2, taking the supporting area rate as an independent variable and the corresponding surface height as a dependent variable, and acquiring a supporting area rate curve.
The step S1 specifically includes:
step S11, inputting the parameter as the bearing area ratio of the surfacep i To be solved for the resultThe fruit has corresponding surface heightc i ；
Step S12, setting the surface heightc i Variable for left end point of search space of (2)clRepresenting that the initial value is the maximum value of the height of the surface sampling point;
step S13, setting the surface heightc i Variable for right end point of search space (2)crExpressing that the initial value is the minimum value of the height of the surface sampling point;
step S14, setting the height of the surface to be obtainedc i By means of variablessmcShow, and set its initial value to (cl + cr) / 2；
Step S15, entering a circulation search, if the surface height issmcCalculated value of area ratio of bearing andp i are equal to each other, thensmcWhen the bearing area ratio isp i Height of surface of timec i Ending the calculation; otherwise, executing step S16;
step S16, if in step S15, the surface height issmcThe calculated value of the area ratio of the bearing is larger thanp i Then set the right end point of the search intervalcrIs equal tosmc(ii) a If in step S15, the surface height issmcThe calculated value of the area ratio of the bearing is less thanp i Then set the left end point of the search intervalclIs equal tosmc(ii) a Order tosmcIs updated tocl + cr) 2; step S15 is repeatedly performed.
The step S2 is specifically realized as follows: obtaining a curve of rough surface bearing area ratio as bearing area ratiopAs independent variable, surface heightcIs a dependent variable; area ratio of supportpEqually spaced or divided at arbitrary intervals, and the division points are marked asp i ，iIs a natural number.
The following is a specific implementation of the present invention.
The method for acquiring the rough surface bearing area rate curve comprises the following specific implementation steps:
(1) Acquiring measurement data of the surface geometric morphology of the mechanical part, wherein the data is the height of the surface sampling by the measurement equipment;
(2) The surface data obtained by measurement are preprocessed, and the purpose of preprocessing is to: filtering out shape errors caused by a measuring platform, a sample and other factors during measurement; filtering high-frequency noise in the measured data; adjusting the average value of the heights of all sampling points on the surface to be 0;
(3) Writing a program for calculating the surface bearing area ratio at the corresponding height by taking the surface height as an input parameter according to the publication [1 ];
(4) Surface area ratio of bearingpIs 0 to 100 percent, is divided at equal intervals or any intervals according to requirements, and the bearing area ratio of each division point is recorded asp i , iIs a natural number;
(5) For eachp i And (5) repeatedly executing the steps (6) to (10) to calculate the corresponding surface heightc i ；
(6) Set surface heightc i Variable for left end point of search space of (2)clRepresenting that the initial value is the maximum value of the height of the surface sampling point;
(7) Set surface heightc i Variable for right end point of search space of (2)crRepresenting that the initial value is the minimum value of the height of the surface sampling point;
(8) Setting the height of the surface to be determinedc i By means of variablessmcIndicates and sets its initial value to (cl + cr) / 2；
(9) Enter a circular search if the surface height issmcCalculated value of area ratio of bearing andp i are equal to each other, thensmcWhen the bearing area ratio isp i The surface height of the time, finish the cyclic search; otherwise, executing the step (10);
(10) If in step (9) the surface height issmcThe calculated value of the area ratio of the bearing is larger thanp i Then set the right end point of the search intervalcrIs equal tosmc(ii) a If in step (9) the surface height issmcArea of bearingCalculated value of rate is less thanp i Then set the left end point of the search intervalclIs equal tosmc(ii) a Order tosmcIs updated with a value of (cl + cr) 2; repeatedly executing the step (9);
(11) From a series of calculated (bearing area ratio)p i Height of surfacec i ) And performing piecewise linear interpolation between the data points to obtain a bearing area rate curve of the rough surface.
The method for obtaining the rough surface bearing area ratio curve according to the present invention is specifically described by taking the polished surface sample shown in fig. 3 as an example.
According to the algorithm flows given by fig. 1 and fig. 2, the specific steps are as follows:
(1) Obtaining original measurement data of the geometric morphology of the rough surface by using a rough surface three-dimensional morphology instrument;
(2) According to ISO 25178 or other relevant documents, filtering out unnecessary low-frequency and high-frequency components in surface measurement data, and adjusting the average value of the heights of the surface sampling points to 0;
(3) According to the publication [1], writing a program for calculating the surface bearing area ratio at the corresponding height by taking the surface height as an input parameter;
(4) The value range of the bearing area ratio is 0-100%, and the bearing area ratio is divided into 20 equal intervals and recorded asp i ，iIs a natural number;
(5) For each obtained by dividing in the above-mentioned stepp i Executing steps (6) - (10);
(6) Set the left end point of the surface height search interval toclThe initial value is the maximum height value of the surface sampling point;
(7) Set the right end point of the surface height search interval tocrThe initial value is the minimum height value of the surface sampling point;
(8) Defining variablessmcThe surface height to be determined is represented and the initial value is (cl+cr)/2；
(9) Judging when the surface height issmcRate of area of support ofWhether the calculated value is ANDp i Equal, if equal, at this timesmcIs the height of the surface to be determinedc i The cycle is exited; otherwise, executing step (10);
(10) If the surface height issmcThe ratio of the bearing area is larger thanp i Let us ordercr=smcOtherwise, tocl=smcUpdatesmc=(cl+cr) Step (9) is repeatedly executed;
(11) Obtaining a series of data points representing (supporting area rate, surface height) according to the calculation result of the steps, and drawing a supporting area rate curve of the polishing surface by adopting a piecewise linear interpolation method (figure 4);
comparing fig. 5 (the curve of the surface supporting area ratio obtained by the prior art method) and fig. 4, the advantage of the present invention is clearly shown, that is, the data points calculated by the present invention are mostly in the core region of the curve of the surface supporting area ratio, and the change of the supporting area ratio and the surface height can be reflected more finely.
Reference:
[1] zhou super, libang, liweifan, etc. the rough surface modeling method [ J ] with specified bearing rate curve, proceedings of Chinese engineering machinery 2016, 14 (2): 114-118.
The above description is only an embodiment and an implementation example of the present invention, but the implementation example is only exemplary and does not limit the scope of the present invention in any way. It will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the spirit and scope of the invention, and that such changes and substitutions are intended to be within the scope of the invention.
Claims (1)
1. A method for acquiring a rough surface bearing area ratio curve is characterized by comprising the following steps:
s1, calculating corresponding surface height when the bearing area ratio of the surface is a preset value by adopting a search algorithm;
s2, taking the supporting area rate as an independent variable and the corresponding surface height as a dependent variable, and obtaining a supporting area rate curve;
the step S1 specifically includes:
step S11, inputting the parameter as the bearing area ratio of the surfacep i The result being the corresponding surface heightc i ；
Step S12, setting the surface heightc i Variable for left end point of search space of (2)clRepresenting that the initial value is the maximum value of the height of the surface sampling point;
step S13, setting the surface heightc i Variable for right end point of search space (2)crRepresenting that the initial value is the minimum value of the height of the surface sampling point;
step S14, setting the height of the surface to be obtainedc i By means of variablessmcShow, and set its initial value to (cl + cr) / 2；
Step S15, entering a circulation search, if the surface height issmcCalculated value of area ratio of bearing andp i are equal to each other, thensmcWhen the bearing area ratio isp i Height of surface of timec i Ending the calculation; otherwise, executing step S16;
step S16, if in step S15, the surface height issmcThe calculated value of the area ratio of the bearing is larger thanp i Then set the right end point of the search intervalcrIs equal tosmc(ii) a If in step S15, the surface height issmcThe calculated value of the area ratio of the bearing is less thanp i Then set the left end point of the search intervalclIs equal tosmc(ii) a Order tosmcIs updated tocl + cr) 2; repeatedly executing the step S15;
the step S2 is specifically realized as follows:
the curve of the rough surface bearing area ratio obtained is used as the bearing area ratiopAs independent variable, surface heightcIs a dependent variable; area ratio of supportpAre equally spaced or divided at arbitrary intervals, and the division points are marked asp i ，iIs a natural number.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810860990.4A CN109063328B (en) | 2018-08-01 | 2018-08-01 | Method for acquiring rough surface bearing area rate curve |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810860990.4A CN109063328B (en) | 2018-08-01 | 2018-08-01 | Method for acquiring rough surface bearing area rate curve |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109063328A CN109063328A (en) | 2018-12-21 |
| CN109063328B true CN109063328B (en) | 2023-02-28 |
Family
ID=64832375
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810860990.4A Active CN109063328B (en) | 2018-08-01 | 2018-08-01 | Method for acquiring rough surface bearing area rate curve |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109063328B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1988128A (en) * | 2005-12-22 | 2007-06-27 | 日本碍子株式会社 | Electrostatic chuck |
| CN101046225A (en) * | 2002-07-15 | 2007-10-03 | 日本精工株式会社 | Wheel-support rolling bearing unit |
| CN103093029A (en) * | 2012-12-19 | 2013-05-08 | 江苏申模数字化制造技术有限公司 | Disperse cramp bar pre-set height calculation method used for dot-matrix flexible tool |
| JP2017090453A (en) * | 2015-11-04 | 2017-05-25 | 学校法人東京理科大学 | Surface roughness evaluation device and surface roughness evaluation method |
| CN108062071A (en) * | 2017-12-20 | 2018-05-22 | 天津大学 | The real-time assay of parameter curve trace track profile errors |
-
2018
- 2018-08-01 CN CN201810860990.4A patent/CN109063328B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101046225A (en) * | 2002-07-15 | 2007-10-03 | 日本精工株式会社 | Wheel-support rolling bearing unit |
| CN1988128A (en) * | 2005-12-22 | 2007-06-27 | 日本碍子株式会社 | Electrostatic chuck |
| CN103093029A (en) * | 2012-12-19 | 2013-05-08 | 江苏申模数字化制造技术有限公司 | Disperse cramp bar pre-set height calculation method used for dot-matrix flexible tool |
| JP2017090453A (en) * | 2015-11-04 | 2017-05-25 | 学校法人東京理科大学 | Surface roughness evaluation device and surface roughness evaluation method |
| CN108062071A (en) * | 2017-12-20 | 2018-05-22 | 天津大学 | The real-time assay of parameter curve trace track profile errors |
Non-Patent Citations (1)
| Title |
|---|
| 机床支撑地脚结合部法向粗糙接触建模;田红亮等;《万方数据期刊库》;20151230;第1-8页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109063328A (en) | 2018-12-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109447380B (en) | Method and device for determining oil well yield | |
| Blateyron | The areal field parameters | |
| CN105205276B (en) | Oscillating bearing wear out failure physical modeling and analysis method | |
| Minet et al. | Analysis and modeling of the topography of mechanical seal faces | |
| US10969321B2 (en) | Method and system for determining contact angle of porous media | |
| CN109063328B (en) | Method for acquiring rough surface bearing area rate curve | |
| US20170175513A1 (en) | Evaluation of Production Performance from a Hydraulically Fractured Well | |
| Pawlus et al. | Two-process random textures: measurement, characterization, modeling and tribological impact: a review | |
| Venner | EHL film thickness computations at low speeds: risk of artificial trends as a result of poor accuracy and implications for mixed lubrication modelling | |
| CN108266184A (en) | A kind of Single Gas Well control reserve evaluation method based on creation data | |
| CN109190081B (en) | Method for acquiring rough surface contour bearing length rate curve | |
| CN110410070B (en) | Method for determining water flooding reservoir flooding rule | |
| CN113552036A (en) | Method and device for determining oil-water relative permeability of medium-low pore permeability reservoir | |
| CN109345041B (en) | Equipment failure rate prediction method combining Weibull distribution and ARMA (autoregressive moving average) | |
| CN108717498B (en) | Oil reservoir facies permeability curve model and facies permeability curve calculation method | |
| Yousfi et al. | The impact of helical slide honing on surface microtexture compared to plateau honing process through relevant characterization methods | |
| CN110084431A (en) | A kind of shale gas well yield analyzing and predicting method and system | |
| Pintaude et al. | Effect of surface finishing on the elastic contact area evaluation of austempered ductile iron | |
| CN109060231A (en) | A kind of measuring device and method for texturing friction surface lubrication fluid pressure | |
| Almeida et al. | Geometrical shape improvement of steel moulds by robot polishing process for polymer optic replication | |
| CN110895628B (en) | Performance degradation model prediction precision verification method | |
| CN113392372A (en) | Gas well full-life-cycle unobstructed flow determination method and system | |
| CN113810025B (en) | Friction coefficient self-adaptive trend filtering processing method | |
| CN111781643B (en) | Seismic interval velocity determination method and device | |
| CN113466938B (en) | Integral method prestack depth migration method and device |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |