Allocation of Radiation Shielding Boards to Protect the Urban Search and Rescue Robots from Malfunctioning in the Radioactive Environments Arising from Decommissioning of the Nuclear Facility

: A group of Urban Search and Rescue (USAR) Unmanned Ground Vehicles (UGVs) works in the city’s dangerous places that are caused by natural disasters or the decommissioning of nuclear facilities such as the nuclear power plant. Consider the multi sensor platform that the USAR UGV is equipped with, protecting the sensors from the danger that dwells in the working environments is highly related to the success of the USAR mission. The radioactive working environments due to the diffusion of the radioactive materials during the decommissioning of the nuclear facilities are proposed in this work. Radiation shielding boards are used to protect the sensors by installing them on the UGV. Generally, the more shielding boards installed on the UGV, the lower the probability of the sensor malfunctioning. With a large number of UGVs that are used in missions, a large number of shielding boards are needed. Considering the cost of the boards, a conﬂict between the ﬁnite budget and the best protection for the UGV swarm occurs. Due to differences in environments and working times, an asymmetric allocation of shielding boards might protect the UGV swarm better than that by allocating evenly. To minimize the absorption radiation doses, two methods including competitive equilibrium analysis and integer linear programming are proposed. The performances that result from the allocation of the shielding boards by two methods with different budgets are studied.


Introduction
Urban Search and Rescue (USAR) is defined as "the strategy, tactics, and operations for locating, providing medical treatment, and extrication of entrapped victims." [1] The robots, which are adopted by the USAR team, are a group of autonomous mobile robots equipped with sensors [2]. The USAR robots work in the dangerous environments that are caused by natural disasters [1,3] or the decommissioning of the nuclear facilities. Such as floods [4], fires [5], earthquakes [6], nuclear accidents and the decommissioning of the nuclear facilities [7][8][9][10][11][12]. Their duties include monitoring, information sharing, searching and rescuing [13,14].
On March 11, 2011, a massive earthquake and accompanying tsunami hit Japan. The Fukushima Daiichi Nuclear Power Plants sustained extensive damage in this disaster [8]. Three buildings that the increase in the payload. In addition, the working time would be decreased if more shields were installed. Generally, the RICA robot without hardening by shielding is a result of operations research considering that the working time is less than the minimum time the components could tolerate that. Although no specific shield was used when carrying out the 20 missions, the research group [43] considers a robot shield that can be used in more irradiation environments.
Shields not only protect robots from malfunctioning, but also decrease detector noises that arise from the irradiation background [44][45][46]. Lead is used to shield the photomultiplier tube (PMT)) mounted on the UGV that is used to detect the energy and the direction of gamma sources [47]. A tracked-wheel robot that incorporates stair-climbing capability is developed for nuclear power plant operations [48,49]. In this robot, cadmium selenium radiation resistant material shell is used to protect the camera and tungsten/silicone resin shielding strap is used to protect the internal electric components of the body.
Payload changes in different works due to different components that are equipped by the UGV, different missions that the UGV carries out and different budgets for the UGV. To carry tools for emergency stop valve manipulation in a nuclear power plant, a robot that could withstand a 100 kg payload is developed [50]. In response to the needs of the Chernobyl accident, a robot that is able to walk and climb could drag loads of 100 kg horizontally [51]. A robot for intervention, sampling and situation awareness in chemical, biological, radiological and nuclear (CBRN) missions is developed [52]. It is able to carry up to 40 kg of additional payload.
Considering the shield drawbacks including being a payload, increasing UGV weight, reducing mobility and the shield advantages including protecting electric components from damage, decreasing detector noises, operations research problems such as "should shields be used" and "how many shields would be used" are generated.
In this work, we build a digital city filled with radioactive materials which generate radioactive environments. We propose a USAR mission that is located in the digital city, covers a small area of 500 × 500 m 2 and is carried out by a small team of UGVs. There are in total 314 deployment points, which are allocated to 3 deployment zones. Due to a small area covered by the mission, the point that is not located in the area would not be covered by UGVs. The team consists of 6 UGVs. Each UGV is equipped with sensors [12,22,30], such as optical cameras, binoculars, laser range sensors and radiation detectors. Radiation shields are installed on the UGV to protect sensors and UGV itself. Due to the finite budget, we cannot allow every UGV to equip radiation shields as much as possible. We must balance the finite budget and the radioactive damage to all UGVs. Thus, we propose an assumption that the budget to purchase shields for all UGVs is finite and shields cannot satisfy every UGV. Our goal in this work is to find an allocation with which the the UGV swarm could reach a feasible average radiation damage and could survive in the environment. Due to differences in working environments and working times, allocating evenly is not the best method. On the contrary, an asymmetric allocation method might enable the UGV swarm to reach a lower average damage than allocating evenly. Thus, we propose two asymmetric allocation methods, each of which considers assumptions including a finite budget to purchase shielding boards that cannot satisfy every UGV, different radioactive environments that damage the UGVs and different working times of the UGVs.
Allocation of shielding boards is a subdomain of resource allocation. Similar works are studied, such as robot deployment in disaster city [4], the roadside healthcare facility location problem [53], the ambulance redeployment [54], the electric vehicle charging problem [55], the flow control in wireless sensor networks [56,57], the rate control for conmmunication [58], allocation for access channel [59], the charging and rate control in traffic [60]. Methods for the allocation include the mixed integer linear programming [4,53,54], the nonlinear programming and the competitive analysis [59]. Notice that competitive analysis is a method transferred from microeconomics. It belongs to nonlinear programming. Due to the concavity of the utility function, competitive analysis does not ask for numerical methods. Thus, the calculation for the competitive analysis is much quicker than that for other methods. The integer constraint does not fit well with competitive analysis.
The integer constraint is that the number of shielding boards is integer. The allocation does not allow a non-integer number, such as 2.6 shielding boards for the UGV. This generates a need for integer linear programming, in which mixed integer linear programming is a subdomain. Although the calculation could last for a much longer time, a better result could be found with integer linear programming. Thus, we propose two allocation methods including competitive analysis and integer linear programming (ILP). In addition to the number of shielding boards, the installation sequence of shielding boards also has impacts on radiation protection. Additionally, the installation sequence affects the weight balance of the UGV. To focus the study on radiation protection, we propose an installation sequence with at most 8 boards and we assume that this is the only one installation sequence that exists in our study.
In this paper, we face the radioactive damage from the electric devices of the UGV, as a serious problem for USAR UGVs that works in the radioactive environments. We propose a radiation protection method that is installing the radiation shielding boards on the UGV. γ-ray attenuations of shielding boards are simulated in the Monte-Carlo particle simulation toolkit "Geant4". According to the attenuation, we propose the Radiation Damage Minimization Problem. To solve the problem, we propose two methods, which are competitive equilibrium analysis from the microeconomics theory and integer linear programming. Assumptions that we consider in this work are listed as follows, • a finite budget to purchase shielding boards that cannot satisfy every UGV; • different radioactive environments that damage the UGVs; • different working times of the UGVs; • only γ-rays that would be shielded; • at most 8 shielding boards that each UGV could carry; • only one installation sequence of the boards that exists in our work; • a simple diffusion process that radioactivities of the wall and the ground are same; • the urban area only contains concrete walls and concrete grounds.
Due to the conflict between the finite budget and the best protection for the UGV swarm, we discuss the performance of the two methods. Works are listed below.

•
We propose a UGV structure and installation of the shielding board; • we evaluate the radiation damage on the UGV with different number of the shielding board; • we propose two methods to minimize the radiation damage for the UGV swarm; • we discuss the performance of the two methods and compare the malfunction ratios that result from the two methods in different budget.
The UGV structure, the installation of the shielding board, two radioactive environments, the radiation damage and the radiation protection are discussed in Section 2. The allocation methods including the competitive equilibrium analysis and the integer linear programming are proposed in Section 3. The simulation and analysis of the radiation damage to the UGV, the radiation protection of the shielding board and allocation results with the two methods are described in Section 4. Conclusions are drawn in Section 5.

Radiation Protection
USAR UGVs are the object that we should protect from radiation. We propose a radiation shielding board that is capable of absorbing radiation. In this study, we propose a shielding board of which the installation is shown in Figure 1. Each shielding board is 500-mm long, 50-mm wide and 3-mm thick. The more shielding boards installed, the lighter the radiation damage would be. Each UGV is able to equip 8 boards at most. Notice that the size of the board is an example to study the problem of shield allocation. Considering the payload of the UGV that is used in different purpose, shielding board should be altered in to a proper size. So dose the material. These also belong to the allocation problem, in which, shielding boards with different sizes and different materials are commodities. Similar methods that use equilibrium analysis or integer linear programming are feasible. The units are divided into 3 deployment zones. Each deployment zone manages a recharging site which is located at the perimeter of the radioactive environment. As a result of the finite battery capacity of the USAR robot, shift works of UGVs is used to realize sustainable working, such as the 3D reconstruction of the environment, the mapping of the radiation doses and the localization of the radioactive sources. Due to different working times, the radiation doses that UGVs absorb are different. This generates a problem of shielding board allocation to the UGV team. At a finite budget, the number of shielding boards is finite. If the number is not large enough to satisfy every UGV, we must find an allocation to minimize the average radiation damage. Thus, in this work, we study the allocation of shielding boards from both the perspective of a single UGV and the whole UGV team.

Damage Evaluation for the UGVs
The radiation dose is used to evaluate the radiation damage on the component of the UGV [7,8,29,30]. The studies show that the electric devices are vulnerable to the radiation. For instance, a laser scanner, a CCD camera, an inertial navigation system, a CAN-USB converter, a DC-DC converter, a matin controller, a wireless router and a wireless AP are found malfunctioning in the radiation tests. The absorption doses vary from 50 Gy to 300 Gy. The highest radiation doses could be up to 10 3 Gy/h [61]. If a UGV is working in this environment, it would malfunction within an hour.
Due to radiation protection methods include redundancy design [62,63] and radiation shielding, a malfunction of an electric device might not disable the UGV. In this work, we focus on radiation shielding. We propose a malfunction dose that might disables the UGV, denoted by D mal = 50 Gy. The parameters used in this work are listed in Table 1.
Due to the complexity of the four-level dominance hierarchy, we propose a two-element vector to index the USAR UGVs. It is denoted by t = [DZ i , k] T ≡ T. The first element is denoted by DZ i , i ∈ {1, 2, 3} = I. It represents the index of the deployment zone that the robot belongs to. The second element is denoted by k ∈ {1, · · · , Size DZ i }, in which, Size DZ i is the number of the UGVs in the deployment zone DZ i . It represents the index of the UGV in its deployment zone.
The absorption dose of all the electric devices in a UGV is denoted by D t . Noticing that the unit of dose is that 1 Gy = 1 J/kg, the radiation dose is not a function of mass. In this work, we consider silicon is the only part that is vulnerable to the radiation in the electric devices. Thus, the absorption dose of all electric devices D t is not the function of the devices' mass. The absorption dose is a function of the geometrical parameters of the devices. For instance, the size, the orientation and the arrangement of the devices. The effects of these geometrical parameters on the absorption dose are not discussed in this work. Thus, we propose a pure-silicon box representing all devices in a UGV. We study absorption doses of this box in the simulation toolkit Geant4.

Radiation Protection of UGVs
The radiation shielding board is able to share a small portion of the radiation doses that are affected by the incident particles. This portion is denoted by η. The more shielding boards, the lighter the radiation damage. To describe the damage on the UGV, we propose the malfunction ratio, denoted by MR t . It is a function of the absorption dose and the malfunction dose, as shown below, in which, MR t ∈ [0, +∞). The higher the malfunction ratio, the worse the UGV. When the ratio is higher than 1, the UGV malfunctions. The radiation shielding boards protect the UGV. The number of the boards and the installation of the boards have effects on the absorption dose. In this work, we propose an installation of the shielding boards, as shown in Figure 2. Thus, the number of the boards, which is denoted by N B,t , is the only parameter that affects the absorption dose, denoted by D t (N B,t ). We propose that the debris of the reactor fuel are the radioactive sources that disperses into the ground and the buildings. The UGV is irradiated when it passes by the contaminated ground or the buildings.
To study the relation of the absorption dose and the number of the boards, we simulate the irradiation from the radioactive contaminated material to the UGV in Geant4, which is a particle simulation toolkit [38]. Two environments are simulated. The first environment is that the UGV is located on the contaminated ground. The second environment is that the UGV sits beside a contaminated wall while the ground is also contaminated. A schematic view of these environments is shown in Figure 3. Due to a typical enrichment of the nuclear fuel that is used in the reactor is 3%, we propose a radioactive material that contains radioactive isotopes 235 U (3%) and 238 U.

Absorption Doses in Different Working Environments
The absorption dose of the UGV is affected by the working environments. We propose two environments in this work. The ground is radioactive in the first environment. Both the ground and the wall are radioactive in the second environment. The distributions of the absorption dose and the number of the shielding boards are different in two environments.
Each UGV has a specific working area. The environment where the working UGV is depends on the distance from the working area to the nearest building, denoted by dis t . We propose a threshold, denoted by th dis , to decide which environment that the UGV in. The index of the radioactive environment is denoted by Idx env ∈ {1, 2}.
We use the integer linear programming to allocate the shielding boards to the UGV. In which, the absorption dose is the object that should be minimized. According to the distribution of the absorption dose and the number of the shielding boards, the linear relation is fitted. We denote the linear relation for the UGV individual by, in which, the parameter D 0,Idx env (t) refers to the absorption dose when the UGV is not shielded. The parameter pDR(Idx env ) is the deceasing ratio of the dose. As a result that pDR(Idx env ) is negative, the more shields, the less radiation dose is received. We simulate the radiation absorption in the Monte-Carlo particle simulation toolkit Geant4. In which, γ-rays, which belong to photons, are emitted out from the concrete ground in the first environment Idx env = 1 or a combination of the concrete ground and the concrete wall in the second environment Idx env = 2. γ-rays would react with these concretes and deposit their kinetic energy to the concretes, as is noted radiation attenuation or self absorption if we consider the concrete as a part of the source. After the γ-ray flying out of the concrete in which it is emitted, it would react with objects in the atmosphere before hitting a UGV. The object includes air and other concrete walls if these walls lie exactly in its way to a UGV. This is an attenuation of the atmosphere. Due to a small cross-section, i.e., probability, of reactions between γ-ray and air, energy that is deposited by the γ-ray and is transferred to the air is relatively smaller than that occurring in the concrete. Then, the γ-ray hits a UGV and deposit its energy to the UGV. If the shielding board lied between the γ-ray and electric devices, energy of the γ-ray would be transferred firstly to the shield. Thus, electric devices are protected. In addition, due to reactions such as Compton scattering, the direction of the γ-ray would be altered after the γ-ray flying out of the shield. Thus, the γ-ray deviates from its original way that would hit and damage devices. Due to these processes that are simulated in Geant4, a UGV with two nearby radiation sources generally is more affected than if it is above one. If an electric device is hit by a γ-ray, reactions between a photon and silicon would be considered in Geant4. After the simulation, we retrieve energy depositions to the device. The energy depositions are used to estimate parameters, such as the absorption doses when the UGV is not shielded D 0,Idx env (t) and the deceasing ratio pDR(Idx env ).

The Competitive Equilibrium Analysis
We use the partial competitive equilibrium analysis in this work. It belongs to the competitive equilibrium analysis. The shielding boards equipped by the UGV are the commodities. Notice that the shielding board of each UGV is regarded as a unique commodity. Thus, the shielding boards in different UGVs belong to different commodities.
The working time and the working environment are the two parameters of the absorption dose. Due to the different delivering distances, the residual battery capacities of the UGVs while they are arriving are different. In order to prevent the UGV from running out of its battery, the longer the delivering distance, the shorter the working time. In the view of the radiation protection, the shorter the working time, the lighter the radiation damage. Thus, the allocation of the shielding boards should tend to the UGV with a long working time.
Two working environments are considered in this work. Due to the additional contaminated wall, the radioactive activity in the second environment is higher. Thus, the UGV working in the second environment would absorb more radiation dose than the UGV working in the first environment. In addition to the activity, an additional absorption dose may be caused by the photons emitted from the wall. This is the result of the Compton scattering of the photon with the shielding board. This phenomenon would be observed in the second environment. Thus, the second environment is more dangerous than the first environment.
To describe the two parameters, we propose the Cobb-Douglas utility function, denoted by, in which, the parameter t is the index of the UGV. The UGV in 3 deployment zones are included. The parameter N B is the commodity vector that consists of the number of the shielding boards of all UGVs, denoted by, The parameter α t considers the working time and the working environment of the UGV. To keep the constraint on α t , we introduce the parameterα t , denoted by, in which, the parameter pWT(t) refers to the working time.
in which, time(t) is the working time of the UGV t through the whole monitoring mission. time m is the time of the monitoring mission. In this work, the monitoring mission lasts for 6 h. The parameter pWE(t) represents the working environment.
Thus, the parameter α t is denoted by, The allocation of the shielding boards maximizes the utility function that we propose. The utility maximization problem (UMP) is as follows.
in which, the parameter p is the price of the shielding board. The parameter W is the budget of the expenditure on the shielding boards for all UGVs. If N * B is a solution to the UMP, then there exists a Lagrange multiplier λ such that, ∇u(N * B ) = λp (11) in which, p is the price vector of the boards, denoted by p = [p, · · · , p] T . The number of the elements in the price vector is The allocation of the shielding boards is denoeted by,

The Integer Linear Programming
We propose the integer linear programming for the shielding board allocation. The object is minimizing the total absorption dose that is the sum of the doses of all UGVs. We refer to the minimization problem as the Radiation Damage Minimization Problem, which is shown below.
In the Radiation Damage Minimization Problem, the first constraint specifies the available number of shielding boards on the UGV. Each UGV is allowed to equip 8 shielding boards at most, denoted by N B,t = 8. The UGV is also allowed to work without shielding, denoted by N B,t = 0. In addition to the minimum and the maximum of the number of the shielding boards, the number must be an integer. The shielding board is a product of the manufacture. If the board were cut into pieces, the properties such as the radiation attenuation would change. The function of the absorption dose would be invalid. In the view of the mechanical property, the shielding board with specific compounds is too hard. It can only be cut with specific tools, of which the preparing are not feasible in the safe zone. Thus, we do not want to cut any shielding board. The second constraint specifies the limitation for the budget. If the budget were not limited, the allocation would be easy. As a result, we would equip every UGV with 8 shielding boards.
The parameter D 0,Idx env is not a function of the number of the shielding boards N B , therefore the Radiation Damage Minimization Problem is equivalently defined as follows. (15) in which, P is the price vector, denoted by P = [p, · · · , p] T . C is the vector of the decreasing ratio pDR(Idx env ), denoted by,

Simulation of Radiation Damages and Shielding
Recall that we simulate radiation shielding for the UGV in two environments. In the first environment, only the ground exists. The ground is made of concrete. We propose that the ground is contaminated by the radioactive isotopes 235 U and 238 U. The enrichment of 235 U is 3%. Distributions of the gamma radiation of the two isotopes are shown in Figure 4. The distribution of the gamma radiation from the mixture of the both isotopes is also shown in the figure. In the second environment, there exist the ground and the wall. Both of them are made of concrete and both of them are contaminated. The enrichment is also 3%. Due to the real diffusion process of the radioactive material is complex and random, we assume a simple diffusion process that the radioactivity of the wall and the ground are the same. Relations between the absorption dose and the number of the shielding boards are shown in Figure 5. In the first environment, only the ground is contaminated. It means that the ground is the only radiation source that emits photons. Shielding boards that are located between the electric devices and the ground are the first 3 boards in the installation method that we propose. Thus, an obvious decreasing of the absorption dose is shown when the number of shielding boards increase from 1 to 3. In addition to the first 3 boards, the shielding boards are not located in the path of the radiation that is going to hit the electric devices. Generally they are useless. Thus, the absorption doses just vibrate in a small range, when the number of the installed boards increase from 4 to 8. A fitting result is shown in Figure 5. Decreasing ratio pDR(Idx env = 1) is −3.6 Gy.
In the second environment, both the ground and the wall are contaminated. An anomaly occurs when the first shielding board is installed. The absorption dose of the devices is larger than the dose with no shielding board installed. The cause of the anomaly is the scattered photon. Due to the Compton scattering, the incident photon that is emitted from the wall reacts with the first shielding board that is located under the electric device. The reaction is the Compton scattering. The scattered photon, which is one of the products of the reaction, is emitted out of the board and flies toward the electric device. After the photon gets into the device, it would transfer its energy to the device. Thus, the absorption dose increases. A fitting result is also shown in Figure 5. Decreasing ratio pDR(Idx env = 2) is −4.5 Gy. The absorption dose rates are 5.7 Gy/h in the first environment and 9.4 Gy/h in the second environment. Recall that the malfunction dose is 50 Gy. If the UGV was not shielded, it would keep working for about 9 h in the first environment until it is disabled. The UGV would keep working for about 5 h in the second environment. The monitoring work that is proposed in this work lasts for 6 h. The total working time of each UGV in this work is less than 5 h in both environments. The distributions of the working time of the UGV in three deployment zones are shown in Figure 6. The mean working times of the UGVs are 1.88 ± 0.37 h in DZ 1 , 1.86 ± 0.34 h in DZ 2 and 1.94 ± 0.36 h in DZ 3 . If the UGVs work in the first environment, the averaged absorption doses would be 10.7 Gy in DZ 1 , 10.6 Gy in DZ 2 and 11.1 Gy in DZ 3 . The malfunction ratio would be 21% in DZ 1 , 21% in DZ 2 and 22% in DZ 3 . The radiation damages on the UGVs in the second working environment are listed in Table 2.

Radiation Damages on the UGV without Shielding
In this work, we propose two working environments where the UGV works. The working environment of a UGV depends on the distance from the UGV to the nearest wall. Considering that the UGV works in the urban area where concrete walls and concrete grounds are dominant components in urban areas, we make an assumption that the urban area only contains concrete walls and concrete grounds. This assumption avoids discussion about radioactive effects on other components such as trees, soils and lakes. Based on this assumption, we provide a simple structure of urban areas. A map of the buildings and the UGVs from 3 deployment zones is shown in Figure 7. By setting the threshold of the distance th dis = 200 m, working environments in the whole digital city are shown in Figure 8. 208 deployment points are in the first environment and 106 deployment points are in the second environment. If all deployment points were covered, we would use 1032 UGVs. In which, 678 UGVs work in the first environment and 354 UGVs work in the second environments. Numbers of the UGVs that work in different environments are listed in Table 3. Recall that the USAR mission just covers a small area of 500 × 500 m 2 and is carried out by a team of 6 UGVs. Only the mean value and the standard deviation in Table 3 are used.  Considering the absorption dose rates, the working times and the malfunction dose D mal , distributions of the malfunction ratio MR t for the UGV without shielding are shown in Figure 9.
In the circumstances of no shielding, all UGVs function well when finishing the 6-h monitoring work, because all the malfunction ratios are less than 100%. Most of the malfunction ratio lies in the range of [20%, 30%]. Specifically, the average malfunction ratios, denoted byMR t , are 25.5% ± 7.9% in DZ 1 , 26.5% ± 8.3% in DZ 2 and 27.0% ± 8.5% in DZ 3 .
Notice that the absorption dose D t for the UGV with only 1 shielding board decreases slightly than that without shielding in the first environment. Even worse, the dose increases in the second environment. We evaluate the malfunction ratio for the UGV that is equipped with only one shielding board. The evaluation results, which are listed in Table 3, show that the average malfunction ratioMR t for the UGV with only one shielding board is larger than that without shielding board. We conclude that the UGV equipped with only 1 shielding board is useless or harmful to the UGV itself. Thus, we remove the installation of only one shielding board from the candidates of the shielding board installation.

Allocation of the Shielding Board
Before the allocating, we determine the parameters including D 0,Idx env , pDR, pWT and pWE. The parameter D 0,Idx env is the absorption dose when the UGV is not shielded. According to the absorption doses that are simulated in this work, as shown in Figure 5, the absorption dose of the UGV without shielding D 0,Idx env =1 is 34.0 Gy and D 0,Idx env =2 is 56.6 Gy. Fitting results shows that the decreasing ratio pDR(Idx env = 1) is −3.6 Gy and pDR(Idx env = 2) is −4.5 Gy. A distribution of the pWT is shown in Figure 6. A schematic view of the working environment pWE is shown in Figure 8. In which, each dot represents a deployment point. Different color refer to different working environment. We suppose that the price of each shielding board p is 100.
Two allocation methods are proposed in this work. They behave quite differently in shielding board allocation. We take an example of shielding board allocation to 6 UGVs. The allocations of the shielding board under different budgets are shown in Figure 10. With the partial competitive equilibrium analysis, which is the first method in this paper, the ratio of the number of the shielding boards between the UGVs is fixed, denoted by (N B,t 1 : · · · : N B,t 6 ) = (α t 1 : · · · : α t 6 ). The integer linear programming method considers the UGV ID = 1 as the most vulnerable UGV in the field, the shielding boards are allocated to the UGV since the budget is 600. It is the only UGV to which the equipment of the shielding boards is completed. Recall that the equipping of only one shielding board just negligibly decrease the absorption dose of the UGV in the first environment, or even worse, it increases the dose in the second environment. We eliminate the installation of only 1 shielding board from the candidate set of the shielding board installation. Due to the eliminating of the installation candidate, the allocation methods need to be adjusted. The allocations are shown in Figure 11. The malfunction ratios MR t are shown in Figure 12. The absorption dose for the UGV decreasing significantly along the increasing of the budget. Due to the fixed ratio of the number of the shielding boards between the UGVs, this phenomenon could be easily recognized from the result of the competitive equilibrium analysis. The result of the integer linear programming shows that the UGVs are protected better at the lower budget than the shielding allocation with the competitive equilibrium analysis. Specifically, the average malfunction ratios at the smallest budget are 31.2% ± 6.9% with the competitive equilibrium analysis and 26.2% ± 5.5% with the integer linear programming. The effects of the radiation protection are inverted at the largest budget. The average malfunction ratios at the smallest budget are 12.9% ± 3.9% with the competitive equilibrium analysis and 13.0% ± 4.3% with the integer linear programming.
Distributions of the average malfunction ratio are shown in Figure 13. Horizontal axis represents an average budget on each robot. Consider the UGV team, total budget ranges from 600 to 4200 and the average budget ranges from 100 to 700. Specifically, the budget W ∈ {116 × l : l = 1, · · · , 6} is studied. Each wealth lies in a budget interval, of which the upper limit satisfies that each UGV could be equipped with a same number of the shielding board. Recall that the total number of the UGVs is 6. The upper limit that every UGV is equipped with l boards is {116 × 1, · · · , 116 × l, · · · , 116 × 8} ≡ L. The budget interval is denoted by BI l = (100 × (l − 1), 100 × l]. The set of the budget interval is denoted by BI = {BI l ; l ∈ L}.
Allocation results, as shown in Figure 13, indicate that the allocation of the shielding boards with the competitive equilibrium analysis has a better performance in radiation protection than that with the integer linear programming when the average budget is larger than 350. In addition to expectations of the ratio, the maximum malfunction ratio that results from the competitive equilibrium analysis is also lower. For instance, the maximum malfunction ratio that results from the competitive equilibrium analysis is 18.6% for the average budget of 580, while the maximum malfunction ratio that results from the integer linear programming is 29.5%. On the contrary, the result with the integer linear programming is better when the average budget is no more than 350. This phenomenon indicates that the average malfunction ratios resulting from competitive equilibrium analysis converge faster than that resulting from integer linear programming. Figure 13. Distributions of the average malfunction ratio. Horizontal axis represents an average budget on each robot. The error bar refers to the standard deviation of the malfunction ratio. The situation allocating only 1 shielding board to the UGV is eliminated. The team consists of 6 UGVs. The term CD refers to the competitive equilibrium analysis and the term ILP refers to the integer linear programming.

Conclusions
In this work, we face radioactive damage to the electric devices in the UGV. This is a serious problem for the Urban Search and Rescue robot that works in a radioactive environment. For instance, the UGVs that work for the decommissioning of the nuclear power plant suffer from radioactive damage. We study radiation protection methods by installing radiation shielding boards on USAR UGVs. Considering the cost of the shielding board, a conflict between the finite budget and the best protection for a UGV swarm occurs. This is a problem of operations research. To allocate the shielding boards, we propose the Radiation Damage Minimization Problem. To solve the problem, we propose two methods, which are competitive equilibrium analysis from microeconomics theory and integer linear programming. The finite budget is considered by the constraint of both methods.
To evaluate the radioactive damage, we simulate the γ-ray attenuation in the Monte-Carlo particle simulation toolkit Geant4. Two radioactive environments are considered. One simulates field work with the ground which is radioactive contaminated. The other one simulates field work with contaminated ground and contaminated walls. We propose that each UGV could be equipped with 8 shielding boards. There are 6 UGVs to be protected in our work.
Results of the shielding board allocation indicate that the competitive equilibrium analysis shows better performances when the budget is larger than 350, and the integer linear programming is better when the budget is no more than 350. The performance consists of the average malfunction ratio and the maximum malfunction ratio of the UGV swarm. In addition to the performance, we find that the average malfunction ratios resulting from the competitive equilibrium analysis converge faster than that resulting from the integer linear programming.
A future research direction is inter-temporal resource allocation. Consider the fact that a UGV does not need radiation protection while it is recharging in a recharge site where there is no threat of radiation damage. A smaller number of shielding boards would satisfy UGVs, if the shielding board on the recharging UGV is retrieved and reallocated to the UGV which is going to work. This allocation concerning times is denoted by inter-temporal resource allocation. We would need an inventory to retrieve, store and reallocate shielding boards. When the inventory is considered, it is feasible to study the depreciation of the boards. Thus, we are going to study how many shielding boards would be saved, how much budget would be saved and what the depreciation would be, if we were using inter-temporal resource allocation.