CN109063328A - A kind of acquisition methods of rough surface bearing area rate curve - Google Patents
A kind of acquisition methods of rough surface bearing area rate curve Download PDFInfo
- Publication number
- CN109063328A CN109063328A CN201810860990.4A CN201810860990A CN109063328A CN 109063328 A CN109063328 A CN 109063328A CN 201810860990 A CN201810860990 A CN 201810860990A CN 109063328 A CN109063328 A CN 109063328A
- Authority
- CN
- China
- Prior art keywords
- bearing area
- area rate
- apparent height
- smc
- height
- 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.)
- Granted
Links
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 present invention relates to a kind of acquisition methods of rough surface bearing area rate curve.It include: that measurement rough surface three-dimensional geometry pattern is gone forward side by side line number Data preprocess；Determine the bearing area rate of apparent height to be calculatedp i ,i=1,2,3,…；To eachp i , corresponding apparent height is calculated using searching algorithmc i ,i=1,2,3,…；According to calculate it is resulting series (p i ,c i ) value, obtain bearing area rate curve.The present invention is by bearing area ratep i Gauging surface heightc i , in the nucleus of the bearing area rate curve of acquisition, under the same conditions, the resolution ratio of bearing area rate and apparent height corresponding relationship is much higher than conventional method to calculation amount；The present invention calculates the apparent height at given bearing area rate using searching algorithm, and computational accuracy is higher.
Description
Technical field
The invention belongs to the analyses of machining surface geometrical morphology and characterization technique field, in particular to a kind of rough surface
The acquisition methods of bearing area rate curve.
Background technique
Machinery part surface has a variety of different surface geometry patterns, various objects of these geometrical morphologies for surface
Reason, chemical property have very important influence, such as friction, abrasion, sealing, lubrication, thermally conductive, conductive, light reflected refraction
Deng.Studies have shown that only three-dimensional geometry pattern could completely embody the geometrical characteristic of rough surface, it is just able to achieve surface property
Prediction and control, therefore with the progress and development of science and technology, the three-dimensional geometry topography measurement of rough surface and it is characterized in section
Learn research and industry in using it is more and more extensive and deeply.
In the surface three dimensions characterization parameter system such as Birminghan14, STANDSURF, ISO 25178, surface is both defined
Bearing area rate curve, the curve and the characterization parameter based on the curve definitions can reflect the branch of rough surface well
Holding property, adjustment characteristic, antiwear property, lubricating oil film Forming ability, lubricating oil oil storage capacity.
Currently, the method for obtaining surface bearing the area ratio curve are as follows: first equidistantly divide on surface along short transverse, so
After calculate corresponding surface bearing the area ratio at each height, obtain a series of (bearing area rate, apparent height) data points, finally
Using bearing area rate as horizontal axis, apparent height is the longitudinal axis, draws surface bearing the area ratio curve.However, part rough surface passes through
After over mechanical processing (especially finishing), the height of most of surface sampled points is often distributed in a relatively narrow section
(core space).But in some places there may be individual particularly pertinent high points, or especially deep deep valley, and it is a other prominent
High point often has no too many substantial effect to surface property out.Therefore it is using the major defect of existing method:
Disadvantage 1: existing method seeks bearing area rate by apparent height.But it is if desired reversed to solve, that is, pass through given bearing
The area ratio seeks apparent height, needs to use interpolation side on the basis of obtaining whole (bearing area rate, apparent height) data points
Method establishes the functional relation of bearing area rate and apparent height, finally can just find out the height on surface.Therefore this method is not clever enough
It is living, and the computational accuracy of apparent height is related to data point amount, it is difficult to adjustment, EQUILIBRIUM CALCULATION FOR PROCESS precision and calculating speed as needed
Between relationship.
Disadvantage 2: surface core space is the important area for characterizing surface property, divides along short transverse, then counts at equal intervals
The method for calculating bearing area rate often makes bearing area rate curve few in the characterization point quantity of core space, characterizes bearing area
The lack of resolution of rate and apparent height variation relation；And the characterization resolution ratio mistake of the prominent position for accidentally occurring on surface
Greatly, computing resource is wasted.
Summary of the invention
The purpose of the present invention is to provide a kind of acquisition methods of rough surface bearing area rate curve, can be with bearing surface
Product rate is independent variable, and apparent height is dependent variable, is obtained to realize and carry out equidistant or any spacing division to bearing area rate
Take bearing area rate curve.
To achieve the above object, the technical scheme is that a kind of acquisition side of rough surface bearing area rate curve
Method includes the following steps:
Step S1, using searching algorithm, the corresponding apparent height when the bearing area rate on surface is preset value is calculated；
Step S2, using bearing area rate as independent variable, corresponding apparent height is dependent variable, obtains bearing area rate curve.
In an embodiment of the present invention, the step S1 is specifically included:
Step S11, input parameter is the bearing area rate on surfacep i , result to be asked is corresponding apparent heightc i ；
Step S12, apparent height is setc i The region of search left end point variableclIt indicates, initial value is surface sampled point
The maximum value of height；
Step S13, apparent height is setc i The region of search right endpoint variablecrIt indicates, initial value is surface sampled point
The minimum value of height；
Step S14, apparent height to be asked is setc i Use variablesmcIndicate, and be arranged its initial value be (cl + cr) / 2；
Step S15, into cyclic search, if apparent height issmcLocate bearing area rate calculated value withp i It is equal, thensmcIt is just when bearing area rate isp i When apparent heightc i , terminate to calculate；It is no to then follow the steps S16；
If apparent height is step S16, in step S15smcThe calculated value of place's bearing area rate is greater thanp i , then search is set
The right endpoint in sectioncrIt is equal tosmc；If in step S15, apparent height issmcThe calculated value of place's bearing area rate is less thanp i ,
The left end point of the region of search is then setclIt is equal tosmc；It enablessmcValue be updated to (cl + cr) / 2；Repeat step S15.
In an embodiment of the present invention, the step S2 is implemented as follows:
The rough surface bearing area rate curve of acquisition is with bearing area ratepFor independent variable, apparent heightcFor dependent variable；Bearing
The area ratiopEquidistantly or with the division of arbitrary interval, branch is denoted asp i ,iFor natural number.
Compared to the prior art, the invention has the following advantages:
(1) when obtaining surface bearing the area ratio curve, the present invention can carry out at equal intervals bearing area rate, or be appointed
Meaning interval divides；
(2) when according to given surface bearing the area ratio gauging surface height, computational accuracy can be adjusted according to demand, is reached
To the balance of computational accuracy and calculating speed.
Detailed description of the invention
Fig. 1 is general flow chart of the invention.
Fig. 2 is the flow chart in the present invention according to bearing area rate gauging surface height.
Fig. 3 is the sampled point altitude information figure of the polished surface after data prediction.
Fig. 4 is the surface bearing the area ratio curve graph obtained with the method for the invention.
Fig. 5 is the surface bearing the area ratio curve graph obtained with existing method.
Specific embodiment
Technical solution of the present invention is specifically described in 1-5 with reference to the accompanying drawing.
The present invention provides a kind of acquisition methods of rough surface bearing area rate curve, include the following steps:
Step S1, using searching algorithm, the corresponding apparent height when the bearing area rate on surface is preset value is calculated；
Step S2, using bearing area rate as independent variable, corresponding apparent height is dependent variable, obtains bearing area rate curve.
The step S1 is specifically included:
Step S11, input parameter is the bearing area rate on surfacep i , result to be asked is corresponding apparent heightc i ；
Step S12, apparent height is setc i The region of search left end point variableclIt indicates, initial value is surface sampled point
The maximum value of height；
Step S13, apparent height is setc i The region of search right endpoint variablecrIt indicates, initial value is surface sampled point
The minimum value of height；
Step S14, apparent height to be asked is setc i Use variablesmcIndicate, and be arranged its initial value be (cl + cr) / 2；
Step S15, into cyclic search, if apparent height issmcLocate bearing area rate calculated value withp i It is equal, thensmcIt is just when bearing area rate isp i When apparent heightc i , terminate to calculate；It is no to then follow the steps S16；
If apparent height is step S16, in step S15smcThe calculated value of place's bearing area rate is greater thanp i , then search is set
The right endpoint in sectioncrIt is equal tosmc；If in step S15, apparent height issmcThe calculated value of place's bearing area rate is less thanp i ,
The left end point of the region of search is then setclIt is equal tosmc；It enablessmcValue be updated to (cl + cr) / 2；Repeat step S15.
The step S2 is implemented as follows: the rough surface bearing area rate curve of acquisition is with bearing area ratepFor certainly
Variable, apparent heightcFor dependent variable；Bearing area ratepEquidistantly or with the division of arbitrary interval, branch is denoted asp i ,iFor nature
Number.
The following are specific implementation processes of the invention.
The specific implementation steps are as follows for a kind of acquisition methods of rough surface bearing area rate curve of the present invention:
(1) measurement data of machinery part surface geometrical morphology is obtained, which is the height that measuring device samples surface；
(2) surface data obtained to measurement pre-processes, and pretreated purpose is: measuring table, sample when filtering out measurement
Form error caused by itself and other factors；Filter out the high-frequency noise inside measurement data；By all sampling point heights in surface
Average value be adjusted to 0；
(3) it is write according to open source literature [1] with apparent height as input parameter, calculates the surface bearing area at corresponding height
The program of rate；
(4) surface bearing the area ratiopValue range be 0 to 100%, according to demand it is carried out at equal intervals or arbitrary interval
It divides, the bearing area rate of each branch is denoted asp i , iFor natural number；
(5) to eachp i , step (6)-(10) are repeated, to calculate corresponding apparent heightc i ；
(6) apparent height is setc i The region of search left end point variableclIt indicates, initial value is that surface samples point height
Maximum value；
(7) apparent height is setc i The region of search right endpoint variablecrIt indicates, initial value is that surface samples point height
Minimum value；
(8) apparent height to be asked is setc i Use variablesmcIndicate, and be arranged its initial value be (cl + cr) / 2；
(9) enter cyclic search, if apparent height issmcLocate bearing area rate calculated value withp i It is equal, thensmcJust
To be when bearing area ratep i When apparent height, end loop search；It is no to then follow the steps (10)；
(10) if apparent height is in step (9)smcThe calculated value of place's bearing area rate is greater thanp i , then the region of search is set
Right endpointcrIt is equal tosmc；If apparent height is in step (9)smcThe calculated value of place's bearing area rate is less thanp i , then set
Set the left end point of the region of searchclIt is equal tosmc；It enablessmcValue be updated to (cl + cr) / 2；Repeat step (9)；
(11) according to a series of (the bearing area rates being calculatedp i , apparent heightc i ) data pair, between data point into
Row piecewise linear interpolation obtains the bearing area rate curve of rough surface.
By taking polished surface sample shown in Fig. 3 as an example, it is specifically described acquisition rough surface bearing area proposed by the present invention
The method of rate curve.
The algorithm flow provided according to Fig. 1 and Fig. 2, the specific steps are as follows:
(1) rough surface three-dimensional appearance instrument is used, the raw measurement data of rough surface geometrical morphology is obtained；
(2) according to ISO 25178 or other in relation to document, filter out unwanted low frequency in surface measurement data and high frequency at
Point, and the average value of surface sampling point height is adjusted to 0；
(3) it according to open source literature [1], writes with apparent height as input parameter, calculates the surface bearing area at corresponding height
The program of rate；
(4) value range of bearing area rate is 0-100%, is divided into 20 equidistant sections, and be denoted asp i ,iFor certainly
So number；
(5) for divided in above-mentioned steps acquisition eachp i , execute step (6)-(10)；
(6) left end point of the setting apparent height region of search iscl, initial value is the maximum height value of surface sampled point；
(7) right endpoint of the setting apparent height region of search iscr, initial value is the minimum height values of surface sampled point；
(8) defined variablesmcIndicate apparent height to be asked, and enable its initial value be (cl+cr)/2；
(9) judge be when apparent heightsmcWhen bearing area rate calculated value whether withp i It is equal, if equal, at this timesmcExactly need apparent height to be requestedc i , exit circulation；It is no to then follow the steps (10)；
(10) if when apparent height issmcWhen bearing area rate be greater thanp i , enablecr=smc, otherwise enablecl=smc, updatesmc=
(cl+cr)/2 repeat step (9)；
(11) according to the calculated result of above-mentioned steps, the volume of data point for indicating (bearing area rate, apparent height) is obtained,
The bearing area rate curve (Fig. 4) of polished surface is drawn using subsection linearity inser value method；
The surface bearing the area ratio curve that comparison diagram 5(existing method obtains) and Fig. 4 can obviously embody it is of the invention beneficial
Nucleus of the data point that place, the i.e. present invention calculate greatly all in surface bearing the area ratio curve, can more meticulously reflect
The situation of change of bearing area rate and apparent height.
Bibliography:
[1] Zhou Chao, Liu Bingbing, Li Weifan wait to have the rough surface modeling method [J] of specified supporting ratio curve, in
State's engineering machinery journal, 2016,14 (2): 114-118.
The above, only embodiments of the present invention and embodiment, but embodiment be only it is exemplary, it is not right
The scope of the present invention constitutes any restrictions.It will be understood by those skilled in the art that without departing from spirit and model of the invention
Can be with the details and forms of the technical scheme of the invention are modified or replaced under enclosing, but these modifications and replacement each fall within this hair
Bright protection scope.
Claims (3)
1. a kind of acquisition methods of rough surface bearing area rate curve, which comprises the steps of:
Step S1, using searching algorithm, the corresponding apparent height when the bearing area rate on surface is preset value is calculated；
Step S2, using bearing area rate as independent variable, corresponding apparent height is dependent variable, obtains bearing area rate curve.
2. method according to claim 1, which is characterized in that the step S1 is specifically included:
Step S11, input parameter is the bearing area rate on surfacep i , result to be asked is corresponding apparent heightc i ；
Step S12, apparent height is setc i The region of search left end point variableclIt indicates, initial value is surface sampled point
The maximum value of height；
Step S13, apparent height is setc i The region of search right endpoint variablecrIt indicates, initial value is surface sampled point
The minimum value of height；
Step S14, apparent height to be asked is setc i Use variablesmcIndicate, and be arranged its initial value be (cl + cr) / 2；
Step S15, into cyclic search, if apparent height issmcLocate bearing area rate calculated value withp i It is equal, thensmc
It is just when bearing area rate isp i When apparent heightc i , terminate to calculate；It is no to then follow the steps S16；
If apparent height is step S16, in step S15smcThe calculated value of place's bearing area rate is greater thanp i , then the field of search is set
Between right endpointcrIt is equal tosmc；If in step S15, apparent height issmcThe calculated value of place's bearing area rate is less thanp i , then
The left end point of the region of search is setclIt is equal tosmc；It enablessmcValue be updated to (cl + cr) / 2；Repeat step S15.
3. method according to claim 2, which is characterized in that the step S2 is implemented as follows:
The rough surface bearing area rate curve of acquisition is with bearing area ratepFor independent variable, apparent heightcFor dependent variable；Bearing
The area ratiopEquidistantly or with the division of arbitrary interval, branch is denoted asp i ,iFor 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 true CN109063328A (en) | 2018-12-21 |
| CN109063328B 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 |
|---|
| 田红亮等: "机床支撑地脚结合部法向粗糙接触建模", 《万方数据期刊库》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109063328B (en) | 2023-02-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Blateyron | The areal field parameters | |
| Challis | A procedure for the automatic determination of filter cutoff frequency for the processing of biomechanical data | |
| CN110889091B (en) | Machine tool thermal error prediction method and system based on temperature sensitive interval segmentation | |
| CN103292723B (en) | The assay method of furrowfilter rod characteristic parameter | |
| CN109063328A (en) | A kind of acquisition methods of rough surface bearing area rate curve | |
| CN114964219A (en) | Hybrid EMD algorithm based on parameter optimization | |
| CN108254038B (en) | A kind of crude oil storage tank oil-water interfaces data go the calculation method of puppet and level gauging | |
| CN116738353B (en) | Pharmaceutical workshop air filter element performance detection method based on data analysis | |
| CN109190081A (en) | A kind of acquisition methods of Rough Surface Profile bearing length ratio curve | |
| CN109614758B (en) | Method for monitoring circular shape error with spatial correlation | |
| CN111339667A (en) | Processing method for denoising pressure and temperature monitoring data in injection molding production process | |
| CN109426216B (en) | Rectangular plane machining depth identification method for numerical control milling machine | |
| CN114167343B (en) | Crystal oscillator frequency adjustment method based on exhaustion intelligent electric meter | |
| CN105133472B (en) | Road surface roughness measurement method and inertial detector | |
| CN110020000A (en) | Determine the method and apparatus of abnormal air speed data | |
| CN112836166A (en) | Response peak analysis-based first-order differential algorithm for monitoring data of equal sampling experiment | |
| CN109029334B (en) | Rock mass structural plane roughness coefficient size effect global search measurement method | |
| Wright et al. | Wear rate calculation in the four-ball wear test | |
| CN109387222B (en) | Gyro circuit dynamic parameter measuring device and method | |
| Gogolewski et al. | Application of wavelet transform to determine surface texture constituents | |
| CN108154292B (en) | Method and system for establishing prediction model of road surface performance | |
| CN111160489A (en) | Multi-dimensional benchmarking analysis server, system and method based on big data and electronic terminal | |
| CN111121615B (en) | Phase shift interference fringe pattern batch selection method based on Hilbert transform | |
| CN117408913B (en) | Method, system and device for denoising point cloud of object to be measured | |
| CN117315002B (en) | Method, system and device for measuring three-dimensional phenotype of object to be measured |
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 |