JP5916400B2 - Energy analysis apparatus, energy analysis method, and energy analysis program - Google Patents

Description

  The present invention relates to a technique for analyzing energy transferred between a radiation device and an absorption device, for example.

  2. Description of the Related Art Conventionally, a technique for analyzing an illumination state by light rays emitted from a light source, such as an illumination simulator, has been proposed. For example, Patent Document 1 describes that the surface shape of the illumination space is divided and the illuminance of the divided unit surface is calculated based on illumination conditions such as the position of the light source, the luminous flux, and the light distribution.

JP 7-36360 A

  However, Patent Document 1 only describes that the illuminance for each unit surface is calculated based on the illumination conditions such as the position of the light source, the luminous flux, and the light distribution. It is not described at all whether the illuminance is calculated by.

  In the first place, the energy on the irradiated surface cannot be accurately analyzed only by illumination conditions such as the position of the light source, the luminous flux, and the light distribution. Moreover, in order to measure energy such as illuminance at an actual site, it is very expensive in view of necessary facilities and labor.

  The present invention has been proposed in order to solve the above-described problems of the prior art. The purpose of the present invention is to actually measure the energy radiated from the radiation device and absorbed into the absorption device using the actual device. It is to avoid the cost of measurement and to analyze accurately by using a mathematical model.

In order to achieve the above object, the energy analysis device of the present invention has the following technical features.
(1) It has a propagation space model storage unit for storing a propagation space model, which is information obtained by modeling an energy propagation space using a three-dimensional orthogonal coordinate system.
(2) A radiation device model storage unit that stores a radiation device model that is information modeling a radiation device having a radiation element that radiates the energy to the outside.
(3) It has an absorption device model storage unit that stores an absorption device model, which is information modeling the absorption device having an absorption element that absorbs the energy.
(4) The radiation device model includes information on the position of the radiation element in the three-dimensional orthogonal coordinate system, a radiation center that is a central point of energy radiation, and energy radiation intensity.
(5) The absorber model includes information on the position, absorption surface, and absorption efficiency of the absorption element in the three-dimensional orthogonal coordinate system.
(6) It has a dividing unit that divides the absorption surface into triangles based on the position of the absorption element and information on the absorption surface.
(7) A spherical triangle calculation unit that obtains a spherical triangle center-projected on a sphere centered on the radiation center by the triangle divided by the division unit based on the position of the radiation element and the radiation center. .
(8) A solid angle calculation unit for obtaining a solid angle of the spherical triangle.
(9) having an integration unit for integrating the solid angles;
(10) It has a power calculation part which calculates | requires the power of the energy absorbed by the said absorption surface based on the integrated said solid angle, the said radiation intensity, and the said absorption efficiency.

  In the invention as described above, the energy radiated from the radiation device and absorbed by the absorption device can be analyzed using mathematical models such as a propagation space model, a radiation device model, and an absorption device model. The analysis is performed by dividing the absorption surface into triangles, obtaining a spherical triangle obtained by centrally projecting the divided triangles, and integrating the solid angles of the spherical triangles.

  In another aspect, the division unit includes a triangle generation unit that generates a triangle obtained by dividing the absorption surface, and a determination point extraction unit that extracts both ends of the longest side of the triangle generated by the triangle generation unit as determination points; The tolerance determination unit for determining whether at least one of the difference in radiation intensity or the difference in absorption efficiency between the two determination points extracted by the determination point extraction unit is within a preset allowable error range; and And a subdividing unit that generates a triangle obtained by subdividing the triangle when the allowable error determining unit determines that it is outside the allowable error range.

  In the above aspect, even when the radiation of energy does not have isotropic directivity, it is possible to analyze the radiation intensity or the absorption efficiency within the allowable error range.

  In another aspect, the apparatus includes a radiation intensity calculation unit that calculates the radiation intensity based on a total radiation flux and directivity characteristics of the radiation device.

  In the aspect as described above, even if the radiation intensity of the radiating element is undefined, the radiation intensity can be obtained and analyzed based on the total radiation flux and the directivity.

  The invention as described above can also be understood from the viewpoint of an energy analysis method and an energy analysis program.

  As described above, according to the present invention, the energy radiated from the radiation device and absorbed by the absorption device can be accurately determined by using the mathematical model while avoiding the cost of actually measuring the energy using the actual device. Can be analyzed.

It is a schematic diagram which shows the radiation | emission apparatus, absorption device, and energy propagation space in an energy analysis system. It is a schematic diagram which shows an example of a radiation device. It is a simplified perspective view which shows an example of an absorber. It is a functional block diagram showing one embodiment of an energy analysis system. It is a schematic diagram which shows an example of a spherical polygon. It is a schematic diagram which shows the internal angle of the spherical polygon of FIG. It is a schematic diagram which shows the spherical polygon which center-projected the absorption surface on the spherical surface. It is a schematic diagram which shows the internal angle of the spherical polygon of FIG. It is a flowchart which shows an example of the whole process of an energy analysis system. It is a flowchart which shows an example of a triangulation process. It is a flowchart which shows an example of a power calculation process.

Embodiments of the present invention will be described with reference to the drawings.
[1. Assumed equipment and space]
First, an apparatus and a space that the present embodiment models as an analysis target will be described with reference to FIGS.

[energy]
As shown in FIG. 1, the energy E to be analyzed has directivity and propagates through an arbitrary space. The amount of energy E can be observed and extracted as a desired physical quantity according to its characteristics. Examples of the energy E include electromagnetic waves such as light and radio waves, sound waves, and moving substances. However, as long as it has the said characteristic, it is not limited to these.

[Radiation device]
The radiation device 1 is a device that radiates energy E to be analyzed to the outside. As shown in FIG. 2, the radiation device 1 has one or a plurality of radiation elements 11. Various radiating elements 11 can be applied depending on the type of radiant energy. For example, a light source, an electromagnetic wave generation source such as a radio wave source, and a sound wave generation source such as a sound source are included. The radiation form of energy E from the radiating element 11 is point radiation or surface radiation. The radiating element 11 is modeled, for example, by defining the following elements.

(1) Placement and posture (position information)
(2) Radiation center
(2) Radiant intensity or total radiant flux
(3) Radiation vector
(4) Radiation directivity

  The radiation center is the center point of the radiation of energy E. The radiation vector is a vector representing a reference direction of radiation of energy E. The radiation directivity characteristic is a directivity characteristic normalized with respect to the radiation center and the radiation vector. The angle that is the definition area of this directivity is called the radiation angle.

  When viewed from the center of radiation, the direction in which the radiation intensity is maximum coincides with the radiation vector. In the case of isotropic directional characteristics, the normal line of the plane on which the radiating element is placed may be a radiation vector.

[Absorber]
As shown in FIG. 3, the absorption device 2 is a device that absorbs energy E to be analyzed from the outside. The absorption device 2 has one or a plurality of absorption elements 21. The absorbing element 21 is modeled by defining the following elements.

(1) Placement and posture (position information)
(2) Absorbing surface
(3) Absorption sensitivity
(4) Absorption center
(5) Absorption vector
(6) Absorption efficiency (absorption directional characteristics)

  Here, the absorption surface 22 is a flat surface and is a closed region of a simple polygon. The absorption sensitivity is the absorption sensitivity of energy E in each absorption element 21. The energy E is absorbed by the absorption surface 22 of the absorption element 21, converted into a desired energy form according to the absorption sensitivity, and output to the outside of the apparatus. It is assumed that there is no reflection or transmission of energy that reaches the absorption surface 22. Even if there is reflection or transmission, the influence is included in the directivity.

  The absorption center is the central point of absorption of energy E. The absorption vector is a vector that represents the reference direction of absorption of energy E. The absorption directional characteristic is a directional characteristic normalized with reference to the absorption center and the absorption vector. The angle that is the definition area of this directivity is called the incident angle.

  Note that the absorption vector coincides with the direction in which the absorption efficiency is maximized when viewed from the absorption center. Absorption efficiency is synonymous with absorption directivity. In the case of isotropic directivity, the normal line of the absorption surface may be an absorption vector.

[Energy propagation space]
Further, a space in which the radiation device 1 and the absorption device 2 are arranged and the energy E radiated from the radiation device 1 propagates to the absorption device 2 is assumed as an energy propagation space 3.

  In this energy propagation space 3, a three-dimensional orthogonal coordinate system is introduced, and position information such as arrangement and posture of each device and each element is expressed in the coordinate system.

  Further, the energy propagation space 3 is modeled by defining as follows.

(1) Energy E goes straight in the radial direction
(2) There is no loss on the energy propagation space 3 between radiation and absorption
(3) There is no propagation of energy E from other than the target radiation device 1

[2. Configuration of Embodiment]
[overall structure]
Next, the configuration of the present embodiment will be described with reference to the functional block diagram of FIG. That is, the energy analysis system A of this embodiment includes the energy analysis device 100 and the interface unit 200.

  The energy analysis apparatus 100 includes a coordinate system processing unit 110, a device model processing unit 120, an analysis unit 130, a storage unit 140, and the like. The interface unit 200 includes an input unit 210 and an output unit 220.

  The energy analysis apparatus 100 can be configured by a computer that can realize the functions of the above-described units by, for example, a predetermined program. Note that the functional block diagram shown in FIG. 4 is conceptual, and various specific circuits and software for realizing these functions are conceivable and are not limited to specific ones.

[Energy analyzer]
[Coordinate system processing section]
The coordinate system processing unit 110 is a processing unit that performs coordinate processing in the energy propagation space. The coordinate system processing unit 110 includes a propagation space model generation unit 111, a coordinate conversion unit 112, and the like. The propagation space model generation unit 111 is a processing unit that models the energy propagation space based on information input from the interface unit 200. The generated propagation space model is stored in the storage unit 140 described later.

  The coordinate conversion unit 112 is a processing unit that performs coordinate conversion between two designated coordinate systems. For example, in accordance with a request from the analysis unit 130 described later, coordinate conversion can be performed, and a central projection image can be generated.

[Device model processing section]
The device model processing unit 120 is a processing unit that generates a device model. The device model processing unit 120 includes a radiation device model generation unit 121, an absorption device model generation unit 122, and the like.

  The radiating device model generation unit 121 is a processing unit that models the radiating device 1 based on information such as the specifications of the radiating device 1 input from the interface unit 200. Information necessary for modeling includes, for example, all specifications necessary for analysis such as the position and orientation of the radiation device 1 in the energy propagation space, radiation intensity, and radiation directivity characteristics, as described above. The generated radiation device model is stored in the storage unit 140 described later.

  The absorption device model generation unit 122 is a processing unit that models the absorption device 2 based on information such as the specifications of the absorption device 2 input from the interface unit 200. Information necessary for modeling includes, for example, all specifications necessary for the analysis, such as the position and orientation of the absorption device 2 in the energy propagation space, absorption sensitivity, and absorption directivity characteristics, as described above. The generated absorption device model is stored in the storage unit 140 described later.

[Analysis Department]
The analysis unit 130 is a processing unit that performs energy analysis based on the propagation space model, the radiation device model, and the absorption device model. The analysis unit 130 includes an analysis model generation unit 131, a division unit 132, a calculation unit 133, an abnormality determination unit 134, an analysis result output unit 135, and the like.

  The analysis model generation unit 131 is a processing unit that generates an analysis model in response to an input from the interface unit 200. The analysis model is a model according to an analysis principle described later, and includes information as described above. The dividing unit 132 is a processing unit that divides the absorption surface of the absorption element 21 into triangles for analysis using an analysis model. The calculation unit 133 is a processing unit that calculates an analysis result using an analysis model.

  The abnormality determination unit 134 is a processing unit that interrupts the analysis when an analysis abnormality occurs due to incompleteness of various models. The analysis result output unit 135 is a processing unit that outputs an analysis result. The analysis results include normal analysis results or abnormal contents. The analysis result is stored in the storage unit 140.

  The dividing unit 132 and the calculating unit 133 include a plurality of processing units as illustrated in FIG. Details of the processing of these processing units will be described in the operation described later.

[Storage unit]
The storage unit 140 is a storage unit that stores various types of information (including processing results) necessary for the present embodiment. Each information storage area can be regarded as a storage unit for each information.

  As the storage unit 140, for example, any storage medium that can be used at present or in the future, such as a memory, a hard disk, or an optical disk, can be used. A mode in which the storage contents can be used for various processes by mounting a storage medium in which information has already been stored in a reading device may be employed.

  Furthermore, the storage unit 140 includes a register, a memory, and the like used as a temporary storage area. Therefore, queues, stacks, and the like can also be realized using the storage unit 140.

[Interface part]
The interface unit 200 is an interface that inputs and outputs information of the energy analysis apparatus 100. The interface unit 200 includes an input unit 210, an output unit 220, and the like.

  The input unit 210 is a device that inputs information, instructions, and the like necessary for the processing of this embodiment. For example, any device that can be used at present or in the future, such as a mouse, a keyboard, a remote control, a switch, and a display (touch panel) can be used. The input unit 210 also includes an interface that accepts input of information from an external device via a communication network.

  The output unit 220 is a device that outputs information necessary for the processing of this embodiment. For example, any device that can be used at present or in the future, such as a display or a printer, can be used. The output unit 220 includes an interface for outputting information to an external device via a communication network.

[3. Operation of the embodiment]
The operation of the present embodiment having the above configuration will be described.
[Energy analysis principle]
First, as a premise, the principle of energy analysis of this embodiment will be described.

[Power of energy]
Consider a power absorbed by an absorbing surface 22 when a part of the radiant flux of energy E radiated from a single radiating element 11 is irradiated to a certain absorbing element 21. At this time, the radiation of the energy E from the radiating element 11 is steady and isotropically diverges. The radiation intensity is assumed to be known. Further, the absorption of energy at the absorption surface 22 is isotropic with no directivity, and all irradiated radiant fluxes are absorbed.

  The radiant intensity I is the radiant flux dΦ [W] per minute solid angle dΩ [sr] as shown in the equation (1), and represents the density of the radiant flux. Therefore, in the case of uniform divergence, the radiant flux at that time is determined if the passing solid angle is determined regardless of the size of the solid angle and the radiation direction.

  That is, if the apparent solid angle of the absorption surface viewed from the radiation element 11 is obtained, the power absorbed by the absorption surface 22 can be obtained based on the equation (1).

[Solid angle]
Therefore, a method for calculating the solid angle will be described. First, consider a polygon 51 having N vertices (N = 13 in FIG. 5) on a spherical surface 50 as shown in FIG. The solid angle Ω which is the area of the polygon 51 has apexes n = {0,. . The interior angle of the N-1} when the _{a n,} can be calculated as Equation (2).

  The inner angle of the vertex n will be described with reference to FIG. That is, consider two great circles 52 and 53 (sphere cross-sections) that intersect at a vertex n and respectively include vertices (n−1, n + 1) on both sides adjacent to n. The interior angle of the vertex n is an intersection angle observed from the three vertices (n, n−1, n + 1) formed by the two great circles 52 and 53.

  This polygon 51 is called a spherical N-gon or simply a spherical polygon. In the case of a polygon having a plurality of holes inside instead of a simple polygon, the solid angles (area) of the plurality of holes inside may be subtracted from the solid angle (area) of the simple polygon.

[Introduction of coordinate system]
Next, a coordinate system introduced into the propagation space model will be described with reference to FIG.

That is, consider a three-dimensional orthogonal coordinate system Σ _C with the radiation center of the radiation element 11 as the origin. When three orthogonal components are represented by x, y, and z, the positive direction of the z axis coincides with the radiation vector. An arbitrary position vector on sigma ^_C, represented by ^C p.

With sigma _C, also introduces a spherical coordinate system sigma _S to the radial center as the origin. 3 component of sigma _S is one radius vector r and two argument theta, expressed by phi. The declination θ is an angle away from the radiation vector. Deflection angle φ is the separation angle from the positive x-axis direction on the xy plane in sigma _C.

Note that the right screw around sigma _C positive. An arbitrary position vector on sigma _S represented by ^{S p.} A position vector on a spherical surface having an arbitrary radius centered on the origin is represented by (θ, φ) ^T using the two declination angles.

[Apparent solid angle]
It is assumed that the target absorption surface 22 is a flat surface and is a simple polygonal closed region. If all of the vertices can be observed without overlapping when the absorbing surface 22 is viewed from the origin, a spherical surface having the same number of vertices is formed on the unit spherical surface 50 centered on the origin as shown in FIG. Center projected as the upper polygon 51.

  The area of the polygon 51 which is the central projection image of the absorption surface 22 is an apparent solid angle. And this area can be calculated | required by said Formula (2).

[Inner angle of spherical polygon]
The interior angle of each vertex of the spherical polygon 51 can be obtained using the vertex of the absorption surface 22 as shown in FIG. First, the position vector of sigma _C vertex n on the absorption surface 22 and ^C v _n. Then, consider a vertex ^C v _n corresponding to the inner angle a _n to be obtained, and two vertexes ^C v _n−1 and ^C v _{n + 1} adjacent thereto.

These vertex identifiers take congruent values using the number of vertices of the target polygon as a modulus. Further, the arrangement of vertices identifier watches surface of the absorbent surface 22 to impart a clockwise ^{_{... C v n-1, C}} v n, C v n + 1 ... and.

Normal vector ^is the same direction as C _{v ^n,} the plane in which the C _{v n} rests, and _{S n.} At this time, the intersection of the straight line passing through the origin ^C O and the vertex ^C v _n−1 and the plane _Sn is defined as ^C s _n−1 . Similarly, _{let the} intersection point for the vertex ^C v _{n + 1 be} ^C s _{n + 1} .

From these definitions, it can be seen that ^C s _n−1 and ^C s _{n + 1} are respectively placed on two planes including a great circle forming an inner angle to be obtained. It can also be seen that ^C O and ^C v _n are the intersections of these two planes.

Furthermore, _Sn is a plane orthogonal to these two planes. ^Therefore, the C _{v ^n,} a straight line passing through the C _{s ^n-1,} and the straight line passing through the C _{v n} and ^C s _{n + 1} is the angular ^{_{∠ C s n-1 C v}} n C s n + 1 , which forms on _{S n} is , that is, the interior angle _{a n} seeking.

Thus, the interior angle _{a n} by using the following equation (3) to (6), can be expressed as Equation (7). Here, Formula (3) is a plane function S (). Equation (4) is a linear function L (). Expression (5) is an intersection function I () between a straight line and a plane. Expression (6) is an angle function A () between three points.

[Triangulation]
From the above, it can be said that energy analysis can be performed by accurately obtaining the solid angle involved from the specifications and position and orientation of the radiating element 11 and the absorbing element 21. However, the conditions are limited such that the energy from the radiating element 11 is point radiation or non-directional.

  Next, let us consider energy analysis under general conditions that are more relaxed, such as when surface radiation or each element has directional characteristics. However, in the case of surface radiation, it is assumed that the absorption element is separated to such an extent that the directivity characteristic does not depend on the radiation distance. In any case, it is assumed that directivity characteristics are defined.

  Due to the directional characteristics, the radiation intensity and the absorption efficiency change depending on the propagation direction of energy. However, in the minute region, the difference may be considered as an error that can be ignored for the accuracy of energy analysis. These errors are called intensity tolerance and absorption tolerance, respectively, and are defined as energy analysis specifications. In addition, the storage unit 140 stores the intensity allowable error and the absorption allowable error that are input in advance from the interface unit 200.

  First, the inside of the target area is finely divided so as to satisfy these tolerances. Then, energy analysis is performed for each of the divided partial regions. Finally, the desired energy analysis can be achieved by integrating the analysis results.

  This division is performed by the dividing unit 132 dividing the absorbing surface 22 into triangles as shown in FIG. The triangulated area is called an absorption surface element 23.

[Calculation of power]
Due to the triangulation, the spherical polygon 51 center-projected on the unit sphere shown in FIG. 7 is also divided into spherical triangles. For this reason, if the solid angle of each divided triangular area is obtained and the solid angles are integrated, the solid angle of the original spherical polygon can be obtained. Therefore, if the apparent solid angle (area of the spherical polygon) of the absorption surface element 23 is integrated, the solid angle of the absorption surface 22 can be obtained.

In an energy analysis system, it is assumed that the absorption surface element k has an apparent solid angle ω _k and the radiation intensity I (ξ _k ) and the absorption efficiency η (ζ _k ) are associated with each other. Then, the calculation part 133 can obtain | require the power P which the whole absorption surface used as a target absorbs as Formula (8) from the principle demonstrated above. Note that ξ _k represents the radiation angle with respect to the absorbing surface element k, and ζ _k represents the incident angle with respect to the same surface element.

[Radiation intensity]
So far, calculations have been described on the assumption that the radiation intensity is known. However, depending on the radiating element 11, only the total radiant flux and directivity may be shown. At this time, the maximum radiant intensity Imax can be obtained by dividing the total radiant flux Φ by the value obtained by dividing the directional characteristic by the area of the unit spherical Ω.

2 declination in the directional characteristics sigma _S theta, when a real function D determined by phi (theta, phi), the maximum radiation intensity Imax can be expressed as Equation (9).

  At this time, the radiation intensity I (θ, φ) can be obtained as shown in Expression (10).

Further, when the directivity is a real function D (θ) that is determined only by the declination θ and can be differentiated M times, if the maximum integer less than or equal to the real number x is [x], the maximum radiant intensity Imax is expressed by the equation (11). ). Here, D ^(k) (θ) represents k-th derivative of D (θ), D ⁽⁰⁾ (θ) = D (θ), and s and t are ranges of θ in which directivity characteristics are defined. And s <t.

  At this time, the radiation intensity I (θ) can be obtained as shown in Expression (12).

  Such calculation of the radiation intensity can be performed by the radiation intensity calculation unit 121a in the radiation device model generation unit 121. The calculated radiation intensity is stored in the storage unit 140 as part of the radiation device model.

[4. Energy analysis processing procedure]
The procedure of energy analysis processing by the energy analysis apparatus 100 of the present embodiment based on the above principle will be described with reference to the block diagram of FIG. 4 and the flowcharts of FIGS. 9 to 11.

[4-1. Overall processing]
First, the outline | summary of the operation | movement procedure of the energy analyzer 100 is demonstrated with reference to the flowchart of FIG. The operation of the energy analysis apparatus 100 mainly includes procedures such as initialization, model construction, energy analysis, and analysis result output.

  It is assumed that all specifications used for energy analysis of the radiation device 1 and the absorption device 2 in the energy propagation space 3 are defined as preconditions for the operation of this device. These specifications are input to the apparatus via the interface unit 200 and stored in the storage unit 140.

  In any operation, the energy analysis apparatus 100 accepts a stop event from the input unit 210 of the interface unit 200. When a stop event is input, the energy analysis apparatus 100 stops the operation at that time, and waits for an event that can be accepted at the time of stop, such as an analysis start event or an internal state output event.

[Initialize]
First, a user or an external device inputs an analysis start event via the interface unit 200 (step 01). In response to this, the energy analysis apparatus 100 resets the previous state and performs initialization for executing a new model construction and energy analysis (step 02).

[Model construction]
After the initialization is completed, the energy analysis apparatus 100 waits for data input of the energy propagation space model, the radiation apparatus model, the absorption apparatus model, and the analysis model from the interface unit 200 (step 03 NO).

  When data of various models is input from the interface unit 200 (step 03 YES), the propagation space model generation unit 111, the radiation device model generation unit 121, the absorption device model generation unit 122, and the analysis model generation unit 131 are Model construction is performed (step 04). Each constructed model is stored in the storage unit 140 (step 05).

  As described in the above analysis principle, the analysis model mainly includes an analysis target, an analysis procedure, an allowable error, and the like. In addition, in order to make into analysis object the energy E which the absorption device 2 absorbs among the energy radiated | emitted from the radiation device 1, and the directivity characteristic of the radiation device 1 whole, the procedure which rotates the radiation device 1 can also be included. .

  Models other than the analysis model mainly include geometric information such as position and orientation and specifications related to radiation absorption such as radiation intensity and absorption efficiency. Examples thereof are as described above.

[Energy analysis]
The division unit 132 and the calculation unit 133 in the analysis unit 130 perform energy analysis using each model in response to model construction completion or a model construction completion event input from the interface (step 06). This analysis follows the aforementioned energy analysis principle. The procedure will be described later. When the analysis ends normally (NO at Step 07), the storage unit 140 stores the analysis result (Step 09).

  If the abnormality determination unit 134 determines that an analysis abnormality has occurred due to incompleteness of various models (YES in step 07), the abnormality determination unit 134 interrupts the analysis (step 08). The storage unit 140 stores the abnormality content as an analysis result (step 09).

[Analysis result output]
The analysis result output unit 135 outputs an analysis result including an abnormal content or an analysis completion event to the interface unit 200 when the analysis is normally completed or interrupted (step 10). The output unit 220 displays an analysis result in response to a request input from the input unit 210.

[4-2. Details of energy analysis]
[Triangulation]
First, the procedure in which the dividing unit 132 divides the absorbing surface 22 into triangular absorbing surface elements 23 will be described with reference to the flowchart of FIG. That is, the analysis unit 130 resets the divided triangle stack at the start of processing (step 11). And the polygon production | generation part 132a divides | segments the absorption surface 22 into a convex polygon (step 12). This convex polygon is called a partial absorption surface.

  The triangle generation unit 132b divides each partial absorption surface into triangles (step 13). This division generates a triangle by connecting a line segment from an arbitrary vertex of each partial absorption surface to all the vertices remaining on the same surface. Then, the triangle generation unit 132b puts the generated triangle geometric information into the divided triangle stack (step 14).

  Next, the determination point extraction unit 132c extracts the top triangle in the stack (step 15), and extracts two points that are the end points of the longest side among the three vertices of the triangle (step 16). Based on the radiation device model and the propagation space model, the intensity difference determination unit 132d obtains the radiation angle and the radiation intensity at the two points viewed from the target radiation element 11, and determines the intensity difference between the two points (step 17). ).

  Then, the allowable intensity error determination unit 132e determines whether the intensity difference between the two points is within the range of the allowable intensity error stored in the storage unit 140 (step 18). When the intensity difference is within the allowable error range (YES in Step 19), the efficiency difference determination unit 132f determines the absorption efficiency difference between the two points (Step 20).

  That is, the efficiency difference determination unit 132f determines an absorption efficiency difference between two points after obtaining the incident angle and the absorption efficiency viewed from the target absorption element 21 based on the absorption device model and the propagation space model.

  The allowable absorption error determination unit 132g determines whether the absorption difference between the two points is within the range of the allowable intensity error stored in the storage unit 140 (step 21). If the absorption difference is within the allowable error range (YES in step 22), the triangle output unit 132i outputs the triangle as a division completion triangle. The storage unit 140 stores the division completion triangle (step 23).

  When the intensity difference is outside the allowable error range (step 19 NO) or when the absorption difference is outside the allowable error range (step 22 NO), the subdividing unit 132h generates a triangle obtained by subdividing the target triangle (step 24). . For example, the triangle is bisected with the line segment connecting the bisector of the longest side obtained at the time of determination point extraction and the remaining vertex as a boundary. The subdivision unit 132h puts both of the subdivided triangles on the stack (step 14).

  Then, the processes of steps 15 to 22, 24, and 14 are repeated until the target triangle is subdivided into triangles within the allowable error range. If a triangle still exists in the stack (NO at step 25), the processes of steps 15 to 24 and 14 are repeated. If the stack becomes empty (YES in step 25), the dividing unit 132 ends the dividing process.

  The above-described division method is an example. Other methods may be used as long as the triangulation can be performed to the extent that the defined tolerance can be satisfied. The allowable error may be determined based on at least one of the radiation intensity and the absorption efficiency.

[Energy calculation]
Next, the procedure by which the calculation unit 133 calculates power based on the triangles divided as described above will be described with reference to the flowchart of FIG. First, the spherical triangle calculation unit 133a calculates a centrally projected spherical triangle for each divided triangle with respect to the unit sphere (step 31).

  Then, the solid angle calculation unit 133b calculates solid angles for all spherical triangles (step 32). Furthermore, the integration unit 133c determines the apparent solid angle of the absorption surface 22 by integrating all solid angles (step 33). The power calculation unit 133d obtains the power absorbed by the absorption surface based on the solid angle, the radiation intensity, and the absorption efficiency (Step 34).

[5. Effects of the embodiment]
According to the present embodiment as described above, the energy that is radiated from the radiation device 1 and absorbed by the absorption device 2 is avoided by using the actual model and accurately using the mathematical model. Can be analyzed.

  That is, the energy transferred between the radiation device having the isotropic directivity and the absorption device can be accurately obtained. Furthermore, even in the case of other directivity characteristics, it can be obtained accurately within the allowable error range.

  Even if the radiation intensity of the radiating element 11 is undefined, if the total radiation flux and the directivity are defined, the radiation intensity can be obtained and the energy can be analyzed. Furthermore, the directivity characteristics of the entire radiation device 1 or the entire absorption device 2 can be obtained.

[6. Other Embodiments]
The present invention is not limited to the embodiment as described above. The following aspects are also included in the present invention.

  (1) The absorption surface is not limited to a flat surface, and may be a curved surface, a surface having other irregularities, or the like.

  (2) In the tolerance determination, when compared with the range (threshold value) stored as the tolerance, it is set to “greater than” or “less than” that includes the threshold, or “greater than” that does not include the threshold. "Or" smaller "is free.

  (3) The energy analysis device 100 and the interface unit 200 may be realized by a common computer, or may be realized by a plurality of computers connected by a communication network. For example, each unit in the energy analysis apparatus 100 may be configured to be distributed among a plurality of computers. The storage unit 140 may be realized by a distributed storage medium.

DESCRIPTION OF SYMBOLS 1 ... Radiation device 2 ... Absorption device 3 ... Energy propagation space 11 ... Radiation element 21 ... Absorption element 22 ... Absorption surface 23 ... Absorption surface element 50 ... Unit spherical surface 51 ... Spherical polygon 52, 53 ... Great circle 100 ... Energy analysis device DESCRIPTION OF SYMBOLS 110 ... Coordinate system process part 111 ... Propagation space model production | generation part 112 ... Coordinate conversion part 120 ... Apparatus model process part 121 ... Radiation apparatus model production | generation part 121a ... Radiation intensity calculation part 122 ... Absorption apparatus model production | generation part 130 ... Analysis part 131 ... Analysis model generation unit 132 ... division unit 132a ... polygon generation unit 132b ... triangle generation unit 132c ... determination point extraction unit 132d ... intensity difference determination unit 132e ... strength allowable error determination unit 132f ... efficiency difference determination unit 132g ... absorption allowable error determination Unit 132h ... Subdivision unit 132i ... Triangle output unit 133 ... Calculation unit 133a ... Spherical triangle calculation unit 133b ... Solid angle calculation 133c ... integrating unit 133d ... power calculation section 134 ... abnormality determination unit 135 ... analysis result output unit 140 ... storage unit 200 ... interface unit 210 ... input section 220 ... Output section

Claims (7)

A propagation space model storage unit for storing a propagation space model which is information obtained by modeling a propagation space of energy by a three-dimensional orthogonal coordinate system;
A radiation device model storage unit that stores a radiation device model which is information modeling a radiation device having a radiation element that radiates energy to the outside;
An absorption device model storage unit that stores an absorption device model, which is information modeling an absorption device having an absorption element that absorbs the energy;
Have
The radiation device model includes information on the position of the radiation element in the three-dimensional orthogonal coordinate system, a radiation center that is a central point of energy radiation, and energy radiation intensity,
The absorber model includes information on the position of the absorber, the absorption surface, and the absorption efficiency in the three-dimensional orthogonal coordinate system,
Based on information on the position of the absorption element and the absorption surface, a dividing unit that divides the absorption surface into triangles;
A spherical triangle calculation unit for obtaining a spherical triangle in which the triangle divided by the dividing unit is center-projected on a spherical surface centered on the radiation center based on the position of the radiation element and the radiation center;
A solid angle calculation unit for obtaining a solid angle of the spherical triangle;
An integration unit for integrating the solid angles;
Based on the integrated solid angle, the radiation intensity and the absorption efficiency, a power calculation unit for obtaining the power of energy absorbed by the absorption surface;
An energy analysis apparatus comprising: The dividing unit is
A triangle generation unit for generating a triangle obtained by dividing the absorption surface;
A determination point extraction unit that extracts both ends of the longest side of the triangle generated by the triangle generation unit as determination points;
An allowable error determination unit that determines whether at least one of the difference in radiation intensity and the difference in absorption efficiency between the two determination points extracted by the determination point extraction unit is within a preset allowable error range; and
A subdivision unit that generates a triangle obtained by subdividing the triangle when the permissible error determination unit determines that it is outside the allowable error range;
The energy analysis apparatus according to claim 1, comprising:   The energy analysis apparatus according to claim 1, further comprising a radiation intensity calculation unit that calculates the radiation intensity based on a total radiant flux and directivity characteristics of the radiation apparatus. A computer or electronic circuit
A propagation space model storage process for storing a propagation space model which is information obtained by modeling a propagation space of energy by a three-dimensional orthogonal coordinate system;
A radiation device model storing process for storing a radiation device model, which is information modeling a radiation device having a radiation element that radiates energy to the outside;
An absorption device model storage process for storing an absorption device model which is information modeling an absorption device having an absorption element that absorbs the energy;
Run
The radiation device model includes information on the position of the radiation element in the three-dimensional orthogonal coordinate system, a radiation center that is a central point of energy radiation, and energy radiation intensity,
The absorber model includes information on the position of the absorber, the absorption surface, and the absorption efficiency in the three-dimensional orthogonal coordinate system,
The computer or electronic circuit is
Based on information on the position of the absorption element and the absorption surface, a dividing process for dividing the absorption surface into triangles;
A spherical triangle calculation process for obtaining a spherical triangle in which the triangle divided by the dividing process is projected onto the spherical surface centered on the radiation center based on the position of the radiation element and the radiation center;
A solid angle calculation process for obtaining a solid angle of the spherical triangle;
An integration process for integrating the solid angles;
Based on the integrated solid angle, the radiation intensity and the absorption efficiency, a power calculation process for obtaining the power of energy absorbed by the absorption surface;
The energy analysis method characterized by performing. On the computer,
A propagation space model storage process for storing a propagation space model which is information obtained by modeling a propagation space of energy by a three-dimensional orthogonal coordinate system;
A radiation device model storing process for storing a radiation device model, which is information modeling a radiation device having a radiation element that radiates energy to the outside;
An absorption device model storage process for storing an absorption device model which is information modeling an absorption device having an absorption element that absorbs the energy;
And execute
The radiation device model includes information on the position of the radiation element in the three-dimensional orthogonal coordinate system, a radiation center that is a central point of energy radiation, and energy radiation intensity,
The absorber model includes information on the position of the absorber, the absorption surface, and the absorption efficiency in the three-dimensional orthogonal coordinate system,
In the computer,
Based on information on the position of the absorption element and the absorption surface, a dividing process for dividing the absorption surface into triangles;
A spherical triangle calculation process for obtaining a spherical triangle in which the triangle divided by the dividing process is projected onto the spherical surface centered on the radiation center based on the position of the radiation element and the radiation center;
A solid angle calculation process for obtaining a solid angle of the spherical triangle;
An integration process for integrating the solid angles;
Based on the integrated solid angle, the radiation intensity and the absorption efficiency, a power calculation process for obtaining the power of energy absorbed by the absorption surface;
An energy analysis program characterized by causing The dividing process is
A triangle generation process for generating a triangle obtained by dividing the absorption surface;
Determination point extraction processing for extracting both ends of the longest side of the triangle generated by the triangle generation processing as determination points;
An allowable error determination process for determining whether a difference in radiation intensity or absorption efficiency between the two determination points extracted by the determination point extraction process is within a preset allowable error range; and
When it is determined that the allowable error determination process is outside the allowable error range, a subdivision process for generating a triangle obtained by subdividing the triangle;
The energy analysis program according to claim 5 , comprising: In the computer,
The energy analysis program according to claim 5 or 6 , wherein a radiation intensity calculation process for calculating the radiation intensity is executed based on a total radiant flux and directivity characteristics in the radiation device.

Images