Simultaneous Long-Term Planning of Flexible Electric Vehicle Photovoltaic Charging Stations in Terms of Load Response and Technical and Economic Indicators

: Photovoltaic charging stations (PVCSs) are one of the most important pieces of charging equipment for electric vehicles (EVs). Recently, the process of designing solar charging stations as ﬂexible sources has been growing and developing. This paper presents a relatively complete design of a solar charging station as a ﬂexible economic resource in a 10-year planning horizon based on a genetic algorithm in two scenarios. PVCSs are not considered in the ﬁrst scenario. This scenario is only to conﬁrm the results, and the proposed method is proposed. However, in the second scenario, the effects of PVCSs and the demand response strategy (DR) on this development are considered. Copula probability distribution functions are used to create appropriate scenarios for vehicles during different planning years. The proposed energy management system shows a stable performance in terms of the annual load growth index and electricity price of each level of demand over the time horizon along with minimizing power losses and costs required, which makes PVCS efﬁciency higher and gives them a suitable structure and stability. The modeling results in terms of uncertainties in the system indicate that the use of load management along with PVCS design and ﬂexible electric vehicle charge control strategies improves power quality parameters and optimizes system cost over a period of 10 years. Compared to the obtained results with the traditional case, it is observed that long-term planning in terms of DR and PVCSs and the technical speciﬁcations of the network have been improved. As a result of this proposed long-term planning, PVCSs are more ﬂexible. growth application considering three electric vehicle parking lots.


Introduction
The expansion of PVCSs is an effective measure to reduce our dependence on fossil fuels, and its widespread expansion across all countries will occur soon. Due to the advancement in flexible electric vehicles, it is expected that numerous EVs will soon be connected to the power grid for charging [1]. The energy needed to charge these vehicles can be modeled as an additional load for the primary power grid [2]. In order to promote the proper development of EVs and their charging equipment in the power grid, the planning of the future distribution network of both EV charging demand factors and their discharge capacity should be considered [3][4][5][6]. Therefore, in addition to meeting the demand for power in routine distribution network scheduling, distribution network performance planning, including charging stations, requires solutions to other technical issues such as 2 of 20 charging station capacity planning, charging station installation location, power feeder scheduling and substation scheduling, includes EV charging stations [7,8]. Distribution network development planning is a traditional and large-scale optimization problem that determines how to distribute distribution network equipment over several years to ensure a reliable supply of increasing load. Usually, the objective function of the DEP problem is modeled as a cost function, which includes the cost of upgrading or replacing equipment, constructing or increasing the capacity of substations over different time periods. Recently, useful meta-heuristic algorithms have been proposed by researchers [3][4][5][6][7][8][9]. In recent years, several models and methods have been proposed to solve the DEP problem using DG. In the reference [9], the modified particle aggregation algorithm (MPSO) is introduced to solve the DEP problem by considering DG and sources of electrical energy storage. The proposed model optimally selects the total investment and operating costs of the DG and the distribution network. The authors in [10] have put on their agenda a dynamic model for the development of distribution networks to determine the appropriate location, time and capacity of the DG and the appropriate planning to increase the capacity of the lines in the distribution network. Additionally, the DGEP problem is solved with the help of SST and optimal load distribution (OPF) in the reference [11]. Similarly, the reference [12] has proposed a combined method to solve the DGEP problem to reduce the peak load of scattered products. Additionally, the issue of DGEP in reference [13], with the approach of minimizing the total investment costs and current operation in distribution networks in the presence of DG, has been considered. In [14], improved genetic algorithm and XGBoost classifier for power transformer fault diagnosis is utilized. In [15], an effective model for the DGEP problem is proposed. By solving the proposed model in the form of a complex multivariate optimization problem with a variety of discrete and continuous variables, the optimal location of DG installation and the optimal strategy to increase the capacity of existing lines in the network for a certain period of time are determined. The main goal of network design development planning is to create a reliable service with high reliability and cost-effectiveness for customers, EVs and owners of parking lots so that the quality of voltage and power within the allowable range is guaranteed [16]. Therefore, considering the objectives of reliability and cost to provide a complete model of the problem seems necessary. Additionally, in recent years, with the increasing influence of distributed generation at the distribution level, the planning of distributed generation resources is usually performed in conjunction with network development. In this regard, once PQ or PV are included in the problem, distributed generation sources can be modeled broadly or precisely as wind power or other types of power plants, as shown in [17,18]. In the reference [19], the uncertainty of scattered wind production is also considered for more accurate modeling. In order to model the wind speed uncertainty, rail probability distribution functions have been used, which are created after the discretization of probabilistic scenarios. In [20], a study has been conducted on determining the optimal charging profile for demands in peak hours. It has also investigated the effect of the charge profile on the distribution network. The optimum parking location for electric vehicles connected to the network was specified by respecting energy loss and reliability as economic constraints. The production of electric vehicles that can travel long distances necessitates the expansion of charging stations offering a notable amount of electrical energy to supply vehicles. Therefore, due to the progressing employment of electric vehicles globally, a futuristic plan for raising the number of such stations is apparent [21]. For this reason, extensive research has already been conducted or is underway in various aspects of the modeling, planning, designing and operating of charging stations. In [22], dynamic combined economic emission dispatch integrating plug-in electric vehicles and renewable energy sources is performed. A network with optimal parking locations and distributed sources is improved so that, in addition to providing the desired reliability, the losses also reach the optimal level. The goal of [23] was to maximize profit by modeling smart parking. In the proposed model, several economic and technical indicators were considered, along with the security constraints of the system. The authors in [6] presented a multi-objective algorithm regarding the number of electric vehicle parking lots, the location and capacity of these parking lots and the energy scheduling of power supply sources in the system [24]. Due to the growing adoption of electric vehicles in urban areas, researchers in [25] have examined electric vehicle charging demand models for urban settings. The authors in [26] developed a time-space model combining transport analysis with a power system analysis. The authors in [27] provided a time-based electric vehicle-demand-forecast model with multiple charging stations in an urban area. Through a cell transfer traffic model, the vehicle entry rate at a charging station was predicted. The authors in [28] estimated the time-space character of a charging station by modeling and simulating according to actual traffic data. In [29], an innovative dual model is proposed; a new model for parking electric vehicles at charging stations is presented at the first level. Moreover, a novel model is developed at the second level to minimize overall system costs, considering technical constraints on progressing distribution networks. Additionally, applying renewable sources such as photovoltaic cells is included. The advancement of solar charging equipment is deemed a clean and workable energy supply infrastructure for electric vehicles [25]. There is no restriction concerning the location of photovoltaic charging stations; such stations have been placed in commercial buildings and residential complexes [30]. Significant efforts have been made to study the energy management strategies of photovoltaic charging stations, such as displaying charging methods for network-connected hybrid electric vehicles (PHEVs). Based on [31], a substantial number of grid-connected hybrid electric vehicles are linked to the grid in a coordinated manner. Following the management of an entity, such as an electric vehicle parking lot, they can be set up as a small virtual power plant with a short start-up time and no considerable cost. Accordingly, in [32][33][34][35][36], fast energy management algorithms have been implemented for grid-connected charging parks in industrial/commercial locations. Additionally, in [37], an innovative operation strategy is proposed for micro-grids of commercial buildings, including photovoltaic array systems. Principles of decisionmaking strategies concerning the improvement of consumption of photovoltaic-based energy and a reduction in adverse effects on the power grid are discussed. In [38][39][40][41], classifications of electric vehicles for photovoltaic charging stations (PVCSs) are introduced, which diminishes energy costs associated with power supply. Further, the power system flexibility issues are also discussed with the renewable energy generation. An Energy Management System (EMS) is essential to achieve the maximum capabilities of micro-grids, such as increasing reliability, improving power quality, reducing the cost of energy supply and reducing greenhouse gas emissions. Thus, the optimal performance of the network in the micro-grid will be attainable. Different structures for EMSs employing different optimization algorithms and different shapes for micro-grids are presented in relevant research, some of which aim to optimize the performance of the source in the micro-grid system and minimize operating costs [42]. Economic aspects have been discussed in other proposed EMS schemes, and associated strategies have been reported to achieve maximum profit [43]. According to [44], various methods have been proposed to establish energy balance in micro-grids due to limitations in applying renewable energy sources. In [45], fault locating transmission lines with thyristor-controlled series capacitors by fuzzy logic method is performed. In [46], short-channel effects improvement of carbon nanotube field effect transistors is discussed. In [47], an overview on functional integration of hybrid renewable energy systems in multi-energy buildings is performed.
Although comprehensive studies have been conducted in the long-term planning of simultaneous photovoltaic charging stations of electric vehicles, a deep understanding of terms associated with load response and economic indicators is demanding. There are various resolutions to meet the requirements of supplying the increasing number of electric vehicles. Some people are of the opinion that increasing the number of power plants will improve the final amount of electricity needed, consequently resulting in environmental side-effects and high costs. This paper aims to optimally plan the charging stations and electric vehicles parked in their lots. It is plausible to consider the charging stations as temporary batteries during non-peak hours by saving the extra amount of electricity from the network and returning this energy during peak hours. Accordingly, the same procedure is presumable for electric vehicles since they are manufactured with high-capacity batteries. Electric vehicle owners are free to choose whether they contribute to meeting the high energy demands during peak hours by linking their vehicles' batteries to the grid while they have parked their EVs. In this regard, this article is structured as follows. In Section 2, the formulation of the optimization problem is stated. In Section 3, the proposed method and an expression of data are asserted. Furthermore, in Section 4, according to the load response programs in different scenarios, the results and outputs of the simulation are given before the conclusion. Section 5 contains the conclusions and policy implications.

The General Framework of the Optimization Problem
In recent years, many methods have been used by researchers to model the uncertainty of possible quantities. Among these techniques, we can mention the methods of probability tree, Monte Carlo, fuzzy method, etc. The capillary method is another method that has received more attention in recent years. In summary, the most attractive advantage of capillary functions is their ability to estimate limit functions and the degree of correlation. Examples of random variables are copulas. Considering the dependence between variables, this creates a multivariate distribution function. The following is a summary of capillary performance steps in various issues as shown below: (1) First, the sample amount of the problem variables that should be examined and analyzed with the information about the problem are determined through the available real statistics. (2) Then, these amounts are normalized and mapped to zero and one.
(3) The capillary function is then applied to the normalized information, and the correlation coefficients between the variables are obtained. (4) Then, random samples are produced in large numbers using real normalized information and the obtained correlation coefficients. (5) Finally, the samples are mapped to their true range from zero to one.
In addition, the charging station owner, who is considered the distribution system manager in this study, can purchase batteries from battery manufacturers in bulk and at a reduced price, lowering the cost to EV owners. Charging stations can charge EVs to supply the energy needed to operate an EV and exchange power with the grid, and utilize the EVs' battery storage potential. EV owners can request a charge based on the pricing for charging EVs that has been announced. The charging station and solar power generation unit are designed to boost the distribution system's capability and meet the driving needs of electric vehicle users, according to the criteria indicated. The placement challenge is designed once the charging station and solar power-producing unit have been modeled. The issue of placement must be addressed in such a way that the operating factors are preserved. The optimal position and capacity of the solar power-producing unit and charging station will result from the planning problem. It should be emphasized that Monte Carlo simulation is utilized to answer the uncertainties in this paper, and multi-objective genetic algorithms are employed to solve the optimization problems. The multi-objective genetic algorithm is depicted in the flowchart below.
The proposed method's convergence is dependent on the number of starting solutions, just as other evolutionary algorithms. The algorithm begins by generating a random generation from the parent population because it is demographic. It uses a rapid nondefeat sorting strategy with computing complexity (MN2) to regulate non-defeated answers from the parent population. M is the number of objectives, and N is the population size. According to the layer in which it is positioned, each of the answers within a Pareto layer is assigned a reasonable value. Within each layer, the density distance criterion is utilized to measure and rank the solutions. This initiative has been launched multiple times for this purpose, with varying initial populations of 50, 100, 150 and 200 people. The proposed algorithm's convergence is depicted in Figure 1, indicating that an initial population of 150 people may be adequate to produce an optimal solution. As a result, we set the population value to 150 in the simulations in this paper. Figure 2 presented the flowchart of multiobjective genetic algorithm. In summary, based on Figure 3, the typical production pattern of the nth generation in NSGA-II is as follows: (1) Combine the population of parents (P o = P t ) and children (Q o = Q t ) and form the population (R t = P t ∪ Q t ) with size 2N. (2) Run the unsuccessful sort on R t to specify the different layers F i . i = 1, 2, . . . , l.
(3) Create a population of P i+1 size N by selecting superior solutions from non-dominant layers (F 1 , F 2 , . . . , F i ) in order of priority. (4) Calculate the population Q i+1 of children of size N using RKGA operators. size. According to the layer in which it is positioned, each of the answers within a Pareto layer is assigned a reasonable value. Within each layer, the density distance criterion is utilized to measure and rank the solutions. This initiative has been launched multiple times for this purpose, with varying initial populations of 50, 100, 150 and 200 people. The proposed algorithm's convergence is depicted in Figure 1, indicating that an initial population of 150 people may be adequate to produce an optimal solution. As a result, we set the population value to 150 in the simulations in this paper. Figure 2 presented the flowchart of multi-objective genetic algorithm. In summary, based on Figure 3, the typical production pattern of the nth generation in NSGA-II is as follows: (1) Combine the population of parents ( = ) and children ( = ) and form the population ( ∪ ) with size 2N.
(2) Run the unsuccessful sort on to specify the different layers . i = 1, 2,…, l (3) Create a population of size N by selecting superior solutions from non-dominant layers ( , , … . ) in order of priority.
(4) Calculate the population of children of size N using RKGA operators.    Figure 4 shows the execution time curve of the NSGA-II-based programming algorithm for the sample system based on the number of initial responses. This almost linear curve of program execution time demonstrates the necessity of selecting the number of early populations optimally. It should be mentioned that all simulations were performed on a computer with a memory of 32 GB, and Table 1 lists the major parameters of the proposed method as well as program execution time for core7, ram16G and graphic8G.
Chromosome coding: A string of chromosomes the size of a row of the dataset (number of characteristics in the dataset) is represented by a binary digit D for the task of selecting a subset of attributes. Each binary digit represents an attribute. A perspective of a chromosome used in this method is shown in Table 2.
In general, coding is the mapping from the search space to the space of strings made with length L and an alphabetic set of n characters, indicating the absence of features. Step of algorithm. Figure 4 shows the execution time curve of the NSGA-II-based programming algorithm for the sample system based on the number of initial responses. This almost linear curve of program execution time demonstrates the necessity of selecting the number of early populations optimally. It should be mentioned that all simulations were performed on a computer with a memory of 32 GB, and Table 1 lists the major parameters of the proposed method as well as program execution time for core7, ram16G and graphic8G.     Table 1 lists the major parameters of the proposed method as well as program execution time for core7, ram16G and graphic8G.

Modeling Uncertainties
Due to the wide use of distributed generation resources, different optimization methods have been proposed to find the optimal location and capacity, which are significantly affected by the load model and electricity price. However, since the load curve of subscribers and the cost of electricity are obtained from forecasting methods, their changes have a percentage of error and are not known definitively. In this study, uncertainties in the load continuity curve and electricity price are probabilistically modeled using the normal distribution. It should be noted that the Monte Vehicle log method has been used to simulate a load model and electricity price. Tables 1 and 2, in Appendix A, provided the descriptions of various parameters and indices utilized in below mentioned equations.

Load Uncertainty
According to the varied behavior patterns of distribution system subscribers in electricity consumption, daily load changes in each year of the planning period are modeled using the following parameters. The first parameter is the amount of load consumed in the first year of the planning period. Each year of the planning period is divided into: N dl . The demand level is h and the duration of each demand level (h) is τ. The second parameter of the load model is the demand level factor, which is obtained using the probability density function as follows [32]: λ is a random variable produced using the normal distribution function with a mean value of zero and a standard deviation of one for each level of demand. µ and σ are the predicted values of the load demand and their standard deviation, respectively. Considering the load growth rate for each year, the amount of electricity consumption at each level of demand is calculated [31]. The demand for active and passive power is presented in Equations (2) and (3), and apparent power demand is in Equation (4). The superscript e shows that the parameters belong to the Monte Carlo experiment e. , , , , , , , , ,

Electricity Price Uncertainty
Depending on the amount of electricity consumption and the behavior of electricity market users, the price is determined for the electricity purchased from the primary grid. Based on Equation (5), this amount changes for the cost of electricity at the specific level of demand each year of the planning period. Changes in this quantity are modeled by Considering α the load growth rate for each year, the amount of electricity consumption at each level of demand is calculated [31]. The demand for active and passive power is presented in Equations (2) and (3), and apparent power demand is in Equation (4). The superscript e shows that the parameters belong to the Monte Carlo experiment e.

Electricity Price Uncertainty
Depending on the amount of electricity consumption and the behavior of electricity market users, the price is determined for the electricity purchased from the primary grid. Based on Equation (5), this amount changes for the cost of electricity at the specific level of demand each year of the planning period. Changes in this quantity are modeled by multiplying two parameters, the base price of electricity for each level of demand and the price level index [35]-Equation (5).
ρ is the price level index that shows the behavior of the electricity market in the year t and the demand level of h. It is calculated using the probability density function as follows (Equation (6)):

Charging Station Modeling
Charging station power depends on the used battery technical specifications in the vehicle, such as the type and capacity, the amount of initial charge and the number of vehicles available at the charging station. Such parameters are not fixed and have high variability. Therefore, the location and operating conditions of the stations should be examined, and their characteristics should be considered in the model used. Electric vehicle owners' behavior is also the most influential factor in modeling these uncertain parameters compared to others discussed in this study. First, the modeling used for these parameters is discussed in this paper.

Type and Capacity of Vehicle Batteries
The companies that manufacture these batteries use different technologies, which has led to a variety of products. This diversity represents non-uniform consumers and makes the design of charging stations more complex. Here, the considered parameters are influenced by factors such as the amount of construction and development of batteries and each type's public acceptance rate in the station locations. In this paper, a 25-kWh battery is used, and its life is estimated at ten years.

The Initial Charge of Vehicle Batteries
The ratio of energy stored in a battery to its capacity is defined as the initial charge between 0 and 100. This parameter increases when the EV battery is being charged. SOC is reduced by consuming battery power (power delivery to the grid or vehicle). The battery storage capacity of the EVs is not used while they are charging at the charging station. Using this feature will depend on the initial charge of the EV's battery. Due to the variability of the distance traveled by EVs, their efficiency and the type of EV batteries, the amount of initial energy remaining in EV batteries will not be known. This paper considers the amount of initial charge of vehicles referring to the charging station at three levels in Table 3. Number of vehicles n 1 n 2 n 3

Vehicle Charging Schedule
By each vehicle at charging stations, the initial and required final charge and the departure time from the charging station are asked from the vehicle owners. Then, according to this information, by considering the profit from using the vehicles as storage and the cost of charging the batteries for driving, the optimal charging/discharging program for vehicles in the charging station is determined. The time required for full charge/discharge of vehicle batteries is calculated using the initial charge amount as follows [33]: SOC max and SOC min are the maximum and minimum flux values for the vehicle battery.
ES is the capacity of a vehicle battery and P v is the rate of power transition. In this section, the mathematical modeling of charge station power is discussed. Due to the uncertainties mentioned in the previous sections and the terms of them, the input and output power of the electric vehicle charging station will be following Figure 7 and Equations (12) and (13). and are the maximum and minimum flux values for the vehicle battery. is the capacity of a vehicle battery and is the rate of power transition. In this section, the mathematical modeling of charge station power is discussed. Due to the uncertainties mentioned in the previous sections and the terms of them, the input and output power of the electric vehicle charging station will be following Figure 7 and Equations (12) and (13).

Modeling the Location Problem
To model the placement problem in the first step, the placement goals must be specified. In this project, maximizing profits for the distribution system manager (charging station owner) is considered the optimization problem's goal. Profit from charging/discharging patterns of electric vehicles in the charging station, profit from reducing the active power losses in the distribution network, the cost of power supply through the upstream network and the investment cost for constructing the charging station as a profit function are considered. In the following, the average of the mentioned objective functions is discussed.

Profits from Discharge Programs
The battery capacity of electric vehicles in electrical energy storage enables charging stations to provide part of the network power at a lower price than the upstream network during peak times. It reduces the cost for the distribution system manager. The profit for using the power transmission technology from the vehicle to the network is calculated from Equations (10) and (11). As shown in Equation (12), the difference between revenue from power transmission to the network ( ) and the total cost for supplying discharge service ( ) is equal to the profit of the charging station and grid cooperation.

Modeling the Location Problem
To model the placement problem in the first step, the placement goals must be specified. In this project, maximizing profits for the distribution system manager (charging station owner) is considered the optimization problem's goal. Profit from charging/discharging patterns of electric vehicles in the charging station, profit from reducing the active power losses in the distribution network, the cost of power supply through the upstream network and the investment cost for constructing the charging station as a profit function are considered. In the following, the average of the mentioned objective functions is discussed.

Profits from Discharge Programs
The battery capacity of electric vehicles in electrical energy storage enables charging stations to provide part of the network power at a lower price than the upstream network during peak times. It reduces the cost for the distribution system manager. The profit for using the power transmission technology from the vehicle to the network is calculated from Equations (10) and (11). As shown in Equation (12), the difference between revenue from power transmission to the network (R discharge total ) and the total cost for supplying discharge service (C discharge total ) is equal to the profit of the charging station and grid cooperation.
c d is the cost of equipment depreciation due to the use of power injection technology from the vehicle to the network and τ i,h is the efficiency rate of charging station inverters.

Profits from Recharge Programs
Most electric vehicle drivers go to the charging station for their daily chores and short trips to recharge their batteries. The owners of the charging stations can increase their profit by providing charging services for vehicle drivers. The gain obtained from recharging vehicle batteries is calculated from the following equations. Same as the previous section, Equation (15) represents the profit of the recharge program, while Equations (13) and (14), respectively, show the total cost (C charge total ) and revenue (R charge total ) gained from supplying charging service for whole planning years.

Profits from Reduced Power Purchases from the Upstream Network
Since most of the power of the distribution network is supplied from the upstream network, by infiltrating the charging stations in the distribution system, this amount of power supply by the upper network can be reduced. The charging station reduces the cost of energy purchased, thereby increasing profits for the distribution system manager. The profit from this activity is shown in Equation (18), which equals the difference between cost savings and income, calculated from Equations (16) and (17).
As shown in Equation (19), the grid power, denoted by the uppercase grid, can be divided into three intervals. The first period is the network fertility period and is dl1 (park demand level). Naturally, during this period, the price of electricity is high, and the stations try to charge the batteries before this period arrives. By entering this period, they inject the batteries' energy into the network. The other interval, dl4 (low load demand l), is the low load interval of the network. This interval is the opposite of the previous interval and the decrease in demand reduces the price of electricity, and this is when station attendants generally charge the batteries. Finally, the intervals indicated by dl2 and dl3 (medium demand level) are placed between the previous two periods in terms of price and demand. The regular activity of stations and networks occurs. Usually, the electricity networks are designed to meet the needs of consumers quickly at these intervals.
One of the essential issues in the design of power plants and power generation centers is the ability of the network to respond to the peak interval. By looking at this subject from the electricity network viewpoint, this section can be justified as follows. Assuming stations are massive batteries, energy is stored in them during low load times. To supply energy during peak hours, instead of increasing the pressure on power plants, batteries meet part of the needs of consumers. By doing this, the network also benefits, which can be calculated using Equation (20). By injecting the energy stored in the electric vehicle battery into the grid, the losses are reduced. The profit from the sale of the stored electricity is increased. at each level of demand are with and without the charging station, respectively.

Numerical Simulation of the Proposed Method
Scattered product locating is a large-scale, multi-constrained issue that requires comprehensive search methods to resolve. Therefore, in the case of this type, intelligent techniques are used. Factors such as load uncertainty and electricity prices make this an optimization issue. In this project, locating the units is solved by the Genetic Algorithm, which is a well-known method to solve optimization problems. The random load and the price of electricity are also modeled by the Monte Carlo method. For evaluating the efficiency of the proposed method, planning has been performed on the 9-bus distribution network introduced in the reference. The studied network and its specifications are given in Figure 8 and Table 3. The test network has a reference bus with a voltage of 20 kV and eight load buses, where switches separate the distribution lines.
battery into the grid, the losses are reduced. The profit from the sale of the stored electricity is increased.
, , at each level of demand are with and without the charging station, respectively.

Numerical Simulation of the Proposed Method
Scattered product locating is a large-scale, multi-constrained issue that requires comprehensive search methods to resolve. Therefore, in the case of this type, intelligent techniques are used. Factors such as load uncertainty and electricity prices make this an optimization issue. In this project, locating the units is solved by the Genetic Algorithm, which is a well-known method to solve optimization problems. The random load and the price of electricity are also modeled by the Monte Carlo method. For evaluating the efficiency of the proposed method, planning has been performed on the 9-bus distribution network introduced in the reference. The studied network and its specifications are given in Figure  8 and Table 3. The test network has a reference bus with a voltage of 20 kV and eight load buses, where switches separate the distribution lines. All loads are selected as candidate bases for the installation of scattered production units. It should be noted that the time of entry and exit of vehicles to the charging station is randomly generated according to the distribution function according to the load curve. In Figure 6, the curve of the EV charging and discharging programs is presented.
The results are stated in Tables 4-6 according to the parameters defined in the article. Each table is provided with a specific number of electric vehicle parking spaces. Figure 9 also shows the change in the total cost of the optimization function at the time of optimization execution, which reaches its optimal value after a few periods. According to the mentioned tables, implementing load management programs reduces the cost of purchasing energy from the overhead network and reduces the overall costs of the network. The execution of load programs also reduces losses during peak hours. To compare the different scenarios in Tables 4-6, the effect of the vehicle load management programs and increasing the number of electric vehicles parking on different network costs is presented. Clearly, the increase in parking lots for electric vehicles will decrease the losses and the cost of purchasing from the overhead network. Still, in the case of the cost related to the All loads are selected as candidate bases for the installation of scattered production units. It should be noted that the time of entry and exit of vehicles to the charging station is randomly generated according to the distribution function according to the load curve. In Figure 6, the curve of the EV charging and discharging programs is presented.
The results are stated in Tables 4-6 according to the parameters defined in the article. Each table is provided with a specific number of electric vehicle parking spaces. Figure 9 also shows the change in the total cost of the optimization function at the time of optimization execution, which reaches its optimal value after a few periods. According to the mentioned tables, implementing load management programs reduces the cost of purchasing energy from the overhead network and reduces the overall costs of the network. The execution of load programs also reduces losses during peak hours. To compare the different scenarios in Tables 4-6, the effect of the vehicle load management programs and increasing the number of electric vehicles parking on different network costs is presented. Clearly, the increase in parking lots for electric vehicles will decrease the losses and the cost of purchasing from the overhead network. Still, in the case of the cost related to the construction of the parking lot, the total profit will fall. Additionally, by carefully looking at the tables, implementing cargo management programs makes electric vehicle parking and scattered production resources more economical. Table 4. Simulation results for the location of an electric vehicle park.

Load Information Without DRP With DRP
Bass number Optimal capacity of the power plant (kW) 3.0070 × 10 7 3.0890 × 10 7 Table 6. Simulation results for the location of three electric vehicle parking lots.

Load Information Without DRP With DRP
Bass number Optimal capacity of the power plant (kW) construction of the parking lot, the total profit will fall. Additionally, by carefully looking at the tables, implementing cargo management programs makes electric vehicle parking and scattered production resources more economical.   The location of the charging station is well represented in Tables 4-6 in the second rows according to the different DR scenarios, including maximum load supply profit, charge program benefit, profit from purchasing energy from the overhead network, profit from loss reduction, investment cost and finally total profit, shown in the section 'Numerical Simulation of the Proposed Method'. In addition to the positive effect of the proposed method on the cost of operation and the peak load, this method has a significant impact on improving the voltage profile and noise reduction. The proposed method can also solve the problems caused by load growth in the long run. Figure 10, presented the rate of active losses of the studied network in different scenarios.
Maximum load supply profit (USD) 2.7785 10 6 2.1116 10 6 Charge program benefit (USD) 1.1038 10 5 0.8388 10 5 Profit from purchasing energy from the overhead network (USD) 4.1348 10 7 3.8693 10 7 Profit from loss reduction (USD) 9.8129 10 5 1.4950 10 6 Investment cost (USD) 1.5139 10 7 1.1505 10 7 Total profit (USD) 3.0085 10 7 3.0883 10 7 The location of the charging station is well represented in Tables 4-6 in the second rows according to the different DR scenarios, including maximum load supply profit, charge program benefit, profit from purchasing energy from the overhead network, profit from loss reduction, investment cost and finally total profit, shown in the section 'Numerical Simulation of the Proposed Method'. In addition to the positive effect of the proposed method on the cost of operation and the peak load, this method has a significant impact on improving the voltage profile and noise reduction. The proposed method can also solve the problems caused by load growth in the long run. Figure 10, presented the rate of active losses of the studied network in different scenarios.  Figures 11 and 12 show the effect of the proposed method in reducing network losses with annual load growth for different years. It can be seen that the combined use of load management programs and parking locations together has significantly diminished both the maximum losses and the total energy wasted in the network.   Additionally, the proposed method can decrease the network peak and modify the network load curve in the first year and subsequent years and control the load growth well. According to Figure 12, it can be seen that the proposed method has reduced the  Figure 11. The amount of energy lost by the studied network in different scenarios.
tained for the two solar parking lots that had the best costs. The results show that with the increase in load in all cases, the network losses increase with a sharp slope, but solar parking can significantly inhibit load growth. Moreover, the proposed method utilization in the first year reduces about 13% of the peak load. About 4% of it is for the share of the load management program, and 9% is related to the share of charge/discharge control programs for electric vehicles. Additionally, the proposed method can increase the load more effectively during non-peak hours and improve the load factor.  Additionally, the proposed method can decrease the network peak and modify the network load curve in the first year and subsequent years and control the load growth well. According to Figure 12, it can be seen that the proposed method has reduced the network peak in the best possible way among the scenarios ( Figure 13). Results are obtained for the two solar parking lots that had the best costs. The results show that with the increase in load in all cases, the network losses increase with a sharp slope, but solar parking can significantly inhibit load growth. Moreover, the proposed method utilization in the first year reduces about 13% of the peak load. About 4% of it is for the share of the load management program, and 9% is related to the share of charge/discharge control programs for electric vehicles. Additionally, the proposed method can increase the load more effectively during non-peak hours and improve the load factor.  It is clearly shown in Figures 13 and 14 that regardless of the number of parking spaces of electric vehicles in all modes of the network, after a specific load growth, it will have a non-standard voltage curve in some buses. The existence of electric vehicles and the growth of other loads lead to improper operation. It is clear from these problems that using only load management programs cannot solve voltage drop and only delays it. In these figures, the bus voltage reached the unauthorized range in the last year. If electric vehicle charge and discharge control programs are used singly, it is observed that the voltage will remain within the allowable range, while the use of electric vehicle management programs and load management better prevents voltage drops due to load growth. Therefore, it can be concluded that cargo management programs, regardless of the number of parking spaces of electric vehicles, have the most excellent effect on improving the network voltage in different years. To better understand the effect of the proposed method on load growth, Figures 13 and 14 of the standard year and tenth year voltage curves are shown by considering the load growth in different scenarios. As shown, load management and locating the electric vehicle parking with the help of distributed generation sources have caused the voltage to remain within the allowable range even after ten years of load growth.
It is clearly shown in Figures 13 and 14 that regardless of the number of parking spaces of electric vehicles in all modes of the network, after a specific load growth, it will have a non-standard voltage curve in some buses. The existence of electric vehicles and the growth of other loads lead to improper operation. It is clear from these problems that using only load management programs cannot solve voltage drop and only delays it. In these figures, the bus voltage reached the unauthorized range in the last year. If electric vehicle charge and discharge control programs are used singly, it is observed that the voltage will remain within the allowable range, while the use of electric vehicle management programs and load management better prevents voltage drops due to load growth. Therefore, it can be concluded that cargo management programs, regardless of the number of parking spaces of electric vehicles, have the most excellent effect on improving the network voltage in different years. It is clearly shown in Figures 13 and 14 that regardless of the number of parking spaces of electric vehicles in all modes of the network, after a specific load growth, it will have a non-standard voltage curve in some buses. The existence of electric vehicles and the growth of other loads lead to improper operation. It is clear from these problems that using only load management programs cannot solve voltage drop and only delays it. In these figures, the bus voltage reached the unauthorized range in the last year. If electric vehicle charge and discharge control programs are used singly, it is observed that the voltage will remain within the allowable range, while the use of electric vehicle management programs and load management better prevents voltage drops due to load growth. Therefore, it can be concluded that cargo management programs, regardless of the number of parking spaces of electric vehicles, have the most excellent effect on improving the network voltage in different years.

Conclusions
This paper presents a long-term simultaneous planning period of electric vehicle solar charging stations in terms of load response and technical and economic indicators in the 10-year planning horizon based on a multi-objective genetic algorithm. In addition, Figure 14. Voltage curve of the studied network in different scenarios and 10-year load growth application considering three electric vehicle parking lots.

Conclusions
This paper presents a long-term simultaneous planning period of electric vehicle solar charging stations in terms of load response and technical and economic indicators in the 10-year planning horizon based on a multi-objective genetic algorithm. In addition, several different scenarios of the implementation of load management programs in the long-term horizon studied to reduce the cost of purchasing energy from the upstream network are presented, which reduced the total cost of the network. The implementation of load management programs has also reduced losses during peak network hours. To compare different scenarios, the effect of power management programs and also increasing the number of electric vehicle parking spaces on different network costs have been presented. It is observed that with the increase in the number of EV parking lots, losses and the cost of purchasing from the overhead network will decrease, but instead, the cost of constructing solar panels in the parking lot will reduce the total profit, which can be compensated in the long run. It was also clearly demonstrated from the numerical results that implementing the proposed load management programs makes the operation of PVCSs much more economical and reduces the adverse effects of charging EVs in the system.

Offers
For investors and future work, solar trees in real-time planning for the development of production and intelligent electrical systems and their integration with solar charging stations can be suggested.

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

Appendix A
The explanation of parameters and indices associated with equations: