The Model of Vehicle and Route Selection for Energy Saving

: The World Bank, United Nations, the Organization for Economic Cooperation and Development, and others are in line with the governments of countries that are strongly interested in the sustainable development of countries, regions, and enterprises. One of the aspects that affects the indicators and prospects of sustainable development is the efﬁciency of energy source use. Nationwide reductions in the greenhouse gas emissions of motor vehicles could have a direct effect on ambient temperature and reducing the effects of global warming, which can affect future environmental, societal, and economic development. Signiﬁcant reductions in fuel consumption can be achieved by increasing the efﬁciency of use, and the performance, of current cargo vehicles. This aspect is directly related to cargo delivery systems and supply chain efﬁciency and effectiveness. The article solves the problem of increasing the effectiveness of cargo delivery and proposes a model that would minimize transportation costs that are directly related to fuel consumption, shortening transportation time. The model addresses the problem of a lack of models evaluating the efﬁciency of cargo to Lithuania that is using several different modes of transportation. For the solution to this problem, the article examines the complexity of the rational use of land and water vehicles depending on the type of cargo transported, the technical capabilities of the vehicles (loading, speed, environmental pollution, fuel consumption, etc.), and the type (cars, railways, ships). The novelty of the ﬁndings is based on the availability to select the most appropriate vehicles, on a case-by-case basis, from the available options, depending on their environmental performance and energy efﬁciency. This model, later in this article, is used for calculations of Lithuanian companies for selecting the most rational vehicle by identifying the most appropriate route, as well as assessing the dynamics of the economic and physical indicators. The model allows for creating dependencies between the main indicators characterizing the transport process—the cost, the time of transport, and the safety, taking into account the dynamics of economic and physical indicators, that lead to a very important issue—reducing the amount of energy required to provide products and services.


Introduction
With the growth of production, and the increasing intensity of transportation of goods, the proper organization of cargo transportation occupies a significant place, because every consignor wants to deliver his cargo at the lowest possible rate, in a timely and safe manner. It is necessary to use the logistic possibilities provided for the effective solution of these tasks. Vehicles are one of the key factors in the supply chain, which affect energy source use. Today, cargo to Lithuania is transported using several modes of transport, therefore raising the question of efficiency. There is a lack of models evaluating the efficiency of such transportation to Lithuania. It details the need for models that are solving this specific issue. The model addresses this research gap problem of models by evaluating the efficiency of several different modes of transport existing in Lithuania. The authors of the article propose a model that could solve problems of efficiency in transporting cargo to the country. The integrated vehicle application guarantees the most rational solutions

Literature Review
The aim of the automotive industry-to use energy more efficiently in a vehicle-is directly related to the efficient organization of freight transport [3]. When analyzing the organization of cargo transportation, and the willingness of every consignor to deliver his cargo more efficiently and at the lowest possible rate, it is necessary to consider many related factors [4,5]. Nakandala, Lau, and Zhang [6] describe ways to reduce the cost of planning for three different road vehicles, where weight constraints prevail when transporting containers; when volume limitation dominates; when volume constraints prevail and the modeling of transport cost optimization to maintain fresh product quality. Lightner-Laws et al. [7] propose to calculate The Global Vehicle Route Problem using a genetic algorithm and an automated work distribution process, reducing the total time that all jobs are delivered, using the lowest number of drivers.
The multifunctional problem of optimization is solved by modeling vehicle movement. The specific software has developed an objective modeling model that can simulate the efficient distribution process by changing parameters dynamically and randomly [8][9][10][11].
The linear mathematical system of road transport infrastructure, presented by Andy Chow and Li [12], is based on traffic flow and concentration ratio and development in space and time. Optimal use of transport operations, namely planning of transport services, diversion of vehicles, control of transport, and selection of optimal routes improves customer service, increases the use of transport assets, reduces transport costs and time, and increases the satisfaction of transport workers [13][14][15].
Karimi et al. [16] proposes the bi-objective mixed-integer linear programming model. The model aims to minimize the transportation costs, and the total delivery time, of the supply chain that has a direct impact on energy efficiency. Analytical and simulation models are presented by Lukinskiy et al. [17] which provide a probabilistic assessment of the unimodal and multimodal international transportation JIT implementation. The first model, in which the order of the implementation of operations does not affect the result, is based on the theory of probability: the composition of distribution laws, random theories of digital characteristics of the variables, and the formula of total probability. The second model reviews the impact of the implementation of operations on transport and their interrelationship, based on the modeling of statistical experiments and presented as an appropriate algorithm that allows us to consider technical organizational constraints. Gen and Syarif [18] proposed a spanning, tree-based hybrid genetic algorithm for solving the multi-product, multi-time period production, distribution, and inventory problem with an increase in energy efficiency. Khalili et al. [19] and Devapriya et al. [20] analyzed a practical problem of a perishable product that must be produced and distributed before it becomes unusable, the goal was to reach the minimum cost. Their proposals create the possibility to reduce the amount of energy required to provide products and services. Chung et al. [21] proposed an optimization methodology that increases energy efficiency, named modified genetic algorithm with crossing date heuristic, to maximize the collection of used tertiary packaging for reuse, meanwhile minimizing the total operating cost by taking the advantages of simultaneous optimization of multi-day planning. Vinay and Sridharan [22] have proven that, when the fixed cost is included in transportation costs, solving the model for large-scale instances optimally requires enormous computational time and effort. A hybrid solution methodology that integrates fuzzy set theory and the ε-constraint method is proposed by Balaman et al. [23]. Wang and Chen [24] presented a mathematical model that solves simultaneous pickup and delivery vehicle routing problems. He proved that the model could significantly reduce the logistics distribution cost and increase energy efficiency. The transportation network is designed to specify the most appropriate transportation mode, and related transportation options, under transfer station availability limitations.
Kuo and Wang [25] surveyed the route plan to calculate the total fuel consumption model. The model also considers three factors that have a significant impact on transport costs-transport distance, transport speed, and loading weight. These factors are used to optimize the route plan and experimental evaluation of the proposed method is performed. Li et al. [26] propose the mathematical optimization model that illustrates how to reduce overall costs, including vehicle operating costs, cost of quality loss, product freshness, cost of delay, energy costs, and greenhouse gas emissions. It also calculates the impact of maximum vehicle load changes on costs and greenhouse gas emissions. Li et al. [27] describes a model that implements a reduction in the cost of a building, and the same principles could be used for reducing the cost of transportation. Jakimavičius and Burinskiene [28] examine how sustainability is defined in terms of energy consumption and air pollution emissions. Authors think that, for this case, the best solution may be more efficient and alternative fuel vehicles. They also examine other strategies, such as Simple Additive Weighting, Technique for Order Preference by Similarity to Ideal Solution, Complex Proportional Assessment that helps achieve planning goals, such as congestion reduction, facility cost savings, increased safety, improved mobility for nondrivers, or more efficient land use development. The proposals could be divided into two groups: the technological and the managerial.
Shankar et al. [29] described a risk analysis approach, that is developed by innovatively integrating the intuitionistic fuzzy set theory and D-number theory to quantitatively model the sustainability risks. The authors also describe the changes in the sustainability of transport in different years. Mesa-Arango and Ukkusuri [30] says, that the key variables used by shippers to determine freight mode choice are cost and delivery times. Often used models that rarely capture the effects of cargo weight, running speeds, network capacity, or route characteristics, even though they are key inputs to any logistical analysis. This study gives a mechanistic method to determine variable rail costs on a single corridor.
Concerns about fossil fuel consumption and the environment, coupled with increasing energy demand, have resulted in collective interest in sustainable and renewable energy. Scientific articles have recently looked for ways to reduce the level of greenhouse gas (GHG) emissions [31,32]. Other methods investigated by researchers [33] relate to the management of the global warming process using renewable energy sources. These authors used a multiobjective model proposed to structure a sustainable supply chain for second-generation biofuel under a triple bottom line approach into the supply chain model incorporating the carbon tax and cap. The model applies the improved augmented ε-constraint approach.
According to Vaquero et al. [31] a key aspect of sustainability in the supply chain is compromise and consideration of factors as the choice of the mode of transport affects the economic, environmental, and social dimensions of the Triple Bottom Line. The tripple bottom line method is widely used in various fields. Bataglin and Ferreira [34] proposed a method for the modularization of sustainable product development with an emphasis on the Triple Bottom Line, using as modularization criteria economic, environmental, and social indicators.
Hollos et al. [35] performed a literature review from which they inferred that sustainable business behavior for the consequences of sustainable supplier cooperation is related to a triple bottom line. The authors investigated which research framework is a path analytic model with six latent variables: sustainable supplier co-operation, strategic orientation, social practices, green practices, operational performance, and, finally, the cost reduction.
Burki et al. [36] make a point on the efficacy of innovation as a determinant of sustainable performance in supply chains. Researchers examine green innovation practices as determinants of the triple bottom line performance.
Sun et al. [37] present a novel approach by monitoring the rail corrugation online using a self-contained energy harvesting system. The authors proposed to improve, by a triple-magnet configuration, the mechanics of the magnetic-floating energy harvester and investigated differences between repellent and attractive configurations.

Structure of Indicators of Functional Dependencies between the Main Transport Process
The transport process can be described by several indicators [44,47]. Usually, the speed, safety, and price of transportation are examined. Each of these indicators evaluates the quality of the carriage in an appropriate aspect. These indicators need to be guided to select a particular route.
The aggregate indicator of the transport process is the total cost of transport. These indicators consist of: the technological costs of the carriage, which include time, transport, loading, storage time, and insurance costs. Total shipping costs are calculated using the Formula (1): here: C a -the technological cost of car transportation; C g -the technological cost of rail transport; C j -technological costs of maritime transport.
Technological costs of transportation include fuel (or other energy) costs, transport supplies (e.g., lubricants), vehicle maintenance and repair, road charges (rail or ports), transportation equipment insurance, driver fees, and forwarding costs.
Fuel (or another energy source) costs are estimated in cash, fuel consumption in kg, or liters multiplied by the price per unit. The fuel or energy consumption itself is predicted according to the relevant methodologies adopted by the modes of transport. Fuel consumption (or energy consumption) is proportional to the distance traveled by the vehicle and is therefore estimated at EUR/km.
The cost of transport materials is also calculated by multiplying the cost of the consumables by the cost per unit of the respective consumables. Consumption of consumables is provided in the technical documentation of the vehicles. These costs are simply proportional to the distance traveled by the vehicle (or proportional to the lifetime of the vehicle). Therefore, these costs are usually estimated at EUR/km. Expenditure on vehicle maintenance is the cost of maintaining the vehicle in the required technical condition (e.g., replacement of parts, diagnostics, technical prophylaxis). The cost of repairs is the cost of restoring the vehicle to its state of repair (e.g., after a breakdown, after an accident, or after a scheduled repair). Maintenance and repair costs are proportional to the distance traveled by the vehicles, so in this case, they are valued at EUR/km. The tolls for the different modes of transport are calculated according to different methodologies [48,50].
The toll is collected as part of the fuel tax. The port fee is calculated by taking into account the length of the berth used and the period of its use. The railway infrastructure charge is calculated by taking into account all the costs of the infrastructure (running costs of the railway itself, the costs of organizing and operating the railway, the additional costs of rolling stock for rolling stock or freight). In most cases, they are estimated by dividing the total annual cost of rail infrastructure in the country by the annual freight turnover in tonne-kilometers. This determines the infrastructure cost for the next year (usually adjusted annually). The railway infrastructure charge is calculated in EUR/tkm or EUR/km. Expenses for transportation equipment insurance are costs for vehicles, mobile loading devices, and containers. They are proportional to the amount of transportation equipment that is prohibited. If the transport equipment is operated evenly (without extraordinary downtime), with the increase in cargo turnover, these costs are increasingly dependent on freight turnover or vehicle mileage. Therefore, the cost of transport equipment insurance, in this case, is estimated at EUR/km. Drivers' salary costs depend on the number of hours worked by drivers (also the cost of social insurance, income tax, disability benefits). If drivers work smoothly, their working hours are proportional to the distance traveled by the vehicles. Therefore, the cost of driver salaries is estimated at EUR/km. Freight forwarding costs depend on the terms and conditions of the forwarding agreements (rates and volume of services provided in them). Continually using the same services, costs are proportional to the distances traveled and are valued at EUR/km. Total time consumption is a product of transportation, loading, storage time, and approximate time input. Approximate time costs are expressed in euros per hour. Such a dimension makes it possible to turn the amount of time into hours into the time of the euro. By expressing the time costs in euros, we will be able to use it as a part of a standard model with other cost expressions. Comparative time costs assess vehicle charges, possible fines, delays in delivery, lost revenue due to occupied vehicles, occupied technological equipment, and busy people.
The total transportation and storage time has two components: transportation time and loading and storage time at the terminals. Total shipment and storage time is calculated according to Formula (3): here: T t -transportation time, h; T s -loading and storage time, h. Transport time is the time needed for the transport process itself. Load and storage time is the time spent in warehouses for unloading, loading, and storage. Total shipment and storage time for each route is specific, depending on route structure and stevedoring terminals.
The structure of the described costs will be reviewed in more detail. Expression of total time consumption is expressed by Formula (4): here: l a , l g , l j -distance caring in different transport mode (cars, rail, sea), km; V a , V g , V j -speed in the appropriate mode of transport, km/h. Insurance costs D is the cost of paying insurance premiums. These costs do not include vehicle insurance premiums (they are valued at the cost of technological transportation). Only cargo insurance premiums are valued. They depend on the nature of the cargo, the route of transportation, and the time of transport.

Calculation of Car Transportation Costs
Meixell and Norbis [5], analyzing tendencies in the motor carrier industry, indicates that as fuel prices have risen, carriers must raise prices, or to find an alternative solution for increasing the effectiveness of energy consumption. Carrier selection criteria is ranked in [45], paying particular attention to why shippers and carriers rank selection factors differently. An optimization model for carrier selection is presenting [45]. Their model has an objective to minimize freight cost, while also including issues of energy efficiency. The model by Liao and Rittscher [49] integrating three groups of decision variables: dynamic procurement lot sizing, supplier selection decision, and carrier selection.
Technological cost of transportation is estimated by determining the carriage rate, EUR/tkm and is calculated according to the formula [46]: here: l a -distance when transporting with cars transport, km; m a -cargo weight transported with cars transport, t; k a -one tkm cost of car transportation, EUR/tkm. The product of cargo mass and distance sometimes is called transport work and measured in tkm. Multiplying it by the cost of one tonne of the carriage, we get the technological cost of transporting transport cars. The cost of transporting one tkm per car depends on the speed of the car: k a = k av ·k v1 here: k av -the average cost of transport by car, EUR/tkm; k v1 -coefficient of car speed. The coefficient of car speed is calculated using the Formula (7): here: a a , b a and c a -equation coefficients for vehicle rolling resistance [37].
Having set the coefficient values, we get the equation: Putting Equation (8) to (6), we get: The average cost of transport by car is 0.02107 EUR/tkm. This figure is estimated to be 20 tonnes per car, 28 L of diesel fuel per 100 km, costs of 1 L diesel are EUR 0.87 without VAT [52], and other car operating costs are estimated by multiplying these costs by a factor of 1.75.
Putting Equation (9) to (5), we get: C a = 0.02107 (0.216 + 0.000952·V + 0.0238 × 10 −3 V 2 ) m a l a ; The graphical interpretation of the Formula (10) when the en-route speed of the car is 70 km/h is shown in Figure 1. The developed formula for convenient practical use by approximate estimation of primary costs for freight vehicles. Depending on the dependence shown in the formula, we can calculate the cost of transportation by knowing the weight of the cargo and the transport distance. Fr example, if the weight of the load is 2000 t and the transport distance by car is 300 km, the cost of technological transportation is 12,642 EUR.

Calculation of the Cost of Rail Transport
Technological transport costs are estimated by determining the rate of rail transport, EUR/tkm, and are calculated according to the formula: here: l g -distance by rail, km; m g -the weight of the cargo that is transported by rail, t; k g -the cost of one tkm of transport by rail, EUR/tkm. The cost of transporting one tkm by rail depends on the speed of the railway: here: k gv -the average cost of transport by rail, EUR/tkm; k v2 coefficient for estimating the speed of the train. The average cost of transport by rail depends on the distance traveled. The tariff table shows the discrete dependence of transportation rates on distance. This data is approximated by the least-squares method, and we get a gradual function. According to 2018 tariffs [53], the gradual function coefficient is 0.2047, and the degree is equal to −0.1515: A degree of less than one means that its effect on function decreases as the transport distance increases. Minus stands for inverse proportionality-decreasing the distance, decreasing the transportation tariff. The coefficient for calculating the railway speed is calculated according to the formula: here: a g , b g and c g -equation coefficients for rail rolling resistance [46].
Having set the coefficient values, we get: to (13) and (15) respectively inserted into Equation (12), we get: After setting the Formula (16) to (11) and following the steps, we get: A graphical interpretation of Formula (17) at a rail en-route speed of 40 km/h is shown in Figure 2. Depending on the dependency shown in Equation (17), we can calculate the technological costs of transportation by knowing the weight of the load and the transport distance. For example, if the weight of the load is 500 tons and the distance of transport by rail is 200 km, then the technological transportation costs are 9173 EUR.

Calculation of Sea Transport Costs
Technological transport costs are estimated by determining the shipping price, the EUR/tkm, the relevant insurance costs and are calculated according to the formula: here: m j -cargo weight transported by sea, t; k j -shipping costs of one tonne (t) by sea, EUR/t. The cost of transporting one tkm by sea depends on the speed of the vessels: here: k jv -the average cost of transport by sea, EUR/tkm; k v3 -coefficient of the speed of the ship. We calculate the average cost of shipping by sea based on the distance-based approximation of the discrete transport tariff dependence on the smallest squares method: here: l j -distance when transporting by sea, km. The coefficient for calculating the railway speed is calculated according to the formula: here: a j , b j , c j -equation coefficients to estimate the motion resistance of ships. After setting the coefficient values to (21), we get: After setting in (20) and (22), we get the Equation (23): After setting Equation (23) in Equation (18), we get: Based on 2019 data [48] function coefficients are equal to −6 × 10 −6 and 0.0677, respectively. The minus sign at the square member coefficient means that the extremity of the function is the maximum. The maximum feature is at a transport distance of about 5000 km. For transportation in the European Union region up to 5000 km, we pay for the maintenance of ships, fuel, and other consumables at European prices, i.e., relatively expensive. For longer distances, maintenance can be carried out in South American or African countries, where prices are lower, leading to lower transportation rates. The ship sails between 4 and 6 thousand kilometers, fuel and maintenance costs increase, but taking advantage of cheap third-party resources does not allow such a sailing distance. This is indicated by the extreme in the appropriate graph spot. A graphical interpretation of the Formula (24) when the ship's en route speed is 25 km/h is shown in Figure 3. Depending on the dependency shown in Figure 3 we can calculate the technological costs of transportation by knowing the weight of the cargo and the distance of trans-portation. For example, in the case of sea transport, the cost of transporting 1000 t of cargo, 4000 km, will be equivalent to EUR 174,000. The same mass will cost 157,000 EUR at 8000 km. With twice the transport distance, transportation costs have not increased but have fallen. The cost reduction was due to much cheaper fuel in Asia, Africa, and South America.

Calculation of Total Technological Transportation Costs
We insert Equations (10), (17), (24) and (4) into the Equation (1). Then the total cost of transporting the goods in different modes of transport (car, rail, and sea) will be: here: V a , V g , V j -route speed for the transport mode concerned, km/h; τ a , τ g , τ j -comparative cost by time for the relevant mode of transport, EUR/h; ∆V a , ∆V g , ∆V j -Differences between planned V pl and actual V f speed V pl − V f = ∆V i , then i = a, g, j, km/h.
They evaluate unplanned vehicle downtime on the road between terminals (customs, road repair, accidents, etc.); T s -storage time, h, which is different at each terminal, because it depends on its technological capacity (loading device, methods), management structure and is calculated or selected individually; τ s -comparing costs by the time of storage in terminals, taking into account the costs of unloading and loading, EUR/h. This indicator may be determinative or variable depending on the duration of storage, type of cargo, and nomenclature of works.
As we can see from (25), in a common model for determining the cost of technology transportation, we must estimate the cost of the cargo at terminals and insurance costs. As mentioned above, these costs are quite different and depend heavily on specific stevedoring, warehousing, insurance, and specific agreements between the parties involved. Therefore, to determine the impact of purely used vehicles on the quality of the transport process, we will not evaluate the last two members of the expression (25). In the model, they will have to be assessed on a case-by-case basis.

Setting the Goal Function
For setting the goal, it used the pre-emptive priority goal programming modeling where a deviation from a higher priority level goal is considered to be infinitely more important than a deviation from a lower priority goal [50]. The quality of the freight process is determined by three key indicators: cost, time, and reliability, which also includes security. The technical safety of cargo in this model is not assessed separately, as it is quite difficult to determine it based on statistical data, and it is also different on separate routes. However, this factor is partly assessed by insurance costs. Since the factors are two, the priority queue will have two options.
For the first option, the goal function is to minimize costs when the transport time is limited. This is a good option when a customer wants to get a service with minimal costs and freight is not urgent. In the case of regular freight, the shipment is on the way and no matter how long it takes on the road. For example, when transporting fuel to the power plant, it is important that the power plant periodically receives a certain amount of fuel with the lowest transportation costs, and for how long the fuel will be transported and on which routes it is not important for the power plant.
Option II: Goal function-T = >min; Restriction-Z < Zmax. For the second option, the purpose of the goal is to minimize the time limit for transport costs. This option is suitable when a customer wants quick delivery of goods and is willing to pay a minimum amount. For example, when delivering parts to collection workshops, they must be delivered at the right time (parts are not accumulated in modern collection workshops but delivered exactly before collection). Therefore, the delivery of parts is a priority.

An Example of Model Realization, According to Two Routes
Two different theoretical routes were modeled using real transport cost data based on data collected from Lithuanian cargo companies. The first of these consists of the railway line from the starting point A to the congestion point 1, the section of the car road between the overload points 1 and 2, the section of the sea transport between the overload points 2 and 3, the car road section from the overload point 3 to the route endpoint B. The second route consists of the car road section of the road from the starting point A to the congestion point 4, the section of the route by sea transport between the congestion points 4 and 3, the car road section from the overload point 3 to the route endpoint B. Examples of the considered routes are schematically shown in Figure 4.  Figure 4 shows: A, B-start and end of routes; 1, 2, 3, 4-congestion points (terminals); a, g, j-car, rail and sea routes. Values of dimensions in the example: expected route speeds: Car transport: V a = 70 km/h; rail transport: V g = 40 km/h; for sea transport: V j = 25 km/h; lengths of carriages are also covered: car transport, route I: l a1 = 150 km; rail transport: l g1 = 600 km; for sea transport: l j1 = 150 km; lengths of carriages, route II: car transport: l a2 = 150 km; transporting with rail transport: l g2 = 0 km; for sea transport: l j2 = 900 km; the weight of cargo carried is m = 2000 t. Unforeseen vehicle downtime on the road will not be assessed yet, τ a , τ g , τ j = 0 (we will only count the technological costs).
We assume that insurance costs are already included in the rate, so we will assume in the formula that insurance costs are zero (D = 0). However, if this is not the case, they must be evaluated. The dependence of technological transportation costs on the amount of transported cargo (according to Formula (2)) at a vehicle speed of 70 km/h is shown in Figure 5. Taking into account the time of transportation, we assume that an entire load of cars is loaded immediately, the number of cars does not limit the volume of transportation [41,42]. If not, i.e., the cargo would be transported with a smaller number of cars, then the road and the transport time would increase accordingly.

Fluctuations in Transport Time and Cost with Vehicle Speed
If we consider the speed of one of the modes (cars) to be variable, then the solutions of the Equation (25) will take the form below (time consumption depends on the type of terminal, its stevedoring facilities, management, organization, etc.), so in each terminal, it may vary, and at the same time, the hourly price may change. In the example to be realized, we take a determinate 1-h price τ = 3 EUR, based on the example of the Klaipeda port terminal.
The calculation results (see Table 2) suggest that the shortest delivery time is 16.16 h on route I, with a car speed of 90 km/h, but the lowest cost will be 110,378 EUR on Route II, with an average car speed of 50 km/h. According to the calculation results, when transporting the first route (road, rail, and sea transport) we get a graphical interpretation, which is shown in Figure 6.  Graphical interpretation according to the second route (when transported by road and sea), changing the route speed of cars is shown in Figure 7. We can see in the graphs that as car speeds increase, transportation time decreases and costs increase as vehicles use more fuel at higher speeds.

Transport Time Dependence on the Number of Cars Used
If the car's load capacity is 20 t, then the transport time dependence on the number of cars used will be as follows: here: n-the number of cars. The graphical expression of the Formula (26) is shown in Figure 8. The greatest impact on the number of cars have up to 50 units. Their total payload is 1000 tonnes or half of the total load. When the number of cars increases, its impact on the transport time decreases accordingly.

The General Solution of the Mathematical Model
For a set of routes according to the given methodology, we calculate the time and cost of transportation for each of them [28]. The general solution for routes I and II, where the vehicle speed is 70 km/h, is calculated based on the following data, see Table 3   Table 3. Estimation of routes I-III, general solution.
Route I: When the speed of the car is 70 km/h we get a shorter time (17.14 h) on route I (road, rail and sea transport), but lower costs (110,814 EUR) on route II (road and sea transport).
A set of points whose coordinates are transportation costs and transportation time make up the overall solution. There were two constraints on the target function setting, one on cost and one on time constraint (see Figure 9). The overall solution was obtained by saying that the transport speed in the modes of transport is constant. In general, the most rational use of vehicles is on routes that do not go beyond time and cost constraints, but the customer can always choose the most appropriate.

Expression of the Cost of Transport Dependency, Taking into Account the Evolution of the Railway Tariff
In different European countries, rail freight rates vary, even several times. So far, we have relied on Lithuanian rail freight rates. In order to apply the calculation methodology to the transport systems of other countries, it is necessary to evaluate the development of the transportation rates of the railway, if necessary other modes [43]. A coefficient of b is used for this purpose. In the expression of the total cost of carriage, we multiply the rail freight rate by factor b. Then the expression of total costs will be as follows: Z = 0.02107 − (0.216 + 0.000952 V a + 0.0238 × 10 −3 V 2 a )·m a l a + b 0.2047·l g 0.8485 m g 0.739 + 0.00325 V g + 0.0814 × 10 −3 V 2 g ) + (−6 × 10 −6 l 2 j + 0.0677 l j )·(0.8483 + 0.00373 +0.09338 × 10 −3 V 2 j ) m j + (l a /∆V a )·τ a + (l g /∆V g )·τ g ·+ (l j /∆V j ) ·τ j ·+ T s τ j + ∑D (27) here: b-the coefficient of the change in the railway tariff. The coefficient b shows the number of times the rail freight rate in the respective region is higher than the Lithuanian rail freight rate. The rail fare may also change for other reasons, as it may be the subject of negotiations. The general solution of the mathematical model estimating the coefficient of variation of railway tariffs b is shown in Figure 10.
In a general solution, instead of set points, we have a set of sections. This means that on each route the cost of transportation will change, with the change of rail tariffs. The railway tariff varies for a variety of reasons: different transport operators, different transport times (moments), the composition of the transport service package. Correspondingly, the total cost of transportation will vary differently, depending on the length of the sections shown in the graph.

Discussion
In the case of Lithuania, carriers of different types of cargo compete with each other for cargo transportation. The question arises as to how to reconcile the economic interests of all types of carriers by bringing them together for the common goal of reducing energy. The need for this research is based on a demand to promote collaboration activities of different types of cargo and to increase the effectiveness of this process, which impacts the sustainability of a particular country.
Route modeling is widely used in the public transport sector, where a passenger can reach his destination using different modes of transport. Different computer programs for routing have been developed and offered to users. In Lithuania, public transport belongs to the state and municipalities, so route modeling is an easy problem to solve. The situation is somewhat different in the field of freight transport, where many private carriers operate, transporting goods of different sizes and weights. The developed model can contribute to the development of computer programs in the field of freight transportation in Lithuania, as well as the logic of its solutions for solving transportation problems can be adapted in other countries.
One of the limitations of the model is that the ideal case is considered. The transportation process involves different private companies competing with each other for cargo that often has different purposes. In this case, the factors that reduce energy consumption are not always an advantage when choosing transportation services. The final cost of transportation also affects the country of origin of carriers who have different rates of remuneration of the drivers. Results of the research are beneficial for the optimization of the entire transport process, minimization of time or cost by selecting the most appropriate route for available roads, terminals, vehicles of a particular country, and creates a possibility to organize this process in a more effective way, creating the added value for implementation goals of Sustainable Development Agenda.
Further studies are needed to assess the factors that would help reconcile the interests of different transport companies in the field of freight transportation. The next step in the study would be to justify the economic benefits for each mode of transport individually, which could encourage the use of such a model in practice.

Conclusions
As a result of the research, the mathematical model was designed to optimize the entire transport process, minimizing time or cost by selecting the most appropriate route for available roads, terminals, and vehicles. The mathematical model for selecting complex vehicles in order to optimize the multimodal cargo transportation process increases energy efficiency, creates an opportunity to reduce the amount of energy required to provide products and services.
A computer experiment has been performed showing the advantages of the mathematical model created by selecting rational types and types of vehicles for cargo transportation. The model allows creating dependencies between the main indicators characterizing the transport process-the cost, the time of transport, the safety, considering the dynamics of economic and physical indicators. In the case of dependencies between the key indicators describing the transport process, it is possible to describe each of the routes in question in a specific case with two main aggregated indicators-cost Z and time T. The model enables the description of the route in a two-dimensional coordinate system, based on the constraint conditions. In the model, the linear programming problem is solved. After solving the linear programming problem, the most suitable routes and vehicles are selected according to the intended conditions. The created model could be used for rational vehicle selection by identifying the most appropriate route, minimizing transportation time and cost for the business companies.
The main advantage of the model is that the mathematical model allows to minimize transportation costs, shorten transportation time, and to select the most appropriate vehicles on a case-by-case basis from the available options depending on their environmental performance and energy saving.
There are two major limitations in this study that could be addressed in future research. First, the study focused on situation analysis in Lithuania and the investigation is based only on the data collected in this country, so this mean research could be extended by analyzing situations in others, firstly neighboring countries. However, it is possible to adapt the model for other countries, taking into consideration peculiarities. Second, the research could be extended by analyzing additional components that affect the fuel savings, as technological factors and increasing efficiency using other methods. Further on, other directions could be involved in the research.

Conflicts of Interest:
The authors declare no conflict of interest.