CN108897904A - Load based on gear structure damage number uses regularity of distribution acquisition methods - Google Patents
Load based on gear structure damage number uses regularity of distribution acquisition methods Download PDFInfo
- Publication number
- CN108897904A CN108897904A CN201810446212.0A CN201810446212A CN108897904A CN 108897904 A CN108897904 A CN 108897904A CN 201810446212 A CN201810446212 A CN 201810446212A CN 108897904 A CN108897904 A CN 108897904A
- Authority
- CN
- China
- Prior art keywords
- load
- gear
- damage
- torque
- distribution
- 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.)
- Pending
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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
The present invention relates to a kind of load based on gear structure damage number to use regularity of distribution acquisition methods, is related to gear destruction analysis technical field.The present invention by calculating gear damage suffered by tired dangerous position under random loading first, then damage number is calculated, finally obtaining load using Non-parameter density estimation method is realized using the process of the regularity of distribution sample size is larger, in the lesser situation of window width, obtains the essential characteristic of load data distribution.
Description
Technical field
The present invention relates to gear destruction analysis technical fields, and in particular to it is a kind of based on gear structure damage number load make
With regularity of distribution acquisition methods.
Background technique
It is when the load cycle frequency is enough, i.e., special when load-time history long enough of the particular elements under typical condition
Determine stochastic variable sample size it is sufficiently large when, seem rambling load-time history will form certain load using point
Cloth rule, this regularity of distribution are that establishment loading spectrum has to the information known.
Summary of the invention
(1) technical problems to be solved
The technical problem to be solved by the present invention is to：The use regularity of distribution of load how is obtained, is provided for establishment loading spectrum
Support.
(2) technical solution
In order to solve the above-mentioned technical problems, the present invention provides a kind of load based on gear structure damage number to use distribution
Regular acquisition methods, include the following steps：
Step 1: calculating gear damage suffered by tired dangerous position under random loading；
Step 2: the calculated result based on step 1 calculates damage number；
Step 3: the result based on step 2, which obtains load using Non-parameter density estimation method, uses the regularity of distribution.
Preferably, step 1 is specially：
Assuming that the SN curve for manufacturing gear material is NSb=C, then gear fatigue under random loading is dangerous
Damage D suffered by position is：
Wherein, N is circulation cycle, and S is load amplitude, and b is the slope of curve, and C is that curve intercept, b and C join for gear material
Number, SiIndicate load change amplitude suffered by gear, TiIndicate the axis torque of corresponding gear load change；
In formula, it is assumed that in a linear relationship between the load amplitude S and axis torque T of tooth root fatigue dangerous point is S=KT, and K is
Preset constant, niTo be T corresponding to axis torque valueiRotation cycle；
From formula (1) as can be seen that since K and C is constant, so authentic representative load makes tooth root fatigue dangerous position
Magnitude at degree of injury is ∑iniTi b, which is known as to damage number.
Preferably, step 2 is specially：Obtain the information of each gear revolving speed and torque, wherein be cut into torque pre-
If sufficiently small section, the revolving speed fallen into respective bins integrates the time, acquires corresponding axis torque TiUnder effect
The cycle n of gear rotationi, obtain (ni, Ti) data pair, to calculate damage number.
Preferably, step 2 is seeking axis torque TiDuring, it is believed that negative torque and positive-torque are dangerous to tooth root fatigue
Damage is in same point caused by position.
Preferably, step 3 is specially：The equal amplitude Two dimensional Distribution of load is divided using the method for Non-parameter density estimation
Cloth fitting, the Non-parameter density estimation refer to, if stochastic variable x obeys a unknown probability density function f (x), benefit
With the sample x therefrom extracted1,...,xnProbability density function f (x) is estimated, wherein n indicates total number of samples.
Preferably, the method for Non-parameter density estimation is specially in step 3：It is assumed that haveing damage data x1,x2,...,xn
∈ [a, b), to section [a, b) make it is following divide, i.e. a=a0<a1<a2<...<ak=b, then enable Ii=[ai-1,ai), i=1,
2 ..., k, therefore haveI ≠ j enables mi=# { xi∈IiIndicate to fall in each section IiIn
Damage number data amount check, then Non-parameter density estimation function be：
Preferably, further include the steps that working out loading spectrum using the regularity of distribution using the load after step 3.
(3) beneficial effect
Then the present invention is counted by calculating gear damage suffered by tired dangerous position under random loading first
Damage number is calculated, finally obtaining load using Non-parameter density estimation method is realized using the process of the regularity of distribution in sample size
In larger, the lesser situation of window width, the essential characteristic of load data distribution is obtained.
Detailed description of the invention
Fig. 1 is rotation histogram schematic illustration of the invention；
Fig. 2 is the damage position schematic diagram of positive-torque and negative torque to gear structure.
Specific embodiment
To keep the purpose of the present invention, content and advantage clearer, with reference to the accompanying drawings and examples, to of the invention
Specific embodiment is described in further detail.
The present invention provides a kind of load based on gear structure damage number to use regularity of distribution acquisition methods, including following
Step：
Step 1: calculating gear damage suffered by tired dangerous position under random loading
Assuming that the SN curve for manufacturing associated gear material is NSb(in formula, N is circulation cycle to=C, and S is load width
Value, b are the slope of curve, and C is curve intercept, and b and C are gear material parameter) then gear fatigue danger under random loading
Damage D suffered by position is：
SiIndicate load change amplitude suffered by gear, TiIndicate the axis torque of corresponding gear load change.
In formula, it is assumed that linearly closed between stress variation amplitude (load amplitude) S and axis torque T of tooth root fatigue dangerous point
System is S=KT (K is preset constant), and it is T that nx, which is corresponding to axis torque value,iRotation cycle (namely gear is carried
Lotus changes corresponding axis torque TiEffect number).
From formula (1) as can be seen that since K and C is constant, so authentic representative load makes tooth root fatigue dangerous position
Magnitude at degree of injury is ∑iniTi b(unrelated with constant K, C), this extremely important numerical value are known as " damage number ".It is making
When determining bench test scheme, exactly rack will be arrived with the damage number " removing " for meeting design object with statistical representativeness
On, to complete relevant examination verifying.It can be seen that realize that this target only needs material parameter b can.
Step 2: the calculated result based on step 1 calculates damage number
In order to calculate gear destruction number, as shown in Figure 1, needing the information of each gear revolving speed and torque.Torque is cut
At sufficiently small section (bin), the revolving speed fallen into respective bins integrates the time, in the hope of corresponding approximately constant water
Flat axis torque TiAct on the cycle n of lower gear rotationi, to obtain sufficiently accurate (ni, Ti) data pair, to calculate damage
Number.
Load of the invention is known as to rotate histogram method using regularity of distribution acquisition methods.With rotation histogram meter
Number principle carries out counting statistics to revolving speed, the torque data of acquisition, calculates damage number, and the use for further obtaining load is distributed rule
Rule provides support for establishment loading spectrum.Wherein, since during gear drive, structure is in transmitting positive-torque and negative torque
When tooth root position tired dangerous point be different, as shown in Figure 2.Therefore, it generally requires and separates positive-torque with negative torque, point
It is not pocessed.The present invention is in the conservative consideration on acceptable degree, in steady-state rotation rack solution formulation, by negative torque
Caused by damage be added to caused by positive-torque in damage, that is, seeking axis torque TiDuring, it is believed that negative torque and positive twist
Square is damaged caused by tooth root fatigue dangerous position in same point.
Step 3: the result based on step 2, which obtains load using Non-parameter density estimation method, uses the regularity of distribution
The use regularity of distribution that load is extracted from actual measurement load data, needs to divide the equal amplitude Two dimensional Distribution of load
Cloth fitting.According to mathematical statistics, parametric method or nonparametric method is can be used in distribution statistics deduction.When the distribution nothing of sample data
When method hypothesis or distributional pattern complexity, parametric method will fail, and must use nonparametric method, that is, use Non-parameter density estimation method obtain
Load is taken to use the regularity of distribution.
Non-parameter density estimation refers to, if one unknown probability density function f (x) of stochastic variable x obedience is not (any
Specific form), to utilize the sample x therefrom extracted1,...,xnProbability density function f (x) is estimated, n indicates gross sample
This number.
Non-parameter density estimation method is the basis of histogram counts.It is assumed that haveing damage data x1,x2,...,xn∈[a,
B), to section [a, b) make it is following divide, i.e. a=a0<a1<a2<...<ak=b, then enable Ii=[ai-1,ai), i=1,2 ..., k.
Therefore havei≠j.Enable mi=# { xi∈IiIndicate to fall in each section IiIn damage number
Data amount check, then Non-parameter density estimation function is as follows：
The above is only a preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art
For member, without departing from the technical principles of the invention, several improvement and deformations can also be made, these improvement and deformations
Also it should be regarded as protection scope of the present invention.
Claims (7)
1. a kind of load based on gear structure damage number uses regularity of distribution acquisition methods, which is characterized in that including following step
Suddenly：
Step 1: calculating gear damage suffered by tired dangerous position under random loading；
Step 2: the calculated result based on step 1 calculates damage number；
Step 3: the result based on step 2, which obtains load using Non-parameter density estimation method, uses the regularity of distribution.
2. the method as described in claim 1, which is characterized in that step 1 is specially：
Assuming that the SN curve for manufacturing gear material is NSb=C, then gear tired dangerous position under random loading
Suffered damage D is：
Wherein, N is circulation cycle, and S is load amplitude, and b is the slope of curve, and C is curve intercept, and b and C are gear material parameter, Si
Indicate load change amplitude suffered by gear, TiIndicate the axis torque of corresponding gear load change；
In formula, it is assumed that in a linear relationship between the load amplitude S and axis torque T of tooth root fatigue dangerous point is S=KT, and K is default
Constant, niTo be T corresponding to axis torque valueiRotation cycle；
From formula (1) as can be seen that since K and C is constant, so authentic representative load causes to damage to tooth root fatigue dangerous position
The magnitude for hurting degree is ∑i niTi b, which is known as to damage number.
3. method according to claim 2, which is characterized in that step 2 is specially：Obtain each gear revolving speed and torque
Information, wherein torque is cut into preset sufficiently small section, the revolving speed fallen into respective bins accumulates the time
Point, acquire corresponding axis torque TiAct on the cycle n of lower gear rotationi, obtain (ni, Ti) data pair, to calculate damage number.
4. method as claimed in claim 3, which is characterized in that step 2 is seeking axis torque TiDuring, it is believed that negative torque
It is damaged caused by tooth root fatigue dangerous position in same point with positive-torque.
5. method as claimed in claim 3, which is characterized in that step 3 is specially：Using the method for Non-parameter density estimation
Fitting of distribution is done to the equal amplitude Two dimensional Distribution of load, the Non-parameter density estimation refers to, if stochastic variable x obedience one is unknown
Probability density function f (x), be using the sample x that therefrom extracts1,...,xnProbability density function f (x) is estimated,
Middle n indicates total number of samples.
6. method as claimed in claim 5, which is characterized in that the method for Non-parameter density estimation is specially in step 3：It is false
Surely it has damage data x1,x2,...,xn∈ [a, b), to section [a, b) make it is following divide, i.e. a=a0<a1<a2<...<ak=
B then enables Ii=[ai-1,ai), i=1,2 ..., k, therefore haveI ≠ j enables mi=# { xi
∈IiIndicate to fall in each section IiIn damage number data amount check, then Non-parameter density estimation function be：
7. such as method described in any one of claims 1 to 6, which is characterized in that after step 3 further include using described
The step of load is using regularity of distribution establishment loading spectrum.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810446212.0A CN108897904A (en) | 2018-05-11 | 2018-05-11 | Load based on gear structure damage number uses regularity of distribution acquisition methods |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810446212.0A CN108897904A (en) | 2018-05-11 | 2018-05-11 | Load based on gear structure damage number uses regularity of distribution acquisition methods |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108897904A true CN108897904A (en) | 2018-11-27 |
Family
ID=64342703
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810446212.0A Pending CN108897904A (en) | 2018-05-11 | 2018-05-11 | Load based on gear structure damage number uses regularity of distribution acquisition methods |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108897904A (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103868692A (en) * | 2014-03-18 | 2014-06-18 | 电子科技大学 | Rotary machine fault diagnosis method based on kernel density estimation and K-L divergence |
| CN106777732A (en) * | 2016-12-27 | 2017-05-31 | 电子科技大学 | A kind of Analysis Method of Fatigue Reliability for considering random load and strength degradation |
| CN106886638A (en) * | 2017-01-22 | 2017-06-23 | 北京理工大学 | A kind of endless-track vehicle power transmission shaft loading spectrum preparation method based on Density Estimator |
-
2018
- 2018-05-11 CN CN201810446212.0A patent/CN108897904A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103868692A (en) * | 2014-03-18 | 2014-06-18 | 电子科技大学 | Rotary machine fault diagnosis method based on kernel density estimation and K-L divergence |
| CN106777732A (en) * | 2016-12-27 | 2017-05-31 | 电子科技大学 | A kind of Analysis Method of Fatigue Reliability for considering random load and strength degradation |
| CN106886638A (en) * | 2017-01-22 | 2017-06-23 | 北京理工大学 | A kind of endless-track vehicle power transmission shaft loading spectrum preparation method based on Density Estimator |
Non-Patent Citations (6)
| Title |
|---|
| DAN XU 等: "Fatigue Damage Mechanism-Based Dependent Modeling With Stochastic Degradation and Random Shocks", 《IEEE ACCESS》 * |
| 刘彦龙: "汽车传动系动态载荷谱提取与台架试验载荷谱编制", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
| 刘海鸥等: "基于核密度估计的履带车辆传动轴载荷谱编制", 《兵工学报》 * |
| 周勃等: "改进二阶统计量算法在风力机主轴承裂纹识别中的应用", 《振动工程学报》 * |
| 梁佳 等: "装载机驱动桥疲劳试验扭矩加载谱编制方法研究", 《中国机械工程》 * |
| 陈道云等: "一种新的高速列车动应力谱分布估计方法", 《机械工程学报》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112735094B (en) | Geological disaster prediction method and device based on machine learning and electronic equipment | |
| CN108319664A (en) | A kind of dam and the recognition methods of project security monitoring data error and system | |
| CN104200113A (en) | Internet of Things data uncertainty measurement, prediction and outlier-removing method based on Gaussian process | |
| CN112381476A (en) | Method and device for determining electric energy meter with abnormal state | |
| CN115600119B (en) | Data processing method and system suitable for wind power generation | |
| CN112288597A (en) | Energy consumption online anomaly detection method based on hierarchical clustering and histogram algorithm | |
| CN102706757B (en) | Multi-axle fatigue analyzing method and application thereof | |
| CN107967571A (en) | Index system generation method based on entropy assessment | |
| CN107209502A (en) | In the case where reducing control based on the method for manufacturing part to the analysis of statistical indicator | |
| CN110226140A (en) | State monitoring method and state monitoring apparatus | |
| CN117113236B (en) | Smart city monitoring system and data processing method | |
| CN102539433B (en) | Nondestructive testing method for quality of leechee seed | |
| CN115936532A (en) | Saline-alkali soil stability assessment method and system based on BP neural network | |
| CN108897904A (en) | Load based on gear structure damage number uses regularity of distribution acquisition methods | |
| CN106407488B (en) | Power failure methods of risk assessment and system | |
| CN108182306B (en) | Method for determining degradation failure threshold of abrasive particle characteristic parameters of vehicle power transmission device | |
| CN110020000B (en) | Method and device for judging abnormal wind speed data | |
| CN104778709B (en) | A kind of construction method of the electronic blackboard based on yarn sequence image | |
| CN106021855A (en) | A reactor period calculating method | |
| CN115496143A (en) | Method and device for detecting abnormity of blast furnace thermocouple temperature measurement data and storage medium | |
| CN108874732B (en) | Method for judging and detecting pure phase similarity of dendrobium aphyllum | |
| CN110010202B (en) | Establishment and judgment standard and judgment method for judging similarity of pure varieties of dendrobium fimbriatum | |
| CN111474293B (en) | Method and system for determining bacterial wilt solution | |
| CN115169751A (en) | Method and system for improving maintenance quality of traditional Chinese medicine decoction pieces | |
| CN118076929A (en) | Monitoring production of material boards, in particular engineered wood boards, in particular using self-organizing maps |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181127 |