The present invention relates to agricultural balers, and, more particularly, to methods of operating agricultural balers.
Agricultural harvesting machines, such as balers, are used to consolidate and package crop material so as to facilitate the storage and handling of the crop material for later use. In the case of hay, a mower-conditioner is typically used to cut and condition the crop material for windrow drying in the sun. In the case of straw, an agricultural combine discharges non-grain crop material from the rear of the combine defining the straw (such as wheat or oat straw) which is to be picked up by the baler. The cut crop material is typically raked and dried, and a baler, such as a large square baler or round baler, straddles the windrows and travels along the windrows to pick up the crop material and form it into bales.
On a large square baler, a pickup unit at the front of the baler gathers the cut and windrowed crop material from the ground. The pickup unit includes a pickup roll, and optionally may include other components such as side shields, stub augers, wind guard, etc.
A packer unit is used to move the crop material from the pickup unit to a duct or pre-compression chamber. The packer unit forms a wad of crop within the pre-compression chamber which is then transferred to a main bale chamber. (for purposes of discussion, the charge of crop material within the pre-compression chamber will be termed a “wad”, and the charge of crop material after being compressed within the main bale chamber will be termed a “flake”). Typically such a packer unit includes packer tines or forks to move the crop material from the pickup unit into the pre-compression chamber. Instead of a packer unit it is also known to use a rotor cutter unit which chops the crop material into smaller pieces.
A stuffer unit transfers the wad of crop material in charges from the pre-compression chamber to the main bale chamber. Typically such a stuffer unit includes stuffer forks which are used to move the wad of crop material from the pre-compression chamber to the main bale chamber, in sequence with the reciprocating action of a plunger within the main bale chamber.
In the main bale chamber, the plunger compresses the wad of crop material into flakes to form a bale and, at the same time, gradually advances the bale toward the outlet of the bale chamber. When enough flakes have been added and the bale reaches a full (or other predetermined) size, a number of knotters are actuated which wrap and tie twine, cord or the like around the bale while it is still in the main bale chamber. The twine is cut and the formed baled is ejected out the back of the baler as a new bale is formed.
Compressing crop residue materials using a large square baler, as described above, makes them easier to handle and reduces labor and transport costs. Since a high density for the material is needed at high capacities, extrusion is preferred over a fixed wall compression. However, due to frictional heating of the main bale chamber this process requires more power resulting in a higher fuel consumption. The power consumption is thus related to the frictional resistance in the bale chamber. Lowering this friction will lower the power consumption but can also cause the final compressed density to decrease. If the final compression density is too high, then excess power was used to generate the bale. There is thus a trade-off between power consumption and bale density during the compression process within the bale chamber
Compression modeling algorithms are known which describe the compression profile of the crop material. However, such compression modeling algorithms are only used on lab type experiments with simulated bale chambers, or they use oversimplified descriptions of the material behavior. The simulated bale chambers may be much smaller than an actual bale chamber, and the crop material is typically cut into small pieces. Such compression modeling algorithms are useful to at least get an idea as to what is actually happening to the crop material within a bale chamber, but lack accuracy in describing what is really happening to the crop material on a dynamic basis within a bale chamber of a baler during operation.
What is needed in the art is a compression modeling method which can be used on-the-fly during baling, such that the baler operation can be optimized.
The present invention provides a method of operating a baler in which input parameters related to the compression profile of crop material in the bale chamber are used in a compression modeling algorithm to adjust physical actuators, on the output side, to thereby optimize operation of the baler.
The invention in one form is directed to a method of operating an agricultural baler, including the steps of: measuring at least one physical input parameter associated with a compression profile of crop material in a bale chamber of the baler; analyzing the compression profile of the crop material, dependent upon the measured at least one physical input parameter; and adjusting at least one physical output parameter of the baler which affects the compression profile of the crop material in the bale chamber.
The invention in yet another form is directed to an agricultural baler, including at least one sensor for measuring one or more physical input parameters associated with a compression profile of crop material in a bale chamber of the baler. Each sensor provides an output signal corresponding to a sensed physical input parameter. An electrical processing circuit is coupled with each sensor and is configured for analyzing the compression profile of the crop material, dependent upon the measured one or more physical input parameters. At least one actuator, under control of the electrical processing circuit, is configured for adjusting one or more physical output parameters of the baler which affect the compression profile of the crop material in the bale chamber.
An advantage of the present invention is that the compression profile of the crop material within the baler can be measured “on-the-fly”, and adjusted to optimize the operation of the baler.
Another advantage is that the compression profile can be measured using a non-linear analysis, which is a more accurate way of determining the compression profile of the crop material.
Yet another advantage is that the measured input parameters relating to the compression profile can be either directly or indirectly related to physical attributes associated with the main bale chamber.
A further advantage is that input parameters such as temperature of the crop material and friction against the crop material can be accounted for with the analysis techniques.
A still further advantage is that the relaxation of the crop material between strokes of the plunger can be lessened through the use of holding structures, such as clamping side walls and/or hay dogs.
Another advantage is that it is possible to determine on an incremental basis which part of the compressed crop material actually undergoes compression with each new stroke of the plunger.
Yet another advantage is that multiple output parameters can be adjusted to affect the compression profile of the crop material, such as the position of the side walls, length of the bale chamber and/or cross sectional shape of the bale chamber.
A further advantage is that the measuring, analysis and adjustment occurs in a closed loop manner, allowing continual optimization of the baler.
The above-mentioned and other features and advantages of this invention, and the manner of attaining them, will become more apparent and the invention will be better understood by reference to the following description of embodiments of the invention taken in conjunction with the accompanying drawings, wherein:
Corresponding reference characters indicate corresponding parts throughout the several views. The exemplifications set out herein illustrate embodiments of the invention, and such exemplifications are not to be construed as limiting the scope of the invention in any manner.
Referring now to the drawings, and more particularly to
Baler 10 operates on a two stage feeding system. Crop material is lifted from windrows into the baler 10 using a pickup unit 12. The pickup unit 12 includes a rotating pickup roll 14 with tines 16 which move the crop rearward toward a packer unit 18. An optional pair of stub augers (one of which is shown, but not numbered) are positioned above the pickup roll 14 to move the crop material laterally inward. The packer unit 18 includes packer tines 20 which push the crop into a pre-compression chamber 22 to form a wad of crop material. The packer tines 20 intertwine the crop together and pack the crop within the pre-compression chamber 22. Pre-compression chamber 22 and packer tines 20 function as the first stage for crop compression. Once the pressure in the pre-compression chamber 22 reaches a predetermined sensed value, a stuffer unit 24 moves the wad of crop from the pre-compression chamber 22 to a main bale chamber 26. The stuffer unit 24 includes stuffer forks 28 which thrust the wad of crop directly in front of a plunger 30, which reciprocates within the main bale chamber 26 and compresses the wad of crop into a flake. Stuffer forks 28 return to their original stationary state after the wad of material has been moved into the main bale chamber 26. Plunger 30 compresses the wads of crop into flakes to form a bale and, at the same time, gradually advances the bale toward outlet 32 of main bale chamber 26. Main bale chamber 26 and plunger 30 function as the second stage for crop compression. When enough flakes have been added and the bale reaches a full (or other predetermined) size, knotters 34 are actuated which wrap and tie twine around the bale while it is still in the main bale chamber 26. Needles 36 bring the lower twine up to the knotters 34 and the tying process then takes place. The twine is cut and the formed bale is ejected from a discharge chute 38 as a new bale is formed.
Plunger 30 is connected via a crank arm 40 with a gear box 42. Gear box 42 is driven by a flywheel 44, which in turn is connected via a drive shaft 46 with the PTO coupler 48. The PTO coupler is detachably connected with the PTO spline at the rear of the traction unit, such as a tractor (not shown).
According to an aspect of the present invention, baler 10 is operated in a closed-loop manner to optimize the compression profile of the crop material within the main bale chamber 26. In general, input parameters associated with the compression profile of the crop material in the main bale chamber 26 are sensed and used in a mathematical analysis method to determine the compression profile of the crop material for each compression stroke of the plunger 30. One or more crop models are used in the mathematical analysis method, utilizing the sensor input parameters, to determine the compression profile of the crop material. The compression profile of the crop material is then used to determine if one or more output parameters of the baler 10 should be adjusted to optimize the compression profile. This sensing, analysis and adjustment occurs in a closed-loop manner to optimize the operation of the baler 10. The mathematical analysis method is used for optimization of the baler 10 with respect to energy consumption of the baler and the compression profile of the crop material.
More particularly, and referring to
The sensor(s) used to measure the physical input parameters may measure, for example, but not limiting:
a force on the plunger 30 which reciprocates within the bale chamber 26;
a position of the plunger 30 within the main bale chamber 26;
a position of at least one structural component defining the bale chamber 26; or
a deformation of the crop material along a length of the bale chamber 26.
The force on the plunger 30 can be measured using a strain gauge (not shown). As another alternative, the speed of the crank arm can be used to determine the force on the plunger. This type of plunger force measuring system is further described in a copending patent application BE2013/0327 entitled “A method and system for determining the plunger load of a baler”, which is incorporated herein by reference.
The position of the plunger may be determined by measuring a rotation of a flywheel 44 which is connected with the plunger 30 via the crank arm 40, and a dead point associated with a change in direction of the plunger 30 (
The one or more structural components of the bale chamber which may be measured include a pair of side doors 64 and a top door 66 (
The deformation of the crop material along the length of the bale chamber 26 may be measured using one or more length wheels 70 located within the bale chamber 26 (
The one or more actuators under the control of controller 58 affect the physical output parameters having an influence on the compression profile of the crop material within the bale chamber 26. These physical output parameters are preferably directly associated with the bale chamber 26, and may include:
a position of one or more side walls associated with the bale chamber;
a fixed length of the bale chamber; and
a cross-sectional shape of the bale chamber.
Referring to
Compression Analysis Example
Current compression profile analyses do not allow a proper simulation of the material behavior in the main bale chamber of a large square baler. To correct for these analysis shortfalls, it is necessary to:
1. use non-linear crop model with interpretable parameters, for a proper description of the material;
2. take the bend in the one or more doors into account for proper simulation of the compression behavior; and
3. calculate the deformation profile along the chamber length which characterizes the pulsation zone and the clamping of the material.
The decay of the compression force along a compression box with constant cross section has been calculated in literature. In deriving the force balances, the friction force is assumed to be the same for all doors. However, it is known that the friction forces will be different at different doors in the baling chamber. As a result, this model for the stress-decay in the compression chamber is not applicable to a compression chamber for extrusion. Also the resulting deformation along the length of the compression box is not calculated when using this equation. After considering the inclination of 3 walls, being the top wall and the two side walls, equations 1-3 are derived for the forces in respectively the x, y and z directions. These equations consider the effect of the inclined side walls on the friction force. Then, a crop model (linear or non-linear) was used in calculating the stress decay with the depth in the compression box as a function of the plunger position and the inclination of the side walls. It is then possible to predict the pressure and the deformations when the crop parameters and the plunger position are known.
The stresses in the compression chamber are modelled as described above. To relate the stresses with the resulting deformations, a crop model is necessary. First the material is described with the 3D, linear elastic material model. Then, a non-linear material model is extended to 3 dimensions for using it in the material description. This last relation is further called the Faborode model, described below.
The Pressure and the Deformation Along the Depth of the Chamber
Material compression in the large square baler is divided into different phases. In the first phase, the newly added material and the expanded material from previous cycles, is compressed against the clamped material in the chamber. In the next phase, the bale starts to move.
This study focuses on the compression of the crop while assuming a stationary condition, i.e., before the bale starts moving. Because the compression continues from the front till the back of the chamber, the material is discretized in slices, as shown in
The friction forces in
The balances in y and z (Equations 2-3) are filled into the balance in x (Equation 1) to give the following relation in x:
Now, the stresses along the length of the compression chamber are derived and these stresses have to be linked with the deformation of the material. For comparison, two relations between stress and deformation are proposed. The first assumes linear elastic and isotropic material behavior and is given by Hooke's general law. In the second relation, the stress is assumed to rise exponentially with deformation.
Linear Elastic Material Behaviour
The balance in x, as shown in Eq. (4), gives the relation of forces on every slice in the chamber. Solving this relation with respect to the deformation, implies knowing the inclinations of the side walls (θl, θt) and knowing the constitutive relations. For linear elastic material behaviour, the constitutive equations are (Hooke's general law):
With E the Young's modulus and ν the Poisson coefficient. The engineering strains (or Cauchy strains) (ϵy,ϵz) are written as a function of the height and width of the chamber as:
ϵy=ln(hx)−ln(h0)
ϵz=ln(bx)−ln(b0)
The assumption of small inclinations of side walls and top door gives the height and width of the chamber as: hx=h0−xθt and bx=b0−2θl. The stress in the compression direction x (in Eq. (5)) is then written as:
Where
A=bx(2μ−μ2θt+θt)
B=2(θl+μ)hx
C=ν(θtbx+2θlhx)
This equation is solved by assuming the movement at the end of the fit to be 0, i.e. before the bale starts moving.
Non-Linear Material Behaviour
The Faborode crop model (Faborode, M. & O'Callaghan, J.; “Theoretical analysis of the compression of fibrous agricultural materials”, Journal of Agricultural Engineering Research, Vol. 35, pp. 175-191 (1986)) relates a uniaxial stress σ to the resulting density ρ as follows:
In this relation, bc and K0 are crop parameters identifying the material behavior and ρ0 is the density before compression.
The three-dimensional non-linear constitutive relations are built up from the Faborode crop model given by equation (7). This relation is expanded into 3 dimensions by assuming the lateral forces to be a function of the compression force only (Ferrero, A.; Horabik, J. & Molenda, M., “Density-pressure relationships in compaction of straw”, Canadian Agricultural Engineering, Vol. 33 (1), pp. 107-111 (1990)). This results in the following relations:
σn=−k1σx (8)
σnl=−k2σx (9)
Assuming small inclinations of the side walls allows the approximation of hx and bx by the constants
where L is the length of the side walls. With the force balance in x (Eq. (4)) and the non-linear constitutive relations (Eqs. (8) and (9)) the stress-decay with x is calculated as:
With A=hxbx and B=k1bx(2μ−μ2θt+θt)+2k2(θl+μ)hx. In this relation, the stress on the plunger σx0 is given by the Faborode model (Eq. (7)). Since the density is related to the engineering strain, the deformation can then be calculated as:
With u0 the initial length of the bale in the chamber and u the length during compression by the plunger.
So far, the side walls and top door are considered flat plates with an inclination. In practice, these side walls may have one or more bends, or even a curved shape, to influence the compression behavior of the compressed material. In the following section, this bend will be implemented.
Deformation by Side Walls and Top Door
The shape of the side walls and top door determines the force and density profile in the chamber and should therefore be considered in the equations. So far, the angle in the side walls was assumed constant with x but the re-definition of the inclination will allow to re-use the derived equations. Indeed, the bend in the side walls will be defined as an extra deformation on top of the inclination of that wall.
The angles in the moment equations should then be altered into the expanded angles:
θt+=θt+α1−Hx-l
θl+=θl+α1−Hx-l
θr+=θr+α1−Hx-l
And the height hx and bx are written as:
hx=ho−x(tan θt+tan α1)+Hx-l
bx=bo−x(tan θl+tan θr+2 tan α1)+2Hx-l
In these relations, Hx-l
The necessary equations are derived above, the measurements for validating them are described in the following section.
Measurements in the Baler
In every plunger cycle, the compression force characterizes the behavior of the crop in the chamber. This compression profile was therefore fitted with the Faborode crop model and the crop parameters allow the calculation of the compression in that cycle. The Young modulus was determined by the slope of the stress-strain relation.
The deformation along the chamber length is measured with wheels penetrating the crop in the chamber. When the plunger compresses the material, the measured rotation relates to the deformation of the material. In total, 12 length wheels are installed along the chamber length. It is obvious that more or less length wheels may be installed, depending on the total length of the bale chamber.
During compression by the plunger, more and more material will be deformed. At a certain plunger force, the continued movement of the plunger will move the bale instead of deforming more material. Hence, the measured deformation profile will decrease along the chamber length and will reach zero at a certain position in the chamber. For a proper fit of the compression, only the length wheels that measure a (positive) deformation are considered and the remaining, part of the bale is left out. Because the fitted profile considers several plunger cycles, the standard deviation on this profile determines the position in the chamber at which the measured deformation becomes zero. Therefore, the wheels where the medians minus 1 standard deviation are strictly positive, are used for fitting the compression models. The point of zero deformations indicates the length of the compression zone in the compression chamber. After this compression zone, there is a clamping zone which determines the clamping force in the chamber. In other words, the crop in this zone forms a wall against which the plunger can compress newly added and expanded crop. In this zone, the density does not change significantly.
At the exit of the compression chamber, the material expands outwards. This expansion can be modelled by assuming the force balances as given in Equations (1-3) but then the material model should be altered for modelling expansion. An example model for the expansion originates from the Faborode model and is given below:
The compression and expansion zone both determine the final bale density. The hold zone determines the clamping force and thus the plunger force at which the deformation of more material stops and the bale starts to move.
Fitting the Compression Profile
The measured compression profiles were then fitted with Eqs. (6) and (11). The fitting algorithm minimizes the sum of the squared errors between the model and the measured data by means of the interior-point method.
The linear elastic material model (Eq. (6)) requires the Young modulus E, the Poisson coefficient ν, the friction coefficient μ and an initial condition for solving the differential equation ϵ0. Analogue to this, the non-linear material description requires the porosity index bc, the force ratio the friction coefficient μ and the initial condition u0/u(x=0). The porosity index be and the Young modulus E are derived from the fit of the force profile. Hence, only the Poisson coefficient, force ratio and friction coefficient needed to be estimated by the model. Due to the sensitivity of the fit to the initial deformation and the inaccuracy of measuring it, the initial conditions were also determined by the fit, resulting in 3 fitting parameters per model.
The cost function of the models shows a large dependency between on one hand the Poisson coefficient ν and the friction coefficient μ, and on the other hand the force ratio k1 and the friction coefficient μ. This dependency causes the optimization to find non-feasible parameter sets. The cost functions for the linear and non-linear material models are shown in
Due to the parameter dependency, multiple initial guesses were used in the optimization. A proper spreading of the guesses inside the feasible parameter domain was obtained with a latin hypercube design. The resulting parameter set with the smallest evaluation of the cost function is then used as a solution of the optimization problem. The resulting fits of the linear and non-linear models are respectively shown in
From the foregoing, it is apparent that a linear approximation of the constitutive material behavior is only valid in a narrow pressure range and expanding this range requires adjusting the parameters during compression. Alternatively, the material can be modelled with a non-linear model. The more proper material description then allows for a simple approximation of the lateral forces, as shown with the force ratio k1. However, this parameter shows a large dependency with the friction coefficient. Indeed, both k1 and μ are ratio's of forces relating two perpendicular stresses. The relation between these stresses is then given by the static momentum equations. As a consequence, many solutions with limited change in fitting performance exist.
The full design of the side walls is implemented and the effect of the bend can thus be investigated. The current models assume small inclinations which implies only a small bend can be implemented. Since a large bend would severely damage the bale shape and would make the compression too aggressive, the assumption of small angles can still hold.
Conclusions of Fitting the Compression Profile
The calculation of the deformation profile implies calculating the zone of compression in the compression chamber. The larger this zone, the more material will be deformed in the compression cycle, which means a limited clamping. The current detection of the point of zero compression is based on the measurements of the length wheels. Therefore, this detection is not accurate enough for proper model fitting.
Both linear and non-linear material models fit the deformation profile but only the non-linear model has a proper material description, resulting in feasible parameters. However, due to the dependency between k1 and μ in the non-linear model, multiple solutions exist. The current optimization is influenced from measurement error caused by the parameter estimation (of parameters E and bc) and the calculation of the point of zero deformation.
Up till now, only the compression is considered. However, when crop is compressed until a certain density and it is held for some time, the required compression force will decrease. The decrease of the compression force with time is called relaxation. Since a bale can stay up to 60 sec. in the compression chamber, this effect should not be neglected.
Relaxation Analysis Example
A simple description of the relaxation behavior is given by the Maxwell model (Equation 17). In this model, the relaxation time η*E−1 describes the time until the stress decays to e−1 (or approximately 36.8%) of its initial value (Sitkei, G., “Mechanics of Agricultural Materials”, Elsevier, Amsterdam (1986)).
σ=−Rη(e(−Et/η)−1) (17)
A second relaxation model is given by Peleg (1980) as:
F0t=k1+k2t
F0−F(t) (18)
where:
F0 is the force at the beginning of the relaxation,
t is the time, and
k1, and k2 are empirical parameters.
It is to be noted that k1−1 is the initial decay rate of the force F with time and k2−1 is the decay rate of the relative force (F0−F(t))*F0−1 close to the end of relaxation. The model of Peleg is used to define the asymptotic modulus Ea as (Peleg, 1980):
Ea=F0(1−1/k2)
A (19)
where A is the cross sectional area of the sample and is the applied relative deformation. The asymptotic modulus is a measure for the power of the material to resist the applied stress (Peleg, 1980). Both the relaxation time and the asymptotic modulus characterize the resistance of the crop to the applied deformation.
Relaxation of the crop material can be lessened using structural features associated with the main bale chamber 26 which tend to hold and thereby inhibit relaxation of the crop material between strokes of the plunger 26. For example, the side doors 64 and/or top door 66 can be used to clamp and hold the crop material between strokes of plunger 26. As another example, so called “hay dogs” 76 can be used in a wall, e.g., top door, of the bale chamber 26 which hold the crop material between strokes of the plunger 26 (
While this invention has been described with respect to at least one embodiment, the present invention can be further modified within the spirit and scope of this disclosure. This application is therefore intended to cover any variations, uses, or adaptations of the invention using its general principles. Further, this application is intended to cover such departures from the present disclosure as come within known or customary practice in the art to which this invention pertains and which fall within the limits of the appended claims.
Number | Date | Country | Kind |
---|---|---|---|
2013/0708 | Oct 2013 | BE | national |
This application is the US National Stage filing of International Application Serial No. PCT/EP2014/065355 filed on Jul. 17, 2014 which claims priority to U.S. Provisional Application 61/860,531 filed Jul. 31, 2013, and Belgium Application BE2013/0708 filed Oct. 21, 2013, each of which is incorporated herein by reference in its entirety.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/EP2014/065355 | 7/17/2014 | WO | 00 |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO2015/014625 | 2/5/2015 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
4624180 | Strosser | Nov 1986 | A |
8794135 | Simon | Aug 2014 | B1 |
9642311 | Posselius | May 2017 | B2 |
20080141870 | Verhaeghe | Jun 2008 | A1 |
20110023732 | Herron | Feb 2011 | A1 |
20120067229 | Lanning | Mar 2012 | A1 |
20120179338 | Dresher | Jul 2012 | A1 |
20120186467 | Lanning | Jul 2012 | A1 |
20120192732 | Lanning | Aug 2012 | A1 |
20120326417 | Lanning | Dec 2012 | A1 |
20130000496 | Posselius | Jan 2013 | A1 |
20130000497 | Posselius | Jan 2013 | A1 |
20140202343 | van Amstel et al. | Jul 2014 | A1 |
20170013772 | Kirk | Jan 2017 | A1 |
20170013782 | Kindt | Jan 2017 | A1 |
Number | Date | Country |
---|---|---|
1020897 | Aug 2014 | BE |
2754138 | Apr 1998 | FR |
Entry |
---|
Faborode M O ; O'Callaghan Jr—Theoretical analysis of the compression of fibrous agricultural materials, Journal of Agricultural Engineering Research, 19861101 Academic Press—ISSN 0021-8634, vol. 35, Nr:3, pp. 175-191. |
Peleg M; Peleg M; Normand M D—Comparison of two methods for stress relaxation data presentation of solid foods, Rheologica Acta, 19830101 Springer Berlin Heidelberg—ISSN 0035-4511, vol. 22, Nr:1, pp. 108-113. |
Number | Date | Country | |
---|---|---|---|
20160165802 A1 | Jun 2016 | US |
Number | Date | Country | |
---|---|---|---|
61860531 | Jul 2013 | US |