JP2009025122A - Cable conductor temperature estimation method in consideration of latent heat effect in tunnel, cable conductor temperature estimation system, and cable conductor temperature estimation program - Google Patents

Abstract

<P>PROBLEM TO BE SOLVED: To provide a system or the like capable of estimating a change of a conductor temperature of a power cable laid in an underground tunnel, in consideration of a latent heat effect by evaporation of water in the tunnel. <P>SOLUTION: A heating value of a cable 4, a heat transfer quantity between an air contact surface such as the inner surface of a tunnel constitution wall or the surface of the cable 4 and the air in the tunnel 1, and a latent heat transportation quantity caused by evaporation of water on a portion wetted by water caused by submergence or the like on the air contact surface are calculated on an analysis face crossing the tunnel 1 wherein the cable 4 is laid. Then, the change of the cable conductor temperature is estimated by using a finite element method based on a heat conduction equation, by using the values, a thermal constant at each portion, an evaporation efficiency of water, or the like. <P>COPYRIGHT: (C)2009,JPO&INPIT

Description

  The present invention relates to a method and system for estimating a conductor temperature of a power cable disposed in the ground, particularly a cable conductor temperature by a finite element method in consideration of a latent heat effect in a sinus where a cable is disposed, The present invention relates to a program for executing a system using a computer and belongs to the field of temperature analysis technology using a finite element method.

  When installing a power cable in the ground, it is necessary to set the type and size of the cable so that the cable conductor temperature does not exceed the predetermined allowable temperature when the target current is passed through the cable. As a guideline, a power cable allowable current calculation manual (JCS168E) is provided by the Japan Electrical Wire Manufacturers Association.

  However, in this manual, the thermal constants such as the heat capacity and thermal conductivity of the soil around the cable are fixed values set on the safe side, so the calculated value of the cable conductor temperature is higher than the actual value. As a result, there was a waste such as limiting the energizing current or increasing the size of the cable more than necessary.

  In response to such problems, in recent years, more accurate estimation of the conductor temperature has been attempted for the purpose of efficient operation of the power cable, etc. An invention relating to a new estimation method has been proposed.

  Among these, the invention disclosed in Patent Document 1 calculates the heat generation amount due to the energization current of the cable, the heat generation amount, the heat constant of the preset cable conduit and the surrounding soil, and the heat balance amount on the ground surface. Based on the above, by simulating the cable conductor temperature at a specific point using the finite element method, it is possible to accurately estimate the temporal change of the conductor temperature.

  Further, the invention disclosed in Patent Document 2 simulates the cable conductor temperature in consideration of the cooling effect in the case where the cooling water pipes are arranged in parallel with the cable in addition to that of Patent Document 1. Further, in the invention disclosed in Patent Document 3, when a cable is laid in a cave provided in the ground, a cable conductor is considered in consideration of a change in temperature due to air movement in the cave. The temperature is simulated.

JP 2004-112964 A JP 2006-038631 A JP 2006-250689 A

  As described above, according to the invention disclosed as the background art, it is expected that the cable conductor temperature is estimated with high accuracy under each condition, but in order to further improve the estimation accuracy, In the case where the cable is laid in the cave, it is necessary to consider the influence of water evaporation in the cave on the cable conductor temperature.

  In other words, underground water may ooze from the floor and side walls of the cave, and external water may enter through the cave's entrances and vents, and vaporize from the cave when the water evaporates. Heat will be taken away as latent heat, but this will affect the temperature in the sinus and thus the cable conductor temperature. In addition, when a cooling water pipe is provided in the cave, water in the pipe can be sprayed into the cave and the inside of the cave can be actively cooled by evaporation thereof. It is necessary to estimate the effect.

  Accordingly, an object of the present invention is to further improve the estimation accuracy of the conductor temperature by estimating the cable conductor temperature in consideration of the latent heat effect accompanying water evaporation in the sinus.

  In order to solve the above problems, the present invention is configured as follows.

  First, the invention according to claim 1 of the present application is a method for estimating the conductor temperature of a power cable laid in an underground cave by using a finite element method by a computer, the cable, the air in the cave and Using the thermal constants of the components around the cave, the heat transfer coefficient between the air contact surface in the cave and the air in the cave, and the evaporation related values related to the evaporation of water at the air contact surface in the cave, Calculating the amount of heat transfer between the air contact surface and the air and the amount of latent heat transport in a predetermined region of the air contact surface, and calculating the amount of heat transfer from the heat transfer amount and the amount of latent heat transport and the amount of heat generated by the cable. It is characterized by estimating a change in the conductor temperature of the conductor.

  Here, examples of the components around the cave include a cave wall for forming a cave in the ground, soil around the cave, and the thermal constant, for example, specific heat capacity of each component And evaporation-related values include, for example, latent heat of vaporization taken from the surroundings when water evaporates, water evaporation efficiency in the latent heat transport amount calculation region, and the like.

  Further, the invention according to claim 2 is the method according to claim 1, wherein the temperature and humidity of the air in the cave used for calculating the heat transfer amount and the latent heat transfer amount take air movement in the cave into consideration. And a value reflecting the upstream temperature and humidity is used.

  On the other hand, in the invention described in claim 3, the conductor temperature of the power cable laid in the underground cave is used as an analysis plane, and the finite element method is applied to a predetermined region of the analysis plane. System to estimate the configuration of the inside and outside of the cave, including the cable layout, and an area to consider the latent heat transport due to water evaporation in the air contact surface in the cave Latent heat transport area setting means, thermal constant setting means for setting the thermal constants of the cable, air in the cave and surrounding components, energizing current setting means for setting the energizing current to the cable, and air contact in the cave A heat transfer coefficient setting means between the surface and the air in the cave, an evaporation related value setting means for setting an evaporation related value related to the evaporation of water on the air contact surface in the cave, and the energizing current setting means Based on the set energizing current A heat generation amount calculating means for calculating the heat generation amount of the bull, and an air transfer for calculating a heat transfer amount between the air contact surface in the cave and the air in the cave using the heat transfer coefficient set by the heat transfer coefficient setting means. A calorific value calculating means; and a latent heat transport amount calculating means for calculating a latent heat transport amount due to water evaporation in the latent heat transport area set by the area setting means using the evaporation related value set by the evaporation related value setting means; Based on the configuration set by the sinus configuration setting unit, model creation means for creating an analysis model by dividing a predetermined analysis region of the analysis surface crossing the sinus into finite elements, and recorded in the thermal constant recording unit Matrix generating means for creating a heat capacity matrix and a heat conduction matrix for each node in the analysis region using the thermal constants, and the heat generation amount of the cable calculated by the heat generation amount calculation means, the air transfer Using the heat transfer amount between the air contact surface in the cave and air calculated by the amount calculation means and the latent heat transfer amount calculated by the latent heat transfer amount calculation means, a thermal load for each node in the analysis region Vector creating means for creating a vector, and each nodal temperature in the analysis region at a predetermined time interval from a predetermined initial state using a heat capacity matrix, a heat conduction matrix and a thermal load vector created by the matrix creating means and the vector creating means And a nodal temperature calculating means for calculating.

  Here, the setting of various data by each of the setting means is carried out by reading data recorded in advance in the recording means, inputting data at the time of analysis, using them together, and further using these data. There is a case.

  According to a fourth aspect of the present invention, in the system of the third aspect, after calculating the latent heat transport amount by water evaporation by the latent heat transport amount calculating means, the humidity of the air in the sinus after the evaporation is calculated. Humidity update means for calculating and updating the humidity value is provided, wherein the latent heat transfer amount calculation means calculates the latent heat transfer amount using the humidity updated by the humidity update means. .

  The invention according to claim 5 is the system according to claim 3 or claim 4, wherein the analysis surface setting means for setting a plurality of analysis surfaces for finite element analysis in a predetermined section of the sinus, Air data setting means for setting data relating to the movement of air in the cave, wherein the air heat transfer amount calculation means sets the air temperature set by the air data setting means as the air temperature used in the heat transfer amount calculation. Depending on the state of movement, the air temperature obtained on the upstream analysis surface is used, and the latent heat transport amount calculating means is set by the air data setting means as the air temperature and humidity used in calculating the latent heat transport amount. The air temperature and humidity on the upstream analysis surface are used according to the state of air movement.

  The invention according to claim 6 is set in the system according to claim 3 or claim 4 by a cooling water pipe setting means for setting a cooling water pipe in the sinus and the cooling water pipe setting means. The heat transfer coefficient setting means between the inner surface of the cooling water pipe and the cooling water in the pipe, and the inner surface of the cooling water pipe and the pipe using the heat transfer coefficient set by the heat transfer coefficient setting means Cooling water heat transfer amount calculating means for calculating the heat transfer amount between the cooling water in the cooling water, the thermal constant recording means records the heat constant of the cooling water pipe and the cooling water in the pipe, The cave configuration setting means sets the configuration inside and outside the cave including the configuration of the cooling water pipe, and the vector creating means is configured to calculate the cooling water pipe inner surface and the cooling water calculated by the cooling water heat transfer amount calculation means. A heat load vector is created using the amount of heat transferred between the two.

  Here, the configuration of the cooling water conduit includes the number, diameter, arrangement position, and the like. When the nodal temperature calculating means starts calculation, initial values of the cooling water pipe line and the cooling water temperature are given as the initial state.

  The invention according to claim 7 is the system according to claim 6, wherein the analysis surface setting means for setting a plurality of analysis surfaces for finite element analysis in a predetermined section of the sinus, and the air in the sinus Air data setting means for setting data relating to the movement of the cooling water, and cooling water data setting means for setting data relating to the movement of the cooling water in the cooling water pipe, wherein the air heat transfer amount calculating means is configured to calculate the heat transfer amount. As the air temperature used in the above, the air temperature obtained on the upstream analysis surface is used according to the state of air movement set by the air data setting means, and the latent heat transport amount calculating means calculates the latent heat transport amount. As the air temperature and humidity used at the time, according to the air movement state set by the air data setting means, the air temperature and humidity on the upstream analysis surface are used to calculate the cooling water heat transfer amount. The stage uses the cooling water temperature obtained on the upstream analysis surface according to the state of movement of the cooling water set by the cooling water data setting means as the cooling water temperature used in calculating the heat transfer amount. Features.

  The invention according to claim 8 is the system according to any one of claims 3 to 7, wherein the meteorological data recording means for recording the meteorological data of each period in each region, and the designated analysis target area And a heat balance amount calculating means for reading out the weather data of the season from the weather data recording means and calculating a heat balance amount of the ground surface based on the weather data, and the vector creating means calculates the heat balance amount The thermal load vector is created using the ground surface heat balance calculated by the means.

  On the other hand, in the invention described in claim 9, the conductor temperature of the power cable laid in the underground cave is used as an analysis plane, and the finite element method is applied to a predetermined region of the analysis plane. A path configuration setting means for setting a configuration inside and outside the sinus including the arrangement of the cable, and an area in consideration of latent heat transport due to water evaporation in the air contact surface in the sinus Latent heat transport area setting means for setting the cable, heat constant setting means for setting the thermal constants of the air in the cave and the surrounding components of the cave, energizing current setting means for setting the energizing current to the cable, the air contact surface in the cave The heat transfer coefficient setting means between the air and the air in the cave, the evaporation related value setting means for setting the evaporation related value related to the evaporation of water on the air contact surface in the cave, and the energization current setting means Based on energizing current A heat generation amount calculating means for calculating a heat generation amount of the cable, and an air transfer for calculating a heat transfer amount between the air contact surface in the cave and the air in the cave using the heat transfer coefficient set by the heat transfer coefficient setting means. A latent heat transport amount calculating means for calculating a latent heat transport amount by water evaporation in a latent heat transport region set by the region setting means, using the evaporation related value set by the evaporation related value setting means; On the basis of the configuration set by the road configuration setting means, model creation means for creating an analysis model by dividing a predetermined analysis region of the analysis surface crossing the cave by finite element, heat recorded in the thermal constant recording means Matrix creation means for creating a heat capacity matrix and a heat conduction matrix for each node in the analysis region using constants, the heat generation amount of the cable calculated by the heat generation amount calculation means, and the air heat transfer amount The heat load vector for each node in the analysis region using the heat transfer amount between the air contact surface in the cave calculated by the exit means and the air and the latent heat transport amount calculated by the latent heat transport amount calculation means A vector generating means for generating the temperature, and each nodal temperature in the analysis region at a predetermined time interval from a predetermined initial state using a heat capacity matrix, a heat conduction matrix and a thermal load vector generated by the matrix generating means and the vector generating means. It is made to function as a nodal temperature calculation means for calculating.

  The invention according to claim 10 is the program according to claim 9, wherein the computer calculates the latent heat transport amount by water evaporation by the latent heat transport amount calculating means, and then the air in the sinus after the evaporation is calculated. Function as a humidity update unit that calculates the humidity value and updates the value of the humidity, and when functioning as the latent heat transfer amount calculation unit, the latent heat transfer amount is calculated using the humidity updated by the humidity update unit. It is made to function like this.

  The invention according to claim 11 is the program according to claim 9 or 10, wherein the computer sets an analysis surface for setting a plurality of analysis surfaces for finite element analysis in a predetermined section of the sinus. And air data setting means for setting data relating to air movement in the cave, and when functioning as the air heat transfer amount calculation means, the air data is used as the air temperature used in the heat transfer amount calculation. When functioning to use the air temperature obtained on the upstream analysis surface according to the air movement state set by the setting means and functioning as the latent heat transport amount calculating means, the latent heat transport amount is calculated. As the air temperature and humidity used in the above, the air temperature and humidity on the upstream analysis surface are used according to the air movement state set by the air data setting means. Characterized in that to function so that.

  The invention according to claim 12 is the program according to claim 9 or claim 10, wherein the computer includes a cooling water pipe setting means for setting a cooling water pipe in the sinus and the cooling water pipe setting means. A heat transfer coefficient setting means between the set inner surface of the cooling water pipe and the cooling water in the pipe, and an inner surface of the cooling water pipe using the heat transfer coefficient set by the heat transfer coefficient setting means When functioning as the cooling water heat transfer amount calculating means for calculating the heat transfer amount between the cooling water in the pipe and the heat constant recording means, the cooling water pipe and the cooling water in the pipe are When functioning to record the thermal constant and functioning the sinus configuration setting means, when functioning to set the configuration inside and outside the sinus including the configuration of the cooling water pipeline and functioning as the vector creation means Is the cooling water heat transfer calculation Characterized in that to function so as to create a thermal load vector using the amount of heat transfer between the calculated cooling water lines the inner surface and the cooling water by means.

  The invention according to claim 13 is the program according to claim 12, wherein the computer sets an analysis surface setting means for setting a plurality of analysis surfaces for finite element analysis in a predetermined section of the sinus, When functioning as air data setting means for setting data relating to air movement in the air and as cooling water data setting means for setting data relating to movement of cooling water in the cooling water pipe, and as the air heat transfer amount calculating means The function of using the air temperature obtained on the upstream analysis surface in accordance with the air movement state set by the air data setting means as the air temperature used when calculating the heat transfer amount, When functioning as the transport amount calculating means, the air temperature and humidity used for calculating the latent heat transport amount are set by the air data setting means. According to the state of air movement, when functioning to use the air temperature and humidity in the upstream analysis surface, and functioning as the cooling water heat transfer amount calculating means, as the cooling water temperature used when calculating the heat transfer amount, According to the cooling water movement state set by the cooling water data setting means, the cooling water temperature obtained on the upstream analysis surface is used.

  According to a fourteenth aspect of the present invention, in the program according to any one of the ninth to thirteenth aspects, the computer includes a meteorological data recording means for recording meteorological data for each period in each region, and specified. The meteorological data of the analyzed area and time is read from the meteorological data recording means, and functions as a heat balance amount calculating means for calculating the heat balance amount of the ground surface based on the meteorological data, and also functions as the vector generating means. In some cases, the heat balance vector is made to function so as to create a thermal load vector using the heat balance amount of the ground surface calculated by the heat balance amount calculation means.

  According to the invention described in each claim of the present application, the following effects can be obtained.

  First, according to the method of claim 1 of the present application, when estimating the conductor temperature of the power cable laid in the underground cave, in addition to the heat generation of the cable and the configuration inside and outside the cave, It is estimated in consideration of the latent heat effect due to water evaporation on the air contact surface in the sinus.

  In other words, the floors and side walls in the cave may be wet with water that has exuded from the ground, water may enter from the outside, or may be actively sprinkled for cooling. As the water evaporates, the heat in the cave is taken away as latent heat of vaporization. According to the method of claim 1, the conductor temperature is calculated in consideration of this latent heat effect. The simulation of the temperature change is performed with high accuracy. As a result, it is possible to avoid the waste of suppressing the energization current more than necessary or increasing the size of the cable more than necessary, which contributes to the improvement of the efficiency of power transportation.

  Further, according to the method of claim 2, in the method of claim 1, the temperature and humidity of the air in the cave that influences the amount of heat transfer and latent heat transport between the air contact surface in the cave and the air in the cave Since the calculation is performed in consideration of the movement of the air, the simulation of the conductor temperature including the latent heat effect is performed with higher accuracy.

  On the other hand, according to the system described in claim 3, the method described in claim 1 is specifically executed by using various hardware resources such as a central processing unit, a recording device, and an input / output device of a computer. By estimating the cable conductor temperature in consideration of the amount of latent heat transport due to water evaporation on the air contact surface in the road, the change can be simulated quickly and accurately.

  In that case, the system according to claim 4 pays attention to the fact that the humidity of the air in the cave required for calculation of the latent heat transport amount changes due to the evaporation of water, and the latent heat transport amount is constantly updated. By using humidity, the simulation is performed with higher accuracy.

  According to the system of claim 5, in the system of claim 3 or claim 4, a plurality of analysis surfaces are set in the sinus, and the temperature of the air in the sinus due to the movement of air between adjacent analysis surfaces. Since the simulation of the conductor temperature on each analysis surface is performed while passing the humidity, the effect of the air movement is correctly reflected and the simulation is performed with higher accuracy.

  Further, according to the system of the sixth aspect, the change in the cable conductor temperature in the case where the cooling water pipe is disposed in the cave is simulated using a computer, so that the cooling water pipe can be effectively arranged. Designing can be done easily and quickly.

  According to the system of claim 7, in the system of claim 6, a plurality of analysis surfaces are set in the sinus, and the temperature and humidity of the air in the sinus due to the movement of air between adjacent analysis surfaces are determined. Conducting the simulation of the conductor temperature on each analysis surface while delivering the cooling water temperature by transferring and moving the cooling water, the effect of the movement of air and cooling water is correctly reflected, and the simulation is performed more accurately It will be.

  Furthermore, according to the system of claim 8, since the cable conductor temperature is calculated including the heat balance on the ground surface, the simulation is performed in a state closer to reality.

  On the other hand, according to the program described in claims 9 to 14, a system similar to the system described in claims 3 to 8 is configured by installing these in a computer. The same effect as these systems is realized.

  Embodiments of the present invention will be described below. In addition, although the following embodiment is about the cable conductor temperature estimation system which concerns on this invention, the method used in this system and its program are the cable conductor temperature estimation method and cable conductor temperature estimation which concern on this invention. An embodiment of the program is configured.

  First, the configuration of a cave model to be simulated will be described with reference to FIG. 1. A cave 1 is provided in the ground, its periphery is surrounded by a wall surface 2, and the outside is covered with soil 3. . And the power cable 4 is laid inside. In addition, an air intake device 5a for taking in air from the ground, an exhaust device 5b for discharging air to the ground, and a blower device 5c for forcibly moving the air in the sinus 1 are disposed in the sinus 1. Internal forced ventilation is possible.

  Further, in the sinus model, a cooling water conduit 6 is disposed in the sinus 1. As shown in FIG. 2, the cooling water pipe 6 extends in one direction from the cooling water supply device 7 and then makes a U-turn at a predetermined position and returns to the cooling water supply device 7 again. The cooling water at a constant temperature (inlet temperature) sent from the supply device 7 to the inlet of the forward passage (hereinafter referred to as “Go pipe”) of the cooling water pipe 6 is increased in temperature by the ambient heat while returning to the return path. (Hereinafter referred to as “Re pipe”), the water is returned to the cooling water supply device 7, cooled to the inlet temperature by the device 7, and then supplied again to the Go pipe inlet.

  In the present embodiment, it is possible to simulate the cable conductor temperature even for a model in which the cooling water pipe 6 is not provided.

  Next, the configuration of a system for performing the above-described model setting or cable conductor temperature analysis will be described.

  As shown in FIG. 3, a computer 10 constituting this system includes a central processing unit 11, a central processing unit 11, an input device 12 used for setting various conditions and controlling the system, and a recording medium such as a CD-ROM. 20, a reading device 13 for reading information such as programs and various data, a recording device 14 for recording programs and data read by the device 13, and further calculation results, and a display device for displaying input screens and calculation results 15 and a printing device 16 for printing calculation results and the like. Further, the central processing unit 11 is provided with a memory 17 for temporarily storing data during processing.

  The recording device 14 includes a soil thermal constant database DB1, a cave wall thermal constant database DB2, a pipe wall thermal constant database DB3, a fluid thermal constant database DB4, a pipe dimension database DB5, a cable database DB6, an air physical property database DB7, An evaporation efficiency database DB8, a cooling water property database DB9, a ground temperature database DB10, a weather database DB11, and a load factor database DB12 are recorded.

  If the structure of these databases is demonstrated in order, first, as shown in FIG. 4, soil heat constant database DB1 will record the name, mass density, specific heat, and thermal resistance for every kind of soil of an analysis object area. It is like that.

  Further, as shown in FIG. 5, the name, mass density, specific heat, and thermal resistance are recorded in the cave wall thermal constant database DB2 for each material constituting the cave wall.

  Further, as shown in FIG. 6, the pipe wall thermal constant database DB3 records the name, mass density, specific heat, and thermal resistance of each cooling water pipe for each material.

  In addition, as shown in FIG. 7, the fluid thermal constant database DB4 records mass density, specific heat, and thermal resistance for the air in the cave in the stationary state and the flowing state, and the cooling water in the stationary state and the flowing state, respectively. It is like that.

  Further, as shown in FIG. 8, the pipe dimension database DB5 records the name, outer diameter, and inner diameter of each specification of the cooling water pipe.

  In addition, as shown in FIG. 9, the cable database DB6 has outer diameters for conductors, insulators, and sheaths constituting the line, using the applied voltage, type name, number of cores, and size as indexes for each type of power cable. Each value of RAC, which is a specific heat capacity, a specific heat resistance, and an electric resistance value per unit length, is recorded.

  Further, as shown in FIG. 10, the air property database DB7 records the kinematic viscosity coefficient, the Prandtl number, and the thermal conductivity of the air in the sinus for each temperature.

  Further, as shown in FIG. 11, the evaporation efficiency database DB8 records in advance the evaporation efficiency of water on the air contact surface in the cave, and the floor surface 2a of the cave 1 shown in FIG. 2b, the evaporation efficiency of water on the surface of the cable 4 and the surface of the cooling water pipe 6, that is, the ease of evaporation of water from the surface when these surfaces are wetted by water or the like is 1 or less. The value is recorded.

  Further, as shown in FIG. 12, the cooling water physical property database DB9 records the kinematic viscosity coefficient, the Prandtl number, and the thermal conductivity of the cooling water for each temperature in the same manner as the air physical property database DB7.

  In addition, as shown in FIG. 13 for an example in Tokyo, the ground temperature database DB10 records the ground temperature at each depth from the ground surface by month for each region.

  In addition, as shown in FIG. 14 for example in Tokyo, the weather database DB11 records each hour of temperature, humidity, wind speed, solar radiation, and precipitation data for each area for each region. It has become.

  In addition, as shown in FIG. 15, the load factor database DB12 has each of the energization amounts at the time when the amount of electricity used is highest during the year, for example, at 14:00 on weekdays in August indicated by diagonal lines. Data indicating the ratio of the energization amount at the time point is recorded by month, day of the week (Sunday, weekday, Saturday), and time.

  The program recorded in the recording device 14 activates the central processing unit 11, and based on the data recorded in each of the databases DB1 to DB12 and the data set by the input device 12, etc., the finite element method Is used to calculate the temperature of each part of the analysis area crossing the cable in the specified analysis area at a specified time interval, and simulate the changes in the cable conductor temperature reflecting the effects of air movement in the cave and cooling water. To work.

  Here, the theoretical background or calculation method of the simulation using the finite element method by this program will be described.

This simulation has the following formula:
(1) K _X (∂ ² θ / ∂ _x ² ) + K _Y (∂ ² θ / ∂ _y ² ) + Q = ρC (∂θ / ∂t)
The equation (1) is applied to a predetermined analysis plane that traverses a cave, a cable, a cooling water pipe, and the like including the ground surface.

This equation (1) indicates the temperature θ at each time t at the point (x, y) on the analysis surface, K _X and _KY are the heat conduction coefficients in the X and Y directions, Q is the unit time, The amount of heat derived from the heat balance per unit volume, internal heat generation, etc., ρ is the mass density, and C is the specific heat.

When this equation (1) is expressed in a matrix for application of the finite element method,
(2) [C] · [dθ / dt] + [K] · [θ] = [Q]
It becomes.

Here, [C] is a heat capacity matrix, [K] is a heat conduction matrix, [θ] is a nodal temperature vector, [dθ / dt] is a time differential vector of the nodal temperature, and [Q] is a thermal load vector. As shown in FIG. 16, when a predetermined analysis region on the analysis surface is divided into a number of elements by finite element division and each node is given a number 1... I ... j, a term C _ij constituting the heat capacity matrix [C]. Indicates the product of the specific heat capacity (specific heat × mass density: J / cm ³ · K) between the nodes i and j and the volume related value therebetween, and the term K _ij constituting the heat conduction matrix [K] is , The product of the thermal conductivity (J / sec · cm ³ · ° K) between the nodes i and j and the value related to the volume between them.

In addition, in the case of a node on the ground surface, the term Q _j constituting the thermal load vector [Q] is a unit area on the ground surface of a region represented by the node, and a heat balance amount Q1 (J / Sec · cm ² ) and the area of the region (length on the analysis surface × unit length).

Among nodes in the ground, for nodes having internal heat generation such as the inside of a conductor, the unit volume in the region represented by the node, the internal heat generation amount Q2 (J / sec · cm ³ ) per unit time, and the region Is the product of the volume (area on the analysis surface × unit length).

Further, the term Q _j of the thermal load vector [Q] for the air contact surface in the cave, that is, the inner wall of the cave, the nodes disposed in the cave and the surface of the surface of the cooling water pipe, is represented by the node. A product of a unit area in a region, a heat transfer amount Q3 (J / sec · cm ² ) per unit time and an area of the region (length on the analysis surface × unit length) is shown.

In addition, a single node representing the air in the cave is set for the inside of the cave, and for that node, the unit area of the air in the cave, the heat transfer amount per unit time Q4 (J / sec · cm ^2). ) And the space cross-sectional area in the cave (the value obtained by subtracting the cross-sectional area of the cable and the cooling water pipe from the cave cross-sectional area).

Similarly, terms Q _j of the thermal load vector [Q] for the nodes of the cooling water pipe inner surface unit area in the region of the conduit interior surface, represented by the nodes, the amount of heat transfer per unit time Q5 (J / sec * The product of cm ² ) and the area of the region (length on the analysis surface × unit length).

In addition, a single node representative of the cooling water is set for the inside of the cooling water pipe, and the unit area of the cooling water and the heat transfer amount Q6 (J / sec · cm ² ) per unit time are set for the node. And the product of the cross-sectional area of the cooling water pipe.

If a node is not on the ground surface, in a cave or in a cooling water pipe, and is not accompanied by internal heat generation, that is, in the soil layer, inside a cave wall or inside a pipe wall, etc. The value of the term Q _j of the thermal load vector [Q] is zero.

In actual simulation, based on the above equation (2), the temperature θ t ′ at time t ′ (= t + Δt) after Δt time from the temperature θ _{t at} each node at time t is _expressed by the following equation:
(3) (2 [C] _{t ″} / Δt + [K] _{t ″} ) [θ] _{t ′}
= (2 [C] _{t ″} / Δt− [K] _{t ″} ) [θ] _t +2 [Q] _{t ″}
Will be asked according to.

Here, the time t ″ is an intermediate time between the times t and t ′ (= (t + Δt) / 2), and the equation (2) is expressed for the time t ″.
[C] _{t ″} · [dθ / dt] _{t ″} + [K] _{t ″} · [θ] _{t ″} = [Q] _{t ″}
In addition,
[Θ] _{t ″} = ([θ] _t + [θ] _{t ′} ) / 2
[Dθ / dt] _{t ″} = ([θ] _{t ′} − [θ] _t ) / Δt
By substituting the relationship, the above equation (3) is obtained.

In this equation (3), each term (C _ij ) _{t ″} constituting the heat capacity matrix [C] _{t ″} is given a value obtained from the specific heat capacity of each element to form the heat conduction matrix [K] _{t ″} . A value obtained from the thermal conductivity of each element is given to each term (K _ij ) _{t ″,} and a node on the ground surface is used as the value of each term (Q _j ) _{t ″} of the thermal load vector [Q] _{t ″.} In addition, the nodes on the air contact surface in the tunnel and the inner surface of the cooling pipe are given the amount of heat transfer at the time t ″ in the region represented by the node, and the nodes representing the air in the tunnel and the cooling water have the air and The amount of heat transfer at the time t ″ of the cooling water is given, and furthermore, the amount of heat generation at the time t ″ in the region represented by the node is given to the node inside the conductor.

In addition, each term of the initial value [θ] _{t = 0} of the node temperature vector is given a ground temperature or a separately set temperature according to the position (depth) or time of each node, and the time of the node temperature Initial value [dθ / dt] of differential vector 0 is given to each term of _{t = 0} , and the nodal temperature vector [θ] is sequentially calculated every predetermined analysis period Δt from that state for an arbitrary period.

  Note that the initial values for the nodes representing the air and cooling water in the cave are the depth and timing when the analysis is started from a steady state (a sufficient time has passed since the cave air and cooling water stopped). The temperature is set to the same temperature as the ground temperature, but when the analysis is started from the state where the air in the cave or the cooling water is flowing or when it starts to flow, the separately set temperature is set as the initial value.

When the specific heat capacity and thermal conductivity of each element are given as a function of time, the value of each term of the heat capacity matrix [C] _{t ″} and the value of each term of the heat conduction matrix [K] _{t ″} are Although it is obtained from a function relating to each time, a constant value is used when there is no time dependency.

On the other hand, in this simulation, as the value of the term Q _j of the thermal load vector [Q], as described above, the heat balance amount Q1 (J / sec · cm ² ) at the node on the ground surface, the internal at the node inside the conductor calorific value ^{Q2 (J / sec · cm 3} ), the amount of heat transfer Q3 at the node of the sinus tract air contact surface ^{(J / sec · cm 2)} , the amount of heat transfer nodes representative of sinus tract air Q4 (J / sec · cm ² ), The heat transfer amount Q5 (J / sec · cm ² ) at the node on the inner surface of the cooling water pipe, and the heat transfer amount Q6 (J / sec · cm ² ) at the node representative of the cooling water are used. In this embodiment, Each is calculated as follows.

First, the heat balance Q1 at the nodes on the ground surface is the solar radiation absorption amount Q11, the radiation amount from the atmosphere to the ground surface Q12, the radiation amount from the ground surface to the atmosphere Q13, and the heat transfer amount from the ground surface to the atmosphere. As Q14,
(4) Q1 = Q11 + Q12-Q13-Q14
It is defined as

Here, the amount of solar radiation absorption Q11 is
(4-1) Q11 = Solar radiation amount x (1-ground surface reflectance)
It is.

Further, the radiation amount Q12 from the atmosphere to the ground surface and the radiation amount Q13 from the ground surface to the atmosphere are σ is a Stefan-Boltzmann constant (J / sec · cm ² · ° K ⁴ ), and ε12 and ε13 are radiation rates (none). Dimension), T is the atmospheric temperature (° K), T ′ is the temperature of the element on the ground surface (° K),
(4-2) Q12 = σ × ε12 × T ⁴
(4-3) Q13 = σ × ε13 × T ′ ⁴
Indicated by

In that case, for the radiation amount Q12 from the atmosphere to the ground surface, the emissivity ε12 is, for example,
No precipitation, less than 50% humidity, ε12 = 0.650
No precipitation and humidity 50% or more, ε12 = 0.850
With precipitation, ε12 = 0.925
And set.

In addition, for the radiation amount Q13 from the ground surface to the atmosphere, the radiation rate ε13 is independent of precipitation and humidity,
ε13 = 0.965
And set.

Furthermore, the amount of heat transfer Q14 from the ground surface to the atmosphere is
(4-4) Q14 = ρ ′ × Cp × D × (T′−T)
It is defined as

Here, ρ ′ is the density of the atmosphere (g / cm ³ ), Cp is the constant pressure specific heat of the atmosphere (J / g · ° K), T ′ is the temperature of the element on the ground surface (° K), and T is the atmospheric temperature. (° K), D is an external diffusion coefficient (cm / sec), and U is a wind speed, for example,
D = 0.0027 + 0.031 × U
And set.

  The amount of solar radiation, atmospheric temperature (air temperature), humidity, precipitation, and wind speed are read from the weather database DB11. In addition, a constant value incorporated in the program is used for the ground surface reflectance. For example, the reflectance is recorded for each soil constituting the ground surface in the soil thermal constant database DB1, and is read and used. Also good.

Then, the value obtained by multiplying the total heat balance amount Q1 on the ground surface represented by the equation (4) by the area of the ground surface region represented by each node is the thermal load vector [Q in the above equation (3). ] Is used as the value of the term Q _j for the node located on the ground surface.

Also, the internal heating value Q2 at the node in the cable conductor is
(5) Q2 = (I ² × RAC) / S
It is defined as

The numerator of the formula (5) indicates the amount of heat generation (J / sec · cm) per unit length of the cable, and the value obtained by dividing this by the cross-sectional area S of the conductor portion of the cable is the unit time and unit of the conductor. The internal heating value per volume is Q2 (J / sec · cm ³ ).

  Here, I is the current (A), RAC is the resistance per unit length of the cable (Ω / cm), and RAC is obtained from the cable database DB6. Further, the current I is obtained by using a separately set value and the load factor database DB12.

Then, the value obtained by multiplying the internal heating value Q2 represented by the equation (5) by the volume of the region represented by the node in the cable conductor is the value of the thermal load vector [Q] in the above equation (3) for the node. used as the value of the term _{Q j.}

Next, regarding the heat transfer amount Q3 (J / sec · cm ² ) of the node on the air contact surface on the inner surface of the cave, the heat transfer amount Q31 between the node and the air in the cave and the region wetted by water Assuming that there is a heat transfer quantity as latent heat of vaporization due to water evaporation, i.e., a latent heat transport quantity Q32 for the nodes in the inside,
(6) Q3 = Q31 + Q32
It is defined as

Among these, as shown in FIG. 17, the heat transfer amount Q31 between the air in the cave is Ta, the temperature of the node representing the cave air, and the temperature of each node on the air contact surface is Tfx (x = 1, 2). , ..., i, j, ..., n, ...)
(6-1) Q31 = α × (Tfx−Ta)
Indicated by

Here, α in the above formula (6-1) is an average heat transfer coefficient (J / sec · cm ² · ° C.) between each node of the air contact surface in the sinus and the air, and heat of the air The conductivity is λ (J / sec · ° C.), the circular diameter (equivalent diameter) of the area equal to the space cross-sectional area in the sinus is d (cm), and the Nusselt number is Nu.
α = Nu × λ / d
Indicated by

The Nusselt number Nu is obtained from the equation of steady-state laminar heat transfer when the flow velocity is 0, that is, when the air in the cave is stopped, and when the air is moving, the turbulent heat in the tube It is obtained from the equation of transmission and is as follows.
When the air is stopped: Nu = 3.66
At the time of air movement: Nu = 0.023 × Re ^0.8 · × Pr ^0.4
Re is the Reynolds number and Pr is the Prandtl number. Re = ω × d / ν, Pr = ν / a
Ω is the average velocity of the air in the cave (cm / sec), ν is the kinematic viscosity coefficient (cm ² / sec) of the air in the cave, and a is the temperature conductivity (cm ² / sec) of the air in the cave .

Therefore, the average heat transfer coefficient α is as follows.
When the air is stopped: α = 3.66 × λ / d
At the time of air movement: α = 0.023 × (ω × d / ν) ^0.8 × (ν / a) ^0.4 × λ / d
Note that the thermal conductivity λ, kinematic viscosity coefficient ν, and Prandtl number Pr (= ν / a) of the air in the above equation are read from the air property database DB7 according to the temperature of the air in the sinus at that time. The equivalent diameter d of the cave is calculated from preset data, and the average flow velocity ω of the cave air is obtained from the wind speed set at the time of analysis.

FIG. 19 shows Reynolds numbers Re _min (air temperature: 60 ° C.) and Re _max (air temperature: −20 ° C.) when air is flowed into the cave under various conditions. Over 2320, the flow under these conditions is shown to be turbulent.

On the other hand, among the heat transfer amounts Q3 of the nodes on the air contact surface on the inner surface of the cave, the latent heat transport amount Q32 due to water evaporation at the nodes wet with water is defined by the following latent heat transport equation.
(6-2) Q32 = ι × β × ρ × Uch
× {q _SAT −q + (dq _SAT / dTa) × (Tfx−Ta)}

  Here, Ta is the temperature of a node representative of the air in the cave, and Tfx is the temperature of each node on the air contact surface (see FIG. 17).

Ι is the latent heat of vaporization (J / g)
ι = 2.5 × 10 ³ −2.4 × Ta
Indicated by

  Β is the evaporation efficiency of water on the air contact surface, and is registered in advance in the evaporation efficiency database DB8 with a value of 1 or less for each type of air contact surface in the sinus.

Ρ is the density of air in the cave (g / cm ³ ), p is the pressure in the cave (hPa), e is the water vapor pressure in the cave (hPa),
ρ = 1.29 × 10 ⁻³ × {273 / (273 + Ta)} × p
× (1-0.38 × e / p) / 1013
Indicated by Note that the water vapor pressure e in the cave is defined as e _SAT is saturated water vapor pressure, and Ha is cave air humidity.
e = Ha × e _SAT
And the saturated water vapor pressure e _SAT is
e _SAT = 6.11 × 10 ^{7.5 Ta / (237 + Ta)}
Indicated by

  Moreover, Uch of Formula (6-2) is a heat exchange rate (cm / sec), and a value at the time of light wind stabilization, 0.01, is adopted as a fixed value.

Further, q represents the specific humidity of the air in the cave, q _SAT represents the saturation specific humidity, and is expressed as follows using the cave pressure p, the water vapor pressure e, and the saturated water vapor pressure e _SAT , respectively.
q = (0.62 × e / p) / (1−0.38 × e / p)
q _SAT = (0.62 × e _SAT /p)/(1−0.38×e _SAT / p)

Further, the rate of change of saturation specific humidity with respect to temperature, dq _SAT / dTa,
dq _SAT / dTa = (de _SAT /dTa)×0.62×p
/(P−0.38×e _SAT ) ²
Indicated by
de _SAT /dTa=6.11×(2500−2.4Ta)
× 10 ^{7.5 Ta / (237 + Ta) / {0.46} × ^{(273 + Ta} ) ² }
It is.

  Note that since the air density ρ and the specific humidity q in the equation (6-2) are functions of the water vapor pressure e, and the water vapor pressure e is a function of the humidity Ha, the latent heat transport amount Q32 is a function of the humidity Ha. However, when the latent heat transport occurs, the humidity Ha changes due to the evaporation of the water that caused it. Therefore, after calculating the latent heat transport amount Q32, the humidity Ha of the sinusoidal air is updated in preparation for the calculation in the next analysis cycle.

That is, first, the absolute humidity (amount of water vapor per unit volume) h ₀ (g / cm ³ ) in the current air is obtained from the current temperature Ta ₀ and the humidity Ha _0, and then this absolute humidity From h ₀ , the mass m of water vapor in the air passing through the cross section of the sinus during the passage of time Δt, which is one analysis period, from the present time using the cross sectional area A of the sinus and the average flow velocity ω of the air in the sinus _{Find 0} (g).
h ₀ = 0.217 × 10 ⁻³ × e _SAT × Ha ₀ / (273 + Ta ₀ )
m ₀ = h ₀ × ω × A × Δt

Further, the sum Σ () of all the nodes of the product of the latent heat transport amount Q32 at the nodes wetted by water among the nodes on the air contact surface in the sinus and the length l on the analysis surface represented by the nodes. Q32 × l) (J / sec · cm) and latent heat of vaporization ι (J / g) are used to determine the mass Δm (g) of water evaporated from the entire air contact surface during the time Δt.
Δm = Σ (Q32 × l) × Δt × Δt × ω / ι

Since the mass Δm of the evaporated water is equal to the mass of the water vapor increased in the air passing through the sinus section during the time Δt, the mass Δm of the evaporated water is set to the mass m ₀ of the water vapor in the air. By adding, the mass m (= m ₀ + Δm) of water vapor in the air in the next analysis cycle is obtained, and the absolute humidity h in the next analysis cycle is calculated by calculating the above formula based on the mass m of the water vapor. Is obtained by the following equation. That is,
h = m / (ω × A × Δt)
It is. Then, by using this absolute humidity h and the newly calculated air temperature Ta in the tunnel, the humidity Ha in the next analysis cycle is obtained by the following equation.
(7) Ha = h × (273 + Ta) / (0.217 × 10 ⁻³ × e _SAT )

  In this way, the humidity Ha is updated every analysis cycle, and the latent heat transport amount Q32 at each node is calculated using the updated humidity Ha. According to the above calculation, when the humidity Ha becomes 100% or more, the evaporation of water is not performed, so the latent heat transport amount Q32 is set to zero.

Next, a description will be given of the heat transfer amount Q4 (J / sec · cm ² ) of the node representative of the air in the cave. This heat transfer amount Q4 is paired with the heat transfer amount Q31 between the air contact surface and the air in the cave. As shown in FIG. 17, the temperature of a node representing the air in the cave is Ta, and the temperature of each node on the air contact surface is Tfx (x = 1, 2,..., I, j,..., N, …)
(8) Q4 = α × Σ (Tfx−Ta)
Indicated by Here, α is an average heat transfer coefficient (J / sec · cm ² · ° C.) between each node of the air contact surface in the sinus and the air used in the equation (6-1).

Furthermore, the amount of heat transfer Q5 (J / sec · cm ² ) at the node on the inner surface of the cooling water pipe and the amount of heat transfer Q6 (J / sec · cm ² ) at the node representing the cooling water are determined between the inner surface of the pipe and the cooling water. It is calculated as follows due to heat transfer between them.

That is, as shown in FIG. 18, if the temperature of a node representing the cooling water is Tw and the temperature of each node on the inner surface of the pipe is Tgx (x = 1, 2,...), Each node on the inner surface of the cooling water pipe The amount of heat transfer Q5 is
(9) Q5 = α ′ × (Tgx−Tw)
The amount of heat transfer Q6 at the node representative of cooling water is
(10) Q6 = α ′ × Σ (Tgx−Tw)
Indicated by

Α ′ in the equations (9) and (10) is an average heat transfer coefficient between each node on the inner surface of the cooling water pipe and the cooling water, and the air in the equations (6-1) and (8). Like the average heat transfer coefficient α for
α ′ = Nu ′ × λ ′ / d ′ (J / sec · cm ² · ° C.)
Indicated by

Here, λ ′ is the thermal conductivity (J / sec · ° C.) of the cooling water, d ′ is the pipe inner diameter (cm) of the cooling water pipe, Nu ′ is the Nusselt number, and as described above, Re ′ is the Reynolds number, Pr ′ is the Prandtl number,
When cooling water is stopped: Nu ′ = 3.66
During movement of cooling water: Nu ′ = 0.023 × Re ′ ^0.8 × Pr ′ ^0.4
Indicated by Re ′ = ω ′ × d ′ / ν ′, Pr ′ = ν ′ / a ′, ω ′ is the average flow velocity (cm / sec) of cooling water, and ν ′ is the kinematic viscosity coefficient (cm ² / sec), a ′ is the temperature conductivity (cm ² / sec) of the cooling water.

Therefore, the average heat transfer coefficient α ′ is
When cooling water is stopped: α ′ = 3.66 × λ ′ / d ′
During movement of cooling water: α ′ = 0.023 × (ω ′ × d ′ / ν ′) ^0.8
× (ν ′ / a ′) ^0.4 × λ ′ / d ′
It becomes.

  Then, the thermal conductivity λ ′, the kinematic viscosity coefficient ν ′, and the Prandtl number Pr ′ (= ν ′ / a ′) of the cooling water are read from the cooling water physical property database DB 9 according to the temperature of the cooling water at that time. . Further, the inner diameter d ′ of the cooling water pipe is read from the pipe dimension database DB5. Further, the average flow velocity ω ′ is obtained from the inner diameter d ′ and a preset flow rate of the cooling water.

FIG. 20 shows a Reynolds number Re ′ _min (cooling water temperature) when cooling water at 0 ° C. and 60 ° C. is flowed through different diameters of pipes generally used as cooling water pipes. : 0 ° C.) and Re ′ _max (cooling water temperature: 60 ° C.), both of the cooling waters exceed the turbulent flow conditions of 2320, indicating that the flow under these conditions is turbulent. ing.

As described above, the heat transfer amounts Q3 to Q6 of each node for the air contact surface in the cave, the air in the cave, the cooling water pipe inner surface and the cooling water are obtained, and the values of the regions represented by the corresponding nodes are obtained. The value multiplied by the area is used as the value of the term Q _j for each node of the thermal load vector [Q] in the above-described equation (3).

  By the way, the temperature Ta and the humidity Ha of the sinusoidal air in the equations (6-1), (6-2), and (8) and the temperature Tw of the cooling water in the equations (9) and (10) Since it is influenced by the flow of air and cooling water, the details will be described later, but the transfer process from the upstream side to the downstream side is performed by the following method.

  In other words, a plurality of analysis surfaces crossing the cave where the cable and the cooling water pipeline are arranged are set, and initial values of air, humidity, and cooling water in the cave are given, and then each analysis surface has a cable. In consideration of the internal heat generation amount, the latent heat effect due to the evaporation of water, the cooling effect due to the cooling water, etc., based on the above equations (3), (4) to (6), (8) to (10) Finite element analysis is performed at time intervals (analysis period). Thereby, the temperature Ta of the node representing the air in the sinus and the temperature Tw of the node representing the cooling water are calculated simultaneously with all the other nodes in the analysis region of each analysis surface. In parallel with this, the humidity Ha of the air in the cave is updated each time.

  In that case, air movement may occur in the sinus due to forced ventilation or naturally occurring, and the cooling water moves in the pipeline at the start of cooling. In each predetermined analysis cycle, a process of passing the air temperature Ta, the humidity Ha, and the cooling water temperature Tw from the upstream analysis surface to the downstream analysis surface is performed.

  That is, first, when at least one of the air in the cave or the cooling water is moving, the longest time during which the air or the cooling water moves between the adjacent analysis surfaces due to the movement, that is, all the adjacent analysis surfaces. The time during which the movement of the fluid ends (the longest movement time) is calculated in advance.

  Then, in the analysis cycle immediately after the longest movement time has elapsed from the time when the analysis starts or when the fluid starts to move, after calculating all the node temperatures individually for each analysis surface and updating the humidity, In the analysis plane, calculate how the temperature and humidity changed during the time required for the fluid to flow from the analysis plane to the downstream analysis plane, and calculate the calculated temperature on the downstream side. The air temperature Ta, the humidity Ha, and the cooling water temperature Tw on the analysis surface are used.

  In this way, the transfer process is performed at every analysis cycle in which the movement of the fluid between all the analysis surfaces is completed, and the air temperature Ta, the humidity Ha, and the cooling water temperature Tw obtained thereby are expressed by the above equation (6-1). ), (6-2), and (8) by using the air temperature and humidity, and the cooling water temperature in equations (9) and (10), and continuing the analysis, Simulate changes in cable conductor temperature on each analysis surface while reflecting cooling and latent heat effects.

  Specifically, the above simulation method includes a program read from the recording medium 20 and recorded in the recording device 14 of the computer 10 shown in FIG. 4, various databases DB1 to DB12, various data input by the input device 12, and the like. Based on the above, the central processing unit 11 will execute, and hereinafter, these simulation operations will be described using a flowchart and a screen displayed on the display device 15 of the computer 10.

  First, when the system is activated on the computer 10, a model selection screen W1 shown in FIG. 21 is displayed on the display device 15. On this screen W1, the model to be analyzed has a cooling water pipe. Choose between a cooling model or an uncooled model without a cooling water line. Here, it is assumed that the cooling model is selected.

  FIG. 22 shows a flowchart of the main routine of this cooling model. First, in step S1, on the screen W2 of FIG. 23 displayed on the display device 15 of the computer 10, the section where the cooling water pipeline is arranged is shown. A plurality of analysis planes that traverse the sinus including the pipes and cables are set. In this case, as shown in FIG. 24, the first analysis plane A (position 0.00) is set at a position closest to the cooling water supply device 7, and the subsequent analysis planes are sequentially directed toward the turning point of the pipe 6. Enter the B, C, and D positions.

  Next, in step S2 of the main routine, the configuration of the cave, such as the arrangement of the cave, cable, and cooling water conduit on the analysis surface and the configuration around the cave, is set. This setting is performed as follows on the screen W3 shown in FIG.

  First, after inputting the number of the analysis surface and the position (X coordinate) of the sinus reference point with respect to the reference point of the predetermined analysis area, the size, width, height, wall thickness, and depth (analysis) of the sinus Enter the Y-coordinate of the cave reference point with respect to the reference point of the area.

  Next, the width and depth of the soil layer around the cave are set. An area indicated by the width and depth is an analysis area on the analysis surface. In this case, if the automatic setting check box is turned on, the analysis area is automatically set based on the size of the sinus. In addition, when there are a plurality of soil layers having different types of soil, the layer thickness is input for each layer (the lowermost layer pressure is calculated from the upper layer pressure and depth).

  In addition, for cables, the number of lines, the configuration of the one-phase laying type and the triple-row type for each line, the arrangement position (X coordinate, Y coordinate) in the tunnel, and the cables that make up each line Set the voltage, type, size (number of cores), etc. At this time, various data about the cable set from the cable database DB6 is read and used for subsequent finite element division and analysis. In the case of the illustrated example, the number of lines is 3, and each line has a three-strip type (see FIG. 27).

  Further, for the cooling water pipeline, the number of pipelines, the center positions of the Go and Re tubes for each pipeline, and the pipeline type are set. And by setting a pipe line type, the outer diameter is read from the pipe line dimension database DB5 and displayed.

  As a result, the configuration related to the cave in the analysis region, the configuration related to the cave walls and soil layers around the cave, the configuration of cables and cooling water pipes in the cave, etc. are determined. A finite element analysis model obtained by dividing the region into finite elements is automatically created, and the entire model is displayed on the screen W4 as shown in FIG. 26, and its main part is displayed on the screen W5 as shown in FIG. Is done.

  Next, as step S3 of the main routine, soil data is set on the screen W6 of FIG. That is, in the input form displayed on the screen W6, the type of soil read from the soil heat constant database DB1 recorded in the recording device 14 is displayed in the pull-down menu, and from this time, the current form is displayed. Specify the type of soil to be analyzed. At this time, the thermal constants such as the density, specific heat, and thermal resistance of the soil read from the soil thermal constant database DB1 are displayed. When a plurality of soil layers are set, the type is designated for each soil layer.

  Next, as step S4 of the main routine, the sinus road data is set on the screen W7 shown in FIG. That is, in the input form displayed on the screen W7, the type of the sinus wall material read from the sinus wall thermal constant database DB2 recorded in the recording device 14 is displayed in the pull-down menu. Specify the type of cave wall material to be analyzed this time. At this time, the thermal constants such as the density, specific heat, and thermal resistance of the sinus wall material read from the sinus wall thermal constant database DB2 are displayed.

  Furthermore, as step S5 of the main routine, latent heat transport amount calculation data is set on a screen W8 shown in FIG. That is, the wetness ranges for the floor surface and the left and right wall surfaces are input in the input field on the left side of the input form displayed on the screen W8. In that case, if the entire surface of the floor is wet due to flooding or the like, the input of the range can be omitted by selecting the entire surface button. In addition, by selecting the cable and cooling water pipe displayed on the analysis model screen that are wet with water and clicking the Add button, the selected cable or pipe is displayed on the right side of the input form. The In this way, an area that is wet with water, in other words, an area that is subject to calculation of the latent heat transport amount is set.

  As described above, various settings are repeatedly performed on all of the plurality of analysis surfaces set on the screen W2. However, the position, shape, properties, and latent heat of the sinus and its internal and external components are also determined. When the calculation target area of the transportation amount is the same on all analysis planes, the configuration set for the first analysis plane A remains unchanged by turning on the check box of the same shape on all analysis planes on the screen W2. The analysis plane is always set.

  When the settings for all analysis surfaces are completed, next, Steps S6 to S7 are executed to set data relating to air movement such as ventilation conditions. This setting is performed on a screen W9 shown in FIG. 31. First, a time schedule for starting and stopping ventilation is set, and a wind speed at the time of starting ventilation and a wind speed at the time of stopping are set.

  In that case, when the ventilation is stopped, there are cases where the movement of air completely stops and there is a case where the flow naturally occurs, so input the wind speed (value greater than 0 or 0) in those cases, In the latter case, specify the direction of the flow at the start of ventilation or the opposite direction. In addition, as the physical property of the air in the cave, it is specified whether it is flowing air or static air, and further, the initial temperature of the cave air and humidity is set. The initial temperature and initial humidity can be automatically set with reference to the temperature and humidity recorded in the weather database DB11.

  Next, in step S8 of the main routine, data relating to the movement of the cooling water such as the cooling start condition is set. This setting is performed using the screen W10 shown in FIG. 32. When the activation condition of the cooling water is selected from the cable energization current or the cable surface temperature, and the energization current is selected, the current for starting the cooling is selected. When the surface temperature is selected by setting the value and the current value at which the cooling is finished, the temperature at which the cooling is started and the temperature at which the cooling is finished are set. In the case of the illustrated example, the cooling water is activated when the analysis value of the cable surface temperature exceeds 60 ° C., and the cooling water is stopped when the analysis value decreases to 40 ° C.

  Further, in this screen W10, in the cooling water supply temperature (inlet temperature) and the supply flow rate, and in the folded portion of the cooling water pipe line, whether the Go water and the Re pipe on the final analysis surface have the same cooling water temperature, Whether to consider the temperature transfer from the Go pipe to the Re pipe is selected. In the latter case, the arrival time of the cooling water from the Go pipe to the Re pipe on the final analysis surface is input.

  The above-described setting of the cooling water data is performed for each pipe line when there are a plurality of cooling water pipe lines, and is repeatedly executed until it is determined in step S9 that the setting for all the pipe lines is completed.

  Furthermore, in step S10 of the main routine, energization conditions are set using a screen W11 shown in FIG. On this screen W11, the previously set track configuration and the like are displayed for each track number. Therefore, for the line indicated by the number, the maximum current to be applied to the cable is input, and whether or not to consider the yearly increase is selected, and when it is considered, the yearly increase rate is input.

  In step S11 of the main routine, data relating to the analysis target area, that is, the area where the simulation is performed and the month when the simulation is started are set using the screen W12 shown in FIG. This gives the ground temperature of the simulation start month at each depth of the area read from the ground temperature database DB10 as an initial value of the nodal temperature of the entire area of the analysis surface at the start of the simulation, and each of the analysis target areas from the weather database DB11. This is because the weather data of the date and time is read out. In addition, when air or cooling water starts moving at the same time as the start of analysis, or when analysis is started with these moving, the initial temperature and initial temperature of the air in the cave preset in the screens W9 and W10 Humidity and cooling water inlet temperature are used.

  Further, in step S12 of the main routine, analysis conditions are set using a screen W13 shown in FIG. That is, on this screen W13, it is selected whether this simulation is an initial analysis or a restart analysis. In the case of the initial analysis, the analysis period is specified by the number of years or the end time, and the unit of the analysis time interval, for example, time, day, month, and the output time interval of the calculation result are set. Here, as the output time interval, either the output of the results of all analyzes performed for each analysis cycle or the time interval separately input is selected.

  On the other hand, in the case of restart analysis, in addition to the above settings, the name of the project to be the basis of restart and the location in the recording device 14 storing the project are designated, and the continuation of the base project is continued. Or restart from an intermediate time point, and in the latter case, the restart point is designated. When the energization condition is changed at the time of restart, the energization current condition is reset on the screen W11 shown in FIG.

  As described above, when all the settings are completed and the set data is recorded in the memory 17, the central processing unit 11 of the computer 10 performs an analysis for the cable conductor temperature simulation as step S13 of the main routine. Execute the process.

  This processing operation will be described in detail again. When the analysis processing is completed, the result is displayed on the display device 15 of the computer 10 as shown in FIGS. Printed out by the printer 16 on demand.

  A screen W14 in FIG. 36 shows the temperature distribution around the sinus when a predetermined time has elapsed after the start of analysis, and a screen W15 in FIG. 37 shows changes in the air temperature in the sinus from the start of analysis. Yes, a screen W16 in FIG. 38 shows a change in the cable conductor temperature from the start of the analysis.

  Next, the operation at the time of execution of the analysis process in step S13 of the main routine will be described with reference to the flowchart of the subroutine in FIG. Here, in the following description, as shown in FIG. 24, it is assumed that four analysis planes A to D crossing the cave and the cooling water pipe are set.

  First, in step S21 of the flowchart of FIG. 39, necessary data is read from the database recorded in the memory 17 or recorded in the recording device 14, and then in step S22, an initial value 0 is set at the current time t.

  In step S23, it is determined whether or not the ventilation activation condition is established based on the ventilation schedule or the like set in advance on the screen W9 of FIG. 31, and if it is established, the wind speed is set in step S24. The ventilation start value set on the screen W9 is set, and if the ventilation condition is not satisfied, the ventilation stop value similarly set on the screen W9 is set in step S25.

  Next, in step S26, an analysis time interval (analysis period) Δt is set, and the longest longest movement time Δtc among the air movement time and the cooling water movement time between the respective analysis surfaces is calculated. These travel times are calculated based on the distance between the analysis surfaces and the wind speed for air, and based on the distance between the analysis surfaces, the flow rate, and the pipe cross-sectional area for cooling water. .

  In step S27, it is determined whether at least one of the condition that the wind speed of the air in the cave is currently greater than 0 and the cooling start condition set in advance on the screen W10 in FIG. 32 are satisfied. When at least one of the conditions is satisfied, that is, when at least one of the air in the cave or the cooling water is moving, in step S28, the current time t is set to the air temperature, humidity, and cooling water temperature. Is set to the delivery reference time tp. Note that when neither of the above conditions is satisfied, that is, when neither the air in the cave nor the cooling water is moving, the reference time tp is not set.

  In step S29, the analysis period Δt is added to the time t, and this is set as the current time t. In step S30, it is determined whether or not the time t has reached the analysis end time.

  Since the analysis end time has not been reached at the beginning of the analysis, step S31 is executed next, and the temperature of all nodes in the analysis region at the previous analysis time (t-Δt) for each analysis surface. From these, the temperatures of all nodes at the current time t are calculated.

  At the time of the first analysis, the temperature of each node is the initial value, that is, for the node representing the air in the cave, the ground temperature corresponding to the initial temperature set on the screen W9 in FIG. 31 or the depth read from the ground temperature database DB9, etc. For the nodes that represent the cooling water, the ground temperature according to the inlet temperature or depth set on the screen W10 in FIG. 32, and the other temperatures for the other nodes are the ground temperature according to the depth, etc. The temperature will be determined.

The calculation for obtaining the temperature of each node at the current time t after the time Δt from the value (or the initial value) at the time of the previous analysis is based on the theory shown in the above formulas (1) to (3), and the flowchart in FIG. This is done as follows by the subroutine shown. Here, it is assumed that the heat capacity matrix [C] and the heat conduction matrix [K] have no dependency on the temperature and time, and each of the terms [C] _{t ″} and [K] _{t ″} in the above equation (3). The values are all fixed values calculated based on values read from the databases DB1, DB2, DB3, DB4, DB6, etc. in which the thermal constants of each element are recorded.

  First, in step S41, for the analysis model divided into finite elements as shown in the screen W4 and the screen W5, using the data such as the mass density, specific heat, and specific heat capacity of each element read from the database, The heat capacity between the nodes is calculated, and a heat capacity matrix [C] is created with this value as the value of each term.

  In step S42, the thermal conductivity between the nodes over the entire analysis region is calculated using the data of the thermal resistance, intrinsic thermal resistance, etc. of each element read out from the database, and this is used as the value of each term. A heat conduction matrix [K] is created.

  Further, in step S43, a thermal load vector [Q] at the time is generated. Among the terms of this thermal load vector [Q], the value Q1 of the term corresponding to the node located on the ground surface is obtained by using the above-described formula ( It is set based on the heat balance amount calculated according to 4).

  Further, the value Q2 of the term corresponding to the node located in the conductor of the cable is the load factor of the corresponding month, day of the week, and time read from the load factor database DB12 and the maximum current and year set on the screen W11 of FIG. The energization current at that time is calculated from the rate of increase, and the internal heating value calculated from the above equation (5) using the current value and the RAC data of the corresponding cable read from the cable database DB6. Is set based on

  Further, the node Q3 of the air contact surface in the sinus is based on the temperature Tfx (x = 1, 2,...), The air temperature Ta, and the humidity Ha of each node of the air contact surface obtained by the previous analysis. The value Q4 of the term obtained from (6), (6-1), (6-2) and corresponding to the node representative of the air in the sinus is similarly calculated for each of the air contact surfaces obtained by the previous analysis. Based on the contact temperature Tfx (x = 1, 2,...) And the air temperature Ta, the temperature is obtained from the above-described equation (8).

  Similarly, the values Q5 and Q6 of the terms corresponding to the nodes on the inner surface of the cooling water pipe and the nodes representing the cooling water are the temperatures Tgx (x = 1, 2,...) Of the respective contacts on the inner face of the pipe obtained by the previous analysis. ) And the cooling water temperature Tw, respectively, are obtained from the aforementioned equations (9) and (10).

In addition, the value of the term corresponding to the node located in the soil layer other than the above is set to 0, whereby the value of each term of the thermal load vector [Q] at the present time is set. Then, this thermal load vector [Q] is set to [Q] _{t ″} in the above equation (3).

When the heat capacity matrix [C], the heat conduction matrix [K], and the thermal load vector [Q] _{t ″} are created as described above, next, in step S44, these are substituted into the above equation (3). and, from the node temperature vector [theta] _t at time t, determining the nodal temperature vector after time Δt [θ] _{t '.}

That is, as the value of each term of the temperature vector [θ] _t , first, the initial value [θ] _{t = 0} of the nodal temperature vector is created by substituting the initial value described above, and the value of each term is If the initial value [dθ / dt] _{t = 0} of the time differential vector of 0 node temperature is created, the temperature vector [θ] _{t ′} after Δt time is sequentially calculated based on these values. In step S45, the temperature vector [θ] _{t ′} after Δt time is recorded in the memory 17.

  In this way, in step S31 of the flowchart of FIG. 39, for each of the analysis surfaces A to D, the temperatures of all the nodes in the analysis region are calculated for each period Δt.

  Further, after this calculation, in step S32, the humidity of the air in the cave is updated for each analysis surface using the above-described equation (7).

  Next, in step S33 of the flowchart of FIG. 39, it is determined whether or not at least one of the condition that the wind speed of the air in the cave is currently greater than 0 and a preset cooling start condition are satisfied. And when both conditions are not satisfied, that is, when neither the air in the cave and the cooling water are moving, the processing for transferring the temperature and humidity of the cave air and the temperature of the cooling water in step S34 and after are executed. Without repeating, steps S23 to S32 are repeatedly executed. That is, in this case, for each analysis surface A to D, the temperatures of all nodes are sequentially obtained independently from the other analysis surfaces for each analysis period Δt.

  On the other hand, when at least one of the two conditions is satisfied, that is, when at least one of the air in the cave or the cooling water is moving, in step S34, the current time t is the temperature delivery reference set in step S28. It is determined whether or not the time (tp + Δtc) obtained by adding the longest movement time Δtc of air or cooling water calculated in step S26 to the time Tp is exceeded. If this time (tp + Δtc) has not been exceeded, steps S29 to S32 are repeatedly executed, and at each analysis period Δt, in step S31, each analysis surface A to D is calculated from each node temperature at the time of the previous analysis. The node temperatures at time t are sequentially calculated, and the humidity of the air in the sinus is updated in step S32.

  Then, when the current time t exceeds the time (tp + Δtc) obtained by adding the longest movement time Δtc to the temperature delivery reference time tp, the air temperature and humidity, and the cooling water temperature are delivered.

  That is, first, it is determined whether or not the air in the cave is moving in step S35. If it is moving, the temperature and humidity of the air in the cave are transferred from the upstream analysis surface to the downstream analysis surface in step S36. In step S37, it is determined whether or not the cooling water is moving. If the cooling water is moving, the cooling water temperature is transferred from the upstream analysis surface to the downstream analysis surface in step S38. Execute the process. In this case, since at least one of the condition that the wind speed of the air in the cave is greater than 0 or the cooling start condition is satisfied, the air temperature and humidity transfer process in step S36 or the cooling water temperature transfer in step S38. At least one of the processes is always executed.

  As described above, when both the air in the cave and the cooling water are not moving, the node temperature is calculated independently for each analysis surface every analysis cycle, and at least the air in the cave or the cooling water is at least When one side is moving, from the temperature / humidity delivery reference time tp set at the start of the movement until the longest movement time Δtc between the respective analysis surfaces of the moving fluid elapses, Similarly, an independent nodal temperature is calculated for each analysis surface for each analysis period. Then, in the analysis cycle immediately after the longest movement time Δtc has elapsed from the reference time tp, that is, in the cycle in which the moving fluid has finished moving between all analysis surfaces, the temperature and humidity of the air or the temperature of the cooling water At least one delivery process is performed.

  Then, when at least one of the air in the cave or the cooling water is still moving at the time when this delivery process is completed, after setting the time t to the delivery reference time tp again in step S28, The transfer process is performed again in the analysis cycle immediately after the longest movement time Δtc of air or cooling water has elapsed from the time. In this way, while at least one of the air and the cooling water is moving, the air temperature and humidity and the cooling water temperature are transferred at every cycle in which the movement is completed between all analysis surfaces. Thereby, the node temperature in each analysis plane is calculated reflecting the movement of temperature and humidity accompanying the movement of air and cooling water.

  Next, specific operations of the air temperature / humidity transfer process in step S36 and the coolant temperature transfer process in step S38 will be described.

First, when describing the transfer process of the sinus tract air temperature, this process is performed by a subroutine shown flowchart in FIG. 41, the initial firstly, in step S51, a preset air temperature T _A in the analysis plane A on the most upstream Set to temperature Tao.

Then, based on the initial temperature Tao, but will calculate the air temperature T _B of the analysis plane B at step S52, the calculation is performed according to a subroutine shown by a flow chart of FIG. 42, first, at step S61, The movement time Δtab of the air between the analysis planes A and B is read from the memory 17 that has been calculated and recorded in advance. Next, in step S62, the movement time Δtab is larger than iΔt and smaller than (i + 1) Δt. I is obtained. That is, counting from the time tp when the air starts to move, the number of times and how many times the analysis period elapses are obtained.

In step S63, the air temperature T _Ai calculated from the reference time tp and calculated in the i-th analysis cycle (time: tp + iΔt) and the (i + 1) -th analysis cycle (time: tp + (i + 1) Δt) are calculated. the air temperature T _{a (i + 1)} read from the memory 17, as shown in FIG. 43, these values by linear interpolation, to calculate the air temperature at time (tp + Δtab) analysis plane a, analyzed surface this Let B be the air temperature TB of the temperature transfer process execution cycle in _B.

Now, for example, when the air movement time Δtab between the analysis surfaces A and B is displayed as a multiple of the analysis period Δt, if it is 0.6Δt, i = 0, and in this case, the air movement start time (reference time) and the air temperature T _a that is calculated by the analysis plane a with a period of tp, the next cycle (time: a tp + Delta] t) air temperature T _a that is computed by, 6: the values interpolated by interior division ratio of 4 the air temperature T _B of the analysis plane B in this transfer cycle.

Similarly, in step S53 of the flowchart of FIG. 41, the air movement time Δtbc between the analysis surfaces B and C is shown in FIG. 42, and in step S54, the air movement time Δtcd between the analysis surfaces C and D is shown in FIG. Processing similar to the processing shown in the flowchart is performed, and the air temperatures T _C and T _D after the transfer processing of the analysis surfaces C and _D are calculated.

  In this way, for each delivery cycle, calculate how many times the air temperature on the upstream analysis surface changed during the time required for the air to move to the next analysis surface, and the calculation The temperature is replaced with the temperature calculated in the next analysis plane in the delivery cycle to obtain the air temperature of the analysis plane.

  The above explanation is for the case where the air in the cave is moving from the analysis surface A side toward the analysis surface D, but the air temperature delivery process when moving from the analysis surface D side toward the analysis surface A is illustrated in FIG. This is performed according to the flowchart of No. 44.

That is, first, in step S71, is set to a preset initial temperature Tao of air temperature T _D in the analysis plane D on the most upstream side. Based on the initial temperature Tao, the air temperature TC of the analysis surface _C is calculated in step S72.

Calculation of the temperature T _C is carried out according to a subroutine shown by a flow chart of FIG. 45, first, at step S81, it reads from the analysis plane D, the memory 17 are recorded in advance and travel times Δtdc of air between C Then, in step S82, i is calculated for which the movement time Δtdc is larger than iΔt and smaller than (i + 1) Δt.

In step S83, the air temperature T _Di counted from the reference time tp of the analysis surface D and calculated in the i-th analysis cycle (time: tp + iΔt) and the (i + 1) -th analysis cycle (time: tp + (i + 1)) The air temperature TD _{(i + 1)} calculated by Δt) is read from the memory 17 and, similarly to the case described above, these values are linearly interpolated to calculate the air temperature at time (tp + Δtdc) on the analysis plane D, This will be the air temperature T _C of the temperature transfer process execution period in the analysis plane C.

Similarly, in step S73 of the flowchart of FIG. 44, the air movement time Δtcb between the analysis surfaces C and B is shown, and in step S74, the air movement time Δtba between the analysis surfaces B and A is shown. performs the same processing as the processing shown in the flowchart of 45, the analysis plane B, the air temperature after passing process a T _B, calculates the T _a. As a result, the intra-cave air temperature delivery process in step S36 in the flowchart of FIG. 39 is completed.

  Further, in the humidity transfer process of step S36, when the air in the cave is moving from the analysis surface A side toward the analysis surface D, the analysis surface D side is moved toward the analysis surface A according to the flowchart of FIG. If it is moving, it is performed according to the flowchart of FIG. Steps S91 to S94 and S101 to S104 in these flowcharts are exactly the same as steps S51 to S54 and S71 to S74 in the flowcharts of FIGS. 41 and 44 showing the temperature transfer process. The humidity transfer process is performed in the same manner as the temperature transfer subroutine between the analysis surfaces shown in FIGS. 42, 44, 45, etc., and the humidity transfer process is performed in the same operation as the temperature transfer process. Is called.

  The cooling water temperature delivery process is also performed in substantially the same manner. This process is performed by a subroutine shown in the flowchart of FIG. 48. First, in step S111, the inlet temperature Two is substituted as the temperature of the Go pipe on the analysis surface A. To do.

  Next, in step S112, the temperature of the Go tube on the analysis surface B is calculated. This calculation is performed as follows according to the flowchart shown in FIG.

  That is, in step S121, the time Δt′ab during which the cooling water moves from the Go pipe on the analysis surface A to the Go pipe on the analysis surface B, which is calculated and recorded in the memory 17 in advance, is read. Find i for which the movement time Δt′ab is larger than iΔt and smaller than (i + 1) Δt. That is, the number of times and the number of analysis cycles are counted from the time when the cooling water is started to obtain the time Δ′tab.

  Now, for example, if Δt′ab = 1.6Δt, i = 1, and this travel time Δt′ab is equal to the first analysis cycle and the second analysis cycle after the temperature delivery reference time tp. Will pass between.

Then, after obtaining i, in step S123, the temperature T _{AGi at} the i-th analysis cycle (time: tp + iΔt) of the Go tube on the analysis surface A and the (i + 1) -th analysis cycle (time: tp + ( i + 1) The temperature TAG ₍ i + 1) at Δt) is read from the memory 17 and these values are linearly interpolated according to the value of the movement time Δt′ab, and the value obtained thereby is converted into the value obtained by calculation. Instead, the Go tube temperature T _GB of the analysis surface B is used.

The calculation of the Go tube temperature T _CG of the analysis surface C and the calculation of the Go tube temperature T _DG of the analysis surface D in steps S113 and S114 of the flowchart of FIG. 48 are also performed by a subroutine similar to the subroutine shown in the flowchart of FIG. The values obtained by interpolating the Go tube temperatures on the upstream analysis surfaces B and C according to the movement times Δt′bc and Δt′cd between the respective analysis surfaces are the Go tube temperatures on the analysis surface C in the temperature transfer period. Let T _{CG be} the Go tube temperature T _DG of the analysis surface D.

  Further, the temperature transfer at the turning point of the cooling water pipe line in step S115 in the flowchart of FIG. 48, that is, the temperature transfer from the Go pipe to the Re pipe in the analysis plane D is performed by the subroutine shown in the flowchart of FIG. To be done.

First, in step S131, it is determined how the cooling water temperature is transferred from the Go pipe to the Re pipe in the final analysis plane D for the cooling water pipe on the screen W10 in FIG. If it is set to be the same as the Re tube temperature, the Go tube temperature T _DG of the analysis surface D obtained by linear interpolation from the temperature of the analysis surface C is set as the Re tube temperature T of the analysis surface D in step S132. Also adopted as _DR .

  On the other hand, if it is set that a certain amount of time Δt′dd is required for the movement of the cooling water from the Go pipe to the Re pipe on the final analysis surface D, the movement time Δt′dd is read in step S133, In S134, i is obtained for which the time Δt′dd is larger than iΔt and smaller than (i + 1) Δt.

Then, after obtaining i, in step S135, the temperature T _{DGi at} the i-th analysis cycle (time: tp + iΔt) of the Go tube of the analysis surface D and the (i + 1) -th analysis cycle (time: tp + ( The temperature T _{DG (i + 1) at} i + 1) Δt) is read from the memory 17, and these values are linearly interpolated according to the value of the movement time Δt′dd, and the obtained value is used as the Re tube temperature on the analysis plane D. _TDR .

  Further, in step S116 of the flowchart of FIG. 48, the temperature of the Re tube on the analysis surface C is performed as follows according to the subroutine shown in the flowchart of FIG.

  That is, in step S141, the time Δt′dc in which the cooling water moves from the Re pipe on the analysis surface D to the Re pipe on the analysis surface C is read from the memory 17, and then in step S142, the movement time Δt′dc is i which is larger than iΔt and smaller than (i + 1) Δt is obtained.

Then, in step S143, the temperature T _DRi in the i-th analysis cycle (time: tp + iΔt) of the Re tube of the analysis surface D and the temperature T in the (i + 1) -th analysis cycle (time: tp + (i + 1) Δt). _{DR (i + 1)} is read from the memory 17 and these values are linearly interpolated according to the value of the movement time Δt′cd to calculate the temperature at the time (tp + Δt′dc) of the Re tube on the analysis surface D. It is referred to as the cooling water temperature T _CR temperature transfer processing execution cycle of Re tube analysis plane C.

Further, the calculation of the Re tube temperature T _{BR on} the analysis surface B and the calculation of the Re tube temperature T _{AR on} the analysis surface A in steps S117 and S118 in the flowchart of FIG. 48 are performed by the same subroutine as the subroutine shown in the flowchart of FIG. The values obtained by interpolating the Re tube temperatures on the upstream analysis planes C and B according to the movement times Δt′cb and Δt′ba between the respective analysis planes are obtained on the analysis plane B in the temperature transfer cycle. The Re tube temperature T _BR and the Re tube temperature T _AR of the analysis surface A are assumed.

  As described above, the cooling water temperature delivery process in step S38 in the flowchart of FIG. 39 is completed.

  Here, the entire simulation operation including the air temperature and humidity delivery process and the cooling water temperature delivery process will be described in detail.

  Now, as shown in FIG. 24, a cooling water pipe 6 extending from the cooling water supply device 7 into the sinus 1 where the cable 4 is arranged and making a U-turn and returning to the cooling water supply device 7 is provided. It is assumed that four analysis surfaces A to D crossing the sinus 1 are set in the section where the cooling water pipe 6 is disposed.

  First, air movement times Δtab, Δtbc, Δtcd between adjacent analysis surfaces and cooling water movement times Δt′ab, Δt′bc, Δt′cd, Δt′dd, Δt′dc between the analysis surfaces. , Δt′cb, Δt′ba, the longest longest travel time Δtc, based on the distance between the analysis surfaces and the wind speed for air, and for the cooling water, the distance and flow rate between the analysis surfaces. Each is calculated based on the pipe cross-sectional area.

  Here, the air in the cave is assumed to flow from the analysis surface A toward the analysis surface D. Furthermore, of the cooling water movement time, Δt'ab, Δt'bc, and Δt'cd indicate the movement time of the Go tube, and Δt'dc, Δt'cb, and Δt'ba indicate the movement time of the Re tube. The movement time between the surfaces is the same for the Go tube and the Re tube. Δt′dd is the time from when the Go tube passes through the analysis surface D to when it makes a U-turn and passes through the analysis surface D through the Re tube.

Then, when these movement times are shown in the form of how many times the analysis period Δt, it is assumed that, for example, the following occurs.
(Air movement time when ventilation starts)
Δtab = 0.4Δt
Δtbc = 0.2Δt
Δtcd = 0.6Δt
(Air movement time when ventilation is stopped)
Δtab = 3.2Δt
Δtbc = 1.6Δt
Δtcd = 4.8Δt
(Cooling water travel time at the start of cooling)
Δt′ab = Δt′ba = 1.6Δt
Δt′bc = Δt′cb = 0.8Δt
Δt′cd = Δt′dc = 2.4Δt
Δt′dd = 1.2Δt

  In addition, when the wind speed at the time of a ventilation stop is set to 0, calculation of the air movement time at the time of a ventilation stop is not performed.

In the case of the above example, the longest movement time Δtc is as follows in each state.
At ventilation start-up and cooling start-up: 2.4 Δt (movement time between analysis surfaces C and D of cooling water)
When ventilation is started and when cooling is stopped: 0.6 Δt (travel time between air analysis surfaces C and D)
At the time of ventilation stop and cooling start: 4.8Δt (movement time between air analysis surfaces C and D)
At the time of ventilation stop and cooling stop: 4.8 Δt (movement time between air analysis surfaces C and D)

  When the wind speed when the ventilation is stopped is 0, the longest movement time Δtc at the start of cooling becomes the longest movement time 2.4 Δt of the cooling water, and the temperature transfer process is not performed when the cooling is stopped.

  Further, it is assumed that ventilation and cooling are started and stopped according to the schedule shown in the time chart of FIG. 52, for example, after the start of energization of the cable due to the establishment of the ventilation schedule, the cooling start condition, and the stop condition.

Furthermore, at the start of the analysis, the air in the cave is in a ventilation start-up state, and the influence of energization on the cable has not yet occurred, so the initial temperature Tao is set to the air temperatures T _{A0 to} T _D0 on the analysis surfaces A to D, The initial humidity Hao is set to each of the humidity temperatures H _{A0 to} H _D0 . In this case, as the initial temperature Tao, a value set on the screen W9, a ground temperature corresponding to the depth of the point read from the ground temperature database DB10, or an air temperature read from the weather database DB11 is set. Further, as the initial humidity Hao, the value set on the screen W9 or the atmospheric humidity read from the weather database DB11 is set. Note that when setting the initial temperature and the initial humidity on the screen W9, different values can be set for each analysis plane.

Also, the cooling water temperature _{T AG0 ~T DG0, T DR0 ~T} AR0 of Go tubes and Re tubes each analysis surface to D, since the cooling water is not moved, the soil temperature in accordance with the depth or the like Each equal value is set.

  The calculation results of the air temperature, humidity, and cooling water temperature of each analysis surface by the delivery process under the above conditions are as shown in FIG. 53, but the air temperature and humidity are delivered at exactly the same timing and with the same calculation contents. Since the process is performed, only the temperature transfer process for air will be described below.

  First, after the initial value is set as described above, in the state where only ventilation is started from the analysis start time (time t = 0) to the third analysis cycle (time t = 3Δt), the longest movement time Δtc is Since the time for air to move between the analysis surfaces C and D is Δtcd = 0.6Δt, the air temperature (and humidity) is transferred every analysis cycle.

That is, as shown in FIG. 53, in the first analysis cycle Δt, the initial temperature Tao is set to the air temperature T _A1 of the analysis surface A, and the calculated temperature T _B1 of the analysis surface B is Instead, since the first temperature T _A0 (Tao) on the analysis surface A and the temperature T _A1 obtained at the time Δt obtained by calculation are Δtab = 0.4Δt, they are interpolated at an internal ratio of 4: 6. The value is _TA0-1 .

Similarly, the temperature T _C1 of the analysis plane _C is the ratio of the initial temperature T _{B0 on} the analysis plane B and the temperature T _B1 obtained at the time Δt obtained by the calculation because Δtbc = 0.2Δt. The value T _B0-1 interpolated by the fractional ratio is _set, and the temperature T _D1 of the analysis surface _D is obtained by setting the first temperature T _{C0 on} the analysis surface C and the temperature T _C1 obtained by the calculation at time Δt to Δtcd = 0. Since 6Δt, the value _TC0-1 interpolated with an internal division ratio of 6: 4 is obtained.

In the next analysis cycle 2Δt, the initial temperature Tao is set again to the air temperature T _A2 of the analysis surface A, and the temperature T _B2 of the analysis surface B is changed to the calculated temperature and analyzed. The temperature T _A1 calculated in the first analysis cycle on the surface A and the temperature T _A2 calculated in the current analysis cycle are interpolated at an internal division ratio of 4: 6, which is a value T _A1-2 .

Similarly, the temperature T _C1 of the analysis surface C is changed to the calculated temperature, and the temperature T _B1 calculated in the first analysis cycle on the analysis surface B and the temperature T _B2 calculated in the current analysis cycle are obtained. , 2: 8 interpolated value for at interior division ratio of the a T _B1-2, the temperature T _D2 are present analysis period the temperature T _C1 calculated in the analysis period of the first in the analysis plane C of the analysis plane D The temperature T _C2 calculated in step _S2 is interpolated with an internal division ratio of 6: 4 to obtain a value T _C1-2 .

Further, in the third analysis cycle 3Δt, the initial temperature Tao is similarly set to the air temperature T _A2 of the analysis surface A, and the temperatures T _B3 , T _C3 , and T _D3 of the analysis surfaces B to D are Instead of the calculated temperatures, the values T _{A2-3 are obtained} by interpolating the temperatures calculated in the second and third analysis periods in the upstream analysis planes A to C with the internal ratios corresponding to the respective moving times. , T _B2-3 , and T _C2-3 .

Since the cooling water does not move during the period from the start of the analysis to the third analysis cycle, as shown in FIG. 53, the Go tube temperatures T _{AG1 to} T _DG1 and T _{AG2 to} T on each analysis surface _{DG2, T} AG3 _{through T DG3,} and Re pipe temperature _{T DR1 ~T AR1, T DR2 ~T} AR2, T DR3 ~T AR3 , the value that is calculated independently for each analysis surface is set.

Next, from the third analysis cycle 3Δt to the ninth analysis cycle 9Δt, cooling is started in a state where ventilation is started, and both the air in the cave and the cooling water move. In this case, first, at the time 3? T, it is set inlet temperature Two of the cooling water in the Go pipe temperature _{T AG3} of the cooling water in the analysis plane A.

  In this case, the longest movement time Δtc of air or cooling water is the movement time Δt′cd (Δt′dc) = 2.4Δt of the cooling water between the analysis surfaces C and D. The delivery process is performed.

  Therefore, as shown in FIG. 53, in the fourth analysis cycle 4Δt, the fifth analysis cycle 5Δt, the seventh analysis cycle 7Δt, and the eighth analysis cycle 8Δt, both the air in the cave and the cooling water are analyzed. The temperature is calculated independently.

That is, the air temperature T _{A4 to} T _D4 of each analysis surface in the fourth analysis cycle, the Go tube temperature T _{AG4 to} T _{DG4 of} the cooling water, the Re tube temperature T _{DR4 to} T _AR4 , and each analysis surface in the fifth analysis cycle. Air temperature T _{A5 to} T _D5 , cooling water Go tube temperature T _{AG5 to} T _DG5 , Re tube temperature T _{DR5 to} T _AR5 , air temperature T _{A7 to} T _D7 of each analysis surface in the seventh analysis cycle, cooling water of Go pipe temperature _T AG7 _~T DG7, Re pipe temperature _T DR7 _{through T AR7,} and eighth analyzes air temperature _T of each analysis plane in period A8 _{through T D8,} Go pipe temperature of the cooling water _T AG8 _{through T DG8,} As for the Re tube temperatures _{TDR8 to TAR8} , values independently calculated on each analysis surface are adopted as they are.

  Then, in the sixth cycle and the ninth analysis cycle, temperature transfer processing is performed on each analysis plane.

That is, first, in the sixth analysis cycle, the initial temperature Tao is set for the air temperature T _A6 of the analysis surface A, and the air temperatures T _B6 , T _C6 , and T _D6 of the analysis surfaces B to _D are Since the travel times Δtab, Δtbc, and Δtcd from the upstream analysis surface are all less than Δt, the temperature after the transfer process obtained in the third analysis cycle 3Δt on the upstream analysis surface and the fourth analysis are performed. values of temperature and calculated in each analysis surface in a cycle and interpolated by interior division ratio corresponding to each of the movement time _{T A3-4, T B3-4,} are _{T C3-4.}

Further, the cooling water, the Go pipe temperature _{T AG6} analysis plane A, with the inlet temperature Two is set, Go pipe temperature _{T BG6} analysis plane B, since a Δt'ab = 1.6Δt, calculated for analyzing surface a in the analysis period of the fourth temperature T _AG4 and fifth temperature T _AG5 calculated in the analysis period, 6: value T _AG4-5 next interpolated by interior division ratio of 4, the analysis surface Since the Go tube temperature T _CG6 of C is Δt′bc = 0.8Δt, the temperature T _BG3 calculated for the analysis surface B in the third analysis cycle and the temperature T _BG4 calculated in the fourth analysis cycle are , 8: 2 interpolated value T _BG3-4 , and Go tube temperature T _DG6 of analysis surface D is Δt′cd = 2.4Δt, so analysis surface C is analyzed in the fifth analysis cycle. calculated temperature T _{CG 5} and the sixth analysis cycle The calculated temperature _{T CG 6} was 4: a value _{T CG5-6} interpolated with 6 internal ratio of.

The Re tube temperature is calculated at the fourth analysis cycle for the Go tube on the upstream analysis surface D because Δt′dd = 1.2Δt for the Re tube temperature T _DR6 on the analysis surface D. value _{T DG4} a fifth and an analysis period in calculated values _{T DG5} 2: a value _{T DG4-5} interpolated with 8 internal division ratio. Note that Re pipe temperature T _DR6 analysis plane D, by previously setting, in some cases set to the same temperature as the Go pipe temperature T _DG6 in the analysis plane D of the upstream side.

Hereinafter, the analysis plane Re pipe temperature _T CR6 of _{C~A, T BR6,} T _AR6 likewise, the value interpolated according Analysis surface of the upstream side to each movement time _{T DR5-6 (Δt'dc} = 2.4Δt), _TCR3-4 (Δt′cb = 0.8Δt), and _TBR4-5 (Δt′ba = 1.6Δt).

Even with the analysis period of the 9th, just like the sixth analysis cycle, the initial temperature Tao, inlet temperature Two are set respectively to the air temperature T _A9 and coolant Go pipe temperature T _AG9 analysis surface A At the same time, for other analysis surfaces, the temperature on the upstream analysis surface is replaced with a value interpolated according to the movement time.

  Furthermore, from the ninth analysis period 9Δt to the fifteenth analysis period 15Δt, cooling starts and ventilation is stopped, but it is assumed that a natural air flow is generated in the sinus.

  In this case, the longest movement time Δtc of air or cooling water is air movement time movement time Δtcd = 4.8Δt when the ventilation between the analysis planes C and D is stopped. Processing will be performed.

Therefore, as shown in FIGS. 53 and 54, in the 10th analysis cycle 10Δt, the 11th analysis cycle 11Δt, the 12th analysis cycle 12Δt, and the 13th analysis cycle 13Δt, The temperature is calculated independently for each analysis plane. In other words, the air temperature _{T A10~ T D10, T A11~ T} D11, T A12~ T D12, T A13~ T D13 of each analysis period, Go pipe temperature _T of the cooling water _{AG10 ~T DG10, T AG11 ~T DG11} , _{T AG12 ~T DG12, T AG13 ~T} DG13, and Re pipe temperature of the cooling water _{T DR10 ~T AR10, T DR11 ~T} AR11, T DR12 ~T AR12, T DR13 ~T AR13 are all in each analysis surface The value calculated independently is adopted as it is.

  Then, in the fifth cycle after the ventilation stops, that is, in the 14th cycle 14Δt from the beginning, a temperature transfer process is performed on each analysis surface.

That is, first, the initial temperature Tao is set for the air temperature T _A14 on the analysis surface A, and the air temperatures T _B14 , T _C14 , and T _D14 on the analysis surfaces B to _D are from the upstream analysis surface. Values T _A12-13 (Δtab = 3.2Δt), T _B10-11 (Δtbc = 1.6Δt), T B interpolated from the temperature calculated on the upstream analysis surface according to the movement times Δtab, Δtbc, Δtcd _C13-14 (Δtcd = 4.8Δt).

Further, the cooling water, the Go pipe temperature _{T AG14} analysis plane A, with the inlet temperature Two is set, Go pipe temperature _T BG14 analysis surface _{B~D, T CG14,} T _DG14 is upstream The values T _AG12-13 (Δt′ab = 3.2Δt), T _BG10-11 (Δt′bc = 1.6Δt) _obtained by interpolating the temperature calculated on the upstream analysis surface according to the moving time from the analysis surface. _), is a _{T DG13-14 (Δt'cd} = _4.8Δt), Similarly, Re pipe temperature _{T DR14, T CR14, T BR14} , T AR14 _{is, T DG10-11} (or _{T DG14), T DR13 −14} , _TDCR10-11 , and _TBR12-13 .

  And cooling is also stopped at the 15th analysis cycle 15Δt. In this case, since the air is moving at the wind speed when the ventilation is stopped, the air temperature is changed every 5 analysis cycles after time 14Δt. A delivery process similar to that at time 14Δt is performed. In this case, the cooling water temperature delivery process is not performed.

  When the wind speed when the ventilation is stopped is 0, the longest travel time Δtc (cooling water travel time Δt′cd between the analysis surfaces C and D) between the time 9Δt and the time 15Δt when the cooling stops is Based on 2.4 Δt), the cooling water delivery process is performed every three cycles. In this case, the air temperature delivery process is not performed. Then, after the time 15 Δt when the cooling is stopped, neither the air temperature nor the cooling water temperature is transferred, and a value calculated independently on each analysis plane is adopted for each analysis period.

  Further, in a model in which a vent is provided in a cave, the temperature Ta and humidity Ha of the cave air are always forced to the temperature and humidity of the atmosphere read from the weather database DB11 in the analysis plane closest to the vent. Set to

  The above is a description of the case where the cooling model with the cooling water pipe is selected on the screen W1 in FIG. 21, but if the non-cooling model without the cooling water pipe is selected, the processing is performed according to the flowchart of the main routine shown in FIG. However, as compared with the main routine of the cooling model in FIG. 22, the flowchart for the non-cooling model does not include the steps for setting the cooling conditions in steps S8 and S9 of the flowchart of the cooling model.

  With this difference, in the setting of the cave configuration in step S152, as shown in the screen W17 in FIG. 55, the cooling water pipeline is not set, and the cooling water data is set on the cooling model screen W10. Not done.

  Analysis plane setting, soil data setting, sinus road data setting, latent heat data setting, air data setting, energization condition setting, regional data setting, analysis of other steps S151, S153 to S155, S157 to S160 Conditions are set using the same screens as the cooling model screens W5 to W9 and W11 to W13 in the same manner as in the cooling model.

  Further, although the output of the analysis result in step S162 is performed in the same manner, the analysis processing in step S161 is performed by a subroutine whose flowchart is shown in FIG.

  When this subroutine is compared with the subroutine shown in FIG. 39 for the cooling model, the subroutine for the non-cooling model is different from the cooling model in that only air movement is determined in steps S177 and S183. The difference is that there is no cooling water temperature transfer process in steps S37 and S38 of the subroutine.

  On the other hand, the calculation of the node temperature for each analysis surface in step S181 is performed by a subroutine similar to the subroutine shown in the flowchart of FIG. 40 of the cooling model, and the air temperature / humidity transfer processing in step S186 is also performed in the cooling process. The model is executed by a subroutine similar to the subroutine shown in the flowchart in FIGS. 41, 42, and 44 to 47.

  Therefore, the specific calculation example of the delivery process shown in FIG. 53 is exactly the same for the air temperature and humidity, and is the same as the cooling model except that the calculation for the cooling water temperature is not performed. It will be.

  As described above, according to the present invention, the change in the conductor temperature of the power cable disposed in the underground cave can be simulated in consideration of the latent heat effect due to water evaporation on the air contact surface in the cave. Is possible. This makes it possible to efficiently design and operate power cables by avoiding excessive cable energization current suppression and cable size increase when designing power cable lines, etc. Will contribute.

It is explanatory drawing of the structural example of the analysis model used as the object of this invention. It is explanatory drawing of the cooling water pipe line arrange | positioned in a cave. 1 is a system diagram showing an overall configuration of an embodiment of the present invention. It is a block diagram of the soil thermal constant database used in the embodiment. It is a block diagram of a cave wall thermal constant database. It is a block diagram of a pipeline wall thermal constant database. It is a block diagram of a fluid thermal constant database. It is a block diagram of a pipeline dimension database. It is a block diagram of a cable database. It is a block diagram of an air physical property database. It is a block diagram of an evaporation efficiency database. It is a block diagram of a cooling water physical property database. It is a block diagram of a ground temperature database. It is a block diagram of a weather database. It is a block diagram of a load factor database. It is explanatory drawing of a heat capacity matrix and a heat conductive matrix. It is explanatory drawing of the finite element analysis model inside a cave. It is explanatory drawing of the finite element model of a cooling water pipe line. It is explanatory drawing of the turbulent flow determination data of the air in a cave. It is explanatory drawing of the turbulent flow determination data of a cooling water. It is explanatory drawing of an analysis model selection screen. It is a flowchart of the analysis main routine of a cooling model. It is explanatory drawing of an analysis surface setting screen. It is explanatory drawing of an analysis surface. It is explanatory drawing of a cave structure setting screen. It is explanatory drawing of the screen which shows a whole analysis model. It is explanatory drawing of the screen which shows the principal part of an analysis model. It is explanatory drawing of a soil data setting screen. It is explanatory drawing of a cave data setting screen. It is explanatory drawing of a latent heat data setting screen. It is explanatory drawing of an air data setting screen. It is explanatory drawing of a cooling water data setting screen. It is explanatory drawing of an energization condition setting screen. It is explanatory drawing of a regional data setting screen. It is explanatory drawing of an analysis condition setting screen. It is explanatory drawing of the screen which shows the temperature distribution around a cave as an analysis result. It is explanatory drawing of the screen which similarly shows the change of the air temperature in a cave. It is explanatory drawing of the screen which similarly shows the change of cable conductor temperature. It is a flowchart which shows the analysis processing operation by a computer. It is a flowchart of a node temperature calculation subroutine. It is a flowchart of an air temperature delivery processing subroutine. 10 is a flowchart of an air temperature calculation subroutine for analysis surface B. It is explanatory drawing of the interpolation method which calculates | requires delivery temperature. It is a flowchart of a reverse direction air temperature delivery process subroutine. 4 is a flowchart of an air temperature calculation subroutine for analysis plane C. It is a flowchart of a humidity delivery processing subroutine. It is a flowchart of a humidity transfer process subroutine in the reverse direction. It is a flowchart of a cooling water temperature delivery processing subroutine. It is a flowchart of a Go tube temperature calculation subroutine for analysis surface B. 5 is a flowchart of a Re tube temperature calculation subroutine for analysis surface D. 10 is a flowchart of a Re tube temperature calculation subroutine for analysis plane C. It is a time chart which shows an example of the start-stop schedule of ventilation and cooling. It is a table | surface which shows an example of the temperature and humidity delivery in each analysis surface. It is a flowchart of the analysis main routine of a non-cooling model. It is explanatory drawing of a cave structure setting screen. It is a flowchart which shows the analysis processing operation by a computer.

Explanation of symbols

DESCRIPTION OF SYMBOLS 1 Cave 4 Cable 6 Cooling water pipe 10 Computer 11 Central processing unit 12 Input device 13 Reading device 14 Recording device 15 Display device 16 Printing device 17 Memory DB1-DB12 Database AD Analysis surface

Claims (14)

A method for estimating a conductor temperature of a power cable laid in a tunnel in the ground using a finite element method by a computer,
Evaporation related to thermal constants of cables, sinus air and surrounding components, heat transfer coefficient between air interface in the sinus and the sinus air, and evaporation of water at the air interface in the sinus Value and
Calculating the amount of heat transfer between the air contact surface and air and the amount of latent heat transport in a predetermined region of the air contact surface;
A cable conductor temperature estimation method in consideration of a latent heat effect in a tunnel, wherein a change in conductor temperature of the cable is estimated from the heat transfer amount and latent heat transport amount and a heat generation amount of the cable. The method of claim 1, wherein:
The temperature and humidity of the air in the cave used for the calculation of the heat transfer amount and the latent heat transport amount take into account the movement of the air in the cave and use values reflecting the upstream temperature and humidity. Cable conductor temperature estimation method considering latent heat effects in the road. A system for estimating the conductor temperature of a power cable laid in the underground cave by using a cross section crossing the cave as an analysis surface and applying a finite element method to a predetermined region of the analysis surface,
A cave configuration setting means for setting a configuration inside and outside the cave including the arrangement of cables; and
Latent heat transport area setting means for setting an area in consideration of latent heat transport due to water evaporation in the air contact surface in the cave,
Thermal constant setting means for setting the thermal constants of the cable, the air in the cave and the surrounding components of the cave;
Energizing current setting means for setting the energizing current to the cable;
Heat transfer coefficient setting means between the air contact surface in the sinus and the air in the sinus;
An evaporation related value setting means for setting an evaporation related value related to the evaporation of water at the air contact surface in the sinus;
A calorific value calculating means for calculating a calorific value of the cable based on the energizing current set by the energizing current setting means;
An air heat transfer amount calculating means for calculating a heat transfer amount between the air contact surface in the cave and the air in the cave using the heat transfer coefficient set by the heat transfer coefficient setting means;
A latent heat transport amount calculating means for calculating a latent heat transport amount by water evaporation in the latent heat transport region set by the region setting means using the evaporation related value set by the evaporation related value setting means;
Model creation means for creating an analysis model by dividing a predetermined analysis region of the analysis surface crossing the sinus finite element based on the configuration set by the sinus configuration setting means,
Matrix creation means for creating a heat capacity matrix and a heat conduction matrix for each node in the analysis region using the thermal constant recorded in the thermal constant recording means;
The calorific value of the cable calculated by the calorific value calculation means, the heat transfer quantity between the air contact surface and the air in the cave calculated by the air heat transfer quantity calculation means, and the latent heat transport quantity calculation means A vector creating means for creating a thermal load vector for each node in the analysis region using a latent heat transport amount;
Node temperature calculation means for calculating each node temperature in the analysis region at a predetermined time interval from a predetermined initial state using the heat capacity matrix, the heat conduction matrix and the thermal load vector generated by the matrix generation means and the vector generation means; A cable conductor temperature estimation system in consideration of the latent heat effect in the cave characterized by having. The system of claim 3, wherein
After calculating the latent heat transport amount due to the evaporation of water by the latent heat transport amount calculating means, the humidity update means for calculating the humidity of the air in the cave after the evaporation and updating the humidity value is provided,
The latent heat transport amount calculation means calculates the latent heat transport amount using the humidity updated by the humidity update means, and the cable conductor temperature estimation system in consideration of the latent heat effect in the cave. In the system according to claim 3 or claim 4,
Analysis plane setting means for setting a plurality of analysis planes for finite element analysis in a predetermined section of the tunnel,
Air data setting means for setting data relating to air movement in the cave,
The air heat transfer amount calculation means uses the air temperature obtained on the upstream analysis surface according to the state of air movement set by the air data setting means as the air temperature used in the heat transfer amount calculation,
The latent heat transport amount calculating means calculates the air temperature and humidity on the upstream analysis surface as the air temperature and humidity used in calculating the latent heat transport amount according to the air movement state set by the air data setting means. A cable conductor temperature estimation system that takes into account the latent heat effect in the cave characterized by its use. In the system according to claim 3 or claim 4,
Cooling water pipe setting means for setting a cooling water pipe in the cave,
A heat transfer coefficient setting means between the inner surface of the cooling water pipe set by the cooling water pipe setting means and the cooling water in the pipe;
Cooling water heat transfer amount calculating means for calculating the heat transfer amount between the inner surface of the cooling water pipe line and the cooling water in the pipe line using the heat transfer coefficient set by the heat transfer coefficient setting means,
The thermal constant recording means records the thermal constant of the cooling water pipe and the cooling water in the pipe,
The cave configuration setting means sets the configuration inside and outside the cave including the configuration of the cooling water pipeline,
The vector creating means creates a thermal load vector using the heat transfer amount between the cooling water pipe inner surface and the cooling water calculated by the cooling water heat transfer amount calculating means, and the latent heat effect in the cave is obtained. Cable conductor temperature estimation system in consideration. The system of claim 6, wherein:
Analysis plane setting means for setting a plurality of analysis planes for finite element analysis in a predetermined section of the tunnel,
Air data setting means for setting data relating to air movement in the cave;
Cooling water data setting means for setting data relating to the movement of the cooling water in the cooling water pipeline,
The air heat transfer amount calculation means uses the air temperature obtained on the upstream analysis surface according to the state of air movement set by the air data setting means as the air temperature used in the heat transfer amount calculation,
The latent heat transport amount calculating means calculates the air temperature and humidity on the upstream analysis surface as the air temperature and humidity used in calculating the latent heat transport amount according to the air movement state set by the air data setting means. Use
The cooling water heat transfer amount calculation means is the cooling water temperature used in the heat transfer amount calculation, and the cooling obtained on the upstream analysis surface according to the state of movement of the cooling water set by the cooling water data setting means. A cable conductor temperature estimation system that takes into account the latent heat effect in a cave characterized by using water temperature. The system according to any one of claims 3 to 7,
Meteorological data recording means for recording the weather data for each period in each region;
A heat balance amount calculating means for reading out the weather data of the designated analysis target area and time from the weather data recording means, and calculating a heat balance amount of the ground surface based on the weather data; and
A cable conductor temperature estimation system in consideration of a latent heat effect in a cave, wherein the vector creation means creates a thermal load vector using the ground surface heat balance calculated by the heat balance calculation means. A program for estimating the conductor temperature of a power cable laid in a underground cave by using a cross section crossing the cave as an analysis surface and applying a finite element method to a predetermined region of the analysis surface,
Computer
A cave configuration setting means for setting a configuration inside and outside the cave including the arrangement of cables,
Latent heat transport area setting means for setting an area in consideration of latent heat transport by evaporation of water in the air contact surface in the cave,
Thermal constant setting means for setting thermal constants of cables, air in the cave and surrounding components of the cave,
Energizing current setting means for setting the energizing current to the cable,
A heat transfer coefficient setting means between the air contact surface in the sinus and the air in the sinus;
An evaporation related value setting means for setting an evaporation related value related to the evaporation of water at the air contact surface in the sinus;
A calorific value calculating means for calculating a calorific value of the cable based on the energizing current set by the energizing current setting means;
Air heat transfer amount calculating means for calculating the heat transfer amount between the air contact surface in the cave and the air in the cave using the heat transfer coefficient set by the heat transfer rate setting means,
A latent heat transport amount calculating means for calculating a latent heat transport amount by evaporation of water in the latent heat transport region set by the region setting means using the evaporation related value set by the evaporation related value setting means;
Model creation means for creating an analysis model by dividing a predetermined analysis region of the analysis surface crossing the sinus finite element based on the configuration set by the sinus configuration setting means,
Matrix creation means for creating a heat capacity matrix and a heat conduction matrix for each node in the analysis region using the thermal constant recorded in the thermal constant recording means;
The calorific value of the cable calculated by the calorific value calculation means, the heat transfer quantity between the air contact surface and the air in the cave calculated by the air heat transfer quantity calculation means, and the latent heat transport quantity calculation means A vector creating means for creating a thermal load vector for each node in the analysis region using the latent heat transport amount, and
Using the heat capacity matrix, the heat conduction matrix and the thermal load vector created by the matrix creating means and the vector creating means, function as a node temperature calculating means for calculating each node temperature in the analysis region at a predetermined time interval from a predetermined initial state. A cable conductor temperature estimation program that takes into account the latent heat effect in the tunnel. In the program according to claim 9,
Computer
After calculating the latent heat transport amount due to the evaporation of water by the latent heat transport amount calculating means, calculate the humidity of the air in the cave after the evaporation and function as a humidity update means for updating the humidity value,
When functioning as the latent heat transport amount calculating means, the cable conductor temperature in consideration of the latent heat effect in the cave characterized by functioning to calculate the latent heat transport amount using the humidity updated by the humidity update means Estimation program. In the program according to claim 9 or 10,
Computer
Analysis plane setting means for setting a plurality of analysis planes for finite element analysis in a predetermined section of the cave, and
While functioning as air data setting means for setting data relating to air movement in the cave,
When functioning as the air heat transfer amount calculation means, the air temperature used in the heat transfer amount calculation is the air temperature obtained on the upstream analysis surface according to the air movement state set by the air data setting means. Function to use temperature,
When functioning as the latent heat transport amount calculation means, the air temperature and humidity used in the calculation of the latent heat transport amount as the air temperature and humidity on the upstream analysis surface according to the air movement state set by the air data setting means. A cable conductor temperature estimation program considering a latent heat effect in a cave characterized by functioning to use temperature and humidity. In the program according to claim 9 or 10,
Computer
Cooling water pipe setting means for setting the cooling water pipe in the cave,
A heat transfer coefficient setting means between the inner surface of the cooling water pipe set by the cooling water pipe setting means and the cooling water in the pipe; and
While functioning as a cooling water heat transfer amount calculation means for calculating the heat transfer amount between the inner surface of the cooling water pipe line and the cooling water in the pipe line using the heat transfer coefficient set by the heat transfer coefficient setting means,
When functioning as the thermal constant recording means, function to record the thermal constant of the cooling water pipe and the cooling water in the pipe,
When functioning the cave configuration setting means, function to set the configuration inside and outside the cave, including the configuration of the cooling water pipeline,
When functioning as the vector creating means, the thermal load vector is made to function using the heat transfer amount between the cooling water pipe inner surface and the cooling water calculated by the cooling water heat transfer amount calculating means. A cable conductor temperature estimation program that takes into account the latent heat effect in the tunnel. In the program according to claim 12,
Computer
Analysis surface setting means for setting a plurality of analysis surfaces for finite element analysis in a predetermined section of the cave,
Air data setting means for setting data relating to air movement in the cave, and
While functioning as a cooling water data setting means for setting data relating to the movement of the cooling water in the cooling water pipeline,
When functioning as the air heat transfer amount calculation means, the air temperature used in the heat transfer amount calculation is the air temperature obtained on the upstream analysis surface according to the air movement state set by the air data setting means. Function to use temperature,
When functioning as the latent heat transport amount calculation means, the air temperature and humidity used in the calculation of the latent heat transport amount as the air temperature and humidity on the upstream analysis surface according to the air movement state set by the air data setting means. Function to use temperature and humidity,
When functioning as the cooling water heat transfer amount calculation means, the cooling water temperature used in calculating the heat transfer amount is determined on the upstream analysis surface according to the state of movement of the cooling water set by the cooling water data setting means. A cable conductor temperature estimation program that takes into account the latent heat effect in the cave, which functions to use the obtained cooling water temperature. In the program according to any one of claims 9 to 13,
Computer
Meteorological data recording means for recording weather data for each region and period, and
Read the meteorological data of the specified analysis target area and time from the meteorological data recording means, function as a heat balance amount calculating means for calculating the heat balance amount of the ground surface based on the meteorological data,
When functioning as the vector creating means, the latent heat effect in the cave is made to function so as to create a thermal load vector using the heat balance amount of the ground surface calculated by the heat balance amount calculating means. A cable conductor temperature estimation program that takes into account.

Images