PLACE IN RETURN BOX to remove this checkout from your record. TO AVOID FINES mum on or before date due. DATE DUE, DATE DUE DATE DUE # MSU Is An Affirmative Action/Equal Opportunity Institution “nod—ll... COMPUTER SIMULATION OF COMBINE HARVESTING AND HANDLING OF SUGAR CANE IN BARBADOS By Winston O'Neale Harvey A DISSERTATION Submitted to Michigan State University in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY Department of Agricultural Engineering 1982 ABSTRACT COMPUTER SIMULATION OF COMBINE HARVESTING AND HANDLING OF SUGAR CANE IN BARBADOS BY Winston O'Neale Harvey The broad objective of this study was to improve the efficiency of combine harvesting of sugar cane in Barbados. The harvesting process was broken down into two subsystems: a field subsystem and a factory yard subsystem. Two computer simulation models, structured in GASP IV simu- lation language, were developed to model the operations involved in these systems. Mbdel FIELDOP simulated the activities involved in the harvest- ing and loading of cane in the field, and in its transportation to the factory for processing. Model FACYARD simulated the weighing and unload- ing activities performed on cane transport units at the factory. Following validation of the models, four different factory yard configurations and eight different field equipment combinations were simulated. Model parameters varied were the number of factory yard scales and the numberscnffield tractors and cane transport wagons assigned to a harvester. Output from the models included utilization factors for the various component machines, daily cane delivery from the field system, and daily amounts of cane handled by the factory yard system. This output was fed into a cost program which calculated unit harvesting costs and total annual cane delivery for the equipment combinations simulated. Results indicated that a second scale at the factory can reduce the humory residence time of transport units by 88 percent, increase combine harvester utilization efficiency by 50-60 percent, increase daily cane receipts at the factory by more than 30 percent, and eliminate milling lost time due to lack of cane. Harvesting systems using two field tractors, rather than one, were also shown to be capable of consistently delivering 15-25 percent more cane per day. The economic analysis demonstrated that harvesting cost per tonne can be significantly reduced by either adding a second field tractor, increasing the number of cane wagons assigned to a harvester, installing a second weigh scale at the factory, or a combination of these. A sensitivity analysis revealed that, as a single measure, adding a second scale to the factory would be two to three times more effective in reducing costs then would either of the other measures. Approved Major Professor Approved Department Chairman ACKNOWLEDGEMENTS The author wishes to express sincere gratitude to his Major Pro- fessor, Dr. Merle L. Esmay, for his support and assistance throughout this study, and to Dr. Ajit K. Srivastava, Dr. Robert H. Wilkinson and Dr. J. Roy Black for serving on his guidance committee. Special thanks are also expressed to the following: Mr. R.A. Baynes of Barclays Bank Caribbean Head Office, who was instrumental in the negotiationg of a loan to finance the research phase of the study in Barbados. Mr. G.B. Hagelberg of the Ministry of Finance and Planning, Barbados, who negotiated some additional funds to help defray research expenses. Dr. John R. Ogilvie, Chairman, School of Engineering, University of Guelph, Canada, who offered some valuable technical assistance. Mr. Michael Gill of Ashbury Plantation, Mr. Geoffrey Skeete of Edgecumbe Estates Ltd. and Mr. Geoffrey Armstrong of Sunbury, whose harvesting operations were the main ones monitored during the data collection exercise. The Managements of Carrington, Foursquare and Andrews Sugar Factories, who readily accommodated the monitors on their prOperties. Ms. Joan Antrobus of the office at Carrington Factory for making weekly factory statistics available. Mr. Ishmael Roett of Barbados '0' Level Institute for his assist- ance in selecting five students who turned out to be excellent monitors. ii The organizers of the Latin American Scholarship Plan for American Universities, for their financial assistance in the form of a LASPAU Fellowship. The University of the West Indies, for granting the author leave to take up the Fellowship. Ms. Lolly Carrington for her friendship and encouragement while she was also a student at Michigan State University. Finally, the author wishes to express his deepest thanks and appreciation to his wife, Ena, for her unfailing moral support and encouragement and for typing a major section of the first draft of the dissertation. iii 2.4 4.3 TABLE OF CONTENTS INTRODUCTION. . . . . . . . . . . . Background. . . . . . . . . . . . . . . . Objective of the Study. . . . . . . . . Approach to the Study . . . . . . . . . OVERVIEW OF THE BARBADOS SUGAR INDUSTRY . . . Introduction. . . . . . . . . . . . . . The Role of the Sugar Industry in the . Barbadian Economy . . . . . . . . . . . . . . Structure of the Industry . . . . . . . . . . .3 3 2 .1 The Production Sector . . . . . . . . 2. 2 The Processing Sector . . . . . . Markets for Barbados' Sugar . . . . . . 2.4.1 The European Economic Community . . Market. . . . . . . . . . . . . . 2.4.2 The World Market. . . . . . . . . . . 2.4.3 The Domestic Market . . . . . . . . . Review of Mechanical Harvesting of Sugar Cane in Barbados. . . . . . . . . . . . . . . LITERATURE REVIEW . . . . . . . . . . . . . . Introduction. . . . . . . . . . . . . . . . . Symbolic and Mathematical Models. . . . . . . Computer Simulation Models. . . . . . . . . . Sugar Cane Simulation Models. . . . . . . . . Computer Simulation Languages . . . . . . . METHODOLOGY AND MODEL DEVELOPMENT . . . . . . Introduction. . . . . . . . . . . . . . . . . Mbdel Development Approach. . . . . . . . . 4.2.1 An Overview of the GASP IV Simulation Language . . . . . . . . . . . . . . The Field Operations. . . . . . . . . . . 4.3.1 Objectives. . . . . . . . . . . . . . 4.3.2 Activities and Events . . iv Page LDC-9H . 10 . l4 . 16 17 . 19 20 . 21 27 27 28 37 . 42 48 52 52 53 55 . 60 . 60 62 4.4 4.5 CHAPTER 5 5.3 5.4 5.5 CHAPTER 6 6.1 6.2 TABLE OF CONTENTS (cont'd.) The Factory Yard Operations. . . . . . 4.4.1 4.4.2 Obj actives O O O O O O O O O O O O O O 0 Activities and Events. . . . . . . . . . . Data Acquisition and Analysis. 4.5.1 Data Acquisition . . . . . . . . . . . . . 4.5.2. Data Analysis. . . . . . . . . 4.5.3 The Data Analysis Program. . . . . . FIELD AND FACTORY YARD COMPUTER SIMULATION mnu s O O C O 0 O O O O C C O O 0 Introduction . . . . . . . . . Field Operations Simulation Model. . . . . . The System Environment . . . . . . . . . . Entities . . . . . . Q . . . . . . . . Endogenous Variables . . . . . . . . . . Output Variables . . . . . . . . . . . The Field Operations Simulation Program (FIELDOP). 5.3.1 General Flow through Program FIELDOP . 5.3.2 Input Data used for FIELDOP The Factory Yard Simulation Model. . . . . . . . . The System Environment . . . . . . . . . . Entities . . . . . . . . . . . . . . . . . Endogenous Variables . . . . . . . . ._. . Output Variables . . . . . . . . . . . . The Factory Yard Simulation Program (FACYARD). . 5.5.1 General Flow through Program FACYARD . . . 5 5 2 Input Data Used for FACYARD. . . . . . . . RESULTS AND DISCUSSION . . . . . . . . . . . . Introduction . . . . . . . . . . . . . . . . . . . Results from the Factory Yard Simulation Model . 6.2.1 FACYARD Results for Current Yard Configuration . . . . . . . . . . 6.2.2 FACYARD Results for Extended Work Period . . . . . . . . . . . . . Page . 64 . 64 . 66 . 68 . 68 . 70 . 75 . 75 . 76 . 76 . 80 . 80 . 81 . 81 .104 .104 O 104 .106 .106 .110 .110 .115 .121 .121 .121 .124 .125 TABLE OF CONTENTS (cont'd.) 6.2.3 FACYARD Results for Configuration with 1 Weigh Queue and 2 Weigh Scales. . . . . . 6.2.4 FACYARD Results for Configuration with 2 Weigh Queues and 2 Weigh Scales . . . . . . 6.3 Results from the Field Operations Simulation Model. . 6.3.1 Validation of model FIELDOP . . . . . . . . . 6.3.2 Output from FIELDOP for Different Equip- ment Combinations . . . . . . . . . . . . . . 6.4 General Discussion. . . . . . . . . . . . . . . . . . 6.4.1 Implications Associated with an Extended work Per 10d 0 C O O O O O O O O O O O I O O 0 6.4.2 Implications Associated with Two-Scale Configurations. . . . . . . . . . . . . . . . 6.4.3 Proposals for Re-organization of Carring- ton Factory Yard. . . . . . . . . . . . . 6.5 Economic Analysis . . . . . . . . . . . . . . . . . 6.5.1 Calculation of Harvesting Machinery Costs . 6.5.2 Sensitivity of System Cost and Output to Various System Parameters . . . . . . . . . . 6.6 The Projected Requirement of Cane Combine Harvesters. 6.6.1 The Theoretical Minimum Harvester Requirements. . . . . . . . . . . . . . . . 6.6.2 Required Number of Harvesters for Various Simulated Annual Output Levels. . CHAPTER 7 CONCLUSIONS AND RECOMMENDATIONS . . . . . . . . . 7 O 1 COnClusionS O O O C O O C O O O O O O 0 O O O O 0 O I 7.2 Recommendations . . . . . . . . . . . . . . . . . . . 7.3 Suggestions for Future Work . . . . . . . . . . . APPENDICES Appendix I Data Collection Manual for Monitors . . . . . Appendix II Fortran Listing of Data Analysis Program DATANAL O O O O O O C O O O O O O O O O 0 Appendix III Fortran Listing of Factory Yard Operations Program FACYARD . . . . . . . . . . . . . . . vi Page .125 .128 .128 .130 .134 .142 .142 .145 .146 .151 .151 .153 .162 .162 .163 .165 .165 .166 .167 O 168 .174 .196 Appendix IV LIST OF REFERENCES . TABLE OF CONTENTS (cont'd.) Fortran Listing of Field Operations Program FIELDOP. . . . . . . . . . vii LIST OF TABLES TABLE PAGE 3.1 Comparison of some contemporary simulation languages. . . . . . . . . . . . . . . . . . . . . . . 51 4.1 Events monitored for the field operations. . . . . . . 65 4.2 Events monitored for the factory yard activities . . . 67 5.1 Entities and their associated attributes (FIELDOP) . . 82 5.2 Files used in model FIELDOP. . . . . . . . . . . . . . 83 5.3 Non-GASP user input variables used in FIELDOP. . . . . 85 5.4 Statistical distributions used for activity times in model FIELDOP . . . . . . . . . . . . . . . . . . . 90 5.5 Input parameters for distributions used in FIELDOP . . 91 5.6 FIELDOP variables monitored by GASP IV subroutine COLCT. O O O O O I O O O O O 0 O O O O I O O O O I I O .103 5.7 FIELDOP variables monitored by GASP IV subroutine TIMSTO O O O O O O O O I O O I O O I O O O O O O O O O 103 5.8 List of non-GASP variables used in FACYARD . . . . . . 107 5.9 Files used in model FACYARD. . . . . . . . . . . . . . 108 5.10 Statistical distributions used for activity times in in model FACYARD . . . . . . . . . . . . . . . . . . . 108 5.11 Attributes used in model FACYARD . . . . . . . . . . . 109 5.12 FACYARD variables monitored by GASP IV subroutine TIMST. I O O O O O O O I O O O O O O O I O O O O O O O 117 5.13 FACYARD variables monitored by GASP IV subroutine COLCT. o o o o o o o o o e e o o o o o o o o o o o o o 118 5.14 Input parameters and distributions used in FACYARD . . 119 6.1 Output from model FACYARD for current yard config- uration at Carrington Factory (Simulation time = 700 minutes) . . . . . . . . . . . . . . . . . . . . . 123 viii TABLE 6.2 6.3 6.4 6.8 6.9 6.10 6.11 6.12 6.13 6.14 LIST OF TABLES (cont'd.) Comparison of observed and simulated values for selected system performance indicators for factory yard operations. . . . . . . . . . . . . Output from model FACYARD for extended work period at Carrington Factory (Simulation time = 960 minutes) I I I I I I I I I I I I I I I I I I I I Output from model FACYARD for configuration with l weigh queue and 2 scales (Simulation time 700 minutes) . . . . . . . . . . . . . . . . . . . Output from model FACYARD for configuration with 2 weigh queues and 2 scales (Simulation time -= 700 minutes) . . . . . . . . . . . . . . . . . . . . Output from model FIELDOP for 10:1:2:2: equipment combination. . . . . . . . . . . . . . . . . . . . . Comparison of observed and simulated values of selected system performance indicators for field operations I I I I I I I I I I I I I I I I I I I I Equipment combinations simulated for field Opera— tions I I I I I I I I I I I I I I I I I I I I I I I I Machinery specifications and cost structure assumed for combine harvesting of sugar cane based on 1982 figures I I I I I I I I I I I I I I I I I I I I I I I Estimated cost of combine harvesting of cane for a field system with 2 field tractors.(One factory yard scale). . . . . . . . . . . . . . . . . . . . . Estimated cost of combine harvesting of cane for a field system with 1 field tractor.(0ne factory yard scale) I I I I I I I I I I I I I I I I I I I I I Estimated cost of combine harvesting of cane for a field system with 2 field tractors. (Two factory yard scales) . . . . . . . . . . . . . . . . . . . Estimated cost of combine harvesting of cane for a field system with 1 field tractor. (Two factory yard seal e8) I I I I I I I I I I I I I I I I I I Response of total system cost/tonne and total annual output to changes in the number of wagons . . ix PAGE 124 . 126 127 129 . 131 . 132 . 135 154 156 . 157 . 158 159 160 TABLE 6.15 6.16 6.17 LIST OF TABLES (cont'd.) PAGE Response of cost per tonne of cane harvested to an additional field tractor or factory scale I I I I I I I I I I I I I I I I I I I I I I I I I I 161 Response of annual cane delivery from the field system to an additional field tractor of factory scale I I I I I I I I I I I I I I I I I I I I I I I I I I 161 Cane delivery rates and the number of harvesters required for various equipment combinations . . . . . . . 164 FIGURE 5.9 5.10 5.11 5.12 5.13 5.14 5.15 6.1 6.2 6.3 6.4 6.5 6.6 6.7 6.8 6.9 LIST OF FIGURES (cont'd.) Flowchart through subroutine SHTDN. . . . . . . . . Flowchart through subroutine STRTUP . . . Causal loop diagram for the factory yard system Flowchart through FACYARD , . . , . Flowchart through subroutine ARIVAL Flowchart through subroutine ENDWGH , , , , Flowchart through subroutine ENDULD , . . . . . . . Simulated daily cane delivery vs. number of wagons in SyStem I I I I I I I I I I I I I I I I I I I I I Simulated daily cane delivery vs. travel distance to factory for mean factory residence time of 9.0 minUteS I I I I I I I I I I I I I I I I I I I I Simulated daily cane delivery vs. travel distance to factory for mean factory residence time of 73'98 m1nu£es o o e o o o o o o o o o o o e o o o o o Harvester utilization vs. travel distance to factory for mean factory residence time of 9.0 minutes,-, Harvester utilization vs. travel distance to factory for mean factory residence time of 73.98 minutes, , Road tractor utilization vs. travel distance to facrory I I I I I I I I I I I I I I I I I I I I I I I Current layout of, and current path of cane transport vehicles through Carrington Factory Yard, , , , , , , Proposal #1 for re-organization of Carrington Factory Yard I I I I I I I I I I I I I I I I I I I I I I I I I Proposal #2 for re-organization of Carrington Factory Yard I I I I I I I I I I I I I I I I I I I I I I I I I xi I I I I I I I I PAGE 99 . 100 . 105 I 111 113 . 114 . 116 . 136 . 138 . 139 . 140 . 141 . 143 I 147 148 150 FIGURE 1.1 2.1 2.2 3.1 3.2 3I3 4.1 5.3 5.4 5.5 5.6 5.7 5.8 The geographic location of Barbados. . The location of estate and small holder land . Percentage of small holders producing various LIST OF FIGURES amounts of cane. Symbolic model of closed circuit queuing system. Generalized models representing the behavior of systems. . Regions in the (B1, BZ) plane for various distribu- tions . . Flow process chart of Operations performed on each transport unit . Basic modes of GASP IV control . Functional flowchart of a GASP IV program. Activities involved in the field subsystem . The 2- and 2- parameter forms of the Gamma distribution . . Flowchart through data analysis program DATANAL. Causal loop diagram for the field system . Flowchart Flowchart Flowchart Flowchart Flowchart Flowchart Flowchart through through through through through through through FIELDOP. . subroutine subroutine subroutine subroutine subroutine subroutine xii I I I ARFLD STLD. NDLD. ARHTR . ARFCT . DPFCT . PAGE . 12 . 13 . 29 . 30 . 31 . 54 . 59 . 61 . 63 . 72 . 73 . 78 . 84 . 87 . 89 . 92 I 95 . 95 . 96 CHAPTER 1 INTRODUCTION 1.1. Background Barbados, located 59° 37' east longitude and 13° 4' north latitude, is the most easterly of the archipelago of islands in the Caribbean Sea. It is a small island of 430 square kilometres, 43,000 hectares of land and a population density of approximately 580 persons per square kilo— metre. (Figure 1.1). Colonized by the British for over 300 years, the island has an economic history which is divisible into three main phases. The first phase, 1627-1650, was characterized by a peasant economy in which an inexcessive number of settlers pursued a relatively self-sufficient and diversified economy based on the production of tobacco, cotton and indigo. The second phase, 1651—1950, saw the island transformed into a rigid and 'lop-sided' export oriented plantation economy dependent almost exclusively on sugar cane production. During this period scarcely any effort was made to exploit what limited opportunities existed for diversification, even within the dominant agricultural sector. The third and current phase, which began in the 19503, has so far been characterized by the political decline of the white planter class and a concomitant increase in the number of somewhat more .moomoumm mo nowumooa owemmumoow one ~.~ ounwfim I. I'llo’l‘... 138,320»; «$4.. .4550 .. . . . ., . «3233 I I IIIIIII - .. II, c w '0 I\ - ...—.....- .Ool o . .I".l .a . I ...tl .aIlO.’ . l.- .O. . o I at. e. . o. n o .o O u u l ... ..N . .. . . a . .t . I ’ I0. I .I,. 0 . I. a .... . | a. I . a o I . I o .H I. c . I I. u "‘2 8922-». 523:: .5 2 . mafiozoouo mumoaos Hamam mo mmwmucmoumm .~.N muswwm .3263 3:29am 2me «3.33m on. .0 3.303838 0cm... .0 mecca» o. c on o~ . - 1 - X wound 14 2.3.2. The Processinngector The manufacture of cane sugar in Barbados began in the 16403 following the introduction of a model of a sugar mill and some skilled artisans from Brazil (Barbados Sugar Review, No. 47, 1981). The early factories were extremely inefficient, the output per factory in the ciecade 1700-1709 averaging just 15 tonnes sugar per year with a tonnes (:ane/tonnes sugar ratio of 18 to 20. In 1709, there were 485 sugar mills, 409 of which were driven by twind and 76 by animal power. Thereafter the processing of cane remained t:echnologically stagnant until the 18408 when the first steam plant was iaitroduced (Barbados Sugar Review, No. 47, 1981). This switch to steam- powered plants generated a steady decline in the number of factories and by 1958 only 26 mills remained. Sugar factories in Barbados are still relatively small and much of ttmdr equipment is outdated. Maintenance of the sugar mills has, there- fcue, posed increasingly difficult problems over the last 20 years and, 618 a result, the number of mills kept in operation has been progressively ‘reduced from 26 in 1958 to 17 in 1969. Prior to 1970, sugar factories were operated largely by individual Owners or private individual companies and closures of factories were based on individual financial decisions, rather than on the interests of the industry as a whole. In June 1970, however, a single company, Barbados Sugar Factories Limited (BSFL), was incorporated to collectively own and operate all existing factories. With the birth of this company, a rigorous rationalisation programme (aimed at planned restructuring of the processing sector of the industry) ensued, leading to the further La 15 closure of the smaller, marginal and uneconomic mills and the moderniza- tion and improvement of those mills that were retained. During the rationalisation process, the island was divided into three zones; North, ICentral and South and two factories were allocated to serve each zone. '10 date five of the original sugar factories are in operation; two each in.the South and Central zones and one in the North. In addition a new tnodern sugar factory has been constructed in the North zone and was com- tnissioned during the latter half of the 1982 harvesting season. The (erection of this new factory represents a major capital injection into tihe processing sector since the last factory established was built as IJDng ago as the year 1920. With these six factories in full operation, tine theoretical throughput capacity of the processing sector is in the rwegion of 480 tonnes per hour or 11,500 tonnes per 24 hour working day. The Barbados Sugar Factories Limited holds majority shares in three siervice companies associated with the processing activity: the Barbados Idolasses Terminal Limited (97%), the Sugar Terminal Limited (56%) and t1m28ugar Transport Limited (100%). The BSFL is primarily a grower- risideration of all relevant economic factors (Barbados Sugar Review, No - 41, 1979). Given normal market conditions, this price tends to be substantially higher than free World market prices. However, since Lome' prices are indexed to EEC price ranges, annual fluctuations in the price received by ACP producers are more strongly influenced by the EEC Supply and demand schedules than by production costs in the ACP coun- tries themselves. The net result of this is that the guaranteed Lome’ Price has not risen as rapidly as production costs have in the Barbados Inch-ISI:ry. The indication is, therefore, that it is essential for the BarbEldos Sugar Industry to keep its production cost trend as closely paralleled as possible to the movement of production costs in the EEC, in oI‘Cier to maintain present profit margins. This implies increased productivity of existing sugar lands and decreased production costs and, Since harveSt labour costs are a major component of production costs, inc )3 eased harvest mechanization. The actual economic return to Barbados from sugar sales to the Eu 1:9 Dean Common Market reflect not only the annually negotiated Protocol Pr iqe but also any premia obtained from individual purchasers within the Community (McGregor et a1., 1979). The largest purchasers of 19 Barbadian sugar within the EEC have been the United Kingdom and Ireland. Sugar supplied to refineries in these two countries commands a premium over the Lome’ negotiated price. 2.4.2. The World Market In global terms, the free (World) market represents a residual out- let for the relatively small proportion of the world's sugar which is not consumed in the producing countries or traded under pre-negotiated Inarket arrangements such as the Sugar Protocol of the Lome’ Convention. Ihie to the residual nature of the world market, cyclical imbalances bet— Ween supply and demand (mainly as a result of weather stimulated supply \rairiations) often generate wide fluctations in the price of 'free' Sugar. Sale of sugar on the world market is governed largely by the llriternational Sugar Agreement (ISA). Under the terms of reference of tillis agreement, Barbados is regarded as a small exporter and is granted E111 annual export entitlement of 70,000 tonnes of raw sugar. This en- t:i‘tlement is not subject to quota adjustments, neither is Barbados ‘Piéquired to observe the ISA's stock provisions with which larger export- ing countries must comply (Barbados Sugar Review, No. 41, 1979). In recent years, the main world market outlets for Barbadian sugar have been the Canadian and U.S. markets. Not all of the sugar exported to North America is in the form of 'tlle crystalline final raw product. A significant portion of the exports tx) the U.S. and the bulk of the exports to Canada are in the form of fancy molasses, (a specialized intermediate product of the processing 20 activity) which normally commands a premium price over raw sugar (Ministry of Finance, 1975-81). The 1977 International Sugar Agreement established a target price range of 24 - 47 U.S. cents per kilogram of raw sugar. The average ISA free market spot price for raw sugar for the period 1977-1979 was 1 7 .90 U.S. cents per kilogram, a price which was below the average f inancial costs of production of most producers in the world (McGregor et a1., 1979). Owing to the large excess supply of world sugar over demand, however, even the lower end of this target price range was not reached until 1981 when the price reached 31.25 U.S. cents per kilogram. 2 - 4 . 3. The Domestic Market The domestic consumption of sugar in Barbados averages 15,000 tonnes annually (McGregor et al., 1979) and, with a population of about 250,000, this works out to be an annual per capita consumption of 60 kilograms. This figure is high by world standards (OECD/OCDE, 1975, 19 79) and, with the demand for sugar in Barbados practically income in- elastic, one cannot envisage any growth in domestic demand for sugar, except a gradual response to total population increase. Before 1976, factories received a fixed price of BDS $376.17 per tonne 96° Pol, or 19 U.S. cents per kilogram, for all grades of sugar. Since then, however, ex-factory prices for domestic sales have risen BUbStantially. In 1977, price differentials between all grades of Sugar (Browns, Yellows and Straws) were established to reflect product (1qu ity and the general level of ex-factory prices even rose above the Dr iQeS obtained for foreign sales other than to the European Economic (3Q tun‘llnity. ¥ 21 Currently confirmed potential sales of the Barbados Sugar Industry total some 139,000 tonnes, raw value; (54,000 tonnes to the EEC; 70,000 tonnes through the ISA and 15,000 tonnes for domestic consumption). McGregor et a1. (1979) have estimated that, if the world sugar prices were to exceed 33 U.S. cents per kilogram, Barbadian exports to the free market could well surpass the 70,000 tonnes allowed under the ISA agrenent. Unfortunately, however, the Industry has not been able to fully exploit the available sales opportunities over the last few years. In 1979, 1981 and 1982 sugar production was only 112,000, 96,000 and 96 ,000 tonnes representing shortfalls in potential sales of 20 percent, 31 percent and 31 percent, respectively. 2.5. History of Mechanical Harvestig of Sugar Cane in Barbados Mechanical harvesting of sugar cane was first introduced to Barbados in 1968 when two machines, a Toft J-150 and a Crichton, were imported frOm Australia. Both of these machines were wholestick harvesters which cut the cane at the base and deposited it in piles for subsequent pickrup either manually or by a mechanical loader. These machines were operated in both burnt and green cane but, even in burnt cane (for Which they were primarily designed) their performance was unsatisfactory. Both machines were set on 1.5 m wheelbases which were not only narrower that). the inter-row spacing of the cane but were inadequate for the rolling terrain and stony conditions that existed in most of the fields at the time. The machines were reportedly very little used (McGregor at 1979) and, since no significant efforts were made to adopt a1- ’ f i Qld conditions and inter—row Spacing to the design constraints of ¥ ‘. turd 5‘! H: r. (I) u 'A . ’_'-? c) 22 the harvesters, it is doubtful whether they were ever really subjected to a serious field evaluation. In 1970, a Cameco wholestick harvester and a push-pile loader were imported from the U.S.A. The harvester was large and proved to be too cumbersome for operation in the small, irregularly shaped Barbadian cane fields. The push-pile loader, however, generated considerable interest among the larger cane growers. In 1971, a decision was taken to switch from wholestick to chopper- type cane harvesters. That same year, two chopper harvesters, a Toft Cid-364 and a Don Mizzi 741, were imported from Australia. The follow- ing year, a second Don Mizzi 741 was brought in and in 1973, a Don Solo chopper harvester was imported. All of these machines were designed to operate in burnt cane and controlled burning of cane was, therefore, a Prerequisite for their use. The Toft CH-364 was a rather large, self-propelled combined har- vester, well furnished with engine power. It was a quick cutter but, under the stony field conditions that existed at the time, numerous problems were experienced, particularly with the machine's cane eleva- tion system. Its wheelbase, though wider than that of the wholestick haI“"esters, was still too narrow and even on relatively gentle lepes, 1El':et‘al stability was a problem. The Don Mizzi harvester was essentially a cane base-cutting, e1 e\’a.ting and cleaning mechanism attached to a reversed L—shaped frame, on to which a standard 60 kw agricultural tractor was mounted. The drive train of the tractor was linked to the main drive shaft of the ha rVesting mechanism by means of a duplex chain and sprocket coupling. ¥ Q. .\ (Q 'h 23 Tractive power, as well as mechanical and hydraulic power for the moving parts of the harvesting mechanism, were derived from the power unit of the tractor. The resulting combination was a relatively compact, self- propelled unit set on an extended wheelbase 3.7 m in width. Like the Toft CH-364, the Don Mizzi 741 machine delivered cane to the left side only, making two-directional cutting impractical in the smaller cane fields. The Don Solo chopper harvester imported in 1973 was a more compact, integrally assembled, self-propelled version of the Don Mizzi with the capability of delivering cane to either the right or left of its direc— tion of travel. On account of this, it was theoretically more suited for work in the smaller cane fields. However, because this machine was set on a narrow wheelbase of just 1.5 m, and because most of the smaller fields are located on sloping and undulating terrain, machine stability during operation became a major problem. Among the machines mentioned above, the Don Mizzi 741 was undoubt— edly the most satisfactory, largely because of the extended width of its wheelbase and because concerted efforts were made to prepare the fields for mechanical harvesting according to the guidelines suggested by Baxter, 1969. By the 1973 harvesting season, many of the teething problems associated with this machine had been overcome, and a well organized and technically efficient mechanical harvesting system had been built around the harvester. Just then, however, two independent St‘lciies conducted by Chase and Eavis (1972) and Hudson (1972) showed that burning of canes at harvest time reduced the organic matter content O f the soil to alarmingly low levels and caused a significant reduction k 24 in subsequent crop yields. In addition, reports from the Entomological Division of the Ministry of Agriculture indicated that the populations of predators, which were heavily relied upon for biological control of economically important pests of sugar cane, were likely to be signifi- cantly depleted if burning continued. Based on these findings, burning of sugar cane prior to harvesting was outlawed. Following this, efforts were made to adapt the Don Mizzi system to green cane harvesting but the harvester's unsuitability to this task prevailed and, in 1974, the Don Mizzi system was finally abandoned. Nonetheless, the experience with the Don Mizzi program proved to be very valuable. Not only did it demonstrate the applicability of chopper lizalrvesting to Barbados, but it also provided a foundation of knowledge firwom which more recent chopper harvester operators were able to benefit. In 1975, two Toft 300 chopper harvesters, reported to be capable cane Of handling green cane, were brought into Barbados by individual Producers. Based on the performance of these, two more chopper har— vesters, one Toft 4000 and one Massey-Ferguson 205, were imported by agricultural machinery distributors in 1979 for contract harvesting operations. Between 1979 and 1982, 17 additional green cane chopper harvesters, 1 Massey-Ferguson 305 and 16 Toft 6000 machines, have been bI.Q"-Jl-ght into the country. In terms of field operating procedures, the Mag Sey-Ferguson and Toft cane harvesters are quite similar but, in terms of design, the Massey-Ferguson machines rely heavily on mechanical drives and linkages while the Toft machines are almost exclusively hydraulically driven. The rapid increase in the number of chopper harvesters in Barbados Q . Vet the last three years has occurred without a corresponding rapid ¥ V a u\u 25 development of machine management capabilities and without adequate re- organization of cane receiving facilities at the factory yards. This has resulted in rather inefficient use of most of the individual mechan- ical harvesting units. In the late 19603, the Barbados Sugar Producers Association, in- spired by the need for a harvesting system specifically suited to Barbados conditions, initiated development of a wholestick machine system in collaboration with a British firm, McConnel Engineering. The first prototype cutter began working in 1972. This machine, which sub- sequently became known as the BSPA/McConnel Stage I, cut the cane at the base and left it in a swath along the row for manual cleaning and loading. Seven additional Stage I machines were introduced in 1973 and by 1976 some 37 of these units were in private ownership. Owing to the large number of workers required to work behind these machines as re- trievers, machine productivity was very low and probably for this reason, most of the machines were little used. In an effort to eliminate the heavy labour requirement behind the Stage I, a BSPA/McConnel Stage II machine was developed, the first one being used in 1975. This was essentially a cleaner-piler which picked up the swath formed by the Stage I cutter, removed most of the cane tops and trash, and deposited the cane in piles for final pick-up by a mechanical grab loader. By 1979, seven of these machines were privately own ed , mainly by the larger estates. They, too, were used very little a “(1 most of the Stage II owners and potential owners have now invested 1 t1 QEine combine harvesters. The Stage I machine has been retained and further refined by a 1 Q in1 firm called CARIB Enterprises. The current model, now known as k .3” av r. Si 26 the CARIB cane cutter, is a vast improvement on the original Stage I prototype. The most major modification was mounting the cutting mechan- ism on a reversed, 4-wheel drive tractor so that considerably more *weight is now on the front end of the machine, and traction and stability .are greatly improved. Based on observations during the 1982 harvesting season, the CARIB machine seems to have a lot of potential for work in ‘the more hilly central and eastern areas of the country that are not zaccessible to the larger chopper harvesters. CHAPTER 3 LITERATURE REVIEW 3.1. Introduction Tomlinson (1973), in an assessment of the problems of the Caribbean ssugar industry, highlighted the need for greater efficiency in organiz- ing and planning the use of harvesting equipment and proposed that the techniques of Operations Research could be found useful in this respect. Operations Research, as described by the Operations Research Society of America, 'is concerned with scientifically deciding how to best design and operate man-machine systems, usually under conditions requiring the allocation of scarce resources' (Phillips et a1., 1976). Broadly speaking, the application of Operations Research techniques requires the construction of a model which incorporates all the impor- tElnt characteristics of the system (Shamblin and Stevens, 1974). Boyce (1972a) classifies three types of models; iconic, analogue and symbolic. IcOnic models physically resemble the subject of inquiry and are char- ac terized by some scaling effect. Analogue models are characterized by the use of a convenient transformation of one set of properties for ant) ther. Symbolic models are characterized by the representation of a. system by mathematical or logical symbols and expressions. 27 28 3.2. Symbolic and Mathematical Models Most harvesting systems can be represented by the symbolic model, as illustrated in Figure 3.1, and can be generally grouped as closed circuit transport syst-s (Boyce, 1972b). In general, a closed circuit tranSport system can be broken down into two separate sub-systems, one located at the field and the other at the installation (factory, mill (Dr other facility). Each sub-system can be represented by a queuing model, in which transport units queue at the field and factory awaiting service. The length of time each arriving unit must wait in line (queue) for service depends on the number of other units already in the <1ueue, the service rate (or service time) and the rate of arrival at the queue (Saaty, 1957). In a real life system, arrival and service rates are not constant but subject to variation about a particular mean. The successful use of mathematical models to represent real queuing Systems, therefore, depends primarily on the identification of the form 0f the probability distributions which best represent the actual arrival and service rate distributions (Page, 1972; Phillips, 1976; Hillier and IJileberman, 1974) . In a general discussion on the use of mathematical models to de- Sc1':ibe phenomena in agricultural systems, Demichele (1976) proposed tI'lli‘ee general forms which could be used, depending on the state of know- ledge (Figure 3.2). Based on the system being studied, these forms may be represented by different statistical distributions. Hahn and Shapiro ( 1967) suggested that, in the engineering context, there are often in— sufficient grounds for choosing a Specific model. They indicated that I“Zirgmre 3.3 (adapted from E.S. Pearson, University College, Brandon) may 29 QUEUING SERVICES Loading -‘ ' f"""""1l ___% |.- - .. ,______T ------- 1‘ - l ___..| -- --..-.4 l ——-—— |r--'""--1' u n ‘_—_d‘. - ‘r v ISON QUEUING SERVICES NON QUEUING SERVICES Travel loaded ' Travel empty . I .-t. J Hi COMMON QUEUING SERVICE Unloading Figure 3.1 Symbolic model of closed circuit queuing system 3O Probability of an outcome 3 all possible x A 1) The distribution function of a system where nothing is known Probability of an outcome x all possible x (L2) The distribution function of a system where some things are now known Probability of an outcome x all possible x (3 ) The distribution function of a system where most things are known JPTlgure 3.2 Generalised models representing the behaviour of systems §3 F igure 3 .3 45 0| OI q 31 ‘1 h . .. 1 aw ".9 ‘ u wrn‘tm| agi‘ fit?“ ‘3 5".“ .&# is“ 4* W“ .nnM" 1.2.3,. .. “W IMPOSSIBLE muss” “53"" '~.a- n.- l “W‘*W‘fi1"'blwi‘ , mwmmaf‘ .‘u t- ‘r.. 2 iv“: ...": ”1.12""“J‘HIE. ‘fi‘tflflwI I 21; a?” {NET-L: L¥* ... " “I I ’ r 1‘» mo ‘1‘ a» “Mes! ...... 3% \s.» ...... ..w v.6“ ° 11"7'3'u‘ asses; ,th ..- ....31ugI\.I‘I'I\1h “1‘31" .I‘ J!" I".|.-.H|I"I"1"'I‘:h"“‘l “.IJI' "um. I W‘! 'IiI'!‘ “M“; 3'” :4‘;\ H " “L ‘J HHIJII {-T‘ ‘1'1'11‘ Ikéxw“\$_'¢ "‘ ‘1‘." 2;: .‘Is‘fl ‘4‘, 41M "3“; éu' WI". "-I “WWW“? 3" «‘9‘ “51"”; -.'- fins-‘33,; 09’!“ . 3’13“ 1 1- Pt}? I. ' . Til? Lfef- -10‘ ‘&1|\\%.' “LN...“ ‘M “Jr“? M?” Var-1‘15 A“? 1 r'.‘ ‘%Id '18 1x ‘ 41.1%.,“ Ilfl‘iw 1" ... w““‘hh1‘|§fi"c'$“h. EXPONENTIAL DISTRIBUTION I 2 3 4 Bl Regions in the (81,82) plane for various distributions 32 be used for such a selection. The values of B1 and 82, the square of the standardized measure of skewedness and the standardized measure of peakedness, respectively, can be estimated using equations (1) and (2) below. The estimates so obtained are, however, very sensitive to a few extreme observations and should, therefore, be used with caution, particularly in situations where the number of observations is less than 200. On account of these factors, it is generally desirable to set up a frequency table to enable the fitted distribution to be com- pared with the observed data. 3/2l2 Bl=b1=IM3/M2) (1) 2 32 = b2 = Mg/(Mz) . . . . . . . . . . . . . . . ... . . (2) mahere N _ 2 M2 = l/N 2 (X1 - x) . . . . . . . . . . . . . . . . . (3) =1 N _ 3 M3 = l/N 2 (X1 - X) . . . . . . . . . . . o . . . . . (4) 1: N _ u M. = l/N 2 (X1 - x) . . . . . . . . . . . . . . . . . (5) i=1 xi = 1th value of the variable x x = sample mean Hahn and Shapiro (1967) further pointed out that the initial selec- tlzirbn of a model should be based on an understanding of the underlying 33 physical phenomena and that a distributional test should be used to evaluate the adequacy of the physical interpretation. Boyce (1976) indicated that the process of selecting suitable statistical distributions and estimating the appropriate parameters could, in some cases, involve a virtually impossible sampling task if the real system were operating in an unsteady and/or erratic way. He further suggested that the use of an optimistic, a pessimistic and a Imost likely estimate of the particular quantity of interest, based on .field samples, might provide a basis for estimating the parameters of as selected distribution. Dumont and Boyce (1972) described three methods of fitting three ssuitable distributions to observed work time data for two agricultural tmnit operations; tranSport loading and combine harvesting: These are as onllows: (1) the method of maximum likelihood for the Gamma distribution, (2) the method of matching moments for the Beta distribution, and (3) the method of matching percentiles for the Johnson SB distribution. The Chi-squared test statistic was used to test the goodness of fit 13‘:>2:r each set of parameters determined. Their general conclusion was that, in practice, any one of the three distributions would be acceptable ‘3‘:’?t? use in simulation modelling to represent the distribution of transport :I‘CDHEiding times, but that the Gamma distribution seemed incapable of pro- xifii~cling an acceptable fit for the combine harvesting times data. Bouland (1967) conducted a study of truck queues which formed at (:CDIIntry grain elevators in the U.S.A. and found that transport unit 34 arrival rates could be described by Poisson distributions for rates of less than 35 arrivals per hour, and by Uniform distributions for rates exceeding 35 arrivals per hour. Service times of weigh scales and unloaders were found to fit Erlang distributions. The probability density functions of these distributions are as follows: The Poisson distribution f(x) e‘M M3 x = 0,1,2. . . . . . . . . . . . . . . . . . (6) where value of discrete variable x M 8 expected value = variance Source: Bhattacharya and Johnson, 1977. the Discrete Uniform distribution f(x) l/N x = l,2,.....,N . . . . . . . . . . . . . . . . (7) where value of discrete variable K Source: Bhattacharya and Johnson, 1977. L393§2:e Erlangidistribution (Anx(n-l)e-Ax)/(n-l)|,x > 0 . . . . . . . . . . . (8) f(x:n,A) = “°qb1£ere x = value of continuous variable A = scale factor n = shape factor (restricted to integers) Source: Manetsch and Park, 1980. 35 Dumont and Boyce (1972) found that the Gamma distribution, the more general form of the Erlang distribution, gave an adequate repre- sentation of transport loading times for various field conditions. The Gamma distribution has the following probability density function. f(x:n,A,u) = An(x - u)n-1 e-A(x-U). . . . . . . . . . . . . . (9) where n > O; A > O; u < x < co; -m < u < w and G(n) is the Gamma function given by G(n)= 2°(x-u)”'1e'(x'“)dx................(10) where x = value of continuous variable A = scale factor n = shape factor n = location parameter Brooks and Shaffer (1971) reportedly developed a mathematical model 11C> predict the output for tipper trucks hauling dirt from an excavation site to a landfill area (Ogilvie et a1., 1978). Their approach consisted of six basic steps: (1) determination of all possible combinations of shovel and truck, each combination being called a state; (2) measurement of the ability of the shovel to change each state defining the rate of transition between states due to the shovel by the mean service rate; (3) measurement of the ability of the truck to change each state defining the rate of transition between states 36 due to the trucks by the arrival rate of the trucks at the shovel; (4) determination of the net effect of changes by shovel and by trucks; (5) determination for steady—state conditions, of the per- cent of time on average that an Operation was in any particular state; (6) calculation of the steady state output of the system by multiplying the service rate by the Production Index (PI), where PI = sum of the percentage of time that the system was in each productive state. The main objective of this model was to improve the opportunity for the contractor to apply queuing theory in practice by providing a pro- (zedure which did not require knowledge of the mathematical rigour of caueuing theory, but at the same time provided a method of predicting C>l1tput, which could be rapidly and easily applied. Based on their model, Brooks and Shaffer (1971) reportedly developed a set of curves for var- iOUS combinations of the mean service and arrival rates and the number of trucks for simple single server and multiple server systems, as well as systems in which the shovel had a hopper. Audsley and Boyce (1973) used a mathematical approach of queuing t".‘Il‘iory to develop models of cyclic transportation systems. In their ‘iICDIEIC, they assumed that service times were independent and identically distributed and could, therefore, be represented by some form of Erlang distribution. They obtained similar results to the production index ‘flilich Brooks and Shaffer obtained in 1971 (Ogilvie et a1., 1978). 37 3.3. Computer Simulation Models An approach to the study of large and/or complex systems which is rapidly gaining popularity is computer simulation. Phillips et al. (1975) stated that simulation has become one of the most widely used and accepted tools of systems analysis. Broadly speaking, simulation is a problem-solving technique for defining and analysing a model of a system (Dent and Blackie, 1979). To simulate means 'to duplicate the essence of a system without actually attaining reality' (Rockwell, 1965). Underlying any simulation is a mathematical abstraction of the system which relates various system functions. In a paper on the use of simu- lation methodology in agriculture, Rockwell (1965) suggested several reasons for using simulation analysis. Most of these reasons have been re-iterated in a more recent work by Naylor (1971). They are as follows: (1) The system may be too complex for analytical solution in which case, simulation can yield valuable insight into which variables are more important in the system and how these variables interact. (2) With simulation, it is possible to build in time delays, non-linearities, irregular distributions and discontinuities into the system. (3) Systems in which time is a critical factor are well suited to computer simulation. For certain types of stochastic problems, the sequence of events may be of particular importance. Information about expected values and moments may not be sufficient to describe (4) (5) (6) (7) 38 the process. In these cases, simulation methods may be the only satisfactory way of providing the required information. Simulation can be used to experiment with new situa- tions about which we have little or no information. The simulator permits later decisions to be based on earlier system output as in a dynamic programming framework. The process of designing a computer simulation model forces the researcher or designer to explicitly describe the system processes and the required data. The know- ledge obtained in the design activity frequently suggests changes in the system being studied. The effects of these changes can then be tested, via simula- tion, before implementing the changes on the actual system. The process of simulation design is in itself a valu- able educational tool and it has been used by many companies as a pedagogic device to help management understand the characteristics of the systems which they control. Simulation permits us to experiment with systems that, in reality, it would not be possible to experiment with. Through simulation, one can study the effects of certain informational, organizational and environmental changes on the operation of a system by making alterations in the model of the system and observing the effects of (8) (9) (10) (ll) (12) ‘ (13) (14) 39 these alterations on the system's behaviour. Simulation can serve as a 'preservice test' to try out new policies and decision rules for operating a system, before running the risk of experimenting on the real system. Simulation can often lead to a training device such as a management game. Since the simulation process gives no valuable insights into the qualitative aspects of human decision making, and because of its natural realism, it often becomes possible to convert the simulator into a simulation training device. Simulation models usually provide better user accept- ance than analytic models. The systems user can see the reality of the simulation and thus will usually have more faith in the conclusions from the simulation output. The interpretation of simulation results does not usually demand a mathematical background on the part of the user. Mbnte Carlo simulations can be performed to verify analytical solutions. Simulation is well suited to sensitivity analysis in which key system parameters are selectively altered to test their contribution to overall system performance. When new elements are introduced into a system, simula- tion can be used to anticipate bottlenecks and other *‘Ifl/ ’5 40 problems that may arise in the behavior of the system. Computer simulation, as defined by Pritsker (1974), is the estab- lishment of a mathematical logical model of a system and the experimental manipulation of the model on a computer. Development of the model re- quires an in-depth analysis of the system in order to identify its important characteristics and components. Such components include entities, their attributes and events, which may be defined as follows: Entities - objects within the boundaries of the system being studied Attributes - characteristics of entities within the system Events - occurrences which cause change in the status of the system. Simulation can be analogue or digital, discrete or continuous, with or without a computer, with or without real time (fast or slow) and with Or without a human decision maker in the simulated process (Rockwell, 31-5?€559. Digital simulation requires that for each instant of model time, a Series of calculations be performed to produce a set of discrete out- ‘:>‘11:s. The reaction of the model to any input function is, therefore, (l‘iilermined by repeating the series of calculations for each instant of time for which the reSponse is required. Analogue simulation requires that model variables simultaneously assume their appropriate values so tjlat parallel recording of these values can reproduce all the important \r a”)... 41 aspects of the system's performance (Link and Splinter, 1970). Generally, the rapid development of high speed digital computer technology, and of simulation languages based on such technology, has made digital simula- tion preferable. Hillier and Lieberman (1974) stated emphatically that simulation models need not be completely realistic representations of the real systems, since representing all of the minute details of real systems often leads to excessive programming and use of excessive amounts of computer time for the benefit of a small amount of additional informa- tion. These authors further suggested that the behavior of system elements is best represented by theoretical distributions which best fit observed data from which random samples can be drawn. In light of this, an important aspect of simulation modelling is the validation of the model to show that it adequately represents the 'real-world' situa- tion. Model validation can be achieved by comparing observed 'real—world' £3)rstem performance data with data generated by the model (Hillier and Lieberman, 1974; Manetsch et a1., 1974). Standard statistical tests E311(2has the Student's 't' test for the comparison of two means and the X2 and F tests for inferences about variances, can sometimes be used to Cletermine whether the two sets of data are statistically different. These procedures are designed to make inferences about the values of 't:11€3 parameters 'u' and 's' that appear in the prescription for the mathematical curve of the normal distribution and are collectively 1‘nown as 'normal-theory' parametric inference tests. Another useful statistical test which may be used for the compari- £3Onof two sets of data is the Wilcoxon Rank-Sum Test originally proposed 42 by F. Wilcoxon (1945). An equivalent alternative version of this test was independently proposed by H. Mann and D. Whitney (1947) and is now known as the Mann-Whitney U Test (Bhattacharrya and Johnson, 1977). Neither the Rank-Sum nor the U Test is restricted by the assumption that the data are normally distributed and they both test specifically whether two samples drawn from different populations have the same dis- tribution. These tests are non—parametric tests and they can be used for both small and large sample sizes, although some test power is lost as the sample size decreases (Bhattacharrya and Johnson, 1977). 3.4. Sugar Cane Simulation Models Simulation models have been used for about 15 years to study farm machinery systems. Specific areas of application have been: predicting expected returns (Sowell et a1., 1967); analysis of specific cropping systems such as cotton production (Stapleton, 1967); forage harvesting (XZOupland and Halyk, 1969); the performance of field machines and trans- I>C>rt units for a row crop planting system (Von Bargen and Peart, 1969); tillea effect of the different harvesting system configurations on closed Circuit cyclic transportation systems (Boyce, 1972); silage harvesting <:Iinssel et a1., 1977) and sugar cane harvesting. The literature on a‘F'Plication of simulation modelling to sugar cane harvesting is examined th‘ some detail below. Sorenson and Gilheany (1970) developed a model for testing different StIrategies and decision rules governing the deployment of equipment on a (lane plantation in the Caribbean. For this model, time for loading in tile field, cane transport travel time for a given distance and unloading 43 time at the factory, were represented as constants, and statistical models were developed for crop rotation, rain and trash, rained-out field conditions, mill capacity and transport unit breakdowns. These authors suggested that, with minor modifications, their program could be used to optimize the length of the cane crop. Farquhar (1972) discussed the potential of adapting the systems approach in general and simulation in particular, to the analysis of sugar cane production systems, and used a generalized model to highlight some of the complexities of such systems. Shukla, Chisolm and Phillips (1973) developed a computer program for analyzing harvesting, loading and transportation of sugar cane. The programme reportedly gave good results for the systems studied but was somewhat locale-specific and depended heavily on the coefficients derived from a time and motion study which had been on-going for several years. Early (1974) used an inventory model to simulate harvesting condi- tiions in the Philippines. The specific objective of this work was to g337nchronize field and factory operations under conditions of random rainfall. He identified three limitations to his model based on the Inetzhodology used in constructing the model, the data used in modelling tlllee system response and the sequencing of rainfall events. Despite these limitations, he concluded, after simulating 10 years of operations, that 'tillea existing policy of allocating a daily quota to all farmers, each in I’lnaportion to their share of the milling capacity, could be improved by adepting a system of reaping zones which showed better comparative yields Siach month of the cropping period. He also estimated that this policy ‘Ehange, depending on the utilization and availability of irrigation, 44 could increase yields by between 12.5 and 16.5 percent. A linear programming model was used by Tonsman (1974) to predict, for different periods, which fields should be harvested, and the amount and type of machinery that should be used at each location. This model included in the optimization process, such agricultural characteristics as seeding type and variety, soil preparation, irrigation regime, fertil- ization, climatic conditions, water availability, machine field capaci- ties, transportation equipment and the factory's production programme, among others. Based on recommendations which resulted from this work, improvements of 14 percent and 30 percent in the ratios of sugar pro- duced to cane milled and sugar produced to harvested area, respectively, were observed. Hoekstra (1973, 1974, 1975) constructed three simulation models for a single cane harvesting and handling system in South Africa. His first model examined the effect of mill stoppages on transport units which off-loaded directly into the factory. The results indicated that the total number of transport unit hours lost over the cropping period, for mill stoppages ranging from less than 10 minutes to 5 hours, was negligible. In his second model, he examined mill yard operations. In this model, four vehicle types were processed through one of two unload- ing and each vehicle had to be weighed in and out. In general, the simulation results showed that improvement in the cane flow through the system could be achieved by improving the communications between the mill yard and the different suppliers, and by increasing the amount of cane delivered directly into the mill. In the third simulation model, the influence of mill operations on the delay time between cutting and milling was examined. In general, it was found that the delay between 45 cutting and milling was dependent only on the respective timing of the cutting. Grinding over six days reduced both the mean and the Spread of delay times. Hoekstra also indicated that the spread of delay times about the mean was not due to the irregular grinding of deliveries but rather to variances from a strict first-in first-out policy of milling cane deliveries. Cochran and Whitney (1975) examined the effect of different numbers of transport units on field loader utilization. The overall results of their work were combined with a theoretical analysis adopted from the work of Melissa (1966) to develop a nomograph which permits graphical prediction of delivery rates for given values of transport unit capabil- ities, mean loading rate, total trip time and the number of transport units. A cost model was also developed to facilitate selection of the optimum number of transport units for any given system configuration. Loader transport system costs were broken down into three categories. (1) Labour cost - defined by equation 11 (2) Equipment fixed cost - defined by equation 12 (3) Equipment operating costs - defined by equation 13 The total loader transport system cost (TC) is obtained by summing the three equations 11-13 to give equation 14. The equations are as follows: Systsn labour cost = Wlo + (Nt)(Wtr)/Del . . . . . . . . . . (11) to System fixed cost = Lfc + (Nco)(Cfc) + Z chi/Thui) i=1 (l/Del) . . . . . . . . . . . . . . . . (12) 46 System operating cost = Loc/R + (Coc + Nt/Nc*Toc) ' (tP/(S)(Cc)) . . . . . . (13) to Total cost, TC = l/Del((Lfc + (Nco)(Cfc/Shu + Z chi/Thui i=1 + W10 + (Nt)(Wtr)) + (Loc/R + (Coc + Nc/Nc- Toc-tP/((S)(Cc)). . . . . . . . . . . (14) where TC = Total transport cost $/tonne Lfc = Loader fixed cost $/year Cfc - Cart fixed cost $/year chi = ith tractor fixed cost $/year Thui = ith tractor hours of use per year hours/year Shu = Season hours of use hours/year Wlo = Labour cost for loader operator $/hour er 8 Labour cost for tractor operator $/hour Loc = Loader operating cost $/hour Coc = Cart Operating cost $/hour Toc = Tractor operating cost $/hour tp = Round trip distance (field to mill) kilometres Nt = Number of transport tractors in use Nc = Number of carts in use R = Loader rate for given conditions tonnes/hour S 8 Average speed of transport units kilometers/hour Cc = Cart capacity tonnes Del = Transport system production rate to mill tonnes/hour Nco Number of carts owned 47 Nto = Number of tractors owned and reserved for transport system use only during harvesting season One of the more recent applications of simulation modelling to sugar cane handling was reported by Ogilvie et a1. (1978). These researchers investigated the handling and transport of hand-cut whole stick sugar cane at the Frame Cooperative Farm in Jamaica. Two computer simulation models were constructed; one for the field operations and one for factory yard operations. They reported that the models developed permitted the testing of current and modified equipment arrangements and management policies, which in turn resulted in higher throughput and optimum use of transport and handling equipment. They also reported that the GASP IV simulation language used for constructing the models not only proved to be most adaptable to the system studied, but also provided the user with excellent intrinsic filing and report formatting routines and cap- abilities. Loewer et a1. (1979) developed a computer model to optimally select sets of equipment for 56 possible alternatives through the sugar cane harvesting network in Brazil, and to compute the labour requirements and fixed, variable, indirect and total costs. The results of the simulation indicated that transportation costs accounted for 60 percent of the total cost of harvesting sugar cane under Brazilian conditions. It was concluded, therefore, that modification of the transportation activity offered the greatest potential for reducing total harvesting and handling costs. 48 3.5. Computer Simulation Languages In implementing a systems model on a computer, the user has the option of using a general purpose programming language or one of several special purpose simulation languages. General purpose programming lan- guages that may be of interest to the simulation modeller are FORTRAN, ALGOL, BASIC, and PL/l. On account of their generality, these languages may be used to construct simulation models of any type of system. How- ever, the modeller is required to develop his own input-output routines, set up his own time clock and switches within the model, and write his own Special purpose routines such as normal pseudorandom number generators (Dent and Blackie, 1979). Modelling using these languages requires con- siderable programming ability and a reasonable knowledge of the computer and its associated systems (Manetsch and Park, 1980). Special purpose simulation languages, on the other hand, have evolved in response to a need to reduce the programming skill and effort required to program and use computer simulation models (Krasnow and Merikallio, 1964). These languages contain specialised facilities, such as automatic time-keeping routines and sophisticated output formatting routines, which are convenient for modelling particular types of problems (Teichrow et a1., 1966). Some of these languages are really supersets of general purpose languages (for example CSMP is a superset of FORTRAN) whereas others, such as GPSS, are self-contained (Tocher, 1965). In any case, most of these languages were originally developed to satisfy the requirements of specific problems and they, therefore, differ in the type and range of their possible applications. Manetsch and Park (1980) described a number of criteria by which simulation languages may be 49 evaluated. These are as follows: Language type - Simulation languages are classified as either 'discrete event' or 'continuous flow' or both of these. Discrete event languages are designed to simulate real world systems when a microscoPic viewpoint, which considers individual objects or events as entities, is appropriate. In contrast, continuous flow languages are designed to simulate systans which are best characterized by aggregates or flows of discrete entities. Models of such systems are normally structured using differential and/or difference equations and continuous flow languages usually contain integration packages designed to efficiently obtain particular solutions for large sets of such equations. Universality - This criterion describes the generality of a language with respect to the range of computers with which it is compatible. Some languages are very machine dependent and are useable only on limited com- puter makes, models and sizes. Higher Order Modularity - Situations are sometimes encountered in which a complex sub-system structure occurs repeatedly in the structure of the larger total system (for example, similar machines in a production process). In such cases, it is often possible in simulation to design one basic model 'higher order building block' which can be used repeatedly in the overall model (with different inputs and structural parameter values where necessary) to model the replicated sub—system whenever it occurs. This capability of a language is termed 'higher order modularity' since it involves modularity vis-a-vis a model component constructed from more basic components. 50 Programmability - This criterion relates to the programming effort required to program and operate models in a given language. The concept is somewhat arbitrary since it is often highly dependent upon the partic- ular problem at hand. Nevertheless, for some types of problems, some languages are easier to work with than others. Optimization Capability - Some simulation languages offer the facility of running simulation models linked with optimization routines (such as Complex or Powell's routines) which use search or other techni- ques to optimize the performance of the model with respect to some criterion. Output Formats - Most special purpose simulation languages include 'canned' output formats which provide selected model outputs in tabular, graphical (versus time), histogram or other forms. Such facilities can significantly reduce model programming requirements. The following is a tabulated comparison of a number of contemporary simulation languages. The table was originally compiled by Manetsch and Park (1979) for six languages, based on the six criteria described above, but has been extended by the author to include an additional simulation language 'GASP IV' and an additional descriptive criterion 'Error Diag- nostics'. Table 3.1 Comparison of some contemporary simulation languages w____. 0 3 >. ~4 u o u m a GH 14 U (0 ° a r: .1 s g 2 15 m .4 u >. .o 'H >. no . .3 '3 a: g 2:: .2 :2 no u :5 m 5H NH 9 a o s u a a -H-H u a u #1 m .4 u .§.a : u a: u u > .c a u: a a. o “ 4'3 5 :1. 39'“ 8 ”“ “ t .3 o 8 => =£ a. 8'3 5 n: CSMP X Limited Yes Good Yes Good Good number of large computers DYNAMO X Moderate for Yes Good No Good Good large scale machines FORDYN X Excellent for Yes Fair Yes Fair Fair most machines FORTRAN XL X Excellent for Yes Poor Yes Poor Poor most machines GASP II X Excellent for No Fair Yes Good Good most machines GASP IV X X Excellent for Yes Fair Yes Very Good most machines Good GPSS X Limited number No Good No Good Good of large machines SIMSCRIPT .X X Good for large ? Fair Yes Good Fair machines only FORTRAN & X .1 Excellent for Yes Fair Yes Fair Fair FORDYN 5 to GASP 11 Good i i___L_L_= CHAPTER 4 METHODOLOGY AND MODEL DEVELOPMENT 4.1. Introduction Most agricultural crop harvesting and tranSportation systems can be represented by some form of closed circuit queuing system as shown in Figure 3.1. Essentially, such a system consists of cutting the crop and loading it into transport units which are then taken to a process- ing, storage or other facility where unloading takes place (often after the loaded tranSport units are weighed). Empty transport units are sub— sequently returned to the field. Since a service facility (e.g., har- vester, scale or unloader) is not always free on arrival of a transport unit, service queues often build up within the system. The calculation of the capacity of cyclic queuing systems can be rather complex, especially in an agricultural context, due to the wide variability of operating conditions that may prevail. One relatively simple approach to finding approximate solutions for such systems is through the use of computer simulation programs. In such programs, the variability of operating conditions can be handled by using appropriate statistical distributions to represent the operating times of the dif- ferent pieces of equipment involved in the system (Dumont and Boyce, 1972). This chapter describes the simulation methodology used to study mechanical harvesting and handling of sugar cane under Barbadian condi- tions. 52 53 4.2. Model Development Approach There are two distinct systems of mechanical harvesting of sugar cane in use in Barbados: chopper harvesting and whole-stick harvesting. In the chopper harvesting system, the cane harvested by a combine is delivered directly into a wagon pulled alongside the harvester by a field tractor. When a wagon is filled, it is moved from under the delivery chute of the harvester and an empty wagon is drawn up in its place in such a manner as to ensure continuous harvester operation. In the case of the whole-stick system, the cane is cut by a mechanical cutter and then loaded by a mechanical grab loader (in a separate opera- tion) into a transport trailer pulled alongside the loader in the field. . In either system, filled transport units are deposited in a queue at the edge of the field and, when a complement of two full wagons or trailers and a road tractor becomes available, a trip is diapatched to the factory. The actual flow of individual operations performed on each transport unit during the entire harvesting and handling process is shown in Figure 4.1. Operations above the 'broken line' on the Figure are associated with the field subsystem, while those below the line are involved in the factory subsystem. For the purposes of this study, the overall cane harvesting and handling system was broken down into two subsystems: a field subsystem and a factory subsystem -- a division not at all unrealistic since in Barbados, the cane production and processing sectors are separately owned and managed entities (see Chapter 2). Each subsystem was then decomposed into its Specific activities, and events identified such that each activity was bounded by a pair of sequential time events, one 54 WAGON WAITS TO BE LOADED WAGON IS MOVED TO HARVESTER WAGON IS LOADED WAGON IS MOVED AWAY FROM HARVESTER WAGON WATTS TO BE MOVED TO FACTORY UNIT (2-WAGON TRAIN) IS MOVED TO FACTORY UNIT WAITS TO BE WEIGHED EACH WAGON IS WEIGHED SEPARATELY UNIT IS MOVED TO UNLOADER UNIT WAITS TO BE UNLOADED EACH WAGON IS UNLOADED SEPARATELY UNIT 18 RETURNED TO FIELD Figure 4.1. Flow process chart of operations performed on each transport unit 55 signalling its start and the other signalling its end. The time to perform a given activity was calculated by finding the lapse time between the relevant pair of sequential events. Two computer simulation models were then developed: a field opera- tions model (FIELDOP) to simulate the activities involved in the field subsystem; and a factory yard operations model (FACYARD) to simulate the activities involved in the factory subsystem. Based on a review of the literature on simulation models and computer simulation languages, the GASP IV simulation language was selected for constructing both models. Details of this language are presented in the following section. 4.2.1. An Overview of the GASP IV Simulation Language (Source: Pritsker, 1974) Simulation is a problem solving procedure for defining and analys- ing a model of a system. Digital computer simulation, as defined by Pritsker (1974), is "the establishment of a mathematical-logical model of a system and the experimental manipulation of it on a digital com- puter." Simulation languages typically provide the structure and the term- inology to facilitate the building of simulations and relieve the user of a considerable amount of personal programming effort. GASP IV is such a language. It helps the user to build computer simulation programs that can be both the model of the system and the vehicle for analysing the system. The philosophical basis for the design of GASP IV is the concept of modelling a system in two dimensions: the time dimension and the state-space dimension. Fundamental to building a GASP IV simulation 56 model is the decomposition of time and state space into manageable ele- ments. The decomposition in the time dimension requires the user to define events and potential system changes generated by these events, to specify the causal mechanisms by which events can occur, and to de— fine the mathematical-logical relations that tranSpire when an event occurs. In the state-space dimension, GASP IV presumes that a system model can be decomposed into its entities, which are described by attributes. Attributes may be discrete or continuous. A discrete attribute is one whose value remains constant between event times, while a continuous attribute is one whose value may change between event times according to a prescribed code of dynamic system behaviour.' Continuous attributes are referred to as "state variables". In essence then, GASP IV provides a formalized view of the world which specifies that the status of a system be described in terms of a set of entities and their associated attributes and state variables. In GASP IV, an event occurs at any point in time beyond which the status of a system cannot be projected with certainty. Events are described in terms of the mechanism by which they are scheduled. Those that occur at a specified projected point in time are referred to as "time events", while those that occur when the system reaches a partic- ular state are called "state events". The behavior of a system model is simulated by computing the values of the state variables at small time steps and the values of attributes at event times. GASP IV automatically decomposes the time axis into points at which events occur, based on the equation form of the state 57 variables, the time of the next event and accuracy and output require- ments. The user is, therefore, relieved of the task of sequencing events. When an event occurs, it can change the status of the system in three ways: by altering the values of the state variables or the attri- butes of the entities; by altering relationships that exist among enti- ties or state variables; and/or by changing the number of entities pre- sent. Methods are available in GASP IV for accomplishing each type of change. At each time step, the state variables are evaluated to determine if the conditions prescribing a state event have occurred. If a state event was passed, the step size was too large and is reduced. If a state event occurs, the model status is updated according to the user's state event subroutine. Step size is automatically set so that no time event will occur within a step, by setting the step size so that the time event ends the step. Since time events are scheduled happenings, certain attributes are associated with them. At the minimum, a time event must have attributes that define its time of occurrence and its type. If the time event is associated with an entity, then either the attributes of that entity must be associated with it or the event must be able to refer to the attri— butes of the entity. For example, if there is an end-of-service event for an item, the attributes for that item must,:hnsome way, be associated with the event. A filing system is provided for storing entities and their associated attributes. System monitoring procedures, which can be pre-scheduled or called as required, are also provided. GASP IV is designed to provide eight specific functional capabili- ties: 58 (1) Event control (2) State variable updating using integration, if necessary (3) System state initialization (4) Program monitoring and event reporting (5) Information storage and retrieval (6) System performance data collection (7) Statistical computations and report generation (8) Random deviate generation. The first four of these capabilities are primary functions which constitute the basic modes of the language as shown in Figure 4.2. The remaining four are support oriented. GASP IV has two forms of program control. One, the executive function, directs the system program into its various modes: initializa- tion, state variable updating, monitoring, etc. The other, the event selection function, operates within the simulation model and sequences the execution of event routines. The modular structure of GASP IV allows event logic to be changed relatively simply to investigate the effect of changes in a system on selected measures of system perform- ance, since each event is, in fact, a separate subprogram. A data pool allows changes made in data inputs to be communicated throughout an entire simulation model. The preparation of reports summarizing the results of simulation runs is simplified by utilizing standard report programs that obtain their information from the common data pool. Model debugging is also facilitated by the provision of access points at which program results can be sampled without interfering with the logic of particular events. 59 l l l . 8:53.39...- . ,2 - 7w: .4 .8»... .....é- I E 3 II s l l r 95.933... 62.2.. :35 8:. ...-=8 cos-22. 25326 L a!” Figure 4.2. Basic modes of GASP IV control. 60 In summary, GASP IV has five distinct features that make it partic- ularly attractive as a simulation language: (1) (2) (3) (4) (5) The language is FORTRAN based and requires no separate compiling system. It is easily maintained and can be implemented on new computing systems and on the comput- ing systems of different manufacturers. GASP IV is modular and can be made to fit on all machines that have a FORTRAN IV compiler. GASP IV is easy to learn since the host programming language is usually known and only the simulation concepts need to be mastered. The implementation of these concepts is apparent to the user. GASP IV can be used for discrete, continuous, and com- bined simulation modelling, and is the only well-docu- mented simulation language with this capability (see Table 3.1) GASP IV can be easily modified and extended to meet the needs of particular applications. A functional flowchart of a GASP IV program is presented in Figure 4.3. 4.3 The Field Operations 4.3.1. Objectives The main objective of simulating the field operations was to identify changes in the cane harvesting and transportation procedures that could lead to increased efficiency of utilization of the expensive chopper CM manual ‘ STA!" ) l 61 "with! M UM" IV “on wound. “amino 7 III m cum-6A9 vaults maul MO Swimm- INTLC Immin- non-GASP I Mud-a GASP MUM I canola. «on warm and tun-mu I ii, I /|\ VII M 1 I 00"“ A I u. I .. M a V“ tulle-m I m? 1 No m Im- l m: "on Swot-um “AV! L L Rm m I Mum STATE I than om. vacuum V“ ESuI-thm SCOND ' We non-m Nahu- WI-um u... an * I l l — fl ; l m swan | SUI-u WI. 7 L m I 1 twin (an: . . Em Two | Type 2 Tu. N I i I Saul-twan- 0"!" I Stimu- SUNNY W mod-m 7 but W m a um“ mm 6! M MO” I _—_——i. Mum to START. RESIN". or 510' Figure 4.3. Functional flowchart of a GASP IV program. combine were as l. 62 harvesters and their associated equipment. Specific objectives follows: Determination of the current levels of utilization of the chopper harvesters and of the cane transporta— tion equipment within the field subsystem. Construction of a computer simulation model to accu- rately reflect the Operational characteristics of the field subsystem. . Use of this model to determine the effect on potential subsystem output of variations in the characteristics, specifications and combinations of different subsystem components. . Use of the simulation model to test recommendations which may be implemented to improve and/or optimize the overall performance of the field harvesting and cane transport operations. 4.3.2. Activities and Events The principal activities involved in the field subsystem are the loading of empty cane transport vehicles and the dispatch of full trans- port units to the factory for processing (Figure 4.4). Inevitably, road tractors are not always available whenever a complement of two full wagons is ready for transportation to the factory, neither are harvesters always free when empty transport wagons return to the field. Queues of full and empty wagons, therefore, form an integral part of the field subsystem. 63 I IFROM FACTORY . r--——-——————- wuoeoumm nuwsoa Soumzw pamfiw mnu “Om Emuwowp aOOH Humnmo .H.m ounwwm amozuma mocmumfi. zuouumw -- / +v\\maowwwz . huaaw . . mo wanna H ,5" mafia \ I, \ \ "- ‘ ..-FH\ \\%66 MW \ avoumu>umsx muwumm>u. memo mo .o: muouomuu I + s I _ moon \ 06.3 / muouomuu +, 053 7 mo .0: “ H03 , mamfim I; Omanoyx .. x L mo .0c 1!! - x . 3.33. + / I \ + K / mumwcnom amumam mmanmfinm> msocmwovcm IIII wuoav %HHMp hamwum> mama 79 therefore, be represnted by the same probability distributions. 7. A transport unit always consists of a road haulage tractor and a complement of two wagons. 8. Delays in field operation due to rainfall are infrequent in occurrence and negligible in duration. 9. Transport units and individual wagons are treated on a first-come, first-served basis throughout the system. 10. Agronomic practice dictates that, for a given field, only one cane variety is used. On the basis Of the above assumptions, the daily cane quota, the transport unit capacity, the length of a shift, the length of the break between shifts, and the operational policies of management are all con- sidered to be constants in the system model. 5.2.2. Entities The entities considered for the model are: l. The mechanical (chopper) harvesters 2. The field tractors 3. The cane transport wagons 4. The road tractors The operating personnel associated with the equipment (e.g., tractor operators, the harvester operator) are considered to be integral parts of the component equipment and do not affect the system as individuals. They are, therefore, not treated as separate entities in the model. 80 5.2.3. Endogenous Variables The endogenous variables generated within the system boundary are as follows: Time in the queue of empty wagons Time in the queue of full wagons Time to complete a round-trip Tonnes of cane harvested per day . Available work time. 5.2.4. Output Variables The output variables of interest for the field subsystem are: 1. 2. 10. The average total round-trip time The average time an empty wagon waits at the field to be loaded . The average time a full wagon waits at the field to be transported to the factory The utilization of the harvester The utilization of a field tractor . The utilization of a road tractor The number of trips completed per Shift . The number of trips in transit at the end of a Shift The number of full wagons in the field The total tonnes of cane delivered to the factory from the field subsystem. 81 5.3. Field Operations Simulation Program (FIELDOP) As previously stated, in GASP IV each entity has one or more attri— butes associated with it. The entities and their associated attributes for the FIELDOP program are listed in Table 5.1. Nine files are used in the program for storing information on entities and combinations of entities at different points in time during the simulation. These files, along with the entities stored in each and the attributes used to rank entries in files, are shown in Table 5.2. Throughout the simulation, the identity of each piece of equipment is maintained by assigning it to the appropriate file. For example, if at the end of a shift, a field tractor and empty wagon combination remains in the field, the com- bination is dissembled and the field tractor put in the file for field tractors (File 5) while the empty wagon is placed in the file for empty wagons (File 2). 5.3.1. General Flow Through Program FIELDOP The general flowchart for program FIELDOP is presented in Figure 5.2. The main program assigns appropriate values to the card reader and printer units, reads in the initial values for non—GASP (user-generated) variables (Table 5.3), prints out a list of these variables, and initi- alizes the GASP IV variables declared in the Dimension, Common and Equivalence statements before making a call to the GASP IV executive subroutine, GASP. The GASP IV subroutines DATIN and INTLC (user-written) are then called from subroutine GASP. DATIN enters input data read in by the main program and subroutine INTLC initializes the start, end, and operating work times, daily cane delivery quotas, report arrays and 82 Table 5.1. Entities and their associated attributes (FIELDOP) Attribute Entity Number Description Harvester Field Wagon Road tractor tractor 1 Event time X X X X 2 Event code X X X X 3 Time of arrival in queue of empty wagon X 4 Time of arrival in queue of full wagons X 5 Time of departure from field X X 6 Time of arrival of a transport unit at factory X X 7 Time of departure of a transport unit from factory X X 8 Harvester number X 9 Field tractor number X 10 Road tractor number X 11 Number of full wagons for one trip to factory X 12 Trip completion indicator X X 13 Identification Of field location (i.e., NEST) X X X X 14 Time of entry into file 8 X X 83 Table 5.2. Files used in model FIELDOP File no. (I) Entries KKRNK (I) IINN (I) 1 Future events 1 3 2 Queue of empty wagons at field 3 3 3 Queue of full wagons at field 4 3 4 Queue of road tractors l3 3 5 Queue of field tractors 13 3 6 Harvesters 13 3 7 Information on road tractors processed at the factory 6 1 8 Field tractors and empty wagons coordinates waiting to be loaded 14 1 9 Transport units in transit at end of the morning shift cameo L INIT RDR PRNT GASP ** USER SUBROUTINE f b EVNTS ‘-____‘ 84 COMMON DIMENSION JEQUIVALENCE OTHER SUB-PROGRAMS USED FUNCTIONS GAMA ERLNG RNORM FILEM INTLC ARFLD STLD NDLD ARHTR ARFCT DPFCT SHTDN STRTUP * GASP IV SUBROUTINE (USER WRITTEN) Figure 5.2. Flowchart through FIELDOP SUBROUTINES DATIN TIMST COLCT FILEM RMOVE ERROR DSTART** REPRT** OTPUT* 85 Table 5.3. NON-GASP user input variables used in FIELDOP Symbol Definition ISC System condition: ISC = 1: simulate ISC = 0: end of simulation NWG Number of wagons NHT Number of harvesters NFT Number of field tractors NRT Number of road tractors NEST Identifies estate simulated QUOTA Estate delivery quota of cane for one day CAP Amount of cane delivered by a transport unit per trip (tonnes) STM Time of starting field operations on first day ENT Time of ending field operations on a work day AMOUNT Amount of quota remaining DWN Status of system: DWN = 0: system Operational DWN = 1: system down BREAK One hour lunch break for which harvesting stOps between shifts 86 equipment combinations, and establishes performance criteria for the simulation. A successful return to the executive routine GASP from DATIN is indicated by a print-out of the initial values established for GASP IV variables. The executive routine GASP then calls subroutine EVNTS (IX) which passes control to individual user-written event sub- routines ARFLD, STLD, NDLD, ARHTR, ARFCT, DPFCT, SHTDN and STRTUP, in accordance with the value of the argument IX. Subroutine ARFLD (Figure 5.3) is called when a transport unit, consisting of 2 wagons and a road tractor, arrives at the field, having delivered a load of cane to the factory. The road tractor is placed in File 4 and the empty wagons in File 2. The amount of the day's quota remaining is then calculated, as well as the round-trip time for the delivery. If the quota is satisfied, a shut-down event is scheduled and the road tractor is set idle. If not, the status of the harvesters is checked and a start of loading event is scheduled for each available harvester. If there is a complement of full wagons at the field, it is combined with the road tractor and a trip is diSpatched to the factory. Otherwise, the road tractor is set idle. Subroutine STLD (Figure 5.4) handled the start of a loading event. The harvester is combined with a field tractor and empty wagon combina- tion and set busy. The end of loading is then predicted by generating a deviate from the Gamma distribution of loading time using parameter set #1 (see Tables 5.4 and 5.5). Subroutine NDLD (Figure 5.5) simulates the activities in the field at the end of loading. The wagon just loaded is placed in the queue of full wagons (File 3) and the field tractor in File 5 (see Table 5.2). 87 C m) Is transport unit Yes Reduce quota arriving from by capacity of facfi::Z//””’r transport unit Calculate ‘ No round-trip time J, Park empty wagons lt in queue 2 1 Park road tractor in queue 4 Is system down Yes { RETURN ) Y Schedule shutdown es . of operations for the day quota fulfilled Check no. harvesters ( ) available J RETURN Figure 5.3. Flowchart through subroutine ARFLD 88 Is NO no. field tractors CT. 0 Yes Is Set no. hgrvesters no. field tractors Yes to be used .EQ. GT. ’ no. empty wagons no. empty wagons available - N Set no. harvesters to be used .EQ. no. field tractors available - N Is Yes N.EQ.O No Set up N field tractor and empty wagon combinations and schedule arrival at harvester Are all wagon complements available for dispatch to factory Set road tractors idle, collect stats. 1 C W} Figure 5.3. (continued) Dispatch as many trips to factory as there are road tractors/full wagon combinations available 1 Schedule arrival time at factory for all units dispatched J 89 C 'sm 3 Record identity of_] field tractor Remove a field tractor/empty wagon combination from file 6 collect statistics Predict end of loading for wagon c m) Figure 5.4 Flowchart through subroutine STLD ‘ Set harvester busy and 9O axon any wanOOH mo uumum mnu pan comma H onmeH wo meHu OHOH m wcmHum mouncHz one consume oaHu omemH amumm>umm mcowm3 HHnw some: H a HHnm o wOMHum wocoomm mo mnmnc Ou oEHu Hm>mua HOOOMHO pHmHm Ho>muu pHmHm comma H a woman m wOMHum mpcooom umumm>umn Ou mEHu Hm>mu9 nouomau pHOHm Hm>muu pHOHm Mao: use moowns N a xuaam a Hmsuoz monuOEOHHM pHmHm on women Hm>muH uouomuu omom Ho>muu vmom huOuomm :H msowma N w wsHmeHco m maamo mounst %uOuumm um oaHu musmpHmom HOuomwu wood can wanwHoz anon you muowms N w HHaw N Hassoz mmuuoEOHHM muouomu Ou women Ho>mue wouomuu pmom Hw>muu vmom H assoc mmuost sowm3 awe oaHu wSvaOH umumm>unm onvaH Hanan: mom mama muHsa paumHSOHmO vo>Ho>sH mu > uo umuoamumm OOHunnHuuan maHu >uH>HuO< >uHusm H H ¢ mOQHmHm Hmvoa CH mmaHu auH>HuOm pom pom: maOHunnHuumHv HOOHumHumum .¢.m mHan 91 Table 5.5. Input parameters for distributions used in FIELDOP Paramter Dist. Set A Min. Max. B type Entity l 6.10 0.00 15.46 2.00 GAMMA Loading time per wagon 2 12.74_ 6.44 24.62 5.89 NORMAL Travel speed to factory 3 5.21 0.00 143.40 4.00 GAMMA Factory yard residence time 4 23.27 19.31 32.18 4.02 NORMAL Travel speed to field 5 0.02 0.00 429.00 1.88 ERLANG Travel time to harvester 6 0.03 0.00 310.00 3.37 ERLANG Travel time to full queue 7 0.06 0.00 102.30 1.00 ERLANG Harvester idle time 8 7.89 6.05 18.00 3.20 Minimum observed times1 9 3.00 5.00 7.00 Distance to factory lMinimum Observed times reading from left to right refer to minimum loading time, minimum factory residence time, minimum travel time to and from the harvester, and minimum harvester idle time, respectively. 92 C ....) 1 Record identification of field tractor and hold in queue 5 l 7 Predict time of arrival of full wagon in queue 3 Is no. wagon/field tractor combinations in queue 8 GT. 0 NO Is no. empty wagons in Yes Set field tractor ] s. field * idle, collect stat LT. 1 No Remove empty wagon and schedule its arrival at the harvester 1 Schedule start of loading Is a road tractor and two full wagons available J Dispatch a trip RETURN g,/ I to the factory Figure 5.5. Flowchart through subroutine NDLD 93 Since the loaded wagon spends some time travelling across the field before it reaches the queue of full wagons, its time of arrival in the queue is predicted by sampling the Erlang distribution of the travel time to this queue using parameter set #6. This time is recorded as attribute 4 (see Table 5.1). If there is a field tractor and empty wagon combination waiting in File 8, it is removed and set in motion towards the harvester. Its time of arrival at the harvester is scheduled by sampling the Erlang distribution for travel time to the harvester using parameter set #5. If there is no such combination in File 8, the harvester is set idle. The field tractor just released after depositing the full wagon is then combined with an empty wagon (if there is one available) and a start of loading event is initiated by setting the harvester busy. If a road tractor and a complement of full wagons are available, a trip is dispatched to the factory. Its time of arrival at the factory is scheduled by sampling the Normal distribution for travel speed to the factory and combining the deviate obtained with the known distance to the factory to Obtain the travel time to the factory. Based on numerous field observations, it takes an average of 3 minutes for a trip to be prepared for departure from the field, so 3 minutes are added to the current time TNOW when calculating the arrival time Of a transport unit at the factory. Subroutine ARHTR (Figure 5.6) simulates a field tractor and empty wagon combination approaching the harvester while it is busy loading another wagon. This occurs only if the number of field tractors in the system is greater than one. Therefore, subroutine ARHTR is called only 94 if there are two or more field tractors assigned to a single harvester. In this case, holding a field tractor and empty wagon combination in Field 8 represents a situation where a wagon is being loaded and a second wagon is waiting very close behind so that it can start receiv- ing cane as soon as the first wagon becomes full, without the harvester having to stop cutting. Subroutine ARFCT (Figure 5.7) simulates the arrival of a transport unit at the factory, predicts its residence time by sampling the Gamma distribution of factory residence time using parameter set #3, and schedules the departure of the emptied unit from the factory. Subroutine DPFCT (Figure 5.8) predicts the arrival of an empty transport unit at the field by sampling the Normal distribution of the travel speed to the field using parameter set #4, and combining the deviate Obtained with the known distance between field and factory to arrive at a travel time. The returning transport unit is labelled as having delivered its load (ATRIB (12) = 1.0) and information on it is recorded in File 7. Subroutine SHTDN (Figure 5.9) is used to close out all activities at the end of a shift or a work day. All entities are returned to their respective files and the future events file (File 1) is cleared in preparation for a new start at the commencement of the next shift or next day. Subroutine STRTUP (Figure 5.10) is called to clear all of the GASP IV statistical data collection arrays in preparationiknrthe new shift. If full wagons and road tractors are available at this time, trips are dispatched to the factory in accordance with the number of complete transport units available. C T“) Put field tractor/empty wagon combination into file 8 II Cm) Figure 5.6 Flowchart throxgh sibroutine ARHTR ‘ ARFCT ) Record time of arrival at factory Predict transport unit's residence time at factory Store residence time data for histogram « 1 CW 3 Figure 5.7 Flowchart through subroutine ARFCT 96 <:: Di[CT_:) Record time of departure from factory 1 Label unit as having delivered its load I l # [Schedule the arrival of \ unit back at the field Store information on road tractor processed at the factory in file 7 1 (mm) Figure 5.8 Flowchart through subroutine DPFCT Cm) 1 [Set system status markerW] Is no. empty wagons and field tractors LE. 0? Put empty wagons and field tractors into their respective files Clear future events file ' (file 1) and sort out entities [Set no. future events - N] LIX-17 l_Calculate round-trip time] Dark field tractors in file 21 v Set road tractor idle hold in file 4 and collect statistics CD + 0 Figure 5.9 Flowchart through subroutine SHTDN 98 [NR=NK+1 Yes - Yes LNNQ(9).I.E.oz—Jl N O LN - we)? Remove N entries from file 9 l and store in file 1 quota satisfied? Yes __> end-of-shift time? Figure 5.9 (continued) 99 __.I Schedule start T—schedule start of next shift I of shift next day a Schedule and Calculate amount of next shift quota remaining ‘flv of next shift Schedule and on next day [ Schedule end time Fix time of next start up event Fix time of next shut down event c m) Figure 5.9 (continued) 100 C ...... D Set sustem status marker DWN - 0.0 1 Clear GASP IV data storage arrays 1 set no. harvesters = N6 ] Is 9 N6OEQoOo Yes {all ERROR (SD IMO Is . NO. empty wagons = 0. Yes No from queue 2 1 Set time in queue I 0 and collect stats. 1 [ Remove N6 field tractors L Remove N6 empty wagons ] from queue 5 Figure 5.10 Flowchart through subroutine STRTUP 101 and collect stats. 1 [Schedule end of loading l Set harvester busy no.road tractors .LT.O? Set no. full wagons = N 1 Calculate no. 2-wagon trains available and set no. = K no. 2-wagon trains .EQ.0. Is no. road tractors .LT. no. 2-wagon trains? Call ERROR (6):) Klano. road tractors available [Kl - no. 2-wagon trains available Cb Figure 5.10 (continued) 102 Remove Kl full wagons from queue 3, set time in queue 3 = 0.0 and collect statistics l Remove K1 road tractors from queue 4 Fix time of departure of ] a trip from the field [Schedule time of arrival at factory] Set road tractor busy and.collect stats. no. empty wagons LT. 0. field tractors Set no. empty wagons = K [Set no. field tractors-= K e 1 Set K field tractor and empty wagon combinations busy _ 1 [Schedule arrival of the 1 combinations at harvester Cm?“ 3 Figure 5.10 (continued) 103 Table 5.6. FIELDOP variables monitored by GASP IV subroutine COLCT Variable D i ti Symbol escr p on RTRIP The total trip time of the road tractor from the field to the factory and back. TIQE The time spent by an empty wagon at the field waiting to be loaded. TIQF The time Spent by a full wagon at the field waiting to be taken to the factory. TIQW The time spent by a field tractor and wagon combination waiting to be loaded when more than one field tractor is used per harvester. Table 5.7. FIELDOP variables monitored by GASP IV subroutine TIMST Variable . Symbol Entity Description BUSHT (I) Harvester Utilization of harvester number I defined by attribute 8. BUSFT (J) Field tractor Utilization of field tractor number J defined by attribute 9. BUSRT (K) Road tractor Utilization of road tractor number R defined by attribute 10. 104 5.4. The Factory Yard Simulation Model 5.4.1. The System Environment The system environment for this model consists of the following exogenous variables: Number of cane harvesters in operation at any given time Number of road tractors, field tractors and wagons operating in locations from which the cane is diSpatched to the factory The field conditions at these locations - Road conditions along routes from the field to the factory. All of these factors affect the total number of transport units dis— patched to the factory which, in turn, affects the inter-arrival time of transport units at the factory (see Figure 5.11 for causal loop diagram). 5.4.2. Entities The entities considered for this model are: 1. The weigh scale 2. The chopped cane unloader (crane) 3. The whole cane unloader (hoist) 4. The chopped cane transport units (2 wagons and 1 road tractor) 5. The whole cane transport units (2 trailers and 1 road tractor) As was the case in the field operations model, Operating personnel of transport units are considered to be integral parts of those units and, 105 seesaw was» muoOumm one How Emuwmfip eOOH Hmmnmo :Zwm ouanm mumHos mo Hanan: \‘I i \ I . manna / . umHos kHOuunm o vmsoueeme muHso menu .383 mo .os humvsnom seamen mmHmOm mo Hones: \ oSHu / N HN>fiHHm a . InousH . a s I uHcs uouoem o vacuumepr muHsn ammo commonu mo .0: GOGQHU mo Humane 106 as such, are not modelled individually. In addition, since a transport unit is handled intact throughout the factory yard process, individual wagons or trailers are not treated as entities. 5.4.3. Endogenous Variables The endogenous variables generated within the system boundary are as follows: 1. Time Spent in the weigh queue 2. Time spent by a unit in the chopped-cane unloader queue (crane) 3. Time spent by a unit in the whole cane unloader queue (hoist). 5.4.4. Output Variables These include the endogenous variables listed above and the following: 1. Average time a unit Spends in the factory yard subsystem 2. Total number of units through the system 3. Number of units processed by the scale server 4. Number of units processed by the crane unloader 5. Number of units processed by the hoist unloader 6. Maximum number of units in the scale queue 7. Minimum number of units in the chopped cane unloader queue 8. Maximum number of units in the whole cane unloader queue 9. Utilization of the scale server 107 Table 5.8. List of NON-GASP variables used in FACYARD Symbol Definition XISYS Number in system TISYS Time in system XVW Number of units through yard KJ Local assessment of type of unit KL Local assessment of type of unloader TIQ2 Time in queue for scale TIQ3 Time in queue for crane (chopped cane) TIQ4 Time in queue for hoist (whole cane) BU82* Status of scale BUS3* Status of crane (chopped cane unloader) BUS4* Status of hoist (whole cane unloader) *Values of 0.0 = free; 1.0 = busy for these variables. 108 Table 5.9. Files used in model FACYARD File Number (I) Entities KKRNK (I) IINN (I) 1 Future events 4 1 2 Scale queue 4 3 3 Chopped cane unloader queue 4 3 4 Whole cane unloader queue 4 3 Table 5.10. Statistical distributions used for activity times in model FACYARD Distribution Parameter Random number Event type set stream Inter-arrival Exponential 1 1 time Weigh time Erlang 2 2 Unload time Erlang 3 3 (crane) Unload time Erlang 4 3 (hoist) 109 Table 5.11. Attributes used in model FACYARD Attribute number Description Value 1 Event time 2 Event code 1 - Arrival 2 - End of weighing 3 - End of unloading 3 Time into system 4 Time into weigh queue 5 Type of unit 1 - ChOpped cane 2 - Whole cane 6 Unloader type 1 - Chopped cane 2 - Whole cane 7 Unloader used 1 - Crane 2 — Hoist 8 Time into queue 3 (crane) 9 Time into queue 4 (hoist) 110 10. Utilization of the crane server 11. Utilization of the hoist server. 5.5. Factory Yard Simulation Program (FACYARD) The system consists of an unlimited number Of cane transport units which arrive randomly at the factory where they undergo two services: weighing and unloading. Two types of transport units arrive: chopped- cane wagons and whole-cane trailers. A transport unit consists of two wagons or two trailers and a road tractor and has a nominal carrying capacity of 10 tonnes. Current factory yard equipment consists of one crane for tipping chopped-cane units directly onto the feeder table and one hoist for unloading whole-cane and either stacking it or placing it directly onto the feeder table. 5.5.1. General Flow through Model FACYARD The flowchart for model FACYARD is presented in Figure 5.12. The main program (user written) assigns appropriate values to the card reader and printer units and initializes the variables declared in the DIMENSION, COMMON and EQUIVALENCE statements, before making a call to the GASP IV executive subroutine, GASP. Subroutine DATIN and INTLC are then called from GASP to enter input data and establish initial conditions and performance criteria for the simulation. A successful return from DATIN is indicated by a print out of intermediate results in which the input data and initial values established for the variables are echoed. The executive routine GASP then calls subroutine EVNTS (IX) from which calls are made to the individual user-written event subroutines C m) 111 GASP IV SUB-PROGRAMS USED [ COMMON FUNCTIONS DIMENSION 1 EQUIVALENCE ERLNG INIT RDR PRNT INTLC # DATIN / ARIVAL EVNTS ENDWGH <:_ 'STOP ) '\\\\\\\\\\. ENDULD * GASP IV SUBROUTINE ** USER SUBROUTINE Figure 5.12.Flowchart through FACYARD SUBROUTINES DATIN Tan COLCT FILEM RMOVE OTPUT* ERROR ** ** ** (USER WRITTEN) 112 ARIVAL, ENDWGH or ENDULD, based on the value of the argument IX. In order to start the simulation, one arrival is scheduled for a chopped- cane transport unit at simulation time 0.0 and one for a whole-cane unit at time 0.5 minutes. With IX - 1, subroutine ARIVAL is called to process the arriving transport unit (Figure 5.13). The type of unit just arriving is re- corded and the next arrival of this type is scheduled by generating a random deviate from the Exponential distribution of inter-arrival times using the distribution parameters stored in parameter sets #1 and #2. The status of the scale is then checked. If free, it is set busy and an end-Of-weighing event scheduled by generating a deviate from the Erlang distribution of weigh times using parameter set #3. If the scale is busy, the arrival is put into the queue to the weigh scale (queue 2) and its time of entering the queue is recorded. With IX = 2, subroutine ENDWGH is called to handle transport units at the end of weighing (see Figure 5.14 for the flowchart of this rou- tine). At the end Of weighing, the transport unit just weighed is sent to either the crane or hoist unloader, depending on whether it is a chopped- or whole-cane unit. If the unloader is free, an end-of-unload— ing event is schedule by sampling the Erlang distribution using the appropriate parameter set for the unloader service time (Table 5.14). Otherwise, the weighed tranSport unit is put into the appropriate un- loader queue and its time of entering the queue is recorded. The first unit standing in the scale queue is then removed and an end-of—weighing event scheduled for it. The time spent in the scale queue by this unit is then calculated by a call to the GASP IV subroutine COLCT. ‘ SUBROUTINE ARIVAL > Schedule next arrival Add one to number in system Collect statistics on number in system 113 Yes Is scale busy, No Set scale busy Predict end—of-weighing Set time in queue Collect statistics on time in queue Put arrival in proper queue. Record time in queue I. :<:: RETURN ::) Figure S.13.Flowchart Of the event ”arrival at factory" 114 < SUBROUTINE ENDWGH ‘:) _Chopped cane Determine type Of cane 1 Whole cane I Go to correct unloader I Put arrival in Put arrival in chopped cane queue ole cane queue Set crane busy ’[Go to bring in Set hoist busy Predict end-Of-service‘ next unit for Predict end-of-service Collect statistics weighing Collect statistics ERROR .g_____.Yes Queue2.1t.0.0 Remove lst entry Yes_§ Set scale busy Schedule end-of-wghng Collect statistics ‘Queue2.gt.0.0 Set scale free _-_ Collect statistics RETURN Figure 5.L4.Flowchart of the event ”end-of—weighing" 115 With IX = 3, subroutine ENDULD is called to handle units at the end of unloading (see Figure 5.15 for flowchart). The type of tranSport unit just unloaded is recorded and, since unloading is the last activity performed on the unit, the time the unit spent in the system is calcu- lated by a call to the GASP IV subroutine TIMST and the number of trans- port units through the factory yard process is updated by one. The first unit in the queue to the appropriate unloader is then removed and an end-of-unloading event is scheduled for it. Subroutine COLCT is then called to calculate the time spent by this unit in the unloader queue. The variables monitored by the GASP IV subroutines TIMST and COLCT in program FACYARD are listed in Tables 5.12 and 5.13, respectively, and a computer listing of the program is given in Appendix III. 5.5.2. Input Data for FACYARD As mentioned earlier, the input data for FACYARD consists of the statistical distribution parameters determined for the various activi- ties of the Factory Yard system, using the data analysis program (DATANAL) previously described. The input parameters used in the model are given in Table 5.14. For the Erlang distribution, A represents the scale factor and B the shape factor and the mean of the distribution is given by A x B. The EXPONENTIAL distribution is a special form of the ERLANG distribution in which the shape factor B is always equal to 1, so that A is, in fact, the mean of the distribution. DATANAL not only outputs values for the scale and shape factors of the distribution, but also generates a sorted list of the observed 116 <::_ SUBROUTINE ENDULD V Reduce number in system by 1 Collect statistics Calculate time in system V Calculate number of units through yard _, Yes Yes Is D Yes unit leavin 1\\\4¥fgj>//§I ERROR Ranove first entry Schedule and of unloading Nc Set crane busy Collect stats Set hoist busy Collect stats No V I Calculate time in queue 3 Calculate time‘ in queue 4 l Set crane free ‘ Set hoist free Collect stats I I C > Collect stats Figure 5415.Flowchart of the event "End—of—Unloading" 117 Table 5.12. FACYARD variables monitored by GASP IV subroutine TIMST Variable Descri tion Symbol P XISYS The number of transport units in the factory yard system at any given time. BUSZ The utilization of the weigh scale BUS3. The utilization of the chOpper-cane unloader (crane) BUS4 The utilization of the whole-cane unloader (hoist) 118 Table 5.13. FACYARD variables monitored by GASP IV subroutine COLCT Variable . Symbol Description TISYS The total time spent in the factory yard system by a transport unit TIQ2 The time spent by a transport unit in the scale queue waiting to be weighed TIQ3 The time spent by a chopped-cane transport unit in the queue to the chopped-cane unloader (crane) waiting to be unloaded TIQ4 The time spent by a whole-cane tranSport unit in the queue to the whole-crane unloader (hoist) waiting to be unloaded. 119 nuoooHao amHo: ozaooo oooouoaoez aaeezmzooxm a a.mm o.o o~.m N Ho>aooo oooosoooooeo am .wnHauoz n .2.m ..z.< Ne uouomuu noon «0 nOHumNHHHuD I Nam: He wouumau woos mo SOHueuHHHu: u Hem: ooum>HHov mono mo mossoe u amazzoe moowes HHom mo amnanz u uzmz enemas homes no nooanz I 0302 ustm no one an uHmswuu OH meHaH u HHmMH OQHuu unease Heuoa u mymmw Ne nouomua vHOHm mo SOHumuHHHus u Nam: He souumuu OHOHm mo OOHumuHHHua u Hem: umumm>umn mo sOHumNHHHus u an: ozone sowm3 HHom sH eaHH moHH moose momma human :H oaHH MOHH maHu eHuu mmmum>< I Hmmam OMHnm mo maez u mmmz woman: man I w=I <2 mmmz H 53 22 g“ 120- 5' 2 field tractors S (U 80'- 1 field tractor V p3 H g / 40 l I . I 1 I 6 8 10 12 14 Mean factory residence time = 9.0 minutes 200' 2 field tractors 1'5 o 16 '- E S g 1 field tractor I-I a: c: 0 80' S H .< c: 40 I J I 1 I 6 8 10 12 14 NUMBER OF WAGONS Figure 6.1 Simulated daily cane delivery vs. number of wagons in system 137 residence times. Generally, cane delivery increased as the number of wagons increased and outfits with two field tractors delivered 10 to 28 percent more cane than those with one field tractor. For the shorter factory yard residence time, increasing the number of wagons caused a large increase in the daily cane delivery to the factory. Systems with less than 10 wagons show a disproportionately lower output than systems with 10 or more wagons. Figures 6.2 and 6.3 Show daily cane delivery versus travel distance between the field and the factory. Indications are that: 1. Generally, the amount of cane delivered per day decreased as the travel distance increased, and the decrease be- came more pronounced as the distance got beyond 5 kilo- meters and the number of wagons became less than ten. 2. For a mean factory yard retention time of 9 minutes, the amount of cane delivered by systems with 2 field tractors and 12 or more wagons was unaffected by in- creases in travel distance until a distance of 5 kilo- meters was exceeded. The affect of travel distance on harvester utilization is presented in Figures 6.4 and 6.5. These graphs are somewhat Similar in trend to those obtained when cane delivery was plotted against travel distance. The following observations were made: 1. Harvester utilization decreased as travel distance increased. 2. For the mean transport unit factory residence time of 74 minutes, harvester utilization varied from 52 percent (with a 14 : l : 2 : 2 equipment combin- ation) to 38 percent, with a 8 : 1 : 2 : 2 combina— tion. 3. Where the mean factory yard residence time was 9 minutes, harvester utilization varied from 78 to 71 percent when two field tractors were used, and from 65 to 50 percent when one field tractor was used. 138 14:1 32" \ 12:1: - 160- ° :9. ga . hi .- c: g 8:1: ° 5 80F H E l l J L J 402 3 4 5 6 7 ZOOP E 0— ,. 14:1: : a: E g 10:1: : t: :3 8()" 8:1: : H 3 l l J l J 402 a 4 5 6 7 TRAVEL DISTANCE (KILOMETRES) Figure 6.2 Simulated daily cane delivery vs. travel distance to factory for mean factory residence time of 9.0 minutes N 139 A160- In E E: L— 53120 E 14:1:2:2 a 0 12:1:2:2 g 80- U a 10 13232 8.1:2:2 S 2 l l l l I D " 402 3 4 5 6 7 A120F- VJ a: z z o 5 \ 80- E \ ‘ 14:1:1:2 E. 12:1:1:2 51' 10:1:1:2 0 4o 8:1:1:2 N I— i 0 E 3‘ o 1 l l 1 _J 2 3 4 5 6 7 TRAVEL DISTANCE (KILOMETERS) Figure 6.3 Simulated daily cane deli-very vs. travel distance to factory for mean factory residence time of 73.98 minutes. HARVESTER UTILIZATION (PERCENT) HARVESTER UTILIZATION (PERCENT) Figure 6.4 100 80 60 12:1:2: T“ 10:1:2: 4C) —_‘ :2: l I l I J 202 3 4 5 6 7 80 60 4O 20 140 ‘r-———______________‘~_—__——‘~—_—‘————_ 14:1:1: 12:1: ' \\ ‘_ 10:1: ‘“‘ 8:1‘1 l I I J J 3 4 5 6 7 TRAVEL DISTANCE (KILOMETRES) Harvester utilization vs. travel distance to factory for mean factory residence time of 9.0 minutes l—i t—l O. 0 NM 141 r '2 I‘ll __ E E: 12:1:2 2 g \ :51 ~10122 H :4 H — \ S 8122 M N [—1 U) m > g 1 1 l J A TRAVEL DISTANCE (KILOMETRES) Figure 6.5 Harvester utilization vs. travel distance to factory for mean factory residence time Of 73.98 minutes 142 Road tractor utilization versus travel distance is shown in Table ‘6.6. The graphs reveal that for a mean transport unit factory residence time of 74 minutes, all field equipment combinations resulted in high levels of utilization (90 - 100%) for both road tractors. Conversely, for the smaller factory residence time of 9 minutes, road tractor utilization tended to vary with the amount of cane delivered to the factory. Equipment combinations with one field tractor recorded 20 to 25 percent lower road tractor utilization than those with two field tractors. The high level of road tractor utilization recorded for the factory yard configuration with one weigh scale reflect the fact that, with this set up, a transport unit spent an average of 65 minutes in the scale queue awaiting the weighing service (see Table 6.1) during which time the road tractor, though non-productive, was in effect busy. 6.4 General Discussion The results presented above indicate that an increase in the total number of cane transport units serviced by the factory yard system can be achieved by either extending the daily cane receiving period or adding a second weighing facility to the system. There are a number of implications associated with either policy. 6.4.1 Implications of an Extended WOrk Period An extension of the daily cane receiving period at the factory would necessitate an equivalent extension of the daily work period for field operations. This implies harvesting and transportation of cane 143 Mean factory residence time = 73.98 minutes 100t- /"f.' 7 80- 60P- E F: l J l 1 J g 4O2 3 4 5 6 7 E: 2: C) El.“ 51 S H e: :a a: 8 9’ Mean factory residence time = 9.0 minutes 2 100,. ~ 3 M 80- #—"."—'______,._. 60- ’5 4c) 1 41 I I J 2 3 4 5 6 7 TRAVEL DISTANCE (KILOMETRES) Figure 6.6 Road tractor utilization vs. travel distance to factory HI—IH G>c>haa~ HI—H—II—I NNNN 00...... NNNN 14:1:2:2 12:1:2:2 10:1:2:2 8:1:2:2 144 during night-time hours, a change likely to be strongly opposed by workers and their trade union representatives. If accepted, the extension would bring about an immediate and significant increase in costs to the sugar industry, in the form of 'over-time' wages. Eventually, however, this increase should be offset by a contraction in the overall length of the cane harvesting season. Sugar factories in Barbados normally grind cane continuously through- out the day for Six consecutive days before shutting down for cleaning and maintenance. With a rated grinding capacity of 90 tonnes per hour, Carring- ton Factory requires a minimum of 2,160 tonnes of cane per day to avoid losing milling time due to lack of cane. Given the current factory yard configuration of one weigh scale and an 11 1/2 hour work period, a maximum of 210 transport units can be handled in a normal work day. With a nominal transport unit capacity of 10 tonnes, this works out to a daily cane deli- very of 2,100 tonnes, 60 tonnes short of the minimum requirement, and the factory is likely to lose an average of one hour per day due to lack of cane (see Table 6.1). When the work period is extended to 16 hours a day, however, a total of 289 transport units can be processed at the factory yard (Table 6.3) making an additional 790 tonnes of cane available to the mill and eliminating mill lost time due to "grinding out". If only chopped-cane transport units were received during the work extension period, a minimum of 29 additional transport units would be unloaded at the factory, yielding an increase in daily cane delivery of 290 tonnes, which would also be enough to eliminate mill lost time due to lack of cane. With the current arrival frequencies of cane transport units at the factory and a single weighing facility, extension of the daily work 145 period is unlikely to have any positive effect on the efficiency of utilization of the mechanical harvesting equipment in the field. In fact, a comparison of Tables 6.1 and 6.3 reveals that the mean factory yard residence for transport units actually increased from 73.98 to 82.83 minutes when the work period was extended. 6.4.2 Implications Associated with Two-Scale Configurations The results presented in Tables 6.4 and 6.5 indicate that if either of the two-scale configurations were implemented, the factory yard system would be able to handle all the cane transport units currently dispatched to it within the normal daily work period. Furthermore, since the utili- zation factors for the service facilities (scales and unloaders) were all less than 0.65, it is evident that, with two scales, the factory yard sys- tem would be able to handle considerably more units than its current allocation. The existing factory yard layout provides only 30 metres between the scale and the unloaders for queuing of transport units leaving the scale. This restricts unloader queue lengths to two transport units and often causes a back-up of weighed units on to the scale platform, result- ing in stoppages Of the weighing process and increased traffic congestion in the yard. The current queue positions also encroach on the area to the north end of the feeder table, which should be used for stacking whole cane. With two scales in Operation, this problem is likely to be accent- uated since transport units will be coming out of the weighing process at about twice the current rate. The solution to the problem lies in 146 apprOpriate location of the second scale, together with re-routing of the factory yard traffic and/or relocation of some of the existing fac- tory yard facilities. Two alternative proposals aimed at achieving this are presented in the following section. 6.4.3 Proposals for Re-organization of Carrington Factory Yard Figure 6.7 shows a plan of the current layout of Carrington Factory Yard and the path of cane transport units through the yard. Arriving units travel along the western and southern borders of the general fac- tory area and enter the factory yard at its south end through gateway A-A. They then continue through a narrow section (10 metres) of paved yard on to the weigh scale located near the north end of the work shop building and, on leaving the scale, make a tight S-turn to get under the boom of the crane or hoist unloader. One proposed re-organization of the factory yard is presented in Figure 6.8. Specific recommendations are as follows: 1. Remove the existing scale and relocate it, along with the second scale to form a 2-scale parallel facility, at a point B adjacent to the south border of the factory area as shown. 2. Let chopped-cane transport units, on leaving either scale, proceed along the path shown through gateway A—A and on to the crane unloader. 3. Let whole-cane units, on leaving the scale area, travel behind the office and workshop buildings along the path indicated to the hoist unloader. 4. Let empty transport units to be tared travel around the main factory building back to the scale area as Shown. As is currently done, units to be tared should be given priority at the scales. If adopted, the above proposal would provide more than enough space for queues to be formed at all service facilities and would significantly 147 a mH ou so H ”onom mHmoz I / one» zaouomm nouwsHuumo nwnounu mOHOHno> uuoemcmuu menu «0 name assayed use .mo uSOAMH unouuoo N.o ouanm ~AHOUUWH Mandi/8H mUflgI / \ mflfiuwu How mUHGDI.II:.II l...\\ muHsn msnOImHonz moanms muHsn mnmulvoeeonu mop—363' wsHumunm mann poveoqn.I.I.II 148 E mH On Eu H "OHeum chow zuouomm GOumcHuumu mo nOHu Emoz A||+Il \ / HcmeOIOH you He Humoaoum 0.0 ouanm .33 new mnHwaH manal waHumu wow muHsDI.II.I. muHsp mmmOImHonz noanmB man: annulvoaaonu omanw3.IIIIIII wcHuoumu muHap munmoqn.l III 149 reduce traffic congestion in the central area of the factory yard. Implementation of recommendation #3 above would necessitate construction of a new road passing behind the office and workshop buildings. This should be a paved road 12—15 metres in width. A gravel surface bound with coal tar or other low-cost bituminous material should be adequate. Proposal two, presented in Figure 6.9, is an attempt to avoid having to relocate the existing weigh scale. The recommendations are as follows: 1. Locate the additional scale (scale #2) along the southern border of the factory area as in the previous proposal. 2. Separate arrivals into two queues, one for chopped-cane transport units and one for whole-cane units. 3. Weigh chopped-cane units at scale #2 and allow them to proceed, via gateway AHA, to the crane unloader. 4. Let whole-cane units by-pass scale #2 and proceed to scale #1 via a new road passing behind the office and workshop buildings. These units will now make their final approach to the scale in a direction opposing that of the current approach (i.e., from north to south). If adopted, proposal two would increase the yard area available for stacking whole cane (as would proposal one) and would provide additional queuing Space for chopped-cane transport units. However, the area avail- able for formation of the queue to the whole-cane hoist will remain virtually unchanged and some backing up of weighed vehicles on to scale #1 may still occur, particularly if there is a break in operation of the hoist. 150 EA: oquH Home.» .Couumm couwsHuumu mo COHumNHameOIou wow Ne Hemonoum 0.0 muanm museum waH>me muHODI'II wsHHMu now muH:DI.II.| "OHmum / I \ 3ch oneOIOHOLB nostmfi mann undulwomoono 025.3311 3.552 I wnHuouao muHsn Hooves I I I I. 151 6.5 Economic Analysis 6.5.1. Calculation of Harvestingenachinery Costs The cost of mechanical harvesting of sugar cane can be calculated on the basis of the standard machinery cost accounting methods outlined in the Agricultural Engineers Yearbook (ASAE, 1982). Component machin- ery costs are defined as follows: 1. Fixed costs: Costs which do not depend directly on the amount of machine use, such as depreciation, interest on investment, taxes, insurance and storage, and which are charged regardless of the extent of machine use. 2. Operating costs: Costs which depend directly on the amount of machine use; namely, repairs, maintenance, fuel and oil. 3. Labour costs: Costs which are directly associated with machinery operators and any auxillary operating personnel. 4. Total cost: The sum of the fixed, operating and labour COStS . The total system fixed cost (TSFC) is given by summing the fixed costs of the component machines involved in the harvesting operation, as illustrated in the following equation: TSFC = l/HCDR [(NHT x HTFC) + (NFT x FTFC) + (NRT x RTFC) + (NWG x WGFC)]. . . . . . . . . . . . . . . . . . . . (1) where: TSFC total system fixed cost ($/tonne) NHT number of harvesters in system NFT NRT NWG HCDR HTFC FTFC RTFC WGFC 152 number of field tractors in system number of road tractors in system number of wagons in system hourly cane delivery rate from system (tonnes/hr) harvester fixed cost per hour field tractor fixed cost per hour road tractor fixed cost per hour wagon fixed cost per hour The total system operating cost is given by: where: TSOC -- l/HCDR [(UHT x HTOC) +FTOC (um + UFT2) +RTOC(URT1 +URT2) + (ch xWGOC)]. . . . . . . . . . (2) UHT = utilization factor for the harvester UFTl = utilization factor for field tractor #1 UFTZ = utilization factor for field tractor #2 URTl = utilization factor for road tractor #1 URTZ = utilization factor for road tractor #2 HTOC = harvester Operating cost per hour ($/hr) FTOC = field tractor operating cost per hour ($/hr) RTOC = road tractor operating cost per hour ($/hr) TSOC = total system operating cost ($/tonne) Total system labour cost is given by equation (3) as follows: TSLC - [(NHT x HTLC) + (NFT x FTLC) + (NRT x RTLC)]. . . . . . (3) where: TSLC = total system labour cost (S/tonne) HTLC = harvester labour cost ($/hr) 153 FTLC = field tractor labour cost ($/hr) RTLC = road tractor labour cost ($/hr) The above equations formed the basis of a Fortran program which com- puted the component and total costs for the eight field equipment combin- ations simulated. The cost structure assumed for these computations is presented in Table 6.9. Costs are based on 1982 figures and equipment utilization factors were obtained from the simulation results previously presented in Figures 6.1 through 6.6. Output from the cost program in- cludes the total system cost per tonne and the total amount of cane harvested per year for each equipment combination. In general, the cost per tonne decreases as annual output increases. The details are shown in Tables 6.10 through 6.13. The least cost combination consists of 2 field tractors, l4 wagons and 2 scales at the factory. 6.5.2. Sensitivity of System Cost and Output to Various System Parameters The response of cost per tonne and annual cane delivery to successive increments of two in the number of wagons is shown in Table 6.14. In all cases, incrementing the number of wagons by two resulted in a decrease in the cost per tonne and an increase in the annual cane delivery. In addition, the response to the last increment (from 12 to 14 wagons) was smaller than the responses generated by the first two increments. Gener- ally, however, the response was neither linear nor uniform. Table 6.15 shows the response of cost per tonne of cane harvested to the addition of a second field tractor or weigh scale to the system, and Table 6.16 shows the effect of these changes on annual cane delivery from the field system. The figures indicate that the addition of a second 154 mN.H Hm.HH oo.m n¢.moa .mm\mamou nmme A asap :wuufius Nm 9 uwuaonm a moamuamcH mom qu.m omw.a mum.ma maam> caov couuauz so Nu e ammumuaH oow omn.m coo.o coo.mm . Aummm\mv Ammasmmm ouou mo waam> wwm>ammv Awesome mafia unwfimuumv aowumwompnmn memou nmme ooo.w ooo.~o ooo.wq ooo.omq A.mvm my mofiua compound HmHuHcH 0H m w o Amummmv mMHH magnoma voumaaumm com cos com com Amazonv .mp5 0H so mama mauxuoa s as“: mxwmz ea I commmm wcaumm>umn mo nuwcmq ma «a as Aua\mmuuHHV coauaasmaou Haas moccou m 33 cm . 33 on u:\mm:c0u mm mufiomamo vmumm monHmums mcfinaoo pom bosomwm musuoauum umoo can maoaumoHMfioQO >umcfisomz .m.o manmh 155 on.o mama com um>o vmmuam Nw.qw ® mu: N: All Homnum>o umw>umn umoa mmv\mu: m 9 mzmv N + mmmu em you zmv\uao: H om.o mm.N vowuma mono mop coo um>o powwow un\mm.qw um Hamsum>o uww>umn umoa smc\mun w e mama OH + mmmv om How zmv\mu=on q usonma w.owcm£omz II II II cm.m Aun\mw.~m e “mumm>umn Hon NV Hmccomuma wcwumumno humaaaxa< II Hm.m Hw.m mm.m momma m.u0umummo mHmoo abomo umoo Hmfluaow Noav mocmcmuafima mam ufimamm I- ms.o mm.o m~.~ Aun\mv «puss uma mm.~w a ago II o~.m mm.m wq.oq Aua\.mvm my muufia pom um.ow e doom mHmoo ozHHm¢m mz¢o Q<—~»I-—n1 x A‘f—WNHOO“. A I E ET OF NON- RAW DATA, READ FROM NEXT SET T. 0. AND. NTS. LE. NUMB) GO TO 40 ORE COLUMNS OF RAW DATA TO CHECK, RETURN TO SUBROUTINE ."ESO2 ANB. .NTS. LE. NUMB) GO TO 25 .5 M) I, AAA 2T6§L"PAGE..",I2/IHO, ."NUMBER OF VALUES I ". O, I"VALUES) IN ASCENDING ORDER..."//(IX, 9F8. 3)) :ém— SUBROUTINE FITA c**** c**** C**** C**** c**** 6 200 20] CAAAA 8 202 I"DEGREES OF FR IRINT I”ACTUAL NUMBER. 175 SUBROUTINE FITA)(P, L) DIMENSION 9(5. COMMON /BLKI/N, XIID COMMON /BLK2/VALU IO COMMON /BLK3/NYVALU(O ET-P I,L PL-PZ EPS-P(3‘L) BUNBRYI I;- NYVALU I - THNOII)-o. No-zo ND-ND+I BUNDRYEI )-I IFIx(x BUNDRY NB)- -IFIx( FIND CLASS BOUND SIZE-(BUNDRY(NB) DO'BII I-2, ND 3U30RY(I)-BUNDRY(J)+SIZE I-2 J-J+I IF J. GT. Nv) CD To I XIJ). GE. DUNDRY(I NYVALUIl-I)-NYVALU( GO TDA I-I +I 5F(liGT. .NB) CD To 2 \600 O O. O 0.0 O )I I- GOT I)+ O 5 I GO TO“ gngLUINOI'NYVALU(NO)+NV- J+I FIND THEORETICAL NUMBER IN CLASSES CONST=(PL**ET)/GAMMA(ET) DO 6 I-I, NO J-I+I A-BUNDRYEI ) -EPS E-BUNERYM J -EPS THNOEI -CONST*NV*SIZE*(RINT(A, PL,ET)+RINT(B. PL, ET)+A. 0* CONTICUE .PL, ET) 6.0 HRI ITE(6E ZOO) FORMAT(IHI/IHO, 5X "LOWER BOUND. ", X."UPPER BOUND.”,5X, X, "THEORETICAL UMBER."/) DO 7‘I-I, ,NO J- I+ HRITE(6, ZOI)BUNDRY(é), ,BUNDRY(J).NYV AL (I ) ,THNO(I) EORMATééH ,SX, F72 X, F7.2.I3X,IA,I I6X .F7 2) ON T IN FIND CHI SQUARED CHIZ-O. O IDF-O I- -I.NO DI-THNO I)—NYVALU(I) IFITHNO I .LT I.O) GO TO 8 CHIZ-CHI2+(DI*DI)/THNO(|) IDF-IDF+I CONTINUE IDF-IDF-A NRITE(6 202) IDF. CHIZ FORMAT(IH0/IHOEE06”VALUE DF CNI SQUARED FOR", I3, 2x, O M 3" FIG. 6 IPAGEBIPAGE+I CALL GRAPH(NO) RETURN END FUNCTION RINT 176 c**** c**** FUNCTION RINT(X.PL.ET) timid: IF(X.LE.0.0) 50 TO I RINT-EXP((ET-I.O)*ALOG(X)-PL*X) RETURN I RINT-0.0 RETURN END SUBROUTINE GRAPH 177 Ctttt c**** C**** PLOT ACTUAL AND THEORETICAL HISTOGRAM SUBROUTINE GRAPH(NO) Ctttt INTEGER UNITSée) COMMON /BLKI/ .x 1000).NTS,NACT(IS).IGO.IPAGE COMMON /BLK2/VALU IOO.20).DIFF IOO,2O).UNITS COMMON /BLK /NYVALU(50),BUNORY 50),THNO(50) OIMENSION A INE(80) DATA DOT,BLANK,AYE.TEA/"."." "."A"."T”/ BlG-NYVAL (I) no I-2. 0 IF( YVALU IIiLE.BIG) GO TO IO LE.BIG) GO TO 9 ‘-VCAZ C CONTINUE 9 FACT-BIG/G d ,0 r LINE(L).L-IO.7O.IO) ORETICAL NUMBERS IN CLASS ..... ", O row—— M—o' cz—Izoo—I-Izzwzwzqzzmzz I I zzr+w3—I:c-Iz~\- wbmm d x AOANE(§;$L-II,7I,IO) IKLINE%L .L-II.80) .70 I) I Amsa‘f" . o—Aav— x—x—— ..a—ao d. Va. DMDMM —-‘A-4MV—II-'IM IO ' ININI' U'I-IINI ALU I FACT+II. O(I;/%£CT+II.5)S) HHT- O HHA- O 6 80 LANK v H ——mr| M4)-- (I).(ALINE(L).L-II.80) 7OAI) - x-- BUNDRY(J), (ALINE(L),L-II,80) —mmrn m pmmm II—-rnm z—-— _Ou V—nm. dV—lV—ldO dvdnnxxovd 20k 2 AMP ZWE-Ar- 4W,- 1’ '4 'H “M'- zmowI-OOIH‘I—Or'oowt-l-ofi'fiIII—r’o ' oomomro,10§ Amr'lCOM—r-lr'lmmdr'l- ' mxmt>ont>>o>omt>>o——33>>ogcm£m£>w-H'I oax—— SUBROUTINE SIMPLE 178 c**** c**** C**** SUBROUTINE SIMPLE(NP. NPS. P, VFUN. L) COMMON /BLKI/NV. .XEIOOO). .NTS, NACT(IS). IGO,IPAGE COMMON IBLKZ/VALU I00 20), DIFF IOO ,26),UNITS COMMON /BLK /NYVALU(NO)I BUNDRY 50) ,THNO(SO) REMENSION P 5.9 {NTEGER M 8:0. 0 h mm. NPI-NPS+ NP2-NPS+2 %T.NPS+? |5|TP AL sVALUES OF FUNCTION NPS I VFUN(J)-LIK(P, J. NP. NPS) NRITE(6, I9 FORMAT(//, ' SIMPLEx ITERATIONS - (MAX 300 PRINT EACH IOTH)"/) OO IOI NT-I. OO FIND MAXIMUM AND MIMIMUN VALUE... 537VEUN(I) [S-VEUNU)P DO 2 J- IF(VFUN(J). PLE. FB) 60 TO 3 H-J FB-VFUN(H) GO TO 2 IE}VEUN(J).GE.E5) GO TO 2 Es-VFUN(L) 2 CONTINUE c**** Ckkflt Cicicfdc OQN zo: E: 2' d—>>c: cor- dryer-Pm Ctkfit CAAAA c**** C c**** c**** Ctkkk c**** C**** CAAAA C**** C**** C**** c**** cAAAR I Cthtk Cfikfix 51 c**** c**** 52 CAAAA C**** C**** CAAAA B INAPO. NNAPT, NNATR, ,g S SUBROUTINE ENDULD 202 SUBROUTINE ENDULD COMMON/GCOMI/ATRI NT, MFA, MFE(IOO). MLE(IOO). MSTOP, NCRDR,N E éééTéo)’ NNTRY, NPRNT, PPARH(50. 6). TNOH, TTDEC Y .TIQ2. TIQ3.TIQ6.DU52.BUS3.DUS6.XVH HHOLE CANE OR CHOPPED CANE T E 2TTCLR, TTFIN, TTR 5 I 2CONHON/UCOHI/XIS I Uhdzm. ( N ( I THIS IS FOR UNLOADING NUMBER IN SYSTEM BY ONE - —AZA cnnznnua -Nv*flv CALL TIMST(XI ISYS. TNOW ,I) CALCULATE NUMBER OF VEHICLES THROUGH YARD XVW=XVH+I.O CALCULATE TIME IN SYSTEM AND COLLECT STATS. TISYS-TNOH-ATRIB(?; CALL COLCT(TISYS CHECK TYPE OF UNIT JUST LEAVING UNLOADER éKL.E8.2 GO TO 2 ATRI (; GT. I. 0) GO TO 53 US 3 OF UE UE (NNQ(3 3))E0 .31 3-0. S T RANE FREE AND COLLECT STATS. =0. O L R;IMST(BUS3.TNOH,3) —1 U1 cn—UKD REMOVE IST ENTRY IN Q3 AND SCHEDULE END OF UNLOADING CALL RMOVE(MFE( (g) .g) ATRIB Ig=TNOW+E LN (A, 3) ATRIB 2 '30 CALL FILEM(I) SET CRANE BUSY AND COLLECT STATS. GALE TIHST(DUS3. TNDN. 3) CALCULATE TIME IN QUEUE.3 TIQ =TNOH- ATRI 8(8 CAL COLCT(TIQ3. 3 RETURN 2 GO TO 53 c**** c**** 53 CARAA c**** EA CAAAA CAAAA 55 C C**** C**** Cflkflz': Cflkfif: CHECK STATUS OF UEUE 6 IF(NNQ(6))50. 56. 5 ESSA-00 SET HOIST FREE AND COLLECT STATS. LREIMST(BUSA, ,TNOW, A) U LQA. GT.O REMOVE IST ENTRY IN QA AND SCHEDULE END-OF-UNLOADING LL RHOVE( (MFE(A). A) I RIBéI ;:TNOH+ERLN G(5, 3) ATRIB 2 CALL FILEM(I) BUSAHDIST BUSY AND COLLECT STATS. CALL TIMST(BUSA, TNOW, A) ALCULATE TIHE IN QUEUE 6 6=TNOH~ -ATRID( (a; CA TI IQ CALL COLCT(TIQA. RETURN SUBROUTINE ENDULD 203 0 CALL ERROR 5 STOP (93) C END SUBROUT I NE OTPUT 206 **** Cittkt SUBROUTINE OTPUT COHHON/OCOHI/ATRID(25).JEVNT,HEA, MFE(IOO), MLE(IOO), MSTOP, NCRDR, INAPO, NNAPT, NNATR, ,NNFIL.NNg(IOO). .NNTRY, NPRNT, PPARM(50, 6), TNON. TTBEC 2, TTCLR, TTEIN. TTRIB(25).TT ET c COHHON7UCOHI7AISYS,TISYS,TI02,TIQ3.TIQ6,DU52.BUS3.DUS6,xVH WRITE(NPRNT§I 00) PPARM(I,I) PPARM(2,I) ’ ANOHT - PPAR M( I)*PPARM( LL AULDTC - PPARM A.I;*PP ARM L. ; AULDTH - PPARM ?,I *PPARM E,A NRITE(NPRNT, IOI AULDTC, Au DTH.xVN c RETURN IOO FORMAT(/I§X, ,"HEAN INTER- ARRIVAL TIME FOR CHOPPED- CANE VEHICLES . " gig ’MEAN INTER- ARRIVAL TIME FOR NHOLE- CANE VEHICLES - '. IOI F RMAT(/I§X, ,“MEAN UNLOADING TIME FOR CHOPPED- CANE VEHICLES - " I I; "MEAN UNLOADING TIME FOR NHOLE- CANE VEHICLES = “. g??? 2{, ,Is "NUMBER OF CANE TRANSPORT UNITS THROUGH YARD - ", END 205 mvump mO+wOOOh. szuhh .0 uwwmhh O Noaomn w Hawmbfi w0+wOO—n. pO+wOOmv. «0+wCva. OO+wOONm. pO+wCOmu. p0+woovo. «0+wOONN. .O‘wOOmp. n0+w0n0'. .O+w00mm. 00+w0050. OO+mOOvn. 90+w000'. N0+w0mnm. .O u0+wOONm. .O+wOCC.. N0+m0wvm. .0 vO+wOO—h. n \ v .0 I U u .0 u 0 ~ .0 n .0.— .O u .0.— O r044u2 O umOmZZ O NoOwZZ 0009 thmzz Om. n>mh22 n Nahmzz O uhaazz m Nzuazz O Nmulzz thDfinN ZO~muw> >~ am<0 OOOOOOOOOOOOOOOuaDmJJ w no . awmlbz 23a «mm. \m \h uhw>¢¢I >m CO. uwm132.h0wfioma Zo~h<423~m nmpnv w .3407? IIn HI 0 FNOVID n 0 Q Q vamuhmdad anmuhm<44 «mam-hm<44 m>m~xthm~huom<44 v Iauuzz v udhmzz wmnvv Hamma— ' w hwm hum hum bum hum A. Ac “N0? 'NC‘)‘ Q05 Nauudfi UQDhmZ auhw2ba2w mu w4~u ml» .0 Nzuhoo .0 N 302k v awmzaz wank no bDDhZHGQ >pa2m mu wank th .0 "Inhoo .0 u 302» n mwmzzz WJHu no FDOhZHQQ >ha£w mu wank th .0 uzuhoo .0 N 302p N ameDZ mama no haohzmua .0 .0 .0+m000N. p0+w000N. .0 00+w000m. p0+w000.. 00+w000m. N N >ah2w .0 .0 .0+w000.. w0+w000.. .0 .0 p0+m000w. .0 N . >mh2w ththOU wquu .0 Nimhco .0 u 3OZh w awmznz w4~u no h00h2~ma N N IUDoth mh—ZD haOQMZQHh wZCU m0 uwmzzz mO.N mw40~1u> w2 wZ deUINJOIJ not wish 4<>~mm mZ~auduawh2~ Z¢w2 9cmh43mwa wh<~0ut¢uh2~¢o 208 QN. N.N 0.N mmo .O+mooo,. .o .C+uooo.. n0+w00mm. m04<> .aZU «0+wv.0N. —0+w000v. n0+wmn.w. 00+w10'v. thuxdl m w4~u 2— aumEDZ SDI—Kc: anm. 20~h<—>u0 ou<02< n0+u0005. m0~hm—h. u 302» w mwmtaz wa—u no h00h2~ma an N dud< woaubz~ mxuh tat~xwo ohm ocwwam<~u<> FZuhmnmuua-w:~» ecu wo-m~p. .o 00+w.w>n. .O+m.mb.. m0+umvmw. .O+wnmoo. 00+uvan. —O+mhmh.. ~0+mpme. 132—Zn: >0 zwo ohm ..20~h<>&wmmo 20 wa mom mumhmnhapmoc —0+m009n. .0+m000v. ¢0+w000¢. 00+m00Nm. -0+w00mN. w0+w00wm. p0+w00NN. —0+w00m.. N0+wCNOp. .0+w00mm. 00+w0050. 00.w00vn. .0+w000'. N0+w0mmm. .0 p0+w00Nm. p0+u000.. N0+w00vm. .0 «0+m00.h. 00+w0005. u wank wamuDU w no — awmtbz 20a Nam. \m \h whw>m0 00. ummzbz homfioaa zouhdaazum ..haoama >a<523m amdo.a UN 'Nflvm OO+un00m 00+wthv 00+umwmm N0+u.mnu zm~x vouh non» N0~h m>m~h awhwidudn awpw2a»Zw 0 0 00+u000N. ¢0+w000N. p0+w000N. n0+wmmm0. 00+wmmmw. p0+w000p. n0+w0hww. N N. >uh2w .0 .0 .0+w000.. p0+w000.. .C+w000—. 00+wnmmw. n0+wmmmo. ~0+m000.. n0+w00¢®. U u. >ah2w .0 .0 p0+w000N. '0+w000N. .0+u000N. n0+umavw. n0+wmmvw. ¢0+w000.. 00+wmmmw. N 0. >uh2u .0 .0 w0+m000.. .0+w000w. .0+w000.. nO+mpmvw. n0+wvmvw. .0+w000.. n0+wnmmw. n m >ah2m .0 .0 p0+u000N. .0+w000N. .09w000N. 00+wmhv0. n0+umhvw. p0+w000.. 00+wmmvw. U m >c»2w .0 .0 .04w000w. .0+u000.. p0+w000p. 90+wmocw. n0+umocm. .0+w000.. 00+w.mvw. N b >uh2w .0 .0 .0+w000w. .0+w000.. .0+w000.. 00+whwv0. 00+mbwvw. u0+w000.. 00+mmwfiw. u w >mh2w .0 .0 .0+w000N. ¢0+w000N. .va000N. n0+wmnv0. 00+wwovw. v0+w000p. n0+mmhvw. N m >mh2u .0 .0 p0+w000w. ¢0+w000p. u0+w000.. n0+wh¢n0. 00+whmnw. .0+w000¢. n0+whwvw. N v >ah2w .0 .0 .0+w000N. .0+w000N. w0+w000N. 00+wchnw. 00+wv>nw. .0+w000—. 60+wmnvw. N n >uh2w .0 .0 .0+u000N. .0+w000N. ¢0¢w000N. n0+wthw. n0+wthw. .0+w000.. n0+wv>nw. u N >ck2w .0 .0 «0+w000N. .0+w000N. .0+w000N. n0+wm>Nw. 00+wthw. .0+m000—. n0+mwth. a _ >a»2w mFZNhZOU w4~m mm wauu Zn awmzaz SDI~Xwo oa< 00+u0005. muuhm~hthw cog Doumwa wimh n0+mcmom. NS—hoo 60+w0005. n 302» N cwmzaz mama no h30h2~ma .0 .0 w0+w000N. .0+w000N. p0+w000N. .0 m0+wrn0h. .0+u000.. 00+w—n05. n v >mh2m .0 .0 '0‘w000N. .0+w000N. .0+w000N. 00+uwaD. no+wthQ. .0+w000N. n0+me05. u n >awa .0 .0 p0+w000N. .0+w000N. u0+w000N. 00+wmme. 00+wmme. u0+u000n. 00+me05. U N >ab2w .0 .0 w0+w000w. .0+w000¢. u0+w000w. .0 00+uvp05. .0+w000p. n0+w9.0h. a . >¢h2w ththOU wguu 210 .O+wooow. .O+mooo.. .O+woooN. ¢O+wOOON. .O+uoooN. .O+woooN. .O+wooo.. 'O+woooN. ¢O+moooN. 00+uoooN. .O+wooo.. .O+wooo.. 'O+wOOON. .O+wooo.. .ovuooop. —O+wooo.. .O+woooN. .O+woooN. .O+wOOON. .O+u000N. .0+w000p. —O+woooN. .O+wOOON. .O+m000w. .O+mooo.. —0+w000w. .04wOOON. .O+wooo.. >hazw mu m4~u th N wJ—u Zn mmeDZ 232~x<2 050.. 20~h4~>wo DfldDZdhm .mnO. wJ—u Zu aumIDZ wo4 00+w000h. mouhm~howuo v0+w000.” n0+wmmww. n0+wmmwwHM .O+uooo.n ooooooooooooood .O+wooo.. n0+wv.0h. .O+wooo.. n0+woomm. —O+uoooN. 00+w.moh. .O+woooN. nO+wmmmm. .O+uoooN. nO+wmhmw. .O+woooN. n0+u.vmm. .04wooo—. n0+u.0mm. w0+woooN. n0+un.mo. .O+moooN. n0+wmmb0. .O+moooN. mO+uvbhw. .O+mooo.. nO+wombo. .O+wooo.. n0+unm>m. .O+woooN. n0+m>mbo. .O+wooo.. nO+wh.hw. .O+wooo.. ll QN mN NN wN 0N mv m. by 0— m. vw >ub2m >¢h2w >mh2w >uh2w >uh2w >ab2w >ah2w >mh2w >mh2w >mh2w >mh2w >mh2m >mh2w >mh2w 211 mo; « a... . a m5 m8 «mm mm. a 2 003.000» . Non »( »m00 Dummmooua mw2~4 4<»0» 040a muzua »Z~mn wa—a »2—un w»wwIm >»a2m mu wquu wI» N0\»N\00 N0\'N\mo N0\.N\00 N0\.N\00 N0\»N\00 w4~u 2m aw0202 EDS~Xw0 0m<02<»m wduu 2m awmzaz w043w>< m0~»m~»<»m «Om oouuwa wtu» 00+w»o»0. u!~»00 00+0000». N 302» Cwmznz mama no »00»2~¢u .0..»0.»». .0».»0.»». .0..»0.».. .0..»0.»u. .0»..0.w.. Results for Extended work period Simulation time = 960 minutes 212 00.00N a 0¢(> 200001» m»~20 »10nm2 w2<0Iw4013 co; wZH» 02~040420 2402 05.0 a mw40~1w> w2¢0I00au010 cog min» 02—01042: 21w! 0N.m - mw40~Im> deuquOII 105 min» 4(>~au(Imm»2~ 2(w! 0'.» I mmJU—Iw> NZCUIOmAQOIU 80$ 02—» 4(>~¢a¢Ium»2~ 24w: .om»40mw¢ m»<—0w!um»2~cc 213 0». N». .mN moN W00 .0 »0+0000». -0+0000.. N0+00000. w34¢> .uao .O+wv.wu. .O+wwva. n0+uuon.. n0+w»nc.. 102~X¢2 n0+00000. 00+00000. 00+00000. 00+00000. J<>uw»2~ wtu» Iemwa0(~¢ .0 .0 .0 v0+00000. t:x~2~x «Izo~»<>uwm .O+woo.n. .O+uoomn. uO+uOno.. .O+wooo.. .O+uooo.. » to — >w>¢(1 .0+0000.. .0 00+wnmmv. .O+wooo.. .0 00+mwmmv. .O+wooo.. .o .OIw.»mv. «0+w00mc. .o .O+mmnoc. 208~K00 0»m <> »2w»m_mawd-wx~» can mu~»m~»<»moo .o+wmv.u. .OIwmovv. 00+wmmmm. .0+wmn»N. .OIwVOOh. 00+0N.Q». OO+mmn»n. .O+wv»o.. «0+wmmmu. OO+mv0vo. .O+ummo.. N0+w0now. >0 24w: no om >wo o»m . no 20 cum¢m mmam<~c<> can mo~»m~»<»ma. .O+w00ov. .0+wooQ.. OO+mOONm. .0+000.0. »0+000NN. .0+0000.. .O+w00mo. OO+woo»o. OO+moo'n. «0+u0mnm. .o .O+woon. «0+womvm. .o .O+uoo.». 00+00000. I win» »20¢¢00 ¢w0202 20¢ N00. \0 \» w»(0 >0 00, 800202 »0000ca 20~»1402~m Io»coamm >cm~x V0~» nO~» Non» m>m~» aw»w!(uOID— IE SCALE IS BUSY, PUT ARRIVAL IN 02 ANO RECORO TIME IN; ATRIOIA)-TN0N ‘ CALL EILEH(2) RETURN IE WHOLE CANE UNIT. GO TO SCALE #2 IF(BUSS) 111,107,108 CAL TIHST(OUS ,TNOH 5) ATRIBEI g-TNOW+ RLNG(S, 2) ATRIB 2 CALL EILEH( SET TIME IN 05 ANO COLLECT STATS. “‘2'"? gALURN COLCT(TIQS. 5) IE SCALE 2 IS BUSY, PUT ARRIVAL IN 05, RECORO TIHE IN. ATRIB(A)-TNON CALL EILEH(5) RETURN CALL ERR0R(87) STOP ENO Ctttt Ctktt c**** Ct ”A Ctttt cAAAA CRAAA C**** c**** c**** c**** c**** c**** c**** I c**** c**** c**** c**** I02 Ctttt c**** c**** Ctttt C**** c**** c**** I03 Ctttt c**** c**** SUBROUTINE ENDWGH 219 SUBROUTINE To NANOLE VEHICLES AT END OF NEICHINC SUBROUTINE ENONCN COHMON/GCOMI/ATRIB(2?),JEVNT.8 NEA ,NEE(IOO), NLEIIOO) NSTOP+ NCROR,N INAPO NNAPT NNATR NNE L.NN2(IOO ), NNTRY, NPRNT, PPARN (5L A). TNOL TTBEG 2. TTCLR, TTEIN ,TTRIB(25).TT ET I§8°M0N7UCOHI7XISYS.TISYS.TIQ2.TIQ3.TIQN.TIQ5,BUSZ.BUS3,BUSN.BUSS. EHICLE FOR WEIGHING g WHOLE CANE OR CHOPPED CANE EHICLE TO CHOPPED CANE CRANE 0R WHOLE CANE HOIST IE(KL.EQ.2) GO TO 2 CHECK CHOPPED CANE CRANE IF(BUSB.GT.0.0) GO TO I03 IF CRANE NOT BUSY. SET CRANE BUSY AND PREDICT END OF UNLOADI NG ATRIB I-TNON+ERLNC(A. 3) ATRIBZ ATRIB7 CALL EILEN(o ) SET CRANE BUSY ANO COLLECT STATS. CALE.l TIHST(BU53, TNOH, 3) TIQ WE- CAL COLCT(TIQB 3) GO TO BRING IN NExT CHOPPED CANE UNIT EOR VEICNINC CO TO 20) ‘ As BUS3-L O. PUT ARRIVAL IN CHOPPED CANE QUEUE ATRl8283-TNOW ATRIB CALL oEILEN(3) GOT CHECK HHOLE CANE HOIST 2 IF(BUSN.GT.0.0) GO TO 105 c**** CAAAA C**** c**** ION CAAAA c**** c**** C**** c**** CARAA CARAA I05 IF HOIST NOT BUSY. SET HOIST TO BUSY AND PREDICT END OF UNLOADING ATRIB I -TNOH+ERLNG . ATRIBIZI'M (5 3) ATRIB - CALL NILEH(I) 35S HOIST BUSY AND COLLECT STATS. h- I. 0 CALL. TIHST(BUSN, TNOW, A) -o. 0 CALL COLCT(TIQA, A) GO TO BRING IN NEXT WHOLE CANE UNIT FOR WEIGHING GO TO 202 AS BUSh-I. 0. PUT ARRIVAL IN WHOLE CANE QUEUE ATRIBE;;:TNOW ATRIB CALL oEILEN(A) GOT c**** 201 cA*A* c**** 1138 c**** c**** Ctttt Ctttt 11h c**** C**** c**** C 202 c**** 116 C c**** 117 nn 112 SUBROUTINE ENDWGH 220 IF(NNQ(2)) 112,113,11h IF US92-06 SET SCALE FREE ANO COLLECT STATS. CALL TIMST(BUSZ, TNON. 2) RETURN S IN OF N ( FE(2). 2) oN8N+ERLNO(3. .2) 2, REHOVE IST ENTRY AND I ND IGHING H EH (I) E BUSY ANO COLLECT STATS. s T(BUSZ TNON 2) E TINEB IN OUEU ON-ATR (h; LCT(T|Q2, 2 IF(NNQ(§)) II2, II6 II EHPTY, SET 5 ALE 2 FREE ANO COLLECT STATS. aus -O. O CAL RETUR EEHST(BUSS.TNON.5) E CANE UNITS ARE IN Q5. SENO ONE TO SCALE #2 OVE(NEE w(3)' NE) g:TNO WE LN (6. 2) LEH(I) HST(BUS E I NH L I IS 20 L IF L CAL H ATR I ATR 2 gAL I CAL ,TNON, 5) TI H- ATR (h 0 a? F T _T C CT(T|Q5.5 5 -IE N 83° cETURN CALL ERROR(92) STOP ENO Results for Configuration with l queue and 2 scales 221 OO.wnN N ou¢> IUDORI» mh~22 baoam2 w2 w2 w2«UIm4013 no; w2~h 4<>~aadnuuh2— z uZ~uu .ano uO+wmmo.. u0+moo.—. .O+wmv.m. n0+wun.u. 838~x<8 n0+uoooh. n0+uooob. nO+mOOOb. n0+woooh. n0+woooh. 4(>¢mh2~ mluh —O+wwmon. .O+wooo.. .o 00+thnv. .O+wooo.. .o 00+umvmv. .O+wooo.. .o oo.m.mmc. .O+wooo.. .o 00+ww>mv. .O+mooom. .o .0+wowm.. SDI—x4: z::~z~: >uo chm camwum<~u<> FZubmumcua-m:_p no; munpmnpmm.. .O+wmmnu. .O+wohh.. 00+w~mmn. .O+whmon. .O+uomm.. 00+w>.o.. .O+wnsm.. 00+u.nmn. 00+un0nw. .O+wonmn. >0 Zau: no om ‘ >uo ohm 838~2~8 oozonh<>awmmo 20 mecm mw4m<~u<> can monbmuhdhch uO+mO~o.. .O+woo.n. .O+w00mn. «0+w0uo.. .O+wooo.. .O+wooo.. — no w >w>u¢1 —O+w00mm. OO+wOOhm. OO+wOO'n. .O+w000v. .O+wOOm.. OO+wOONm. .O+woovm. ¢O+mOONN. .O+wOOm—. .O+wOOmm. OO+wOOhO. OO+wOOvn. NO+wOmnm. .O .O+wOONm. NO+wowvm. .O —O+mOO.b. nO+wOOOb. A ulnh Fanunu uwmiaz 23a «mm. \m \h whdo >m DO. uwmtaz howfioun 20~h<432~m cohuoawu >a¢823w athcc "Nfl'lnm OO+whwvw 00+w0n>m OO+u.ohv OO+wmmem .O+mn>0n zm~x v0~b non» «Duh m>m~h mwhwzdcda awhw2¢an 33:" CT A RIVAL TIME AT FACTORY (I Fl ATLB {- ATR;D(5)+(PPARM(9. .NNRUN)/RNORM(2. 2))*60. O N c**** pR C c E**** SHUT DOWN OPERATIONS EOR TODAY 25 ATRIBEI}-TNOW+I0.0 ATRIB 2 - .0 CALL FILE (n c NIO-INT(ATRIB(IO)) E**** SET ROAD TRACTOR IDLE BUSRT(NIO)-0 CALL TIMST(BUSRT(NIO). TNOW. (A+NIO)) GO TO 302 O -o §OIL CA fiRROREI 3 0 CA ERROR 2 02 fig SUBROUTINE STLD c**** 230 SUBROUTINE STLD nnnnnnnnn Ctttt C c**** c**** Cttfit THIs SUBROUTINE SIMULATEs THE START OE LOADING I THE IDENTITY OE THE EIELD TRACTOR IS RECORDED As N I THE HARVESTER Is COMBINED HITH THE ENTITIEs INVOLv D THE LOADING PROCESS 3) THE END OE LOADING Is PREDICTED COMMON /GCOMI/ ATRIB(25).J JEVNT MEA MEE(IOO), MLE(IOO).M$TOE, NCRDRN INAPO. NNAPT. ,NNATR,NNFIL.NN2ET (IOOf, NNTRY, NPRNT, ’EEARM(5O.A). TNOH, TTBEG 2TTCLR, TTEIN WRIB(23).TT% 2COMMON’ /uCOMi/ Isc .2 OTA AH ENT ENTM BREAK STTM NRT, NHT NET,NNG. ICAP. BUSHT(Z Bu SET (). Bué RT h; .TIQE. TIQE, RTRIP, éHEND .NNTM RECORD IDEN 3T9 YOF EIE ELD TA N9-ATRIB(9 HARVESTER IS COMBIN CALL RHOVE(HFE(6).6 :é- N((ATR?B(8)) BAL MST(BUSHT(N8). TNON. N8 T 4.8 ) CT END OE LOADING BY SAMPLING THE GAMMA DISTRIBUTION T3 LOAD TI IME E N em D WITH THE ENTITIES IN THE LOADING PROCESS gg-TNOH+GAMA(I.1)+PPARH(8,I) ILER(I) mxn>>ovng 2m>—i 4*!” U-fl‘fl” mt- cr——-!Dr-I—: SUBROUTINE NDLD cAAAA Csnsasusaxx C nan C c**** 5 N2-INT(ATRIB) (9) CAAAA c**** C c**** c**** IO c**** CAAAA C CAAAA CAAAA '5 C c**** ca*** 25 26 C CAAAA c**** 30 CAAAA INAPO.NNAPT, N 2 TTCLR.TTFIN II IbNNDAY .NNPT 231 SUBROUTINE NDLD THIS SUBROUTINE SIMULATES THE ACTIVITIES AT THE FIELD AT THE END OF LOADING COMMON /GCOM ,MFA, MFE(IOO). MLE(IOO). MSTOP, NCRDR, N T O),NNTRY, NPRNT, PPARM(50, A), TNOW, TTBEG N T.ENTM, BREAK, STTM NRT, NHT, NET. NNG. IEE. TISF RTRIP SHEND,H B LR, MMNIT, MMON, ,NNAME(%)$ NNOE NRNS, NNRUN. NNSTR, NNYR SSEED6 133?‘ C V“ cm-‘m—wz’) I N I 9 9T COMMON /GCOM5 4‘ CMMON IUCOMi ' RECORD IDENTI A'fl \Z\W\ I AT CLL EILEM PREDICT TIM OE ARRIVAL OE EULL HAGON IN QUEUE OE EULL NAGONS (EILE 3) ATRIB(A)-TNOH+ERLNG(6. I)/6O. O+PPARH(8. 3)/6O. O A-ATRI BIA) CALL EILEM(3) CHECK NUMBER OF WAGON AND FIELD TRACTOR COMBINATIONS WAITING FOR LOADING. IF NUMBER.GT.O SCHEDULE STLD FOR ONE IF(NN8(8)) go CALL MOVE( E TIQN-TNOH-ATR CALL COLCT(TI ATRIBEI g-TNOW ATRIB 2 CALL EILEM(I CHECK NUMBER) OE EMPTY WAGONS AT EIELD AND SET UP ARRIVAL OE ONE AT HARVESTER IE NUMBER.GT.O IF(NN (2). E 0 TO 30 CAL MOVE( 2) g%7.1)+PPARM(8,h) THE CYCLE TIME IS NOT LESS THAN HE QUEUE OE EULL NAGONS ;+ GO TO IS +A-TNOW (2)) MOVE( 7 I)+PPARM(3.A) O 26 -TNOW IF NUMBER OF EMPTY WAGONS AT FIELD. EQ. 0, SET FIELD TRACTOR IDLE AND COLLECT STATS. BUSFT(N3M CALL TI ST(BUSFT(N§), ,TNOH, (2+N )) Rgco RD ID OE HARVE TER AND HOL IN EILE 6 =INT(ATRIB(8)) SUBROUTINE NOLO 201 c:*** C**** 202 20 ctttz CAAAA c**** c**** 232 BUSHT(NB) mo CALL TIMST BUSHT(NB), TNow, N8) CALL EILEM 6 NII-INT(ATRI (II)) IF THERE IS A COMPLEMENT OF FULL WAGONS AND A ROAD TRACTOR AVAILABLE. DISPATCH A TRIP TO THE FACTORY IF§NNS(3).GE.NII.AND.NNQ(A).GT.O) GO TO 202 Do EOEL TIRMOVEIMEEQMD L COLCT(TIQF, 3) ROAD TRACTOR EROM QUEUE OE ROAO TRACTORS, Qh VE(MFE h) h OE DEPARTURE EROM EIELO IME O ARRIVAL AT FACTOR ATSIB(5)+(PPARM(9. .NNRUNI/RNORM(2. 2))*60. O .... 2 23m 4. (I) TRIB 0)) TRAC OR BUSY ANO COLLECT STATS, Iu#22 v Na#m22 0 nbzafinn 20—mmw> >~ amdu # >w>m41 .3 u2~u## .0 «.va n0.nv n000## 0 ”Dav—u . .0+00000. .040000m. N0+w000.. #0+w0m00. n0+unn0.. .0 n0.w00,n. .0 n0+womnv. .0 n0+mm.wn. m0+w.nm#. m0+uvn... .0 N0+wn0vn. .0+w0vvw. N0+umvm.. .0 0 — . n n 0 v. 0 . n— n— u304II m? u4wUZZ .0 .0 .0 .0 .0 .0 .0 .0 .0 .0 um0w22 0 "00022 0000’ u#4n22 m NSQQZZ # awmxbz 20“ N0 «000?? .0+uOOOn. .O¢w0#m#. n .OIwoowm. u .OIwownn. u .OIuohmp. u ~O+w#wnw.. u .O+wwONm. u n0+wvbnw. u .0+woo'o. n n 0' UCJUCJUCJUCJUCJ H~~~~~H~~~ u#mm22 um—IZZ .o..00000000000ma:m44 \m. \ow mpqo >0 00. ameDZ #0010ua 20~#<402~m _0nvv # nu4011 "NNQLDLDFOO’? n n v m #0u.#auIm<44 .m.»am:mu»m<44 .mvkamamupm444 .v.»um:mu»m<44 .n.#aw:mupm<44 An.»um:mupm<44 ...#um:mu»m<44 .u.pum:muhm<44 ..vpumamu#m<44 .vaxmamupm444 ..vpzmnmu»m¢44 )0—#M00<44 u0~#u00<44 w0~#u00<44 a~a#muum<44 0 N4~u22 0w u<#m22 .mpnm noum__ e ”4:2. # MQO#mI #um aw#w£<&