In this section, we explain the details of the proposed model, consisting of notations, assumptions, an objective function and constraints, and thus, the mathematical formulation of the model is presented based on the nature of the inventory system behavior for growing items.
3.1. Notations and Assumptions
The notation of the proposed mathematical model is expressed as follows:
Index:
: The index of the time.
Parameters:
: The weight of a unit item at time (weight unit);
: The maximum allowable length of the breeding period (day);
: The minimum allowable length of the breeding period (day);
: The constant demand rate per weight unit (weight unit/year);
: The holding cost per weight unit (monetary unit/weight unit/year);
: The purchasing cost per weight unit (monetary unit/year);
: The setup cost per growing cycle (monetary unit/setup);
: The production (feeding) cost per unit item during the growing cycle (monetary unit/unit item);
: The tax of carbon dioxide production (monetary unit/liter/day. weight unit);
: The disposal cost of carcass (monetary unit/unit carcass);
: The polynomial function of the fraction of dead items during the growing cycle (percent);
: The polynomial function of production (feeding) consumption (unit items);
: The polynomial function of carbon dioxide production (liter/day. weight unit);
: The growing rate;
: The asymptotic weight;
: The integration constant of the growing function;
: The shape parameter of the growing function.
Dependent variables:
: The total weight of the inventory (weight unit);
: The breeding period (day);
: The consumption period (year);
: The annual setup cost in a year (monetary unit);
TP: The annual purchasing cost in a year (monetary unit);
: The annual holding cost in a year (monetary unit);
: The annual disposal cost in a year (monetary unit);
: The annual production (feeding) cost in a year (monetary unit);
: The annual carbon dioxide production tax in a year (monetary unit);
: The total cost in a year (monetary unit).
Decision variables:
: The slaughter age (day);
: The total number of growing items ordered at the beginning of a cycle (unit items).
At the beginning of the breeding period (
) in the poultry farm, a lot of newborn birds
with an initial weight
(the weight of newborn birds) are purchased from suppliers and are then raised until they reach the market slaughter weight (
). Once the birds reach the proper size and weight, at the slaughter date (
), they are killed for market consumption annually. During the breeding period, several chickens are dead before reaching the slaughter date, mainly because they have reduced walking ability and due to lameness and reduced access to water and feed, which lead to debilitation and death. The number of dead birds at date
is equal to
, where
is the percentage of the cumulative daily mortality. Therefore, on the slaughter date, the number of live birds
are killed for consumption, and the rest of them,
, are disposed. After the slaughtering process, the consumption period (
) starts with the constant demand rate (
) until the weight inventory level reaches zero. The behavior of the weight inventory system of growing items is illustrated in
Figure 1.
For the growth function, there are several valuable studies that have measured bird growth (Rezaei [
6]); however, the Richards function, proposed by Richards [
33], is one of the most important growth functions for bird weight (see also Goliomytis et al. [
34]). Therefore, this function is used in this study, as follows:
where
is the weight of the body of a live bird at age
.
As mentioned before, many birds are dead during breeding. In the current study, we use the polynomial function of the percentage of cumulative daily mortality, which relates the percentage of daily mortality to the age of birds, based on the work of Xin et al. [
35], as follows:
During the breeding period, birds grow and nurture with a feeding function, as presented by Goliomytis et al. [
34]. This function is a polynomial function, which depends on the birds’ age. This function is fitted based on the collected data and is estimated as follows:
Moreover, for carbon dioxide production function, we choose the commonly polynomial function in the real world, which depends on the age of the birds, and it is fitted and estimated according to Leonard et al. [
36], as follows:
3.2. Objective Function and Constraints
We consider a condition where a poultry farm purchases newborn birds, grows them up to the market slaughter weight, kills them and responds to customer demand. The total cost of the inventory system includes the setup cost, purchasing cost, holding cost, feeding (production) cost, disposal cost and carbon dioxide production tax. Next, each component of the total cost is obtained as follows:
At the beginning of the breeding cycle, some activities and processes, such as cleaning and maintenance, are needed to start the breeding process, and the costs of these activities are imposed as the setup cost (
) on the company for each cycle. Because we should obtain the setup cost for a year, we can divide the setup cost per cycle by
, as follows:
where
is the consumption period, which can be obtained from Equation (6).
Considering
, the annual setup cost is as follows:
- -
Purchasing cost
As
Figure 1 shows, at first, a number of newborn birds (
) with an initial weight (
) are received from the supplier, and subsequently, the purchasing cost per growth cycle becomes equal to
, where
is the purchasing cost per weight unit. Therefore, the annual purchasing cost is computed as follows:
- -
Holding cost
It is easy to derive from
Figure 1 that the average weight of the inventory level during the consumption period is
, and the length of the consumption period is
. The holding cost per cycle is
, where
is the annual inventory cost per weight unit. Therefore, the annual holding cost is calculated as follows:
- -
Disposal cost
Several birds are dead during the breeding process as result of losing their ability to walk, and the percentage of cumulative daily mortality at the slaughter date
is
. Hence, the total number of birds that die during the breeding period is equal to
, and subsequently, the disposal cost per cycle is computed by multiplying
by
, where
is the disposal cost of a dead bird. Finally, the annual disposal cost is obtained as follows:
- -
Feeding (production) cost
The feeding function per weight unit, which depends on the age of the chicken, is stated in Equation (3). Considering
and
, the feeding cost per unit and the length of the breeding period, respectively, the total feeding cost per cycle is
. Thus, the annual feeding cost is calculated as follows:
- -
Carbon Dioxide Production Tax
According to Broucek and Cermák [
37], CO
2 production by animals is relative to their metabolic heat production and consequently to their metabolic body weight, which, in turn, is affected by bird activity and temperature. Carbon dioxide is produced during the breeding process. The carbon dioxide production function, which depends on the birds’ age, is indicated in Equation (4). Considering
and
, the tax cost of carbon dioxide production and the length of the breeding period, respectively, the annual tax of carbon dioxide production is obtained as follows:
- -
Total cost
The annual total cost is formulated as follows:
Substituting
from Equation (1), the annual total cost is as follows:
- -
Constraints
When the company wants to order a number of newborn birds, the number of ordered items must be an integer number, as the company can buy only live birds. Moreover, the slaughter date must be an integer number between the minimum allowable length of the breeding period () and the maximum allowable length of it (), because, in the real world, each domestic animal has a growth period that is determined by the market, breeding process and its nature. Therefore, this constraint is .
- -
Final Model
According to the objective function in Equation (14) and the constraints stated in the above subsection, the final model of the proposed study is as follows: