Wherein, f (j, beta)₁，β₂，β₃) The distribution equation of Weibull of the actually measured powder removing rotating speed is shown, j is the number of the rotating speed of the powder removing device, and j is 1, 2, 3. Beta is a₁Size parameters for the Weibull distribution equation; beta is a₂Is the shape parameter of Weibull distribution equation; beta is a₃Position parameters of Weibull distribution equation; e is the number of fingers.

Then calculating the corresponding value Z of the intermediate density by substituting f (j) into the following formula_jSo as to obtain the corresponding actual rotation speed of the powder removal

Then the rotation speed n between the two and the theoretical preset rotation speed is obtained through calculation₀The ratio of 200r/min is recorded as α

Wherein, f (j) is a Weibull distribution equation of the actually measured powder removing rotating speed; j is a type number for measuring the rotating speed of the powder removing device, and j is 1, 2, 3.. Q; q is the total number of types for measuring the actual rotating speed; c [ f (x)]_maxThe probability maximum value under the equation mode of all the actual rotating speeds is obtained; c [ f (x)]_minThe probability minimum value in the equation mode of all the actual rotating speeds is obtained; z_jIntermediate density value of actual rotation speed;

intermediate density value Z for actual rotation speed_jThe corresponding powder removing rotating speed; n is₀Taking 200r/min as a theoretically preset rotating speed; alpha is the ratio of the actual rotational speed to the theoretical rotational speed.

Further, calculating the deviation degree of the total mass center of the loose flower spike in the pollen removing device: the radius of the roller in the pollen removing device is R according to the measurement of S1, but due to the fact that the shapes of the pine pollen ears actually entering the inside of the roller are different and the mass distribution is not uniform, the coordinate of the ith pine pollen ear in the roller device in the plane is (x) detected by a coordinate measuring machine_i，y_i) N, wherein N is the total number of the picked pine pollen ears entering the roller, and the total mass center is (x) calculated through the mass center₀，y₀) Thereby obtaining the actual rolling radius d after correction

Wherein M is_iN, i is the mass of the i-th picking of the pine spikes measured by the mass sensor, 1, 2, 3.. N; n represents the total number of the picked pine pollen ears entering the roller; (x)₀，y₀) Calculating the position of the total mass center of the pine flower spike in the pollen removing device; (x)_i，y_i) The coordinates of the ith pine pollen ear in the plane in the roller device are obtained through the detection of a coordinate measuring instrument; n ═ 1, 2, 3.; n represents the total number of the picked pine pollen ears entering the roller; d is the actual rolling radius of the pine pollen after deviation correction; r is the inner radius of the drum.

Further, calculating the pollen pini dryness coefficient: the moisture content in the loose flower spikes is detected through a penetration type microwave moisture meter, the phenomenon that the flower spikes are easy to adhere in the rolling pollen removal process is avoided, the dryness of the loose flower spikes can cause certain influence on the pollen removal effect of the loose flower spikes, and the moisture content measured through the penetration type microwave moisture meter is recorded as k₁Then with an empirical standard k₀Comparing to obtain the dryness coefficient of the pine cones, and recording as k

Wherein k is the dryness coefficient of the pine pollen ears; k is a radical of₁Measured moisture content of the pine pollen ear, k₀Is the common standard moisture content of the pine pollen.

Further, determining the optimal rotation speed of pollen pini pollen removal based on deviation correction: according to the parameter data calculated and determined in S1-S4, the optimal rotation speed of pollen pini pollen removal after deviation correction is

Wherein n is the optimal rotation speed of pollen pini pollen removal after deviation correction; f is pollen pini completely removedAn optimal centripetal force; m_iThe mass of the i-th picking of the pine pollen spikes, i is 1, 2, 3.. N, and N is the total number of the picked pine pollen spikes entering the roller; d is the actual rolling radius of the pine pollen after deviation correction; alpha is the ratio of the actual rotating speed to the theoretical rotating speed; k is the pollen pini dryness coefficient.

Has the advantages that: according to the method, the optimal rotating speed of pollen pini pollen shedding is determined under the correction of a plurality of influence deviations, so that accurate picking and vibration screening can be realized for different types and heights of pines in pine forests, full-automatic integrated pollen pini collection pollen shedding work can be realized, the yield of pollen pini pollen shedding is greatly improved, the cost of collecting and producing seeds of pollen pini is reduced, and the pollution and waste of pollen pini are greatly reduced.

Drawings

FIG. 1 is a flow chart of the method of the present invention.

Detailed Description

The method for determining the optimal pine pollen ear pollen-removal rotating speed based on deviation correction comprises the following steps:

s1, determining main working parameters of pollen pini pollen removal:

the quality of the pine cones after each picking is measured by a quality sensor and recorded as M_iNumbering the sequence i of each picking, wherein i is 1, 2, 3.. N, and N is the total number of picked pine pollen ears entering the roller; obtaining the inner radius of the roller as R through actual measurement; and (3) introducing a three-dimensional model of the whole pollen-removing roller into Adams multi-body dynamics simulation software, setting boundary constraint and rolling load, then carrying out simulation calculation by using a post-processing module, and finally extracting the optimal centripetal force F for pollen emergence of the flower spikes.

S2, determining a powder removal rotating speed completion value based on Weibull distribution under an actual working condition:

pollen pini pollen removing operation is carried out according to the known parameters of S1, the actual rotating speed of the pollen pini pollen removing device is measured for multiple times by a rotating speed sensor before operation, the number of each measurement is j, and the actual rotating speed obtained by the jth measurement is recorded as n_j，j＝1，2，Q, Q is expressed as the total number of types for measuring the actual rotational speed thereof, and then compared with the theoretically preset rotational speed n₀Comparing at 200 r/min; after a plurality of times of measurement, all the detected actual rotating speeds of the powder removing device are led into MATLAB software to automatically fit the Weibull distribution curve

Wherein, f (j) is a Weibull distribution equation of the actually measured powder removing rotating speed; j is a type number for measuring the rotating speed of the powder removing device, and j is 1, 2, 3.. Q; q is the total number of types for measuring the actual rotating speed; c [ f (x)]_maxThe probability maximum value under the equation mode of all the actual rotating speeds is obtained; c [ f (x)]_minSummary of equation modes for all actual rotational speedsA minimum value of rate; z_jIntermediate density value of actual rotation speed;

S3, calculating the deviation degree of the total mass center of the loose flower spikes in the pollen removing device:

the radius of the roller in the pollen removing device is R according to the measurement of S1, but due to the fact that the shapes of the pine pollen ears actually entering the inside of the roller are different and the mass distribution is not uniform, the coordinate of the ith pine pollen ear in the roller device in the plane is (x) detected by a coordinate measuring machine_i，y_i) N, wherein N is the total number of the picked pine pollen ears entering the roller, and the total mass center is (x) calculated through the mass center₀，y₀) Thereby obtaining the actual rolling radius d after correction

Wherein M is_iN, i is the mass of the i-th picking of the pine spikes measured by the mass sensor, 1, 2, 3.. N; n represents the total number of the picked pine pollen ears entering the roller; (x)₀，y₀) Calculating the position of the total mass center of the pine flower spike in the pollen removing device; (x)_i，y_i) The coordinates of the ith pine pollen ear in the plane in the roller device are obtained through the detection of a coordinate measuring instrument; n ═ 1, 2, 3.; n is expressed as pluckedThe total number of the loose panicles entering the roller; d is the actual rolling radius of the pine pollen after deviation correction; r is the inner radius of the drum.

S4, calculating the dryness coefficient of the pine pollen ears:

the moisture content in the loose flower spikes is detected through a penetration type microwave moisture meter, the phenomenon that the flower spikes are easy to adhere in the rolling pollen removal process is avoided, the dryness of the loose flower spikes can cause certain influence on the pollen removal effect of the loose flower spikes, and the moisture content measured through the penetration type microwave moisture meter is recorded as k₁Then with an empirical standard k₀Comparing to obtain the dryness coefficient of the pine cones, and recording as k

S5, determining the optimal rotation speed of pollen pini pollen removal based on deviation correction:

according to the parameter data calculated and determined in S1-S4, the optimal rotation speed of pollen pini pollen removal after deviation correction is

Wherein n is the optimal rotation speed of pollen pini pollen removal after deviation correction; f is the optimal centripetal force for completely removing pollen from the pine pollen; m_iThe mass of the i-th picking of the pine pollen spikes, i is 1, 2, 3.. N, and N is the total number of the picked pine pollen spikes entering the roller; d is the actual rolling radius of the pine pollen after deviation correction; alpha is the ratio of the actual rotating speed to the theoretical rotating speed; k is the pollen pini dryness coefficient.

Claims

Translated fromChinese

1.一种基于偏差修正的松花穗脱粉最佳转速确定方法，其特征在于：包括如下步骤：1. a method for determining the optimum rotational speed of pine flower ear de-powdering based on deviation correction, is characterized in that: comprise the steps:

S1、松花穗脱粉主要工作参数的确定；S1. Determination of main working parameters of pine flower ear powder removal;

S2、基于Weibull分布实际工况下脱粉转速完成值的确定；S2. Determination of the completion value of powder removal speed under actual working conditions based on Weibull distribution;

S3、脱粉装置内松花穗总质心偏差程度的计算；S3. Calculation of the deviation degree of the total centroid deviation of the pine flower ear in the de-powdering device;

S4、松花穗干燥度系数的计算；S4. Calculation of dryness coefficient of pine flower ear;

S5、基于偏差修正后的松花穗脱粉最佳转速确定。S5. Determine the optimum rotational speed for de-powdering of the pine flower ear after the deviation correction.

2.根据权利要求1所述的基于偏差修正的松花穗脱粉最佳转速确定方法，其特征在于：所述S1的松花穗脱粉主要工作参数的确定方法如下：松花穗每次采摘完通过质量传感器测量并记录其每次采摘的质量，记为M_i，对每次采摘的顺序i进行标号，i＝1，2，3...N，N表示为采摘的松花穗进入到滚筒中的总数；通过实际测量得到滚筒内部半径为R；在Adams多体动力学仿真软件中导入脱粉滚筒整机的三维模型，通过设置边界约束和滚动载荷，然后利用后处理模块进行仿真计算，最后提取出花穗出粉的最佳向心力为F。2. the method for determining the optimum rotational speed of the pine flower ear dedusting based on deviation correction according to claim 1, it is characterized in that: the determination method of the main working parameter of the pine flower ear de-powdering of the described S1 is as follows: the pine flower ear is picked every time and passes through The quality sensor measures and records the quality of each picking, marked as M_i , and marks the order i of each picking, i=1, 2, 3...N, N means that the picked pine flower ears enter the drum The internal radius of the drum is R by actual measurement; the three-dimensional model of the complete de-powdering drum is imported into the Adams multi-body dynamics simulation software, and the boundary constraints and rolling loads are set, and then the post-processing module is used for simulation calculation, and finally The best centripetal force to extract the flower ear and powder is F.

3.根据权利要求2所述的基于偏差修正的松花穗脱粉最佳转速确定方法，其特征在于：所述S2的基于Weibull分布实际工况下脱粉转速完成值的确定方法：3. the method for determining the optimum rotational speed for de-powdering of pine flower ear based on deviation correction according to claim 2, is characterized in that: the method for determining the completion value of de-powdering rotational speed based on the Weibull distribution actual working condition of described S2:

根据S1的已知参数进行松花粉脱粉工作，在工作前通过转速传感器多次测量松花穗脱粉装置的实际转动速度，对每次的测量进行编号为j，第j次测量得到的实际转速记为n_j，j＝1，2，3...Q，Q表示为测量其实际转动速度的总类型数，然后与理论上预设的转动速度n₀＝200r/min进行比较；经过多次测量后，将所有检测到的脱粉装置的实际转动速度导入到MATLAB软件中自动拟合其Weibull分布曲线，According to the known parameters of S1, the pine pollen dedusting work is carried out. Before the work, the actual rotation speed of the pine flower ear dedusting device is measured for many times through the rotational speed sensor, and each measurement is numbered j, and the actual rotation speed obtained by the jth measurement is The shorthand is n_j , j=1, 2, 3...Q, and Q is expressed as the total number of types whose actual rotational speed is measured, and then compared with the theoretically preset rotational speed n₀ =200r/min; After the first measurement, the actual rotation speed of all the detected powder removal devices was imported into the MATLAB software to automatically fit its Weibull distribution curve.

其中，f(j，β₁，β₂，β₃)为实测脱粉转速的Weibull分布方程，j为测量脱粉装置转速的编号，j＝1，2，3...Q，Q表示为测量其实际转动速度的总类型数；β₁为Weibull分布方程的尺寸参数；β₂为Weibull分布方程的形状参数；β₃为Weibull分布方程的位置参数；e为指数量，Among them, f(j, β₁ , β₂ , β₃ ) is the Weibull distribution equation of the actual measured powder removal speed, j is the number of the measured speed of the powder removal device, j=1, 2, 3...Q, Q is expressed as β₁ is the size parameter of the Weibull distribution equation; β₂ is the shape parameter of the Weibull distribution equation; β₃ is the position parameter of the Weibull distribution equation; e is the number of exponents,

然后通过将f(j)代入下式计算出其中间密度对应值Z_j，从而得到其对应的脱粉实际转速为

然后通过计算得到其与理论预设的转动速度n₀＝200r/min的比值记为αThen, by substituting f(j) into the following formula, the corresponding value Z_j of the intermediate density is calculated, and the corresponding actual speed of powder removal is obtained as:

Then, the ratio between it and the theoretical preset rotational speed n₀ =200 r/min is obtained by calculation and denoted as α

其中，f(j)为实测脱粉转速的Weibull分布方程；j为测量脱粉装置转速的类型编号，j＝1，2，3...Q；Q为测量其实际转动速度的总类型数；C[f(x)]_max为所有实际转动速度的方程模式下的概率最大值；C[f(x)]_min为所有实际转动速度的方程模式下的概率最小值；Z_j为实际转动速度的中间密度值；

为实际转动速度的中间密度值Z_j对应的脱粉转速：n₀为理论上预设的转动速度，取200r/min；α为实际转动速度和理论转动速度的比值。Among them, f(j) is the Weibull distribution equation of the actual measured rotation speed of the powder removal device; j is the type number of the rotation speed of the powder removal device measured, j=1, 2, 3...Q; Q is the total number of types that measure the actual rotation speed of the powder removal device. ; C[f(x)]_max is the probability maximum value under the equation mode of all actual rotation speeds; C[f(x)]_min is the probability minimum value under the equation mode of all actual rotation speeds; Z_j is the actual rotation speed the median density value of the velocity;

is the powder removal speed corresponding to the intermediate density value Z_j of the actual rotational speed: n₀ is the theoretically preset rotational speed, taking 200r/min; α is the ratio of the actual rotational speed to the theoretical rotational speed.

4.根据权利要求3所述的基于偏差修正的松花穗脱粉最佳转速确定方法，其特征在于：所述S3的脱粉装置内松花穗总质心偏差程度的计算方法：4. the method for determining the optimum rotational speed of pine flower ear dedusting based on deviation correction according to claim 3, it is characterized in that: the calculation method of the total centroid deviation degree of pine flower ear in the dedusting device of described S3:

根据S1的测量得到脱粉装置中滚筒的半径为R，但是由于实际进入到滚动内部的松花穗形状不一、质量分布不均等原因，通过坐标测量仪检测得到在滚筒装置内第i个松花穗在平面内的坐标为(x_i，y_i)，i＝1，2，3...N，N表示为采摘的松花穗进入到滚筒中的总数，通过质心计算得到其总质心为(x₀，y₀)，从而得到了修正之后的实际滚动半径dAccording to the measurement of S1, the radius of the drum in the de-powdering device is R, but due to the different shapes and uneven mass distribution of the pine flower ears that actually enter the rolling interior, the i-th pine flower ear in the drum device is detected by the coordinate measuring instrument. The coordinates in the plane are (x_i , y_i ), i=1, 2, 3...N, N represents the total number of pine flower ears picked into the drum, and the total centroid is calculated by the centroid as (x₀ , y₀ ), thus obtaining the corrected actual rolling radius d

其中，M_i为通过质量传感器测量到的第i次采摘松花穗的质量，i＝1，2，3...N；N表示为采摘的松花穗进入到滚筒中的总数；(x₀，y₀)为计算得出的脱粉装置内部松花穗的总质心的位置；(x_i，y_i)为通过坐标测量仪检测得到在滚筒装置内第i个松花穗在平面内的坐标；i＝1，2，3...N；N表示为采摘的松花穗进入到滚筒中的总数；d为进行偏差修正之后的松花穗实际滚动半径；R为滚筒的内部半径。Among them, M_i is the quality of the i-th picking pine flower ears measured by the quality sensor, i=1, 2, 3...N; N represents the total number of picked pine flower ears entering the drum; (x₀ , y₀ ) is the calculated position of the total centroid of the pine flower ear inside the powder removing device; (x_i , y_i ) is the coordinate in the plane of the i-th pine flower ear in the drum device detected by the coordinate measuring instrument; i =1, 2, 3...N; N represents the total number of pine flower ears picked into the drum; d is the actual rolling radius of the pine flower ears after deviation correction; R is the inner radius of the drum.

5.根据权利要求4所述的基于偏差修正的松花穗脱粉最佳转速确定方法，其特征在于：所述S4的松花穗干燥度系数的计算方法：5. the method for determining the optimum rotational speed of pine flower ear de-powdering based on deviation correction according to claim 4, is characterized in that: the calculation method of the pine flower ear dryness coefficient of described S4:

该步骤通过穿透式微波水分仪检测松花穗内部的水分含量，避免在滚动脱粉过程中易造成花穗粘连现象，其干燥度大小会对其脱粉效果造成一定的影响，通过穿透式微波水分仪测量的水分含量，记为k₁，然后与经验标准值k₀做比较，得到了松花穗的干燥度系数，记为k，In this step, the moisture content inside the pine flower ear is detected by a penetrating microwave moisture meter, so as to avoid the sticking phenomenon of the flower ear during the rolling and powder removal process. The moisture content measured by the microwave moisture meter is denoted as k₁ , and then compared with the empirical standard value k₀ to obtain the dryness coefficient of the pine flower ear, denoted as k,

其中，k为松花穗干燥度系数；k₁为测量出的松花穗水分含量，k₀为松花穗的一般标准含水量。Among them, k is the dryness coefficient of the pine flower ear; k₁ is the measured water content of the pine flower ear, and k₀ is the general standard water content of the pine flower ear.

6.根据权利要求5所述的基于偏差修正的松花穗脱粉最佳转速确定方法，其特征在于：所述S5的基于偏差修正后的松花穗脱粉最佳转速确定方法：6. the method for determining the optimum rotational speed of pine flower ear dedusting based on deviation correction according to claim 5, is characterized in that: the method for determining the optimum rotational speed of pine flower ear dedusting based on the deviation correction of described S5:

根据S1-S4中计算和确定的参数数据，基于偏差修正后的松花穗脱粉最佳转速是According to the parameter data calculated and determined in S1-S4, the optimal rotation speed of the pine flower ear after the deviation correction is

其中，n为基于偏差修正后的松花穗脱粉最佳转速；F为松花穗完全脱粉的最佳向心力；M_i为第i次采摘松花穗的质量，i＝1，2，3...N，N表示为采摘的松花穗进入到滚筒中的总数；d为进行偏差修正之后的松花穗实际滚动半径；α为实际转动速度和理论转动速度的比值；k为松花穗干燥度系数。Among them, n is the optimal rotation speed of the pine flower ear after the deviation correction; F is the best centripetal force for the pine flower ear to completely remove the powder; Mi is the quality of the pine flower ear picked at the_i -th time, i=1, 2, 3.. .N, N represents the total number of pine flower ears picked into the drum; d is the actual rolling radius of the pine flower ears after the deviation correction; α is the ratio of the actual rotation speed to the theoretical rotation speed; k is the dryness coefficient of the pine flower ears.