JP7634472B2 - Building health assessment system - Google Patents

Description

The present invention relates to a soundness assessment system that diagnoses and evaluates the soundness of a building.

Various building integrity assessment systems have been proposed that can grasp the structural performance, i.e., integrity, of a building, such as the degree of damage to the building, without directly visually inspecting the building after an earthquake. It is desirable for such building integrity assessment systems to grasp the integrity of buildings with higher accuracy.
For example, Patent Document 1 discloses a technology configured to estimate the presence or absence of damage spreading in a judgment layer using a building damage spread detection model that learns the relationship between acceleration data in the observation layer and the presence or absence of damage spreading in the observation layer based on acceleration data obtained from a sensor installed in the observation layer of a multi-layer building and damage spread information indicating the presence or absence of damage spreading in the observation layer, and acceleration data acquired by a sensor installed in the judgment layer that is the subject of judgment.
Although the configuration disclosed in Patent Document 1 can estimate the presence or absence of damage spreading in the observation floor of the building, it cannot perform a detailed assessment of the soundness of each component constituting each floor of the building. Therefore, there is a limit to how much accuracy can be improved in assessing the soundness.

In response to this, Patent Document 2 discloses a configuration in which vibration sensors are installed at the joints of multiple structural members that form the structural frame of a structure, the joints and the multiple structural members that make up the joints are divided into substructures, and the detection information of the vibration sensors installed on each structural member joined to the joints is input and the vibration sensors at the joints are used as output, and the presence or absence and the extent of damage to the structural members that make up the substructures is detected based on the input/output relationship of the dynamic characteristics of each substructure.
In the configuration disclosed in Patent Document 2, detailed soundness can be determined for each structural member, but it is necessary to install a vibration sensor at each joint of the structural members. This requires the installation of a large number of vibration sensors, which makes the configuration complicated and therefore difficult to realize.

Patent Document 3 discloses a system that includes multiple sensors installed at multiple positions on a building, and measures the impact of an earthquake on the building at each position based on the impact of the earthquake on the building from before the arrival of the main motion to after the arrival of the main motion, which is measured by the multiple sensors. In this configuration, the amount of displacement and the story drift angle are calculated from the measurement results of the multiple sensors, and the soundness of the building is evaluated for each floor of the building and each structural element (member) of the building.
Although the configuration disclosed in Patent Document 3 can also perform detailed evaluation of the soundness of each component, it is necessary to install a sensor for each component in order to do so, which requires the installation of a large number of sensors, making the system configuration complicated.

JP 2020-8332 A JP 2015-4526 A JP 2020-143895 A

The problem that this invention aims to solve is to provide a building integrity assessment system that has a simple configuration and can accurately determine integrity.

In order to solve the above problems, the present invention employs the following means.
That is, the building soundness evaluation system of the present invention is a soundness evaluation system for diagnosing and evaluating the soundness of a building, and includes an analysis unit that performs a static incremental analysis in advance using a three-dimensional frame model configured to correspond to the building, including columns and beams, and calculates an analysis value of the energy absorption amount of each story and an analysis value of the energy absorption amount of each member including the columns and beams for each incremental step, which is a step of increasing the load in the analysis, and records them in an analysis result recording unit; an earthquake information recording unit that records earthquake information obtained from a sensor installed in the building; and an engine that calculates the story shear force and story displacement of each story based on the earthquake information and the weight of each story of the building, and calculates an estimated value of the energy absorption amount of each story based on the story shear force and the story displacement. an incremental step identification unit which selects, from the analysis values of the energy absorption of each layer recorded in the analysis result recording unit, a minimum difference analysis value of the energy absorption of each layer which has the smallest difference from the estimated value of the energy absorption of each layer, and identifies the corresponding incremental step as the incremental step with the smallest difference; and a soundness determination unit which obtains an estimate of the energy absorption of each component from the analysis value of the energy absorption of each component which corresponds to the incremental step with the smallest difference, compares the estimated value of the energy absorption of each component with a damage determination threshold set for each component to calculate the degree of damage, and determines the soundness of the building based on the degree of damage of each component.
According to this configuration, a static incremental analysis is performed using a three-dimensional frame model configured to correspond to a building including columns and beams, and the analysis value of the energy absorption amount of each story and the analysis value of the energy absorption amount of each member including columns and beams, which are calculated for each incremental step, which is a step of increasing the load in this static incremental analysis, are calculated by the analysis unit and recorded in the analysis result recording unit. When an earthquake occurs, the energy absorption amount estimation unit calculates an estimated value of the energy absorption amount of each story of the building based on earthquake information obtained from a sensor installed in the building and the weight of each story of the building. The soundness determination unit selects, from among the analysis values of the energy absorption amounts of each story, a minimum difference analysis value of the energy absorption amount of each story that has a minimum difference from the estimated value of the energy absorption amount of each story, and identifies the incremental step in which the minimum difference analysis value corresponding to the selected minimum difference analysis value of the energy absorption amount of each story is output as a minimum difference incremental step. In the minimum difference incremental step identified in this way, the load is increased so that each story absorbs energy equivalent to the energy absorbed by each story when the earthquake occurred for which earthquake information was obtained. Therefore, the analysis value of the energy absorption amount of each component at the incremental step with the smallest difference in the static incremental analysis is considered to be a value close to the amount of energy absorbed by each component when the earthquake for which earthquake information was obtained occurred. Therefore, by obtaining an estimate of the energy absorption amount of each component from the analysis value of the energy absorption amount of each component corresponding to the incremental step with the smallest difference, and calculating the degree of damage based on the estimated energy absorption amount of each component by comparing it with the damage judgment threshold, it is possible to evaluate the degree of damage to each component that constitutes the building and judge the soundness of the building.
In this way, by estimating the presence or absence of damage to the components that make up a building and the extent of that damage, damage assessment can be performed for each component, allowing for a highly accurate and reliable assessment of the soundness of the building.
Furthermore, in the above-described configuration, although damage assessment can be performed for each component, there is no need to provide a sensor for each component, and therefore a large number of sensors are not required.
Therefore, it is possible to provide a building soundness assessment system that has a simple configuration and is capable of accurately determining soundness.

In one aspect of the present invention, the analysis unit performs static incremental analysis on each of a story collapse model in which parameters are set for the three-dimensional frame model so that the building will undergo story collapse, and a global collapse model in which parameters are set for the three-dimensional frame model so that the building will undergo global collapse, and calculates, for each incremental step, an analytical value of the energy absorption amount of each story and an analytical value of the energy absorption amount of each member, and records them in the analysis result recording unit. The incremental step identification unit, for each combination of the analytical value of the energy absorption amount of each story in the story collapse model and the analytical value of the energy absorption amount of each story in the global collapse model, which are recorded in the analysis result recording unit, calculates, for each story, the analytical value of the energy absorption amount of each story in the story collapse model and the analytical value of the energy absorption amount of each member. the difference between the sum of the analytical value of the energy absorption amount of each layer in the layer collapse type model and the estimated value of the energy absorption amount of each layer is calculated, and a minimum difference combination which is a combination such that the sum of the differences over all layers is minimized is selected as a minimum difference analytical value of the energy absorption amount of each layer, and the incremental step in the story collapse type model and the incremental step in the total collapse type model which correspond to the analytical value of the energy absorption amount of each layer in the story collapse type model and the analytical value of the energy absorption amount of each layer in the total collapse type model in the minimum difference combination are identified as the incremental step with the first smallest difference and the incremental step with the second smallest difference, respectively, and the soundness assessment unit assesses the soundness of the building based on the incremental step with the first smallest difference and the incremental step with the second smallest difference.
Here, the "story collapse model with parameters set to cause story collapse" specifically refers to a three-dimensional frame model in which the strength of each beam and panel is set to be sufficiently large so that a story collapse mechanism is formed in which the columns yield first at the column-beam joints. Also, the "global collapse model with parameters set to cause total collapse" specifically refers to a three-dimensional frame model in which the strength of columns other than the first story column bases and the top story column capitals is set to be sufficiently large so that a total collapse mechanism is formed in which the beams or joint panels yield first at the column-beam joints.
According to this configuration, a static incremental analysis is performed on each of a story collapse model in which parameters are set so that a story collapse occurs for the three-dimensional frame model, and a global collapse model in which parameters are set so that a global collapse occurs for the three-dimensional frame model. For each combination of the analysis value of the energy absorption amount of each story in the story collapse model and the analysis value of the energy absorption amount of each story in the global collapse model recorded in the analysis result recording unit, the difference between the sum of the analysis value of the energy absorption amount of each story in the story collapse model and the analysis value of the energy absorption amount of each story in the global collapse model and the estimated value of the energy absorption amount of each story is calculated for each story, and a minimum difference combination that is a combination that minimizes the sum of these differences in all stories is selected as the minimum difference analysis value of the energy absorption amount of each story. Then, the incremental step in the story collapse model and the incremental step in the global collapse model corresponding to the minimum difference combination selected in this way are identified as the first minimum difference incremental step and the second minimum difference incremental step, respectively.
Here, by combining the states in which the story collapse and the total collapse have progressed to the increment step with the first minimum difference and the increment step with the second minimum difference identified as described above by the static increment analysis, the energy absorption amount of each story is close to the estimated energy absorption amount at the time of the occurrence of the earthquake whose earthquake information was obtained. In other words, the increment step with the first minimum difference and the increment step with the second minimum difference can be regarded as the progress degree of each of the story collapse and the total collapse in the earthquake whose earthquake information was obtained.
By determining the soundness of a building based on the incremental step with the smallest first difference and the incremental step with the smallest second difference identified in this manner, it is possible to estimate a collapse state that is closer to the damage caused to the building by the actual earthquake load, and to evaluate the soundness of the building with greater accuracy.

In one aspect of the present invention, the soundness assessment unit calculates, for each component, the sum of the analytical value of the energy absorption amount of each component in the layer collapse type model corresponding to the incremental step with the smallest first difference and the analytical value of the energy absorption amount of each component in the global collapse type model corresponding to the incremental step with the smallest second difference, and compares this as an estimate of the energy absorption amount of each component with the damage assessment threshold to calculate the degree of damage.
The analytical value of the energy absorption amount of each member in the story collapse model corresponding to the incremental step of the first minimum difference identified as above, and the analytical value of the energy absorption amount of each member in the total collapse model corresponding to the incremental step of the second minimum difference are obtained. As already described, the incremental step of the first minimum difference and the incremental step of the second minimum difference are the progress of each story collapse and total collapse in the earthquake whose earthquake information was acquired, so the sum of the analytical values of the energy absorption amount of each member corresponding to the incremental step of the first minimum difference and the incremental step of the second minimum difference obtained as above is considered to be a value close to the energy absorbed by each member when the earthquake whose earthquake information was acquired occurred. The soundness of the building can be evaluated with higher accuracy by comparing this sum with the damage judgment threshold value for each member as an estimated value of the energy absorption amount of each member.

The present invention makes it possible to accurately determine the health of a device using a simple configuration.

1 is a diagram showing a schematic configuration of a building soundness evaluation system in this embodiment. FIG. FIG. 1 is a diagram showing a schematic diagram of a layer-collapse type collapse mechanism. FIG. 13 is a diagram showing a schematic diagram of a total collapse type collapse mechanism. FIG. 13 is a diagram showing an example of a collapse mechanism during an actual earthquake, and is a schematic diagram showing a collapse mechanism that combines a story collapse type and a total collapse type. FIG. 2 is a diagram for explaining members that constitute a three-dimensional framework model. FIG. 13 is a diagram showing an example of static incremental analysis of a story-collapse type three-dimensional frame model. FIG. 13 is a diagram showing an example of static incremental analysis of a total collapse type three-dimensional frame model. This is a diagram showing how the amount of energy absorbed by each story of a building during an earthquake can be estimated by combining the amount of energy absorbed by the story-collapse model and the total-collapse model. FIG. 13 is a diagram showing the amount of energy absorbed by a member when the member undergoes one cycle of deformation. 2 is a flowchart showing the flow of a building health evaluation method in the building health evaluation system. 1 is a diagram showing an example of verification in which the energy absorption amount of each component constituting a building is estimated for model 1 using the building soundness evaluation system in this embodiment. FIG. 13 is a diagram showing an example of verification in which the energy absorption amount of each component constituting a building is estimated for model 2 using the building health evaluation system in this embodiment. FIG. 13 is a diagram showing an example of verification in which the energy absorption amount of each component constituting a building is estimated for model 3 using the building health assessment system in this embodiment. FIG. 13 is a diagram showing an example of verification in which the energy absorption amount of each component constituting a building is estimated for model 4 using the building health evaluation system in this embodiment. FIG. 13 is a diagram showing an example of verification in which the energy absorption amount of each component constituting a building is estimated for model 5 using the building health evaluation system in this embodiment. FIG. 13 is a diagram showing an example of verification in which the energy absorption amount of each component constituting a building is estimated for model 6 using the building health assessment system in this embodiment. FIG.

Hereinafter, an embodiment of a building soundness evaluation system 1 according to the present invention will be described with reference to the accompanying drawings.
FIG. 1 shows a schematic configuration of a building soundness evaluation system 1 according to this embodiment.
This invention is a building health monitoring system that targets buildings including columns and beams, estimates the amount of energy absorption in each story and each component of a building based on the results of static incremental analysis using a three-dimensional frame model composed of columns and beams, and earthquake information obtained from sensors installed in the building, and then compares the estimated energy absorption amount with a damage judgment threshold set for each component to calculate the degree of damage to each component, thereby diagnosing and evaluating the health of the building.
1, the building health assessment system 1 includes a building 10 and an assessment device 20. The building health assessment system 1 assesses the health of the building 10 after an earthquake occurs.
The building 10 is constructed on the ground and has a plurality of stories 11 in the vertical direction. In this embodiment, the building 10 is constructed as a so-called rigid frame structure. The number of stories and the structure (reinforced concrete structure, steel frame structure, reinforced concrete structure, etc.) of the building 10 are not limited in any way.
The building 10 is provided with sensors 12. In this embodiment, the sensors 12 are installed on each floor 11 of the building 10. When each sensor 12 detects shaking exceeding a preset threshold during an earthquake, it detects acceleration waveform data for each floor 11 of the building 10. Each sensor 12 has a memory (not shown) and a communication unit (not shown). Each sensor 12 stores the detected acceleration waveform data in the memory (not shown) as earthquake information. After the earthquake ends, the sensor 12 transfers the acceleration waveform data stored in the memory from the communication unit to the evaluation device 20 via the external network 100.
Here, the external network 100 is, for example, a public wireless network that can wirelessly communicate with the communication unit. The communication unit transmits earthquake information detected by the sensor 12 to the evaluation device 20 via the external network 100 as shown in FIG.

The evaluation device 20 is connected to an external network 100 via a wired or wireless connection. The evaluation device 20 evaluates the soundness of the building 10 based on earthquake information detected by the sensor 12.
Fig. 2 is a schematic diagram of a collapse mechanism of a story collapse type, Fig. 3 is a schematic diagram of a collapse mechanism of a total collapse type, and Fig. 4 is a schematic diagram of a collapse mechanism that combines a story collapse type and a total collapse type, showing an example of a collapse mechanism during an actual earthquake.
When the rigid frame that constitutes the building 10 receives a strong earthquake motion, the members are deformed and a collapse mechanism is formed. The collapse mechanism is roughly divided into a "story collapse type" in which the columns yield first at the column-beam joints as shown in FIG. 2, and a "total collapse type" in which the beams or the joints between the columns and beams yield first as shown in FIG. 3. The collapse mechanism is assumed at the time of design based on the strength ratio of the columns, beams, and the joints between the columns and beams, but since it is affected by the input direction and the variation in the member strength during an actual strong earthquake, the collapse mechanism does not necessarily form as expected. In addition, unless the member-to-member strength ratio is an extreme value, it is highly likely that the damage distribution will show a mixture of the "story collapse type" and the "total collapse type" during a strong earthquake, as shown in FIG. 4. Therefore, even if the story deformation angle of each story is known by monitoring with the sensor 12, it is considered that the members that are damaged and the degree of damage will differ depending on the collapse mechanism formed at that time.
Therefore, the evaluation device 20 estimates the progress of each of the story collapse type collapse mechanisms and the total collapse type collapse mechanisms that have occurred in the building 10 based on the earthquake information detected by the sensor 12, estimates the degree of damage that has occurred in each component of the building 10 based on the estimation results, and judges the soundness of the building based on the degree of damage to each component.
Here, the "story collapse model with parameters set to cause story collapse" specifically refers to a three-dimensional frame model in which the strength of each beam and panel is set to be sufficiently large so that a story collapse mechanism is formed in which the columns yield first at the column-beam joints. Also, the "global collapse model with parameters set to cause total collapse" specifically refers to a three-dimensional frame model in which the strength of columns other than the first story column bases and the top story column capitals is set to be sufficiently large so that a total collapse mechanism is formed in which the beams or joint panels yield first at the column-beam joints.

As shown in FIG. 1 , the evaluation device 20 mainly includes an analysis unit 21, an analysis result recording unit 27, an earthquake information recording unit 22, an energy absorption amount estimation unit 23, an incremental step identification unit 24, and a soundness determination unit 25.
The analysis unit 21 calculates an analytical value of the energy absorption amount of each story 11 of the building 10 and an analytical value of the energy absorption amount of each member constituting the building 10 by performing an analysis in advance using a three-dimensional frame model M configured to correspond to the building 10 as shown in Figures 2 and 3. The analysis result recording unit 27 is a so-called database, in which the analytical value of the energy absorption amount of each story 11 of the building 10 and the analytical value of the energy absorption amount of each member constituting the building 10 calculated by the analysis unit 21 are recorded.
FIG. 5 is a diagram for explaining members that constitute the three-dimensional framework model.
The analysis unit 21 uses a three-dimensional frame model M including columns 16 and beams 17 for the analysis, as shown in Fig. 5. In the three-dimensional frame model M, the columns 16 and beams 17 are replaced with wire rods as elements that can take into account deformation in the axial, shear, bending, and torsional directions, respectively. Since both ends of the columns 16 and beams 17 become plastic during an earthquake, the energy absorption amount is calculated for both ends 16a, 16b, 17a, and 17b of the columns 16 and beams 17. The joints between the columns 16 and beams 17 are treated as panels 19 composed of rigid beams 19a and diagonal members 19b formed by elastic-plastic spring elements, and the energy absorption amount is calculated.

The analysis unit 21 performs static incremental analysis using a story collapse model M1, in which parameters are set so that the building 10 undergoes story collapse with respect to the three-dimensional frame model M as shown in Fig. 2, and a global collapse model M2, in which parameters are set so that the building 10 undergoes global collapse with respect to the three-dimensional frame model M as shown in Fig. 3. The story collapse model M1 has parameters set so that the strength of beams and panels is sufficiently large so that story collapse progresses during static incremental analysis. The global collapse model M2 has parameters set so that the bending strength of the end of columns is sufficiently large except for the column bases of the lowest story and the column capitals of the highest story so that global collapse progresses during static incremental analysis.
6 and 7 are diagrams showing examples of static incremental analysis of a story-collapse type space frame model, respectively. FIG. 7 is a diagram showing examples of static incremental analysis of a whole-collapse type space frame model.
In the static incremental analysis, the load, i.e., seismic force, acting on the three-dimensional frame models M1 and M2 of the building 10 is gradually increased, and the behavior of the building is analyzed at each incremental step, which is the step of increasing the load. In each of Figures 6 and 7, the left side shows the i1th incremental step, which is a state in which the energy absorption amount of each story and each member is small shortly after the static incremental analysis starts. The right side shows the i2th incremental step, which is a state in which the load is further applied from the i1th incremental step, the energy absorption amount of each story and each member increases, and story collapse or total collapse has progressed.
In the static push-back analysis, the Ai distribution, which is the external force distribution, is determined from the height of the building 10, the weight of each story 11, and the primary design natural period. In the static push-back analysis, calculations are performed so that the external force is gradually increased while maintaining the Ai distribution.

Next, based on the results of the static incremental analysis, the analysis unit 21 calculates, for each incremental step in the static incremental analysis, the analytical values of the energy absorption amount of each story and the analytical values of the energy absorption amount of each component, including the columns 16 and beams 17.
For this purpose, first, based on the results of the static incremental analysis of the story collapse model M1, the analysis unit 21 calculates the energy absorption amount of each story and the energy absorption amount of each component (columns 16, beams 17, panels 19) for each incremental step, which is the stage at which the load is increased in the static incremental analysis, with the X direction, which is a direction located in a horizontal plane, and the Y direction, which is a direction perpendicular to the X direction in the horizontal plane, as the load directions, and records these in the analysis result recording unit 27 as the analysis values of the energy absorption amount of each story and the analysis values of the energy absorption amount of each component.
In practice, the analysis unit 21 calculates the following values.
_sEwcX (i, j): analytical value of the energy absorption amount of the jth story from the bottom at the i-th incremental step when loaded in the X direction in the static incremental analysis result of the story collapse model M1. _sEwcY (i, j): analytical value of the energy absorption amount of the jth story from the bottom at the i-th incremental step when loaded in the Y direction in the static incremental analysis result of the story collapse model M1. _eEwcX (i, k): analytical value of the energy absorption amount of the kth member at the i-th incremental step when loaded in the X direction in the static incremental analysis result of the story collapse model M1. _eEwcY (i, k): analytical value of the energy absorption amount of the kth member at the i-th incremental step when loaded in the Y direction in the static incremental analysis result of the story collapse model M1. As shown in Figure 6, in the story collapse model M1, when the applied load in the load direction (X direction, Y direction) along the horizontal plane is increased step by step through each incremental step, the bottom story finally collapses, and the energy absorption amount in the bottom story is significantly larger than that of the upper stories. Also, the collapse in the bottom story is mainly due to damage to the columns, and the energy absorption amount in the columns of the bottom story is significantly larger than that of the beams, panels, etc. other than the columns.

In addition, based on the results of the static incremental analysis of the total-collapse model M2, the analysis unit 21 calculates the energy absorption amount of each story and the energy absorption amount of each component (columns 16, beams 17, panels 19) for each incremental step in each of the X and Y directions, and records these in the analysis result recording unit 27 as the analysis values of the energy absorption amount of each story and the analysis values of the energy absorption amount of each component.
In practice, the analysis unit 21 calculates the following values.
sE _osX (i, j): analytical value of the energy absorption amount of the j-th layer from the bottom at the i-th incremental step when loaded in the X direction in the static incremental analysis result of the globally collapse type model M2. sE _osY (i, j): analytical value of the energy absorption amount of the j-th layer from the bottom at the i-th incremental step when loaded in the Y direction in the static incremental analysis result of the globally collapse type model M2. eE _osX (i, k): analytical value of the energy absorption amount of the k-th member at the i-th incremental step when loaded in the X direction in the static incremental analysis result of the globally collapse type model M2. eE _osY (i, k): analytical value of the energy absorption amount of the k-th member at the i-th incremental step when loaded in the Y direction in the static incremental analysis result of the globally collapse type model M2. As shown in Fig. 7, in the total collapse model M2, when the applied load in the load direction (X direction, Y direction) along the plane is gradually increased, the energy absorption amount increases in all stories of the building 10. Similarly, the energy absorption amount of each component increases. In each story, the energy absorption amount increases not only for columns but also for beams, panels, etc.

When an earthquake occurs, the earthquake information recording unit 22 records earthquake information transmitted from the communication unit 14 of the sensor 12 in the building 10 via the external network 100. The earthquake information recording unit 22 records acceleration waveform data in each story 11 of the building 10 detected by the sensor 12 when an earthquake occurs.
The energy absorption estimation unit 23 calculates the story shear force and inter-story displacement of each story 11 based on the earthquake information detected by the sensor 12 and the weight of each story 11 of the building 10, and calculates an estimated energy absorption amount of each story based on the story shear force and inter-story displacement.

と表される。なお、ＦＸ（ｊ）、ＦＹ（ｊ）、ＦＺ（ｊ）は、Ｍａｓｓ（ｊ）をｊ番目の層の重量、ａｂＡＣＣＸ（ｊ）、ａｂＡＣＣＹ（ｊ）、ａｂＡＣＣＺ（ｊ）を、それぞれ、ｊ番目の層におけるＸ方向、Ｙ方向、Ｚ方向の観測加速度とすると、次のように表される。
ＦＸ（ｊ）＝－Ｍａｓｓ（ｊ）×ａｂＡＣＣＸ（ｊ）
ＦＹ（ｊ）＝－Ｍａｓｓ（ｊ）×ａｂＡＣＣＹ（ｊ）
ＦＺ（ｊ）＝－Ｍａｓｓ（ｊ）×ａｂＡＣＣＺ（ｊ）
層間変位は、各層１１における加速度波形の２階積分を上下層で差分することにより算出される。 First, the energy absorption estimation unit 23 calculates the story shear force in each story 11 based on the acceleration waveform data of each story 11 obtained by the sensor 12 during an earthquake and the weight of each story 11 .
More specifically, let QX(j), QY(j), and QZ(j) be the story shear forces in the X direction, Y direction, and Z direction perpendicular to the X and Y directions, respectively, on the jth story, and FX(j), FY(j), and FZ(j) be the inertial forces in the X direction, Y direction, and Z direction, respectively, on the jth story. At this time, on the top floor, QX(j), QY(j), and QZ(j) are respectively:
QX(j)=FX(j),
QY(j)=FY(j),
QZ(j)=FZ(j),
In addition, except for the top floor, the story shear force acting on the jth story from the bottom is the sum of the inertial forces in the upper stories, so if the total number of stories is T, then

Note that FX(j), FY(j), and FZ(j) can be expressed as follows, where Mass(j) is the weight of the jth layer, and abACCX(j), abACCY(j), and abACCZ(j) are the observed accelerations in the X, Y, and Z directions, respectively, in the jth layer.
FX(j)=-Mass(j)×abACCX(j)
FY(j)=-Mass(j)×abACCY(j)
FZ(j)=-Mass(j)×abACCZ(j)
The inter-story displacement is calculated by taking the difference between the second integral of the acceleration waveform in each story 11 and the upper and lower stories.

Next, the energy absorption estimation unit 23 calculates an estimated value of the energy absorption amount in each story 11 for each of the two horizontal directions (X direction and Y direction) based on the story shear force and inter-story displacement in each story 11 calculated as described above.
In practice, the energy absorption amount estimating unit 23 calculates the following value.
sE _obsX (j): an estimate of the amount of energy absorbed by the jth layer from the bottom in the X direction, calculated based on earthquake information. sE _obsY (j): an estimate of the amount of energy absorbed by the jth layer from the bottom in the Y direction, calculated based on earthquake information.
Specifically, the energy absorption estimation unit 23 calculates the estimated values sE obsX (j) and sE obsY (j) of the energy absorption in each story 11 based on the story shear force in each story 11 calculated from the earthquake _information and the history area in relation to the inter-story _displacement .

The incremental step identifying unit 24 selects the minimum difference analytical value of the energy absorption amount of each layer which has the smallest difference from the estimated values sE _obsX (j), sE _obsY (j) of the energy absorption amount of each layer from the analytical values of the energy absorption amount of each layer 11 recorded in the analysis result recording unit 27, and identifies the corresponding incremental step as the incremental step with the smallest difference. The incremental step identifying unit 24 will be described in detail below.
Let _nwcX be the total number of incremental steps for increasing the load when loading in the X direction in the static incremental analysis result of the story collapse type model M1, _nwcY be the total number of incremental steps for increasing the load when loading in the Y direction in the static incremental analysis result of the story collapse type model M1, _nosX be the total number of incremental steps for increasing the load when loading in the X direction in the static incremental analysis result of the global collapse type model M2, and _nosY be the total number of incremental steps for increasing the load when loading in the Y direction in the static incremental analysis result of the global collapse type model M2.

The total number of incremental steps when loading in the X-direction in the static incremental analysis result of the story-collapse type model M1 is _nwcX , and the total number of incremental steps when loading in the X-direction in the static incremental analysis result of the global-collapse type model M2 is _nosX . Therefore, there are nwcX ×nosX combinations of the analytical values of the energy absorption amount of each story when loading in the X-direction in the static incremental analysis result of the story-collapse type model M1 and the analytical values of the energy absorption amount of each story when loading in the X-direction in the static incremental analysis result of the _global- collapse type model _M2 .
The incremental step identification unit 24 calculates the layer energy difference e _x at the time of loading in the X direction for all of these n _wcX ×n _osX combinations by the following formula (1), selects the minimum difference combination which is the combination that makes the layer energy difference e _x at the time of loading in the X direction the smallest, as the minimum difference analytical value of the energy absorption amount of each layer, and identifies a combination of the incremental step i wcX (1≦i wcX ≦n wcX ) at the time of loading in the X direction in the static incremental analysis result of the story collapse type model M1 and the incremental step i _osX (1≦i _osX ≦n _osX ) at the time of loading in the X direction in the static incremental analysis result of the global collapse type model M2, which corresponds to the analytical value of the energy absorption amount of each layer in the story collapse type model M1 and the analytical value of the energy absorption amount of _{each layer} in the global _collapse type _{model M2, among the minimum} difference combinations.

Similarly, there are nwcY ×nosY combinations of the analytical values of the energy absorption amount of each layer when loaded in the Y direction in the static incremental analysis results of the story collapse type model M1 and the analytical values of the energy absorption amount of each layer when loaded in the _Y direction in the static incremental analysis results of the global collapse type model _M2 .
The incremental step identification unit 24 calculates the layer energy difference _eY when loaded in the Y direction for all of these _nwcY ×n _osY combinations using the following equation (2), selects the minimum difference combination that results in the smallest layer energy difference _eY when loaded in the Y direction as the minimum difference analytical value of the energy absorption amount of each layer, and identifies, among this minimum difference combination, a combination of the incremental step iwcY (1≦iwcY ≦nwcY ) when loaded in the Y direction in the static incremental analysis result of the story collapse type model M1 and the incremental step i _osY (1≦i _osY ≦n osY ) when loaded in the Y direction in the static incremental analysis result of the global collapse type model M2, which corresponds to the analytical value of the energy absorption amount of each layer in the story collapse type model _M1 and the analytical value of the energy absorption amount of each _layer in the global _collapse type model _M2 .

The operation of the incremental step identification unit 24 and the soundness assessment unit 25 described later is based on the premise that the amount of energy absorbed by each story (and each component) when the building 10 is subjected to earthquake motion after damage has progressed to a certain extent is the same as the amount of energy absorbed by each story (and each component) when the building 10 is subjected to the same earthquake motion in its newly constructed state.
For example, consider a case where the entire building 10 collapses, and then a story collapse progresses. In such a case, the amount of energy Ea that the building 10 ultimately absorbs can be considered to be the sum E1 + E2 of the amount of energy E1 that the building 10 absorbs from when it is newly constructed until the point at which the entire collapse has progressed and the amount of energy E2 that the building 10 absorbs from the point at which the entire collapse has progressed until the point at which the story collapse has progressed.
According to the above assumptions, the amount of energy E2 absorbed from the time when the overall collapse has progressed to the time when the story collapse has progressed is equal to the amount of energy E2' absorbed from when the building 10 was newly constructed to the time when the story collapse has progressed (without going through overall collapse). In other words, the above sum E1+E2 is equal to the sum E1+E2' of the amount of energy E1 absorbed from when the building 10 was newly constructed to the time when the overall collapse has progressed and the amount of energy E2' absorbed from when the building 10 was newly constructed to the time when the story collapse has progressed (without going through overall collapse).
Therefore, based on the above assumptions, from the perspective of the amount of energy absorption of each story and each component, even if story collapse and total collapse occur in building 10 one after the other, or simultaneously, this can be regarded as the sum of the states in which story collapse and total collapse have occurred individually in building 10, each of which has progressed to a certain degree.

FIG. 8 is a diagram conceptually showing how the energy absorption amount in each story of a building during an earthquake is estimated by combining the energy absorption amount in the story collapse model and the energy absorption amount in the total collapse model.
Based on the above considerations, in this embodiment, it is assumed that when an earthquake actually occurs, the earthquake causes story collapse and total collapse of the building 10, and an investigation is conducted into whether the estimated energy absorption amounts sE _obsX (j), sE _obsY (j) of each story calculated from earthquake information can be expressed as the sum of the analytical values sE _wcX (i wcX , j), sE wcY (i _wcY , j) of the energy absorption amounts of each story at any stage (i.e., incremental step i _wcX , i _wcY ) of the static incremental analysis using the story collapse model M1 recorded in the analysis result recording unit 27, and the analytical values sE _osX (i _osX , j), sE _osY (i _osY , j) of the energy absorption amounts of each story at any stage (i.e., incremental _step i _osX , i _osY ) of the static incremental analysis using the total collapse _model M2.
If the estimated values sE _obsX (j) and sE _obsY (j) of the energy absorption amount of each story calculated from the earthquake information coincide with or are close to the sum of the analytical values sE _wcX (i _wcX , j) and sE wcY (i _wcY , j) of the energy absorption amount of each story at some incremental step i wcX , i _wcY in the static incremental analysis using the story collapse model M1 recorded in the analysis result recording unit 27 and the analytical values sE _osX (i _osX , j) and _{sE osY (} i _osY , j) of the energy absorption amount of each story at some incremental step i _osX , i _osY in the static incremental analysis using the global collapse model M2, then it is judged that in the static incremental analysis using the story collapse model M1 _, story collapse has occurred at incremental steps i _wcX , i osY in the earthquake. This can be considered to be a combination of a state in which the collapse has progressed by _wcY and a state in which the collapse has progressed by incremental steps i _osX and i _osY in the static incremental analysis using the global collapse model M2.

Here, the story energy difference e _x at the time of loading in the X direction expressed by the above formula (1) is the sum of the differences between the analytical value sE _wcX (i _wcX , j) of the energy absorption amount at a certain incremental step i _wcX at the time of loading in the X direction in the static incremental analysis result of the story-collapse type model M1 and the analytical value sE _osX (i _osX , j) of the energy absorption amount at a certain incremental step i _osX at the time of loading in the X direction in the static incremental analysis result of the global-collapse type model M2 and the energy absorption amount sE _obsX (j) in the X direction calculated from the earthquake information for each story.
Therefore, by selecting a minimum difference combination, which is a combination of the analytical value _sEwcX ( _iwcX , _j ) of the energy absorption amount when loaded in the X direction in the static incremental analysis result of the story collapse type model M1 and the analytical value _sEosX (iosX, j) of the energy absorption amount when loaded in the X direction in the static incremental analysis result of the global collapse type model M2, which minimizes the story energy difference e x when loaded in the X direction in equation (1), as the minimum analytical value of the energy absorption amount in the X direction, the incremental step _iwcX and incremental step iosX corresponding to the analytical value _sEwcX ( _iwcX , j) of the energy absorption amount of each story in the story collapse type model M1 and the analytical value _sEosX ( _iosX , j) of the energy absorption amount of each story in the global collapse type model M2 in this minimum difference combination are set as the incremental step _{iwcX with} the first minimum difference and the incremental step _iosX with the second minimum difference, respectively. It can be identified as _osX .
Similarly, the layer energy difference _eY during loading in the Y direction expressed as equation (2) as above is the sum of the differences between the analytical value _sEwcY ( _iwcY , j) of the energy absorption amount at a certain incremental step _iwcY during loading in the Y direction in the static incremental analysis results of the layer-collapse type model M1 and the analytical value _sEosY ( _iosY , j) of the energy absorption amount at a certain incremental step _iosY during loading in the Y direction in the static incremental analysis results of the total-collapse type model M2, and the energy absorption amount _sEobsY (j) in the Y direction calculated from the earthquake information for each layer.
Therefore, by selecting a minimum difference combination, which is a combination of the analytical value _sEwcY ( _iwcY , j ₎ of the energy absorption amount when loaded in the Y direction in the result of the static incremental analysis of the story-collapse type model M1 and the analytical value _sEosY (iosY, j) of the energy absorption amount when loaded in the Y direction in the result of the static incremental analysis of the global collapse type model M2, such that the story energy difference eY when loaded in the Y direction in equation ( ₂ ) is minimized, as the minimum analytical difference value of the energy absorption amount of each story in the X direction, the incremental step iwcY and incremental step iosY corresponding to the analytical value _sEwcY ( _iwcY , j) of the energy absorption amount of each story in the story-collapse type model M1 and the analytical value _sEosY ( _iosY , j) of the energy absorption amount of each story in the global collapse type model M2 in this minimum difference combination are set as the first minimum _difference incremental step _iwcY and the second minimum difference incremental step _iosY , respectively. It can be identified as _osY .

Based on the above considerations, the incremental step identification unit 24 calculates, for each layer 11, the difference between the sum of the analytic value _sEwcX (i, j) of the energy absorption amount of each layer in the layer collapse type model M1 and the analytic value _sEosX (i, j) of the energy absorption amount of each layer in the global collapse type model M2 and the estimated value _sEobsX (j) of the energy absorption amount of each layer for each _nwcX × _nosX combination of the analytic value sEwcX (i, j) of the energy absorption amount of each layer in the layer collapse type model M1 and the analytic value _sEosX (i, j ₎ of the energy absorption amount of each layer in the global collapse type model M2 in the X direction, and selects a minimum difference combination, which is a combination of analytic values that minimizes the sum e _x of the differences in all layers, as a minimum difference analytic value of the energy absorption amount of each layer, and selects the analytic value _sEwcX (i, _j) of the energy absorption amount of each layer in the layer collapse type model M1 from the minimum difference combination. The incremental step i _wcX in the layer-collapse type model M1 and the incremental step i osX in the global collapse type model M2 corresponding to the analytical value sE _osX (i _osX , j) of the energy absorption amount of each layer in the global collapse type model M2 are identified as the incremental step i _wcX with the first smallest difference and the incremental step i _osX with the _second smallest difference, respectively.
Similarly, for each of n wcY ×n osY combinations of the analytic value sE _wcY (i, j) of the energy absorption amount of each layer in the layer collapse type model M1 and the analytic value sE osY (i, j) of the energy absorption amount of each layer in the global collapse type model M2 in the Y direction, the _{incremental step} identification unit 24 calculates, for each layer 11, the difference between the sum of the analytic value sE _wcY (i, j) of the _energy absorption amount of each layer in the layer collapse type model M1 and the analytic value sE _osY (i, j) of the energy absorption amount of each layer in the global collapse type _model M2 and the estimated value sE _obsY (j) of the energy absorption amount of each layer, and selects a minimum-difference combination, which is a combination of analytic values that minimizes the sum e _Y of the differences in all layers, as a minimum-difference analytic value of the energy absorption amount of each layer, and selects the analytic value sE _wcY (i, _j) of the energy absorption amount of each layer in the layer collapse type model M1 in the minimum-difference combination. The incremental step iwcY in the layer-collapse type model M1 and the incremental step i _osY in the total collapse type model M2 corresponding to the analytical value sE _osY (i _osY , j) of the energy absorption amount of each layer in the total collapse type model M2 are identified as the incremental step _iwcY with the first smallest difference and the incremental step i _osY with the _second smallest difference, respectively.

このようにして、健全性判定部２５は、部材ごとに、第１差分最小の増分ステップｉ_ｗｃＸ、ｉ_ｗｃＹに対応する、層崩壊型モデルＭ１における各部材のエネルギー吸収量の解析値ｅＥ_ｗｃＸ（ｉ_ｗｃＸ，ｋ）、ｅＥ_ｗｃＹ（ｉ_ｗｃＹ，ｋ）と、第２差分最小の増分ステップｉ_ｏｓＸ、ｉ_ｏｓＹに対応する、全体崩壊型モデルＭ２における各部材のエネルギー吸収量の解析値ｅＥ_ｏｓＸ（ｉ_ｏｓＸ，ｋ）、ｅＥ_ｏｓＹ（ｉ_ｏｓＹ，ｋ）の和を計算することにより、各部材のエネルギー吸収量の推定値ｅＥ_ｏｂｓ（ｋ）を導出する。 The soundness determiner 25 determines the soundness of the building 10 based on the incremental steps _iwcX and _iwcY with the first smallest difference and the incremental steps iosX and iosY with the second smallest difference in the _X and _Y directions, respectively, identified as described above.
As described above, the incremental step identification unit 24 identifies the incremental steps _iwcX , _iwcY with the smallest first difference and the incremental steps iosX, _iosY with the smallest second difference in each of the _X and Y directions, so that, for an earthquake for which earthquake information has been acquired, it is possible to specify up to what stage, i.e., up to what incremental step, the layer collapse and the global collapse have progressed in each of the static incremental analysis using the layer collapse model M1 and the static incremental analysis using the global collapse model M2. If this is the case, then in an actual earthquake, each component can be regarded as having absorbed energy equal to the sum of the analytical values _eEwcX ( _iwcX , k) and _eEwcY ( _iwcY , k) of the energy absorption amount considered to have been absorbed by the component at the incremental steps _iwcX and _iwcY in the static incremental analysis using the story- _collapse model M1, and the analytical values _eEosX ( _iosX , k) and _eEosY (iosY, k) of the energy absorption amount considered to have been absorbed by the component at the incremental steps _iosX and _iosY in the static incremental analysis using the global collapse model M2.
Therefore, specifically, the soundness determiner 25 calculates the estimated value eE _obs (k) of the amount of energy absorbed by each member by the following equation (3).

In this way, the soundness determination unit 25 derives an estimated value eE obs(k) of the energy absorption amount of each component by calculating, for each component, the sum of the analytical values eE _wcX (i _wcX , k ₎ , eE _wcY (i _wcY , k) of the energy absorption amount of each component in the layer collapse type model M1 corresponding to the incremental step i _wcX , i _wcY with the smallest first difference, and the analytical values eE _osX (i _osX , k), eE _osY (i _osY , k) of the energy absorption amount of each component in the total collapse type model M2 corresponding to the incremental step i _osX , i _osY with the smallest second difference.

The soundness determination unit 25 further determines the soundness of the building 10 based on the estimated energy absorption amount eE _obs (k) of each component derived based on the first difference minimum incremental steps i _wcX , i _wcY and the second difference minimum incremental steps i _osX , i _osY . The soundness determination unit 25 compares the estimated energy absorption amount eE _obs (k) of each component with a preset damage determination threshold to calculate the degree of damage. When the estimated energy absorption amount eE _obs (k) of each component is greater than the damage determination threshold, the soundness determination unit 25 determines that there is a possibility that damage has occurred in that component.
Here, as the damage judgment threshold, for example, a damage judgment threshold μ _cr of the plasticity factor can be set according to the plastic deformation performance of the component. As shown in Fig. 9, when a component is deformed for one cycle at the plasticity factor μ, the energy absorption amount E of the component is expressed by the following formula (4).
E=4(μ-1)Myθy...(4)
Therefore, the damage threshold E _cr of the energy absorption of the member is
E _cr =4(μ _cr -1) M _y θ _y
It becomes.
For example, for a square steel pipe with a half length from the end to the center of 1,750 mm, a yield strength of 295 N/ ^mm2 , cross-sectional dimensions of 300 mm x 300 mm, and a thickness of 9 mm (bending moment at yield My = 327 kNm, member angle at yield θy = 0.0065 rad), if μ _cr = 3.0, then E _cr = 17.0 kNm.
Such damage determination thresholds are not limited to the above, but may be appropriately set by other methods, and may be set in different ways for each type of component, such as pillars 16, beams 17, etc.
The soundness assessment unit 25 outputs information indicating the assessment results, such as the calculated degree of damage and whether or not damage has occurred to the components, via a monitor device or another terminal accessible via the external network 100.
The damage level mentioned above is evaluated on a scale of Damage Level I to Damage Level IV. Damage Level I is defined as a state where there is almost no damage to the columns, beams, and bearing walls. Damage Level II is a state where minor cracks have occurred in the columns, beams, and bearing walls. Damage Level III is defined as a state where significant shear cracks are observed in the columns, beams, and bearing walls, and Damage Level IV is defined as a state where large shear cracks are observed in the columns, beams, and bearing walls, resulting in major damage.

In this way, the soundness assessment unit 25 compares the estimated energy absorption amount eE _obs (k) of each component corresponding to the incremental step i _wcX , i _wcY , i _osX , i _osY with the smallest difference with the damage assessment threshold set for each component to calculate the degree of damage, and assesses the soundness of the building 10 based on the degree of damage of each component.
The soundness of the building 10 is evaluated as safe, caution required, or danger. The safety index based on the soundness evaluation corresponds to damage levels I and II. The caution index based on the soundness evaluation corresponds to damage level III, and the danger index corresponds to damage level IV.

(Soundness assessment method)
FIG. 10 is a flowchart showing the flow of a building soundness evaluation method in the building soundness evaluation system 1.
To evaluate the soundness of the building 10 using the evaluation device 20, a static incremental analysis is performed in advance in the analysis unit 21 using the story collapse model M1 and the global collapse model M2 as a preliminary step. The analysis unit 21 calculates an analytical value of the energy absorption amount of each story and an analytical value of the energy absorption amount of each component for each incremental step, which is a stage at which the load is increased in the static incremental analysis. The calculated analytical values of the energy absorption amount of each story and the analytical values of the energy absorption amount of each component are recorded in the analysis result recording unit 27.
Thereafter, when an earthquake occurs, the evaluation device 20 acquires earthquake information detected by the sensor 12 provided in the building 10 from the external network 100 (step S11).
When the earthquake information is acquired, the energy absorption estimation unit 23 of the evaluation device 20 calculates the story shear force and inter-story displacement of each story 11 based on the earthquake information detected by the sensor 12 and the weight of each story 11 of the building 10, and calculates an estimated value of the energy absorption of each story based on the calculated story shear force and inter-story displacement (step S12).

Next, the incremental step identification unit 24 selects the minimum difference analytical value of the energy absorption amount of each layer, which has the smallest difference from the estimated value of the energy absorption amount of each layer, from the analytical values of the energy absorption amount of each layer recorded in the analysis result recording unit 27, and identifies the corresponding incremental step as the incremental step with the smallest difference (step S13).
More specifically, for each combination of the analysis value of the energy absorption amount of each layer in the layer collapse type model M1 and the analysis value of the energy absorption amount of each layer in the global collapse type model M2 recorded in the analysis result recording unit 27 in the X direction and the Y direction, the incremental step identification unit 24 calculates, for each layer, the difference between the sum of the analysis value of the energy absorption amount of each layer in the static incremental analysis result of the layer collapse type model M1 and the analysis value of the energy absorption amount of each layer in the static incremental analysis result of the global collapse type model M2 and the estimated value of the energy absorption amount of each layer, and selects a minimum-difference combination, which is a combination in which the sum of the differences in all layers is minimum, as the minimum-difference analysis value of the energy absorption amount of each layer, and identifies the incremental step in the layer collapse type model M1 and the incremental step in the global collapse type model M2 corresponding to the analysis value of the energy absorption amount of each layer in the layer collapse type model M1 and the analysis value of the energy absorption amount of each layer in the global collapse type model M2 in the minimum-difference combination as the first minimum-difference incremental step and the second minimum-difference incremental step, respectively.

The soundness assessment unit 25 derives an estimate of the energy absorption amount of each component by calculating, for each component, the sum of the analytical value of the energy absorption amount of each component in the layer collapse type model M1 corresponding to the incremental step with the smallest first difference and the analytical value of the energy absorption amount of each component in the total collapse type model M2 corresponding to the incremental step with the smallest second difference (step S14).
Thereafter, the soundness determiner 25 determines the soundness of the building 10 based on the estimated energy absorption amount of each component derived based on the first minimum difference increment step and the second minimum difference increment step (step S15). The soundness determiner 25 compares the estimated energy absorption amount of each component with a preset damage determination threshold to calculate the degree of damage. If the estimated energy absorption amount of each component is greater than the damage determination threshold, the soundness determiner 25 determines that there is a possibility that damage has occurred in that component.
Furthermore, the soundness assessment unit 25 outputs information indicating the assessment results, such as the calculated degree of damage and whether or not damage has occurred to the components, via a monitoring device or another terminal accessible via the external network 100 (step S16).

(Verification of estimation accuracy using soundness evaluation method)
Here, the estimation accuracy of the building soundness evaluation system 1 of the present invention was verified.
Six models, Model 1 to Model 6, were prepared as the framework for realizing the building to be verified.
Model 1 has four stories, the column span in the X direction is 2, the column span in the Y direction is 1, and the column strength ratio γ, which is the strength ratio of the columns to the beams and panels, is 1.0.
Model 2 has four stories, the column span in the X direction is 2, the column span in the Y direction is 1, and the column strength ratio γ is 1.5.
Model 3 had 8 stories, the column span in the X direction was 3, the column span in the Y direction was 2, and the column strength ratio γ was 1.0.
Model 4 had 8 stories, the column span in the X direction was 3, the column span in the Y direction was 2, and the column strength ratio γ was 1.5.
Model 5 has 12 stories, the column span in the X direction is 4, the column span in the Y direction is 3, and the column strength ratio γ is 1.0.
Model 6 has 12 stories, the column span in the X direction is 4, the column span in the Y direction is 3, and the column strength ratio γ is 1.5.
In all of the above models, the structural characteristic coefficient Ds was set to 0.4.
For each of these models, the Hachinohe wave level 3 was input as the input seismic wave from a direction at an angle of 45° to both the X and Y directions. In this case, the energy absorption amount of each component in each of the X and Y directions calculated from the acceleration information detected by the sensor 12 was taken as the correct value, and the estimated value eE _obs (k) of the energy absorption amount of each component estimated by the soundness evaluation system 1 was taken as the estimated value.

11 to 16 are diagrams showing examples of verification in which the building soundness evaluation system 1 in this embodiment estimates the energy absorption amount of each member constituting a building for each of models 1 to 6. In Fig. 11 to 16, the differences that are the comparison results between the correct values and the estimated values are shown for both ends of columns and beams, and for panels.
From each figure, it can be seen that in the model in which the column strength ratio γ is set to a small value of 1.0, the story collapse mechanism of the lowest story becomes dominant, and the static analysis results of the story collapse model M1 are directly reflected, allowing an estimate close to the correct value to be obtained.
Even in a frame where the column strength ratio γ is set to 1.5 and the overall collapse mechanism is dominant, the energy absorption amount of the columns can be estimated with good accuracy, and the points where the energy absorption of the beams and panels is large are also generally consistent. Here, there are correct values where the energy absorption amount at both ends of a single beam differs greatly, but this is because, in the response to the dominant direction, the strength of one end increases due to the composite effect of the slab, and plasticization hardly occurs, and energy absorption is concentrated at the other end, which is thought to be the cause of the difference from the estimated value based on a static analysis that does not consider the composite effect.
In this verification, examples were presented that were based on steel rigid frame structures; however, the same phenomenon occurs in reinforced concrete structures, in that the ends of columns and beams and panels also become plastic. Therefore, it is believed that the validity is ensured at response levels where extremely severe damage does not occur.

The soundness assessment system 1 for the building 10 as described above is a soundness assessment system 1 for diagnosing and assessing the soundness of the building 10, and performs a static incremental analysis in advance using a three-dimensional frame model M (M1, M2) configured to correspond to the building 10 and including the columns 16 and beams 17. In the analysis, for each incremental step, which is a stage of increasing the load, analytical values _sEwcX (i,j), sEwcY (i,j), _sEosX (i,j), and _sEosY (i,j) of the energy absorption amounts of each story, and analytical values _eEwcX (i,k), _{eEwcY (i,k), eEosX (} i,k), and _eEosY (i,j) of the energy absorption amounts of each member including the columns 16 and _beams 17 are calculated. an analysis unit 21 that calculates the story shear force and story displacement of each story 11 based on the earthquake information and the weight of each story 11 of the building 10, and calculates estimated values sE obsX (j), sE obsY (j) of the energy absorption amount of each story based on the story shear force and the story displacement; and an energy absorption estimation unit 23 that calculates the estimated values sE _obsX (j), sE _obsY (j) of the energy absorption amount of each story from the analytical values sE _wcX (i,j), sE _wcY (i,j), sE _{osX (i,j), and sE osY (} i,j) of the energy absorption amount of each story recorded in the _analysis result recording _{unit 27} . an incremental step identification unit 24 which selects a minimum analytical difference value of the energy absorption amount of each layer that minimizes the difference with the minimum analytical difference value of the energy absorption amount of each layer (a combination of the analytical values _sEwcX (i,j), _sEwcY (i,j ₎ of the energy absorption amount of each layer in the story collapse type model M1 and the analytical values _sEosX (i, _j) , _sEosY (i, _{j) of} the energy absorption amount of each layer in the global collapse type _model M2) and identifies the corresponding incremental step as _{the incremental step iwcX} , _iwcY , _iosX , _{iosY with} the minimum difference _; The system is equipped with a soundness assessment unit 25 that obtains an estimated value eE _obs (k) of the energy absorption amount of each component from _osY (i _osY , k), compares the estimated value eE _obs (k) of the energy absorption amount of each component with a damage assessment threshold set for each component to calculate the degree of damage, and assesses the soundness of the building 10 based on the degree of damage of each component.
According to this configuration, a static incremental analysis is performed using a three-dimensional frame model M (M1, M2) configured to correspond to the building 10 and including the columns 16 and beams 17. In this static incremental analysis, analytical values _sEwcX (i, j), sEwcY (i, j), _sEosX (i, j), sEosY (i, j) of the energy absorption amount of each story and analytical values _eEwcX (i, k), _eEwcY (i, k), _eEosX (i, k), _eEosY (i, k) of the energy absorption amount of each member including the columns 16 and beams 17 are calculated for each _incremental step, which is a _stage of increasing the load, by the analysis unit 21 and recorded in the analysis result recording unit 27. When an earthquake occurs, the energy absorption estimation unit 23 calculates estimated values sE obsX (j) and sE obsY (j) of the energy absorption in each story of the building 10 based on the earthquake information obtained from the sensors 12 installed in the _building 10 and the weight of each _story 11 of the building 10. The soundness determiner 25 selects a minimum difference analytical value of the energy absorption amount of each layer which has a minimum difference from the estimated _value sE _obsX (j), sE _obsY (j) of the energy absorption amount of each layer from the analytical values sE _wcX (i, j), sE _wcY (i, j), sE _osX (i, j), sE osY (i, j) of the energy absorption amount of each layer, and identifies the incremental step for which the minimum difference analytical value corresponding to the selected analytic value sE _wcX (i, j), sE _wcY (i, j), sE _osX (i, j), sE _osY (i, j) of the energy absorption amount of each layer is output as the incremental step with the minimum difference i _wcX , i _wcY , i _osX , i _osY . In the incremental steps _iwcX , _iwcY , _iosX , and _iosY with the smallest difference identified in this manner, the load is increased so that each story 11 absorbs the same amount of energy as the energy absorbed by each story 11 when the earthquake for which earthquake information was obtained occurred. Therefore, the analytical value _eEobs (k) of the energy absorption amount of each member at the stage of the incremental steps _iwcX , _iwcY , _iosX , and _iosY with the smallest difference in the static incremental analysis is considered to be a value close to the amount of energy absorbed by each member when the earthquake for which earthquake information was obtained occurred. Therefore, the estimated value eE obs( _k ) of the energy absorption amount of each component is obtained from the analytical values eE _wcX (i _wcX , _k ), eE _wcY (i _wcY , k), eE _osX (i _osX , k), eE _osY (i _osY , k) of the energy absorption amount of each component corresponding to the incremental step i wcX , i _wcY , i _osX , i osY with the smallest difference, and the degree of damage is calculated based on the estimated value _{eE obs (} k) of the energy absorption amount of each component by comparing it with a damage judgment threshold, thereby making it possible to evaluate the degree of damage to each component constituting the building 10 and judge the soundness of the building 10.
In this way, by estimating the presence or absence of damage to the components that make up the building 10 and the extent of that damage, damage assessment can be performed for each component, allowing for a highly accurate and reliable assessment of the soundness of the building 10.
Furthermore, in the above-described configuration, although damage assessment can be performed for each component, there is no need to provide a sensor for each component, and therefore a large number of sensors are not required.
Therefore, it is possible to provide a building healthiness assessment system 1 that has a simple configuration and can accurately determine healthiness.

The analysis unit 21 also performs a static incremental analysis on each of the story collapse model M1 in which parameters are set for the three-dimensional frame model M (M1, M2) so that the building 10 will undergo story collapse, and the global collapse model M2 in which parameters are set for the three-dimensional frame model M (M1, M2) so that the building 10 will undergo global collapse, and calculates, for each incremental step, analytical values _sEwcX (i,j), sEwcY (i,j), _sEosX (i,j), and _sEosY (i,j) of the energy absorption amount of each story and analytical values eEwcX (i,k), _eEwcY (i,k), _eEosX (i,k), and eEosY (i,k) of the energy absorption amount of each member, and records them in the analysis result recording unit 27. The incremental step identification unit 24 calculates the analytical values _sEwcX (i,j), _sEwcY (i,j), _sEosX (i,j), and sEosY (i,j) of the energy absorption amount of each story in the story collapse model M1 recorded in the analysis result recording unit 27. For each combination of _sEwcX (i,j), _sEwcY (i,j) and the analytical values _sEosX (i,j), _sEosY (i,j) of the energy absorption amount of each layer in the global collapse type model M2, the difference between the sum of the analytical values _sEwcX (i,j), _sEwcY (i,j) of the energy absorption amount of each layer in the layer collapse type model M1 and the analytical values _sEosX (i,j), _sEosY (i,j) of the energy absorption amount of each layer in the global collapse type model M2 and the estimated values _sEobsX (j), _sEobsY (j) of the energy absorption amount of each layer is calculated, and a minimum-difference combination which is a combination that minimizes the sums _eX , _eY of the differences in all layers is selected as a minimum-difference analytical value of the energy absorption amount of each layer, and the analytical values sE The incremental steps _iwcX , iwcY in the story collapse model M1 and the incremental steps i _osX , i osY in the total collapse model M2, which correspond to sE _osX (i, j), sE _osY (i _osY , j) and the analytical values sE _{osX (} i, _j) , sE _osY (i _osY , j) of the energy absorption amounts of each story in the total collapse model M2, are identified as the incremental steps _iwcX , _iwcY with the smallest first difference and the incremental steps i _osX , i _osY with the smallest second difference, respectively, and the soundness determination unit 25 determines the soundness of the building 10 based on the incremental steps _iwcX , _iwcY with the smallest first difference and the incremental steps i _osX , i _osY with the smallest second difference.
According to this configuration, a static incremental analysis is performed on each of a story collapse type model M1 in which parameters are set so that the building 10 undergoes story collapse relative to the three-dimensional frame model M (M1, M2), and a total collapse type model M2 in which parameters are set so that the building 10 undergoes total collapse relative to the three-dimensional frame model M (M1, M2). For each combination of the analytical values _sEwcX (i,j), sEwcY(i,j) of the energy absorption amount of each layer in the layer collapse type model M1 and the analytical values sEosX(i,j), _sEosY (i,j) of the energy absorption amount of each layer in the global collapse type model M2 recorded in the analysis result recording unit 27, a difference between the sum of the analytical values _sEwcX (i,j), _sEwcY (i,j) _of the energy absorption amount of each layer in the layer collapse type model M1 and the analytical values _sEosX (i,j), _sEosY (i,j) of the energy absorption amount of each layer in the global collapse type model _M2 and the estimated values _sEobsX (j), _sEobsY (j) of the energy absorption amount of each layer is calculated for each layer 11, and the sums _eX , eE The minimum difference combination that minimizes _Y is selected as the minimum difference analysis value of the energy absorption amount of each layer. Then, the incremental steps _iwcX , _iwcY in the layer collapse type model M1 and the incremental steps iosX, iosY in the whole collapse type model M2 corresponding to the minimum difference combination selected in this way are identified as the first minimum difference incremental steps _iwcX , _iwcY and the second minimum difference incremental _{steps iosX} , _{iosY ,} respectively.
Here, by combining the states in which the story collapse and the total collapse have progressed to the first minimum difference incremental steps _iwcX , _iwcY and the second minimum difference incremental steps _iosX , _iosY identified as above by the static incremental analysis, the energy absorption amount of each story is close to the estimated energy absorption amount _sEosX (i,j), _sEosY (i,j) at the time of the earthquake whose earthquake information was obtained. In other words, the first minimum difference incremental steps _iwcX , _iwcY and the second minimum difference incremental steps _iosX , _iosY can be regarded as the progress degree of each of the story collapse and the total collapse in the earthquake whose earthquake information was obtained.
By judging the soundness of the building 10 based on the incremental steps _iwcX , _iwcY of the first difference minimum and the incremental steps _iosX , _iosY of the second difference minimum thus identified, it is possible to estimate a collapse state that is closer to the damage caused to the building 10 by an actually occurring earthquake load, and to evaluate the soundness of the building 10 with higher accuracy.

In addition, the soundness determination unit 25 calculates, for each component, the sum eE obs(k) of the analytical values eE _wcX (i _wcX , k), eE _wcY (i _wcY , k) of the energy absorption amount of each component in the layer collapse type model _M1 corresponding to the incremental step i _wcX , i _wcY with the smallest first difference, and the analytical values eE _osX (i _osX , k), eE _osY (i _osY , k) of the energy absorption amount of each component in the total collapse type model M2 corresponding to the incremental step i osX , i _osY with _{the smallest second difference, and compares this as the estimated value eE obs (} k) of the energy absorption amount of each component with the damage determination threshold to calculate the degree of damage.
The analytical values _eEwcX (iwcX, k) and eEwcY ( _iwcY , k) of the energy absorption amounts of each member in the story collapse model M1 corresponding to the incremental steps _iwcX , _iwcY with the smallest first difference identified as above, and the analytical values eEosX (iosX, k) and _eEosY ( _iosY , k) of the energy absorption amounts of each member in the total collapse model M2 corresponding to _the incremental _{steps iosX} , _iosY with the _smallest second _difference are respectively determined. As already explained, the incremental step _iwcX , _iwcY of the smallest first difference and the incremental step iosX, _iosY of the smallest second difference represent the respective progress of story collapse and total collapse in the earthquake for which earthquake information was obtained. Therefore, the sum eEobs(k) of the analytical values _eEwcX ( _iwcX , _k) , eEwcY(iwcY, k), _eEosX ( _iosX , k), _eEosY ( _iosY , _k ) of the energy absorption amounts of each component corresponding to the incremental step _iwcX , _iwcY of the smallest first _difference and the incremental step _iosX , _iosY of _the smallest second difference obtained as described above is considered to be a value close to the energy absorbed by each component when the earthquake for which earthquake information was obtained occurred. This sum is used as an estimated value eE _obs (k) of the energy absorption amount of each component, and the soundness of building 10 is determined by comparing it with the damage assessment threshold for each component, thereby enabling the evaluation of the soundness of building 10 with greater accuracy.

(First Modification of the Embodiment)
The soundness assessment system 1 of the present invention is not limited to the above-described embodiment explained with reference to the drawings, and various other modified examples are conceivable within the technical scope thereof.
For example, in the above embodiment, when the direction in which the earthquake force acts is approximately along the X direction, the incremental step identification unit 24 and the soundness determination unit 25 may perform analysis only in the X direction without performing analysis in the Y direction.
In this case, the incremental step identification unit calculates only the layer energy difference e _x when loading in the X direction in formula (1), and does not need to calculate the layer energy difference e _Y when loading in the Y direction in formula (2). In addition, when calculating the estimated value eE _obs (k) of the energy absorption amount of each member in formula (3), the soundness determination unit 25 adds only the analytical values eE _wcX (i, k) and eE _osX (i, k) of the energy absorption amount of each member in the X direction, and does not use the analytical values eE _wcY (i, k) and eE _osY (i, k) of the energy absorption amount of each member in the Y direction.
The same applies when the direction of the earthquake force is approximately along the Y direction.

(Second Modification of the Embodiment)
Alternatively, in cases where it is known in advance that story collapse will dominate over overall collapse due to the structure of the building, the building health assessment system 1 may be configured to operate only with the story collapse model M1 as the target.
In this case, the analysis unit 21 performs a static incremental analysis on the story collapse model M1 in which parameters are set so that the building 10 will undergo story collapse for the three-dimensional frame model M (M1), and calculates, for each incremental step, the analytical values _sEwcX (i,j), _sEwcY (i,j) of the energy absorption amount of each story and the analytical values _eEwcX (i,k), _eEwcY (i,k) of the energy absorption amount of each component, and records them in the analysis result recording unit 27.
Next, the energy absorption amount estimating unit 23 calculates the estimated values sE _obsX (j) and sE _obsY (j) of the energy absorption amounts in each layer 11, similarly to the above embodiment.
Then, for each of the analytical values _sEwcX (i,j), _sEwcY (i,j) of the energy absorption amount of each layer in the layer collapse type model M1 recorded in the analysis result recording unit 27, the incremental step identification unit 24 calculates the difference between the analytical value _sEwcX (i,j), _sEwcY (i,j) of the energy absorption amount of each layer in the layer collapse type model M1 and the estimated value _sEobsX (j), _sEobsY (j) of the energy absorption amount of each layer for each layer, selects the analytical value of the energy absorption amount that minimizes the summation _eX , _eY of the differences in all layers as the minimum analytical value of the difference of the energy absorption amount of each layer, and identifies the corresponding incremental steps _iwcX , _iwcY in the layer collapse type model M1 as the incremental steps _iwcX , _iwcY with the smallest first difference.
Furthermore, the soundness determination unit 25 calculates, for each component, the sum eE obs (k) of the analytical values eE _wcX (i, _k ), eE _wcY (i, k) of the energy absorption amounts of each component in the layer collapse model M1 corresponding to the incremental steps i _wcX , i wcY with the smallest _{first difference, and compares this as the estimated value eE obs (} k) of the energy absorption amount of each component with the damage determination threshold to calculate the degree of damage.

(Third Modification of the Embodiment)
Furthermore, in cases where it is known in advance that overall collapse is dominant over story collapse due to the structure of the building, the building health assessment system 1 may be configured to operate only with the overall collapse model M2 as the target.
In this case, the analysis unit 21 performs a static incremental analysis on the total collapse type model M2 in which parameters are set so that the building 10 will collapse totally in relation to the three-dimensional frame model M (M2), and calculates, for each incremental step, the analytical values sE _osX (i, j), sE _osY (i, j) of the energy absorption amount of each story and the analytical values eE _osX (i, k), eE _osY (i, k) of the energy absorption amount of each component, and records them in the analysis result recording unit 27.
Next, the energy absorption amount estimating unit 23 calculates the estimated values sE _obsX (j) and sE _obsY (j) of the energy absorption amounts in each layer 11, similarly to the above embodiment.
Then, for each of the analytical values sE _osX (i, j), sE _osY (i, j) of the energy absorption amount of each layer in the global collapse type model M2 recorded in the analysis result recording unit 27, the incremental step identification unit 24 calculates, for each layer, the difference between the analytical value sE _osX (i, j), sE _osY (i, j) of the energy absorption amount of each layer in the global collapse type model M2 and the estimated value sE _obsX (j), sE _obsY (j) of the energy absorption amount of each layer, selects the analytical value of the energy absorption amount that minimizes the summation e _X , e _Y of the differences in all layers as the minimum analytical value of the difference of the energy absorption amount of each layer, and identifies the corresponding incremental steps i _osX , i _osY in the global collapse type model M2 as the incremental steps i _osX , i _osY with the second smallest difference.
Furthermore, the soundness assessment unit 25 calculates, for each component, the sum eE _obs (k) of the analytical values eE _osX (i, _k ), eE _osY (i, k) of the energy absorption amounts of each component in the total collapse type model M2 corresponding to the incremental steps i osX , i osY with the smallest _{second difference, and compares this as the estimated value eE obs (} k) of the energy absorption amount of each component with the damage assessment threshold to calculate the degree of damage.

In this manner, in the above embodiment, the static incremental analysis is performed using a plurality of models M1 and M2, and each of the plurality of models has analysis values in a plurality of directions X and Y. However, as shown in the above first to third modified examples, the static incremental analysis may be performed to have two analysis values. Alternatively, for example, the first modified example and the second modified example may be combined to perform the static incremental analysis to have a total of one analysis value for only one direction of the story collapse model, or the first modified example and the third modified example may be combined to perform the static incremental analysis to have a total of one analysis value for only one direction of the global collapse model.
In addition, without departing from the spirit of the present invention, it is possible to select and discard the configurations given in the above embodiments and modifications, or to change them to other configurations as appropriate.

REFERENCE SIGNS LIST 1 Soundness assessment system 22 Earthquake information recording unit 10 Building 23 Energy absorption estimation unit 11 Story 24 Incremental step identification unit 12 Sensor 25 Soundness determination unit 16 Column 27 Analysis result recording unit 17 Beam M Three-dimensional frame model 19 Panel M1 Story collapse model 21 Analysis unit M2 Whole collapse model

Claims (3)

A soundness evaluation system for diagnosing and evaluating the soundness of a building, comprising:
an analysis unit that performs a static incremental analysis in advance using a three-dimensional frame model configured to correspond to the building and including columns and beams, calculates an analysis value of the energy absorption amount of each story and an analysis value of the energy absorption amount of each member including the columns and beams for each incremental step, which is a step of increasing a load in the analysis, and records the calculated values in an analysis result recording unit;
an earthquake information recording unit that records earthquake information obtained from a sensor installed in the building;
an energy absorption estimation unit that calculates a story shear force and a story displacement of each story based on the earthquake information and the weight of each story of the building, and calculates an estimated value of the energy absorption amount of each story based on the story shear force and the story displacement;
an incremental step identification unit that selects a minimum difference analytical value of the energy absorption amount of each layer that has a minimum difference from the estimated value of the energy absorption amount of each layer from the analytical values of the energy absorption amount of each layer recorded in the analysis result recording unit, and identifies the incremental step corresponding to the minimum difference as the incremental step with the smallest difference;
a soundness determination unit that obtains an estimate of the energy absorption amount of each component from the analysis value of the energy absorption amount of each component corresponding to the incremental step with the smallest difference, compares the estimate of the energy absorption amount of each component with a damage determination threshold value set for each component to calculate a degree of damage, and determines the soundness of the building based on the degree of damage of each component;
A building health assessment system comprising: the analysis unit performs a static incremental analysis on each of a story collapse model in which parameters are set so that the building will undergo a story collapse with respect to the three-dimensional frame model, and a total collapse model in which parameters are set so that the building will undergo a total collapse with respect to the three-dimensional frame model, calculates an analysis value of the energy absorption amount of each of the stories and an analysis value of the energy absorption amount of each of the members for each incremental step, and records them in the analysis result recording unit;
The incremental step identification unit calculates, for each combination of the analysis value of the energy absorption amount of each layer in the layer collapse type model and the analysis value of the energy absorption amount of each layer in the global collapse type model recorded in the analysis result recording unit, a difference between the sum of the analysis value of the energy absorption amount of each layer in the layer collapse type model and the analysis value of the energy absorption amount of each layer in the global collapse type model and the estimated value of the energy absorption amount of each layer, and selects a minimum difference combination that is a combination in which the sum of the differences in all layers is minimum as a minimum difference analysis value of the energy absorption amount of each layer, and identifies the incremental step in the layer collapse type model and the incremental step in the global collapse type model corresponding to the analysis value of the energy absorption amount of each layer in the layer collapse type model and the analysis value of the energy absorption amount of each layer in the global collapse type model in the minimum difference combination as a first minimum difference incremental step and a second minimum difference incremental step, respectively;
The building health assessment system according to claim 1, characterized in that the health assessment unit assesses the health of the building based on the increment step of the first minimum difference and the increment step of the second minimum difference. The building health assessment system described in claim 2, characterized in that the health assessment unit calculates, for each component, the sum of the analytical value of the energy absorption amount of each component in the story collapse type model corresponding to the incremental step with the smallest first difference and the analytical value of the energy absorption amount of each component in the total collapse type model corresponding to the incremental step with the smallest second difference, and uses this as an estimate of the energy absorption amount of each component, compares it with the damage assessment threshold, and calculates the degree of damage.

Images