Calculation device and calculation program

Abstract

To enable total monitoring of a building by determining safety of a layer and a degree of damage to posts caused by interlayer displacement.SOLUTION: A calculating device acquires a strain value measured at each time (t) for a predetermined section of an elastic behavior area of posts from a strain sensor installed so as to measure a section of posts in each floor of a building. The calculating devices obtains horizontal displacement of posts themselves on a floor concerned obtained from a time history of stress and deformation and interlayer displacement of the floor concerned by using a rotational angle of upper edges of posts in a floor below the floor concerned for each floor, and determines safety of layers on the basis of whether interlayer displacement is equal to or more than a threshold. The calculating device obtains a degree of damage by using a predetermined evaluation method of the degree of damage on the basis of a time history of stress and deformation, and determines whether posts are damaged.SELECTED DRAWING: Figure 2

Description

[Technical Field]

[0001] The present invention relates to a computing device and a computing program. [Background technology]

[0002] Conventionally, there is known a technology relating to an inter-story displacement measurement system that can determine inter-story displacement with high accuracy (for example, Patent Document 1). In this technology, target markers installed on predetermined floors of a building 2 are photographed, and the difference in the positions of the target markers in each image is determined as the inter-story displacement.

[0003] Furthermore, there is known a technology for improving the accuracy of estimating the damage status of buildings after an earthquake (for example, Patent Document 2). This technology uses a vibration sensor to acquire actual responses based on an earthquake, and uses the actual responses as input to estimate the damage status of buildings through machine learning. [Prior art documents] [Patent documents]

[0004] [Patent Document 1] Patent Publication No. 2021-169983 [Patent Document 2] Japanese Patent Application Publication No. 2020-128951 Summary of the Invention [Problem to be solved by the invention]

[0005] Conventional technology involves estimating inter-story displacements. However, conventional methods for estimating inter-story displacements involve, for example, installing accelerometers only on a few floors (approximately one to four floors) and using information on the building's vibration mode shapes to estimate only the inter-story displacements, and then using these values ​​to determine the soundness of the building. However, due to the influence of secondary members and the influence of member plasticity, it is difficult to say that the estimation accuracy is sufficient. Furthermore, it is difficult to directly grasp the location and extent of damaged members.

[0006] In consideration of the above, the present invention aims to determine the safety of stories and the degree of damage to columns due to story displacement, thereby enabling integrated monitoring of buildings. [Means for solving the problem]

[0007] To achieve the above object, the calculation device of the present invention includes: an acquisition unit that acquires strain values ​​measured at a predetermined cross section of the column's elastic behavior region at each time t from strain sensors installed to measure the cross section of the column on each floor of a building; a calculation unit that calculates stresses at the upper and lower ends of the column for each floor using the acquired strain values ​​and column cross-sectional information at the measurement positions of the column for each floor, calculates the horizontal displacement of the column itself and the rotation angle of the column top end as deformation values ​​obtained by applying the calculated stresses to a column element model, and outputs the time history of the stresses and the deformations; an inter-story determination unit that calculates inter-story displacements for each floor using the horizontal displacement of the column itself at that floor and the rotation angle of the column top end at a floor below that floor obtained from the time history and determines the safety of the story based on whether the inter-story displacements are equal to or greater than a threshold; and a damage determination unit that calculates the damage level using a predetermined damage level evaluation method based on the time history of the stresses and deformations and determines whether the column is damaged. This enables integrated building monitoring by determining the safety of the story and the damage level of the column due to inter-story displacements.

[0008] The computing device of the present invention may further include a frame determination unit that determines the safety of the frame of each floor based on whether the deformation of the column due to the horizontal displacement of the column itself is equal to or greater than a threshold value. This allows monitoring of the safety of the building frame as well.

[0009] In addition, in the calculation device of the present invention, the damage evaluation method can also determine column damage by applying at least one of yield judgment, plasticity factor, and cumulative plastic deformation magnification, thereby making it possible to quantify the damage degree by taking multiple methods into consideration.

[0010] In the calculation device of the present invention, the acquisition unit acquires strain values ​​of each of the cross sections, Cross Section 1 and Cross Section 2, and the calculation unit calculates the stress by finding the stress at the lower end a1 and the upper end b1 in a column element model including the column end and the cross section as elements from a geometric relationship, and then determines the deformation by performing an elasto-plastic analysis on the stress using a predetermined stiffness equation. This allows analysis to be performed using the characteristics of the strain value of the column itself.

[0011] The calculation program of the present invention causes a computer to execute the following processes: acquire strain values ​​measured at each time t for a predetermined cross section of the column's elastic behavior region from strain sensors installed to measure the cross section of the column on each floor of a building; calculate the stress at the upper and lower ends of the column for each floor for each time t using the acquired strain value and column cross-sectional information at the measurement position of the column; calculate the horizontal displacement of the column itself and the rotation angle of the column top end as the deformation value obtained by applying the calculated stress to a column element model; output the time history of the stress and the deformation; calculate the inter-story displacement of each floor using the horizontal displacement of the column itself at that floor and the rotation angle of the column top end on floors below that floor obtained from the time history; determine the safety of the story based on whether the inter-story displacement is above a threshold; calculate the damage level using a predetermined damage level evaluation method based on the time history of the stress and deformation, and determine whether the column is damaged. [Effects of the Invention]

[0012] According to the present invention, it is possible to determine the safety of a story and the degree of damage to columns due to story displacement, thereby enabling comprehensive monitoring of a building. [Brief explanation of the drawings]

[0013] [Figure 1] 10 is an example of installation of a strain sensor on a pillar in this embodiment. [Figure 2] FIG. 1 is a block diagram showing the configuration of a computing device according to an embodiment of the present invention. [Figure 3] FIG. 10 is a diagram showing an example in which a strain sensor is installed in an elastic behavior region. [Figure 4] FIG. 10 is a diagram showing an example in which a beam element with end rotation springs is used as a column element model. [Figure 5] 10 is a flowchart showing a calculation process in the calculation device according to the present embodiment. DETAILED DESCRIPTION OF THE INVENTION

[0014] [Embodiments of the present invention] Hereinafter, an embodiment of the present invention will be described with reference to the drawings.

[0015] Before describing this embodiment, the background on which this embodiment is based will be described.

[0016] Structural health monitoring systems (SHM systems) are becoming increasingly popular. Practical SHM systems install accelerometers, such as MEMS sensors, on the floors of buildings and use their data to estimate inter-story displacements and assess building safety. While it would be ideal to install displacement sensors on each floor of a building to directly measure inter-story displacements, this approach presents several challenges, including the need for jigs to install sensors between the floor and ceiling, the influence of deformation of the installed beams, and the difficulty of securing installation space for maintenance. Therefore, MEMS sensors, which are less accurate than conventional servo accelerometers but are easier to install and relatively inexpensive, have become the norm. While accelerometers are sometimes installed on every floor of a building, cost constraints mean that they are typically limited to a few floors (e.g., one to four). Because these sensors are installed only on a few floors, techniques for estimating inter-story displacements are often used, such as using design models (frame models of the structural framework) or building vibration mode shapes (hereinafter referred to as building mode shapes) based on measurements at the time of completion.

[0017] However, building mode shapes based on design models generally do not take into account the effects of secondary members, etc., and often do not match actual measurement results. Building mode shapes based on measurements at the time of completion mostly show characteristics at the minute amplitude level, since it is not possible to subject an actual building to large vibrations. For this reason, it is difficult to say that estimating the inter-story displacement of a building using building mode shapes will produce appropriate results. To begin with, building mode shapes represent the characteristics of a building in an elastic state, and there are limitations to using them to capture the behavior of the building when it is damaged.

[0018] On the other hand, while story displacement can be an effective indicator for the broader purpose of determining whether a building can remain occupyable after an earthquake, it is difficult to say that it is an accurate indicator on its own for judging the soundness of the frame because it includes the influence of bending deformation of the building frame. Furthermore, it is difficult to determine the location and extent of "damaged components" when a building is damaged using this indicator. To do so, it is necessary to install and monitor additional sensors (such as strain gauges) at or near the predicted damage points of components that are thought to be susceptible to damage. However, there are challenges, such as the need to install sensors in multiple locations because they are installed in components rather than in the layers, and the difficulty of maintaining sensors installed in areas such as beam ends, which are hidden by fire-resistant coatings and finishing materials. As a result, technology for monitoring the location and extent of damaged components is currently not advancing.

[0019] Now that SHM systems are beginning to become widespread, there is a need for new measurement technologies that can measure or estimate story displacement more accurately and identify damaged components and the extent of that damage. Therefore, the method of this embodiment estimates story displacement as well as the degree of damage to columns, enabling integrated monitoring of buildings.

[0020] FIG. 1 shows an example of strain sensor installation on columns in this embodiment. In this embodiment, as shown in FIG. 1, strain sensors 90 for determining stress are installed on column members that are continuous from the lower floor to the upper floor. The strain sensors 90 are installed on the left and right sides of the vertical direction to measure the cross-sectional stress of the columns on each floor of the building, and measure the cross-sectional strain value. Strain gauges, optical fiber sensors, etc. are used as strain sensors 90. Note that the installation locations of the strain sensors 90 are not limited to the left and right sides of the cross section, and they may be installed on, for example, all four sides of the column. While FIG. 1 shows an example of a configuration that measures strain in the column axis direction, a configuration that also measures strain perpendicular to the column axis may be used to directly determine the shear force of the column.

[0021] 2 is a block diagram showing the configuration of a computing device according to this embodiment. The computing device 100 includes an acquisition unit 110, a calculation unit 112, an inter-story determination unit 114, a damage determination unit 116, and a frame determination unit 118. The computing device 100 is realized by a computer including a CPU (Central Processing Unit), a ROM (Read Only Memory) that stores programs for implementing various processing routines, a RAM (Random Access Memory) that temporarily stores data, a memory as a storage means, a network interface, and the like.

[0022] The acquisition unit 110 acquires strain values ​​measured at each time t at a predetermined cross section in the elastic behavior region of the column from the strain sensor 90. The elastic behavior region generally refers to the region near the mid-height of the column, where plastic deformation and buckling do not occur even during a major earthquake. The strain sensors 90 are installed at positions required to determine the stress (axial force, moment, shear force, etc.) of the column in the elastic behavior region. Figure 3 shows an example of strain sensors installed in the elastic behavior region. Strain sensors 90 are installed on the side of the column at two cross-sectional positions (cross section 1 and cross section 2) within the elastic behavior region of the column (D in Figure 3) to measure the strain in the column axial direction.

[0023] Explain the stress of the column. Regarding the strain value, "Cross section 1:1ε 右 , 1ε 左 ”, “Cross section 2:2ε 右 , 2ε 左 " where A is the cross-sectional area of ​​the column, Z is the section modulus of the column, and L is the length between the sections. From these, the axial force, moment, and shear force can be calculated. The axial force is A·(1ε 右 +1ε 左 ) / 2 or A·(2ε 右 +2ε 左 ) / 2. The moment is cross section 1:-Z·(1ε 右 -1ε 左 ) / 2, cross section 2:-Z·(2ε 右 -2ε 左 ) / 2. The shear force is (Z·(1ε 右 -1ε 左 ) / 2-Z·(2ε 右 -2ε 左 ) / 2) / L.

[0024] As described above, the acquisition unit 110 acquires the strain values ​​of the cross sections 1 and 2 for the columns on each floor.

[0025] The calculation unit 112 calculates the stress at the top and bottom ends of the column for each floor at time t using the acquired strain value and the column cross-sectional information at the measurement position of the column. The calculation unit 112 calculates the horizontal displacement of the column itself and the rotation angle of the column top end as the deformation value obtained by applying the calculated stress to the column element model, and outputs the time history of the stress and deformation.

[0026] The principles of the column element model that can be applied to the calculations of the calculation unit 112 will be explained below.

[0027] (Estimation of stress and deformation using the column element model of the i-th floor) In the following, a column element model is assumed. Using the measured column stress, the stresses at both ends a1 and b1 of the column element model, which includes the column end and cross section as elements, are calculated from the geometric relationship.

[0028] The support condition for the displacement at the bottom end of the column element model is fixed, and an elastic-plastic analysis is performed using the stresses at both ends a1 and b1 to determine the behavior of the column. Note that the elastic-plastic characteristics are set by taking into account a design model (for example, a trilinear system model) or the results of separate experiments.

[0029] For example, a beam element with rotational springs at the ends shown in Figure 4 is used as a column element model. From the stress (axial force N, moment M, shear force Q) at the position of cross section 1 obtained from the strain sensor, the stresses at both ends a1 and b1 of the column element model shown in Figure 4 are calculated as follows:

[0030] Axial force: N a1 =-N, N b1 =N Moment: M a1 =-(M+Q L a1 ), M b1 =MQ·L b1 Shear force: Q a1 =-Q, Q b1 =Q

[0031] The distance between a-a1 and b-b1 is 0. L a1 , L b1 is the length La1 from the bottom end a1 of the column element model to the upper section 1, and the length L from the top end b1 to the lower section 1 b1 The stresses at the top and bottom ends of the column are calculated from the above. Next, the stiffness equation for the beam element with end rotation springs is given by {S} t =[K]{u} t where {S} t , {u} t , [K] is expressed as follows:

[0032]

number

[0033] If the support condition for the displacement at the bottom of the column element model is fixed and the displacement at the top of the column is unknown, the stiffness equation for the incremental displacement expression is {ΔS} t =[K]{Δu} t Here, {ΔS} t , {Δu} t is expressed as follows:

[0034]

number

[0035] The stiffness matrix [K] is the tangent stiffness at time t corresponding to the elastic-plastic properties of the member end rotational spring. For the stiffness matrix [K], an element stiffness matrix that relates the member end load and member end displacement expressed in the member coordinate system, as explained in "Element Stiffness Matrix" in Reference 1, can be used. [Reference 1] "Frame Structure Analysis (Matrix Structure Analysis)," http: / / www.arc.hokkai-su.ac.jp / ~kusiyama / MC_1 / structure_analysis.html

[0036] Using the stiffness equation for incremental displacement, elastic-plastic analysis is performed on the stress obtained by measurement to calculate the displacement {u} of the column element model at each moment. t For the sake of convenience, the description is based on the assumption of displacement in one horizontal direction, but in reality, displacement in two horizontal directions is assumed, and δ and θ are added.

[0037] (Deformation of the column itself on the i-th floor) The estimated time history of deformation is the value when the bottom end of the column is fixed. The horizontal displacement δ of the column on the i-th floor obtained b1 of i δ, rotation angle θ of the top end of the column b1 of i The deformation is used in the inter-story determination unit 114 and the frame determination unit 118.

[0038] The inter-story determination unit 114 calculates the inter-story displacement for each floor using the horizontal displacement of the column itself on that floor obtained from the time history of stress and deformation and the rotation angle of the column top end on the floor below that floor.Then, the inter-story determination unit 114 determines the safety of each floor based on whether the inter-story displacement is equal to or greater than a threshold value, and outputs the determination result.

[0039] The method for determining the inter-story displacement will be described. The inter-story determination unit 114 determines the displacement (inter-story displacement) of the i-th floor when viewed as the entire building. i Find Δ. Story drift of i-th floor: i Δ= i δ+ i Θ× i H

[0040] where: i H is the height of the i-th floor, i Θ is the rotation angle of the bottom end of the column on the i-th floor when viewed as the entire building, and is calculated using the following formula. Rotation angle of column bottom end on i-th floor: i Θ= i-1 Θ+ i-1 θ

[0041] For example, if the i-th floor is the third floor, the rotation angle of the top of the column for the second floor, the first floor, and all floors below the third floor can be calculated to calculate the inter-story displacement of the i-th floor. i Δ is calculated. θ simply represents the rotation angle of the upper end when the lower end of the column element model is fixed. In reality, the lower end of the column is tilted by the column on the floor below, causing rotation, so this amount is added from the bottom floor to evaluate the rotation (floor rotation) as a story. Then, the inter-story determination unit 114 calculates the inter-story displacement as a story. i It is determined whether Δ exceeds a certain threshold value. This determination means a rough determination of the safety of the entire frame, or a determination of whether curtain walls or exterior walls have fallen off. The determination result of the safety of a story can be in any format that can be regarded as a determination result, such as the degree of safety or whether safety is present or not. As described above, by estimating and determining the inter-story displacement for each floor, the determination result of the inter-story displacement is output.

[0042] The damage determination unit 116 determines whether or not the pillar is damaged based on the time history of stress and deformation and using a predetermined damage level evaluation method, and outputs the determination result regarding the damage to the pillar.

[0043] In calculating the degree of damage, the elements necessary for evaluating the degree of damage are extracted from the time history of stress and deformation obtained using a column element model, and a damage evaluation method is used to determine whether the column is damaged or not. If the beam element with end rotation springs shown in Figure 4 is used as the column element model, the relationship between the moment and rotation angle of the end rotation springs and the relationship between axial force and axial displacement are the targets. Damage evaluation methods include "yield judgment," "plasticity factor (plastic deformation magnification)," and "cumulative plastic deformation magnification."

[0044] When using the evaluation method for determining stress yield, the axial force and bending moment at the end of the column are used to determine whether the stress state exceeds the yield surface of the column member. In this case, it is not a matter of determining the degree of damage, but rather determining whether damage has occurred.

[0045] When using the evaluation method of the ductility factor, the maximum deformation angle θ of the member u (Reference 2 c θ pmax ) is the elastic deformation angle θ p Plasticity factor (θ u / θ p ) minus the elastic component (≒1) u / θ p The damage level is determined by whether the damage level is equal to or greater than a threshold value. [Reference 2] "Architectural Institute of Japan: Guidelines for Plastic Design of Steel Structures, 3rd Edition, 2017: Figure C1.5.2 Limiting Plastic Deformation During Cyclic Deformation"

[0046] When using the cumulative plastic deformation ratio evaluation method, refer to Reference 2 and c θ pl Calculate θ u The value of R calculated as (=θ u / θ p -1) is the degree of damage.

[0047] The damage assessment method may use at least one of the above methods, or multiple assessment methods may be used to obtain assessment results. As described above, by outputting the damage assessment results for each floor, the damage level of the columns themselves for each floor of the building can be estimated. The assessment results regarding column damage can be in any format that can be regarded as an assessment result, such as the degree of damage or the presence or absence of damage.

[0048] The frame determination unit 118 determines the safety of the frame of each floor based on whether the deformation of the column due to the horizontal displacement of the column itself is equal to or greater than a threshold value, and outputs the determination result. The determination result of the safety of the frame can be in any form that can be regarded as a determination result, such as the degree of safety or the presence or absence of safety. Deformation of the column itself i It is determined whether δ exceeds a certain threshold value. This determination means determining the safety of the frame of the i-th floor.

[0049] Next, the operation of the computing device 100 of this embodiment will be described. Fig. 5 is a flowchart showing the computation process in the computing device 100 according to this embodiment. The CPU reads out and executes programs and various data from the ROM, causing the CPU to function as each part of the computing device 100 and perform computation process.

[0050] In step S100, the acquisition unit 110 acquires, from the strain sensor 90, strain values ​​measured at each time t for a predetermined cross section of the elastic behavior region of the pillar.

[0051] In step S102, the calculation unit 112 calculates the stress at the top and bottom ends of the column for each floor at time t using the acquired strain value and the column cross-sectional information at the measurement position of the column.

[0052] In step S104, the calculation unit 112 calculates the horizontal displacement of the column itself and the rotation angle of the column top end as the deformation value obtained by applying the calculated stress to the column element model for each floor at time t, and outputs the time history of the stress and deformation.

[0053] In step S106, the inter-story determination unit 114 determines the inter-story displacement for each floor using the horizontal displacement of the column itself on that floor obtained from the time history of stress and deformation and the rotation angle of the column top end on the floor below that floor.

[0054] In step S108, the inter-story determination unit 114 determines the safety of each story based on whether the inter-story displacement is equal to or greater than a threshold value, and outputs the determination result of the safety of the story.

[0055] In step S110, the damage assessment unit 116 assesses whether or not the columns are damaged for each floor based on the time history of stress and deformation using a predetermined damage level assessment method, and outputs the assessment result regarding the damage to the columns.

[0056] In step S112, the structure determination unit 118 determines the safety of the structure of each floor based on whether the deformation of the column due to the horizontal displacement of the column itself is greater than or equal to a threshold value, and outputs the determination result of the safety of the structure.

[0057] As described above, the computing device 100 according to this embodiment determines the safety of a story and the degree of damage to columns due to story displacement, enabling integrated monitoring of a building.

[0058] The method for estimating the story displacement and column damage level described in this embodiment has the following advantages.

[0059] In this embodiment, a strain sensor 90 is used, but unlike conventional techniques, the sensor does not need to be installed at the same location as or near the location where damage is expected.

[0060] This embodiment requires only the characteristics of the column itself, and is not affected by secondary members, unlike conventional estimation methods that use building modes. Furthermore, the characteristics of the column itself, to which no secondary members or slabs are attached, can be accurately represented in the design model. Therefore, it is possible to estimate story displacements with higher accuracy than with conventional estimations using accelerometers. Furthermore, it is possible to accurately estimate the deformation of the column itself.

[0061] Since the stress in the column itself is measured directly, it can be used to determine the yield of the column members. Furthermore, the degree of damage to the column can be determined using the stress-displacement relationship in the column element model.

[0062] In the past, when installing a displacement meter, it was necessary to install it within the structural surface of a partition wall, etc., along with a jig, etc., which made installation difficult due to architectural design and planning restrictions. However, in this embodiment, the strain sensor is installed on the column itself, so it is less likely to be subject to architectural design and planning restrictions.

[0063] The present invention is not limited to the above-described embodiment, and various modifications and applications are possible without departing from the spirit and scope of the present invention. Points to note when applying this embodiment are as follows.

[0064] It is desirable that the columns on which the strain sensors are installed are made of steel. When applying to CFT, SRC, or RC structures, it is necessary to keep in mind that it is difficult to measure strain on concrete members and that estimation accuracy using the column element model (affected by accuracy in terms of rigidity and plastic properties) is important.

[0065] It is desirable that the columns on which strain sensors are installed are columns on which no load acts other than at both ends (both joints). For example, in cases where a seismic wall is installed on a column or where a highly rigid brace is connected to the column, it is necessary to note that estimation cannot be made using a simple model such as the beam element with rotational springs at the end shown as an example.

[0066] In the above embodiment, as shown in Fig. 1, strain sensors for determining stress are installed on column members that are continuous from the lower floor to the upper floor. However, this embodiment may also be applied to a specific column. In this case, although it is not possible to grasp the inter-story displacement including the influence of bending deformation, it is possible to estimate the deformation of the column itself and the degree of column damage.

[0067] In addition, any method for determining the degree of damage to a column can be used in this embodiment as long as it is determined from data on all or any of the displacement, rotation angle, and column stress of the column itself. Also, regarding the type of strain sensor, any means that can measure strain can be used in this embodiment. [Explanation of symbols]

[0068] 100 Computing equipment 110 Acquisition Department 112 Calculation Department 114 Interlayer determination section 116 Damage determination section 118 Frame determination section

Claims

1. an acquisition unit that acquires strain values ​​measured at each time t for a predetermined cross section of the elastic behavior region of the column from strain sensors installed to measure the cross section of the column on each floor of the building; a calculation unit that calculates, for each floor at time t, stresses at the upper and lower ends of the column using the acquired strain value and column cross-sectional information at the measurement position of the column, calculates the horizontal displacement of the column itself and the rotation angle of the column upper end as a deformation value obtained by applying the calculated stress to a column element model, and outputs the time history of the stress and the deformation; an inter-story determination unit that determines, for each floor, the inter-story displacement of the floor using the horizontal displacement of the column itself at the floor and the rotation angle of the column top end at a floor below the floor, obtained from the time history, and determines the safety of the floor based on whether the inter-story displacement is equal to or greater than a threshold value; a damage determination unit that determines the degree of damage using a predetermined damage degree evaluation method based on the time history of the stress and deformation, and determines whether the column is damaged; 1. A computing device comprising:

2. The calculation device according to claim 1 , further comprising a frame determination unit that determines the safety of the frame of each floor based on whether the deformation of the column due to the horizontal displacement of the column itself is equal to or greater than a threshold value.

3. The calculation device according to claim 1 or 2, wherein the damage assessment method determines damage to the column by applying at least one of a yield judgment, a plasticity factor, and a cumulative plastic deformation magnification.

4. the acquisition unit acquires strain values ​​of each of the cross sections, cross section 1 and cross section 2, The calculation unit calculates the stress by determining the stress at the lower end a1 and the upper end b1 in a column element model including the column end and cross section as elements from a geometric relationship, and determines the deformation by performing an elastic-plastic analysis on the stress using a predetermined stiffness equation. A computing device according to any one of claims 1 to 3.

5. A strain value measured at each time t for a predetermined cross section of the elastic behavior region of the column is acquired from a strain sensor installed to measure the cross section of the column on each floor of the building, For each floor, for time t, calculate the stress at the upper and lower end portions of the column using the acquired strain value and the column cross-sectional information at the measurement position of the column, apply the calculated stress to a column element model to calculate the horizontal displacement of the column itself and the rotation angle of the column upper end as the deformation value obtained, and output the time history of the stress and the deformation, For each floor, a story displacement of the floor is calculated using the horizontal displacement of the column itself at the floor and the rotation angle of the column top end at a floor below the floor, which are obtained from the time history, and the safety of the floor is determined based on whether the story displacement is equal to or greater than a threshold value; calculating the degree of damage using a predetermined damage degree evaluation method based on the time history of the stress and deformation, and determining whether the column is damaged; A calculation program that causes a computer to execute a process.

Images