CN111597740A - Harvester health monitoring method based on mesoscale ultrasonic teratocardio-band - Google Patents

Harvester health monitoring method based on mesoscale ultrasonic teratocardio-band Download PDF

Info

Publication number
CN111597740A
CN111597740A CN202010292707.XA CN202010292707A CN111597740A CN 111597740 A CN111597740 A CN 111597740A CN 202010292707 A CN202010292707 A CN 202010292707A CN 111597740 A CN111597740 A CN 111597740A
Authority
CN
China
Prior art keywords
harvester
mesoscale
ultrasonic energy
value
damage
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
Application number
CN202010292707.XA
Other languages
Chinese (zh)
Other versions
CN111597740B (en
Inventor
朱林
郭广明
孙进
王鹏
黄嘉铭
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Yangzhou University
Original Assignee
Yangzhou University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Yangzhou University filed Critical Yangzhou University
Priority to CN202010292707.XA priority Critical patent/CN111597740B/en
Publication of CN111597740A publication Critical patent/CN111597740A/en
Application granted granted Critical
Publication of CN111597740B publication Critical patent/CN111597740B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • A—HUMAN NECESSITIES
    • A01—AGRICULTURE; FORESTRY; ANIMAL HUSBANDRY; HUNTING; TRAPPING; FISHING
    • A01D—HARVESTING; MOWING
    • A01D41/00—Combines, i.e. harvesters or mowers combined with threshing devices
    • A01D41/12—Details of combines
    • A01D41/127—Control or measuring arrangements specially adapted for combines
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

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)
  • Life Sciences & Earth Sciences (AREA)
  • Environmental Sciences (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

The invention discloses a harvester health monitoring method based on a mesoscale ultrasonic teratocardio band, which comprises the following steps: s1, determining the damage position of the harvester structure; s2, performing criticality screening on the mesoscale region in the damage position; s3, collecting ultrasonic energy in the mesoscale area; s4, solving the ultrasonic energy of the medium-scale area in a non-damage state to obtain an ultrasonic energy distortion band; and S5, calculating the health state of the harvester based on the mesoscale ultrasonic energy teratocarcinoma band. The monitoring method is real-time, efficient and accurate.

Description

Harvester health monitoring method based on mesoscale ultrasonic teratocardio-band
Technical Field
The invention relates to a harvester health monitoring method, in particular to a harvester health monitoring method based on a mesoscale ultrasonic teratocarcinoma band.
Background
With the deepening of agricultural large-scale production, the continuous operation time and the working condition severity of the harvester continuously refresh the history, and the safety in the operation process is more and more emphasized. Minor damage can be a significant hazard due to the continuity of the work cycle and deterioration of the climatic environment. Therefore, real-time monitoring of the health of the harvester during operation is essential. At present, most of health monitoring on the structure of the harvester is carried out by means of nondestructive inspection under the stop off-line state, and the health state is not monitored by ultrasonic indexes with more detailed damage in a medium-scale damage mode under the real-time working condition, so that the monitoring precision is low. Therefore, it is necessary to monitor the health status of the harvester on-line from the ultrasonic refined index of mesoscale damage.
Disclosure of Invention
The purpose of the invention is as follows: the invention aims to provide a harvester health monitoring method based on a mesoscale ultrasonic teratocardio band.
The technical scheme is as follows: the invention provides a harvester health monitoring method based on a mesoscale ultrasonic teratocarcinoma, which comprises the following steps:
s1, determining the damage position of the harvester structure;
s2, performing criticality screening on the mesoscale region in the damage position;
s3, collecting ultrasonic energy in the mesoscale area;
s4, solving the ultrasonic energy of the medium-scale area in a non-damage state to obtain an ultrasonic energy distortion band;
and S5, calculating the health state of the harvester based on the mesoscale ultrasonic energy teratocarcinoma band.
Further, the method for determining the damage location in S1 includes: establishing a three-dimensional model of the harvester, guiding the three-dimensional model into finite element analysis software, dividing grids of the harvester by adopting an automatic mode, defining and setting corresponding constraint conditions in a pre-processing module, applying corresponding load to the harvester model, solving the static stress and strain of the harvester model in a post-processing module after the constraint conditions are finished, and determining the damage position and the stress distribution condition of the damage position of the harvester structure under the working condition.
Further, the screening method in S2 is as follows: based on the analysis result of S1, cloth is arranged at the position of the injuryThe acoustic emission device is arranged at the geometric center of the damage position; determining an impact load loading area according to the analysis result of S1, wherein the selected position of the area is the area corresponding to the minimum stress value in the finite element analysis result, applying impact excitation load at the selected load point position, and using an acoustic emission device to impact quantity C at the selected load point positioniDetecting, i is the number of the damage position, i is 1, 2, 3.. N, N is the total number of the damage position, calculating the critical impact value J of the mesoscale,
Figure BDA0002449289640000021
in the formula, J is a mesoscale critical impact value in the structure of the harvester to be analyzed; n is the total number of the damage positions; ciThe impact quantity of each damage position in the harvester structure is calculated, the structural load applicability coefficient is calculated and is 1.2 when the borne load is a composite load, β is a welding structure quality coefficient, the harvester is classified according to the service life, β is 1 when the service life is less than 1 year, β is 0.9 when the service life is 1-3 years, β is 0.8 when the service life is 3-5 years, β is 0.7 when the service life is more than 5 years, and delta K is a stress intensity factor corresponding to the damage position in the harvester structure, the impact quantity C is obtained by a finite element analysis method, and the impact quantity C is calculated and is calculated at the ith damage positioniIf the critical impact value J is larger than the critical impact value J of the mesoscale, the damage position with the number i is considered as the mesoscale region.
Further, the ultrasonic energy acquisition method of the S3 meso-scale region: arranging ultrasonic device in the mesoscale region determined in S2, and setting ultrasonic energy value I at position IiCollecting, using the vibration of the harvester in the working process as an input impact load, and calculating the maximum value of the input vibration impact load which can be generated by the harvester under the working condition by using a finite element means, and recording the maximum value as Fmax。
Further, the method for solving the ultrasonic energy of the non-damage state of the medium scale region in S4 to obtain the ultrasonic energy distortion band includes: input vibration impact determined in S3Maximum value of load FmaxOn the basis, the load is led into a finite element analysis model of the harvester, the loading point position is the region corresponding to the minimum stress value determined in S1, and then after the constraint and post-processing conditions are set, the change interval of the ultrasonic energy value of the mesoscale region is solved in an acoustic analysis module of finite element analysis software, namely 1.37Q is takeniTo 2.11QiIn the interval of ultrasonic energy, QiThe maximum ultrasonic energy value of the ith position in a non-damage state is obtained through finite element means analysis.
Further, in the S5, the method for calculating the health status of the harvester based on the mesoscale ultrasonic energy teratocardio band is as follows: based on S3 and S4, the ultrasonic energy value I of the mesoscale region is measured in real timei、1.37Qi、2.11QiSubstituting the following formula into the health state value Z of the ith areaiCalculating, i is the serial number of the damage position,
Figure BDA0002449289640000022
in the formula, ZiIs the health status value of the ith area; qiThe maximum ultrasonic energy value of the ith position in a non-damage state is obtained through finite element means analysis; i isiIs the ultrasonic energy value of the mesoscale region; i is the number of the damage position; and delta K is a stress intensity factor corresponding to a damage position in the harvester structure and is obtained by a finite element analysis method.
Has the advantages that: the invention can realize the real-time monitoring of the health state of the harvester under the working condition, roughly determine the damage area through the ultrasonic detection of the mesoscale damage unit, and further estimate the health state value of the harvester through the abnormal narrow band interval value in the ultrasonic index, thereby being beneficial to accurately acquiring the health state value of the harvester under the working condition in real time and further more effectively ensuring the safety of the harvester in the operation process.
Drawings
FIG. 1 is a flow chart of the method of the present invention.
Detailed Description
As shown in fig. 1, the method of the present embodiment includes the following steps:
s1, determining a structural damage position of a harvester:
establishing a three-dimensional model of the harvester, guiding the established three-dimensional model into finite element analysis software, dividing grids of the three-dimensional model by adopting an automatic mode, defining and setting corresponding constraint conditions in a pre-processing module, applying corresponding load to the harvester model, solving the static stress and strain of the harvester model in a post-processing module after the completion of the operation, and determining the damage position and the stress distribution condition of the damage position of the harvester structure under the working condition.
S2, criticality screening of the medium scale area in the damage position:
on the basis of the analysis result in S1, an acoustic emission device with general characteristics is arranged in the damage position determined by the finite element analysis, and the arrangement position of the acoustic emission sensor is the geometric center of each damage position determined in S1. In the results of the analysis at S1, an impact load loading zone is determined, the selected location of the zone being the zone corresponding to the minimum stress value in the finite element analysis results. Then applying 150N impact excitation load at the selected load point position, and applying the acoustic emission sensor with general characteristics to impact quantity C at the positioniDetection is performed (i is the number of the damaged position, i is 1, 2, 3.. N, N is the total number of damaged positions).
Calculating the critical impact value J of the mesoscale,
Figure BDA0002449289640000031
wherein J is a mesoscale critical impact value in the structure of the harvester to be analyzed; n is the total number of the damage positions; ciThe impact quantity of each damage position in the harvester structure, the structural load applicability coefficient, when the borne load is a composite load, the value is 1.2, β is a welding structure quality coefficient, the harvester is classified according to the service life, when the service life is less than 1 year, the value β is 1, when the service life is 1-3 years, the value β is 0.9,when the service life is 3-5 years, β is 0.8, when the service life is more than 5 years, β is 0.7, and delta K is a stress intensity factor corresponding to a damage position in the harvester structure and is obtained by a finite element analysis method.
After the calculation, the impact quantity C of the ith damage positioniIf the critical impact value J is larger than the critical impact value J of the mesoscale, the damage position with the number i is considered as the mesoscale region.
S3, collecting ultrasonic energy in a mesoscale area:
arranging the transducer of the ultrasonic device in the mesoscale region determined in S2 for the ultrasonic energy I at the I-th positioni(i is a number) was collected. The vibration of the harvester in the working process is used as an input impact load, and the maximum value of the input vibration impact load which can be generated by the harvester under the working condition is calculated by using a finite element method and is recorded as Fmax。
S4, solving the ultrasonic energy of the medium-scale area in the non-damage state:
input vibration impact load maximum value F determined in S3maxOn the basis, the load is introduced into a finite element analysis model of the harvester, and the loading point position is the region corresponding to the minimum stress value determined in the step S1. Then, after the constraint and post-processing conditions are set, solving the change interval of the ultrasonic energy value in the mesoscale area in an acoustic analysis module of finite element analysis software, namely taking 1.37QiTo 2.11QiIn the interval of ultrasonic energy (Q)iMaximum ultrasound energy value in a non-damaged state for the ith position analyzed by finite element means).
S5, calculating the health state of the harvester based on the mesoscale ultrasonic teratocardio band:
based on S3 and S4, the ultrasonic energy value I of the mesoscale region is measured in real timei、1.37Qi、2.11QiSubstituting the following formula into the health state value Z of the ith areaiThe calculation is performed (i is the number of the damaged position),
Figure BDA0002449289640000041
wherein Z isiIs the health status value of the ith area; qiThe maximum ultrasonic energy value of the ith position in a non-damage state is obtained through finite element means analysis; i isiIs the ultrasonic energy value of the mesoscale region; i is the number of the damage position; and delta K is a stress intensity factor corresponding to a damage position in the harvester structure and is obtained by a finite element analysis method.

Claims (6)

1. A harvester health monitoring method based on a mesoscale ultrasonic teratocardio band is characterized in that: the method comprises the following steps:
s1, determining the damage position of the harvester structure;
s2, performing criticality screening on the mesoscale region in the damage position;
s3, collecting ultrasonic energy in the mesoscale area;
s4, solving the ultrasonic energy of the medium-scale area in a non-damage state to obtain an ultrasonic energy distortion band;
and S5, calculating the health state of the harvester based on the mesoscale ultrasonic energy teratocarcinoma band.
2. The harvester health monitoring method based on mesoscale ultrasound teratocarcinoma according to claim 1, characterized in that: the method for determining the damage position in S1 includes: establishing a three-dimensional model of the harvester, guiding the three-dimensional model into finite element analysis software, dividing grids of the harvester by adopting an automatic mode, defining and setting corresponding constraint conditions in a pre-processing module, applying corresponding load to the harvester model, solving the static stress and strain of the harvester model in a post-processing module after the constraint conditions are finished, and determining the damage position and the stress distribution condition of the damage position of the harvester structure under the working condition.
3. The harvester health monitoring method based on mesoscale ultrasound teratocarcinoma according to claim 2, characterized in that: the screening method in the S2 comprises the following steps: based on the results of the S1 analysis, acoustic emission is placed at the site of the lesionThe arrangement position of the shooting device is the geometric center of the damage position; determining an impact load loading area according to the analysis result of S1, wherein the selected position of the area is the area corresponding to the minimum stress value in the finite element analysis result, applying impact excitation load at the selected load point position, and using an acoustic emission device to impact quantity C at the selected load point positioniDetecting, i is the number of the damage position, i is 1, 2, 3.. N, N is the total number of the damage position, calculating the critical impact value J of the mesoscale,
Figure FDA0002449289630000011
in the formula, J is a mesoscale critical impact value in the structure of the harvester to be analyzed; n is the total number of the damage positions; ciThe impact quantity of each damage position in the harvester structure is calculated, the structural load applicability coefficient is calculated and is 1.2 when the borne load is a composite load, β is a welding structure quality coefficient, the harvester is classified according to the service life, β is 1 when the service life is less than 1 year, β is 0.9 when the service life is 1-3 years, β is 0.8 when the service life is 3-5 years, β is 0.7 when the service life is more than 5 years, and delta K is a stress intensity factor corresponding to the damage position in the harvester structure, the impact quantity C is obtained by a finite element analysis method, and the impact quantity C is calculated and is calculated at the ith damage positioniIf the critical impact value J is larger than the critical impact value J of the mesoscale, the damage position with the number i is considered as the mesoscale region.
4. The harvester health monitoring method based on mesoscale ultrasound teratocarcinoma of claim 3, characterized in that: the ultrasonic energy acquisition method of the S3 meso-scale region comprises the following steps: arranging ultrasonic device in the mesoscale region determined in S2, and setting ultrasonic energy value I at position IiCollecting, using the vibration of the harvester in the working process as an input impact load, and calculating the maximum value of the input vibration impact load which can be generated by the harvester under the working condition by using a finite element means, and recording the maximum value as Fmax。
5. The harvester health monitoring method based on mesoscale ultrasound teratocarcinoma of claim 4, characterized in that: the method for solving the ultrasonic energy of the nondestructive state of the medium scale area in the S4 to obtain the ultrasonic energy distortion band comprises the following steps: input vibration impact load maximum value F determined in S3maxOn the basis, the load is led into a finite element analysis model of the harvester, the loading point position is the region corresponding to the minimum stress value determined in S1, and then after the constraint and post-processing conditions are set, the change interval of the ultrasonic energy value of the mesoscale region is solved in an acoustic analysis module of finite element analysis software, namely 1.37Q is takeniTo 2.11QiIn the interval of ultrasonic energy, QiThe maximum ultrasonic energy value of the ith position in a non-damage state is obtained through finite element means analysis.
6. The harvester health monitoring method based on mesoscale ultrasound teratocarcinoma of claim 5, characterized in that: in the S5, the method for calculating the health state of the harvester based on the mesoscale ultrasonic energy teratocardio-band comprises the following steps: based on S3 and S4, the ultrasonic energy value I of the mesoscale region is measured in real timei、1.37Qi、2.11QiSubstituting the following formula into the health state value Z of the ith areaiCalculating, i is the serial number of the damage position,
Figure FDA0002449289630000021
in the formula, ZiIs the health status value of the ith area; qiThe maximum ultrasonic energy value of the ith position in a non-damage state is obtained through finite element means analysis; i isiIs the ultrasonic energy value of the mesoscale region; i is the number of the damage position; and delta K is a stress intensity factor corresponding to a damage position in the harvester structure and is obtained by a finite element analysis method.
CN202010292707.XA 2020-04-14 2020-04-14 Harvester health monitoring method based on mesoscale ultrasonic abnormal narrow band Active CN111597740B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010292707.XA CN111597740B (en) 2020-04-14 2020-04-14 Harvester health monitoring method based on mesoscale ultrasonic abnormal narrow band

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010292707.XA CN111597740B (en) 2020-04-14 2020-04-14 Harvester health monitoring method based on mesoscale ultrasonic abnormal narrow band

Publications (2)

Publication Number Publication Date
CN111597740A true CN111597740A (en) 2020-08-28
CN111597740B CN111597740B (en) 2023-04-25

Family

ID=72184979

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010292707.XA Active CN111597740B (en) 2020-04-14 2020-04-14 Harvester health monitoring method based on mesoscale ultrasonic abnormal narrow band

Country Status (1)

Country Link
CN (1) CN111597740B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112345642A (en) * 2020-10-28 2021-02-09 扬州大学 Harvester health monitoring method considering primary and secondary structure decomposition equivalence
CN113177337A (en) * 2021-04-20 2021-07-27 扬州大学 Reed harvester safety evaluation method based on correlation factor characteristic value fluctuation interval
CN114413804A (en) * 2020-10-28 2022-04-29 扬州大学 Method for determining optimal operation parameters of harvester based on local monitoring strain energy
CN115048731A (en) * 2021-09-09 2022-09-13 海沃机械（中国）有限公司 Health assessment method and device for pull arm of garbage transfer vehicle

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR3050530A1 (en) * 2016-04-20 2017-10-27 Safran METHOD OF CHECKING OR CERTIFYING A GLUE ASSEMBLY
CN108956668A (en) * 2018-07-23 2018-12-07 湖南大学 A kind of crack tip ope ning angle degree measurement method based on SEM in situ
CN109270170A (en) * 2018-11-21 2019-01-25 扬州大学 A kind of sensitivity amendment loading machine Structural Damage Identification considering Jie's scale
CN110008650A (en) * 2019-05-17 2019-07-12 杭州电子科技大学 A kind of modeling method of three-D photon crystal internal flaw molding positioning

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR3050530A1 (en) * 2016-04-20 2017-10-27 Safran METHOD OF CHECKING OR CERTIFYING A GLUE ASSEMBLY
CN108956668A (en) * 2018-07-23 2018-12-07 湖南大学 A kind of crack tip ope ning angle degree measurement method based on SEM in situ
CN109270170A (en) * 2018-11-21 2019-01-25 扬州大学 A kind of sensitivity amendment loading machine Structural Damage Identification considering Jie's scale
CN110008650A (en) * 2019-05-17 2019-07-12 杭州电子科技大学 A kind of modeling method of three-D photon crystal internal flaw molding positioning

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
路腾腾 等: "超声振动对纯钛TA1微圆柱体压缩过程的影响" *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112345642A (en) * 2020-10-28 2021-02-09 扬州大学 Harvester health monitoring method considering primary and secondary structure decomposition equivalence
CN114413804A (en) * 2020-10-28 2022-04-29 扬州大学 Method for determining optimal operation parameters of harvester based on local monitoring strain energy
CN112345642B (en) * 2020-10-28 2023-03-24 扬州大学 Harvester health monitoring method considering primary and secondary structure decomposition equivalence
CN114413804B (en) * 2020-10-28 2023-08-22 扬州大学 Harvester optimal operation parameter determination method based on local monitoring strain energy
CN113177337A (en) * 2021-04-20 2021-07-27 扬州大学 Reed harvester safety evaluation method based on correlation factor characteristic value fluctuation interval
CN113177337B (en) * 2021-04-20 2023-05-26 扬州大学 Reed harvester safety assessment method based on association factor characteristic value fluctuation interval
CN115048731A (en) * 2021-09-09 2022-09-13 海沃机械（中国）有限公司 Health assessment method and device for pull arm of garbage transfer vehicle
CN115048731B (en) * 2021-09-09 2023-10-31 海沃机械（中国）有限公司 Health assessment method and device for pull arm of garbage transfer truck

Also Published As

Publication number Publication date
CN111597740B (en) 2023-04-25

Similar Documents

Publication Publication Date Title
CN111597740A (en) Harvester health monitoring method based on mesoscale ultrasonic teratocardio-band
CN105510038B (en) Wind turbine generator fault monitoring method and device
JP6752866B2 (en) Cutting tool condition inspection system and method
US20140358453A1 (en) Robot reducer predictive trending
CA2764092A1 (en) Spot weld data management and monitoring system
CN115950609B (en) Bridge deflection anomaly detection method combining correlation analysis and neural network
CN111460702B (en) Structural member damage identification method based on forward and reverse damage feature fusion
CN115374959A (en) Health management method and system for aircraft engine
CN111354496B (en) Nuclear power plant accident online diagnosis and state tracking prediction method
CN109059813B (en) Method for detecting corrosion strength of steel structure of hoisting machinery
JP2020176883A (en) Maintenance management system and maintenance management method
CN113487149B (en) Welding spot abnormity identification system and method based on Catboost K-fold cross verification
CN111595939B (en) Harvester damage identification method based on ultrasonic energy spectrum
US11167378B1 (en) Techniques for determining weld quality
JP2965716B2 (en) Piping life diagnosis
CN109979622A (en) Nuclear power plant's voltage-stablizer fatigue life monitors assessment system and method on-line
CN110736671B (en) Method for monitoring abnormal part of pipe fitting hardness
JP4212419B2 (en) Mold crack initiation prediction system
JP2020153938A (en) Tile deterioration diagnosis method, estimation method for building wall surface repair cost, tile deterioration diagnosis system, and estimation system for building wall surface repair cost
CN109470773A (en) A kind of caliberating device and its working method of ultrasonic detecting probe
CN113449396B (en) Off-line inspection-based subcritical boiler drum body state evaluation method
US20220382247A1 (en) Estimation method and estimation system
Kollment et al. Force and acoustic emission measurements for condition monitoring of fine blanking tools
CN118260564A (en) Part life sampling inspection method and system based on Bayesian algorithm
CN113532898A (en) State evaluation method for steam-water separator body of supercritical or ultra-supercritical boiler

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