JP2017198610A - Safety diagnostic device, safety diagnostic method, and safety diagnostic program - Google Patents

Abstract

PROBLEM TO BE SOLVED: To provide a generally applicable method for diagnosing the safety of a building after an earthquake occurrence.SOLUTION: The safety diagnostic device includes: an eigenvalue analysis unit for performing an eigenvalue analysis, using a response value obtained from a sensor and the mass of a predetermined building, and determining the value of the rigidity of the building as well as determining an eigenmode and a participation function based on the combination of the mass and the value of the rigidity; an input value estimation unit for determining an input value of an earthquake wave, using a transfer function obtained from a natural frequency estimated with the response value and from the participation function, and the response value; a response value estimation unit for determining an absolute acceleration response by adding a relative acceleration response obtained by multiplying the difference between the waveform of each eigenmode of the response value obtained using the transfer function and the input value by an eigenmode shape, to the input value; and a safety diagnostic unit for integrating the absolute acceleration response and determining the safety of a structure based on an interlayer deformation angle obtained by an absolute displacement obtained for each floor of the structure and the height of floors of the structure.SELECTED DRAWING: Figure 2

Description

  The present invention relates to a safety diagnostic device, a safety diagnostic method, and a safety diagnostic program.

  Conventionally, when a large-scale earthquake occurs, the owner or user of the building cannot judge the safety and may not know whether the building can be used continuously. In addition, it takes a long time for an expert in a building structure to conduct a field survey to determine safety, and there is a possibility that the accuracy may not be sufficient by, for example, visual judgment alone.

  In addition, a seismometer is installed only in the lower part of the structure to measure seismic waves and input to a simple model to calculate the response value of each layer. If the corresponding damage level is greater than or equal to a predetermined threshold level, a more detailed model A technique of calculating a response value based on the above has been proposed (for example, Patent Document 1).

JP 2012-168008 A

  Conventionally, a system has been proposed in which an input value is measured by a seismometer provided only in the lower part of the structure (for example, the first floor) to determine damage to the structure due to the earthquake. However, it is necessary to prepare a pre-determined analysis model and calculate the response value. Depending on the structure, it is necessary to create a new dynamic analysis model based on the design document, which is not universal. It was. In addition, when the response value is obtained analytically, the validity of the estimated value cannot be evaluated at all, that is, the system is established only under the condition that the estimated value (or dynamic analysis model) is correct. The establishment of the condition remains as an issue.

  Then, an object of this invention is to provide the safety diagnostic method of the building after the earthquake occurrence which can be applied universally.

  The safety diagnostic device diagnoses the safety of the structure after the earthquake occurs. Specifically, eigenvalue analysis is performed using a response value acquired from a sensor for measuring a response value generated in the structure and a predetermined mass of the structure to obtain a rigidity value of the structure, and the mass and An eigenvalue analysis unit for obtaining an eigenmode and a stimulus function based on a combination of the obtained stiffness value, a natural frequency estimated from the response value, and a transfer function obtained using the stimulus function obtained by the eigenvalue analysis unit; An input value estimation unit that obtains an input value of the seismic wave for the structure using the response value, a waveform for each eigenmode of the response value obtained using the transfer function, and an input value obtained by the input value estimation unit A response value estimation unit that obtains an absolute acceleration response or an absolute velocity response by adding a relative acceleration or relative velocity obtained by multiplying the difference between them by an eigenmode shape obtained by eigenvalue analysis and an input value, and a response value Integrate the absolute acceleration response or absolute velocity response obtained by the fixed part, and based on the absolute displacement obtained for each floor of the structure and the interlayer deformation angle obtained using the floor height of the structure, the safety of the structure A safety diagnosis unit for determining

By using the mass of the predetermined structure, the safety of the structure can be diagnosed based on the response value without using an analysis model. That is, the stiffness value can be obtained by eigenvalue analysis, and the input value and the response value of each layer can be estimated. Therefore, it can be said that a safety diagnosis method that can be applied universally can be provided.

  Further, the eigenvalue analysis unit may obtain the value of the stiffness matrix so that the ratio of the secondary or higher natural frequency to the primary natural frequency approaches the actual measurement value. Specifically, the stiffness value can be obtained in this way.

  In addition, a transfer function obtained by multiplying a single-degree-of-freedom transfer function obtained based on each natural frequency estimated from the response value and a predetermined damping constant by a stimulus function and combining them with multiple degrees-of-freedom is obtained. A transfer function setting unit may be further provided. Thus, the transfer function can be obtained by using a predetermined attenuation constant.

  The contents described in the means for solving the problems can be combined as much as possible without departing from the problems and technical ideas of the present invention. The contents of the means for solving the problems can be provided as a device such as a computer or a system including a plurality of devices, a method executed by the computer, or a program executed by the computer. A recording medium that holds the program may be provided.

  ADVANTAGE OF THE INVENTION According to this invention, the safety diagnostic method of the building after the earthquake occurrence which can be applied universally can be provided.

It is a figure which shows an example of the system configuration | structure which concerns on this embodiment. It is a functional block diagram showing an example of a system concerning this embodiment. It is an apparatus block diagram which shows an example of a computer. It is a processing flowchart which shows an example of a safety diagnostic process. It is a graph which shows an example of the value of the acceleration which the acceleration sensor measured. It is a graph which shows an example of the Fourier amplitude spectrum calculated | required from acceleration. It is a figure which shows an example of the setting of mass. It is a figure which shows an example of the initial value of the rigidity ratio used for a search. It is a figure which shows an example of an eigenmode shape. It is a graph which shows an example of the transfer function (amplitude) calculated | required about 1 degree of freedom. It is a graph which shows an example of the transfer function (phase difference) calculated | required about 1 degree of freedom. It is a graph which shows an example of the transfer function (amplitude) developed (combined) with multiple degrees of freedom. It is a graph which shows an example of the transfer function (phase difference) developed by multiple degrees of freedom. It is a figure which shows an example of the estimated value of the input wave in the 1st floor. It is a graph which shows the example which isolate | separated the waveform which the acceleration sensor measured from the primary mode to the tertiary mode. It is a figure which shows an example of the estimated value of the absolute acceleration response in a certain floor. It is a figure which shows an example of the largest interlayer deformation angle calculated about between each floor. It is a figure which shows the example of an output of the diagnostic result of safety.

  Hereinafter, embodiments of the present invention will be described with reference to the drawings. In addition, the following embodiment is an illustration and this invention is not limited to the following structure. In the present embodiment, the safety of each floor of the structure is determined from the response value of the structure measured by at least one acceleration sensor when an earthquake occurs. Safety judgment shall be made on the ground floor.

<System configuration>
FIG. 1 is a diagram illustrating an example of a system configuration according to the present embodiment. The system of FIG. 1 includes an acceleration sensor 1 and a safety diagnostic device 2. In the example of FIG. 1, the present system is provided in the seven-story structure 3 above the ground, but the height, structure, construction method, and the like of the structure 3 are not particularly limited. The acceleration sensor 1 is a triaxial acceleration sensor, for example, and measures the acceleration applied to the device itself. In addition, one acceleration sensor 1 is provided on two or more floors (also referred to as “measurement floors”). When an earthquake occurs, the response on the measurement floor of the structure 3 is not an input value on the first floor (ground). Measure the value. Here, the acceleration measured on the first floor of the structure 3 is called an “input value”, and the acceleration measured on the second floor or more of the structure 3 is called a “response value” with respect to earthquake motion. In the example of FIG. 1, the acceleration sensor 1 is disposed on the roof of the structure 3. The acceleration sensor 1 and the safety diagnostic device 2 are communicably connected, and the measurement value measured by the acceleration sensor 1 is output to the safety diagnostic device 2. The safety diagnosis apparatus 2 is a general computer, and realizes functions to be described later when a processor executes a program according to the present embodiment.

<Functional configuration>
FIG. 2 is a functional block diagram illustrating an example of a system according to the present embodiment. The system of FIG. 2 includes an acceleration sensor 1 and a safety diagnostic device 2, and both are connected. The acceleration sensor 1 continuously measures the acceleration applied to its own device, and transmits the measured value to the safety diagnostic device 2. Further, for example, when the measured acceleration exceeds a predetermined threshold, the acceleration sensor 1 may transmit the value of the acceleration measured in a predetermined period to the safety diagnostic device 2 or the measured acceleration. May be transmitted to the safety diagnostic apparatus 2.

  The safety diagnostic apparatus 2 includes a measurement value acquisition unit 201, a storage unit 202, a natural frequency estimation unit 203, a natural value analysis unit 204, a transfer function setting unit 205, an input value estimation unit 206, and a mode decomposition unit. 207, a response value estimation unit 208, an interlayer deformation angle estimation unit 209, and a safety diagnosis unit 210.

  The measurement value acquisition unit 201 acquires the acceleration measured by the acceleration sensor 1.

  The storage unit 202 is composed of a main storage device or an auxiliary storage device, and stores data input to the system, output data, data for setting operation conditions, data generated in the middle of processing, and the like. . In the present embodiment, it is assumed that information such as the number of floors of the structure 3, the mass of each floor, the initial value used for the stiffness search process described later, and the floor height is stored in the storage unit 202 in advance.

  The natural frequency estimation unit 203 performs frequency analysis on the acceleration acquired from the acceleration sensor 1 to estimate the natural frequency. The natural frequency appears as a local peak (that is, local maximum) in the Fourier transform spectrum of acceleration.

  The eigenvalue analysis unit 204 performs eigenvalue analysis of a mass system model having degrees of freedom (from the second floor to the rooftop) and sets eigenmodes and stimulus functions. Specifically, eigenvalue analysis is performed by setting a predetermined mass and changing the stiffness value. In addition, an appropriate stiffness is searched using a predetermined evaluation function, and an eigenmode and a stimulation function are estimated from a result of eigenvalue analysis performed using a predetermined mass and appropriate stiffness.

  The transfer function setting unit 205 sets a transfer function at the measurement floor where the acceleration sensor 1 is arranged. Specifically, using the estimated natural frequency, stimulation coefficient, and predetermined damping constant, the transfer function of the measurement floor for the first floor is synthesized into a transfer function with multiple degrees of freedom in the complex number region.

Moreover, the input value estimation unit 206 estimates an input wave on the first floor using the measured acceleration and the transfer function. Specifically, the input waveform on the first floor is obtained by dividing the Fourier amplitude spectrum of the acceleration measured by the acceleration sensor 1 by the transfer function and returning to the waveform using inverse Fourier transform.

  The mode decomposing unit 207 multiplies the Fourier amplitude spectrum of the acceleration measured by the acceleration sensor 1 by the ratio of the absolute value of the transfer function of each mode to the transfer function of multiple degrees of freedom, and separates the waveform of each mode.

  The response value estimation unit 208 multiplies the difference between the measured acceleration-specific waveform and the estimated first-order input wave by the eigenmode shape, and calculates the relative acceleration for each mode on each floor. Moreover, the absolute acceleration response in each floor is estimated by adding up the relative acceleration for each mode and the estimated first floor input wave in each floor. Alternatively, the absolute acceleration response in each floor may be estimated by multiplying the estimated first floor input wave by a transfer function set for each layer.

  The interlayer deformation angle estimation unit 209 calculates the absolute displacement by integrating the absolute acceleration response of each floor by the second floor. And the difference of an absolute displacement is calculated | required about each between continuous upper and lower floors, and the amount of interlayer deformation is obtained. Further, the interlayer deformation angle is obtained by dividing the interlayer deformation amount by the floor height.

  The safety diagnosis unit 210 reads a predetermined correspondence between the interlayer deformation angle and the safety from the storage unit 202, and diagnoses safety on each floor based on the interlayer deformation angle obtained by the interlayer deformation angle estimation unit 209. . In addition, the diagnosis result can be output by various methods. For example, it may be displayed on the display of the safety diagnostic apparatus 2, or an e-mail or the like may be transmitted via a network (not shown) such as the Internet and notified to a predetermined destination. Moreover, you may make it print from a printer (not shown), while electric power is supplied by the uninterruptible power supply.

<Device configuration>
FIG. 3 is an apparatus configuration diagram illustrating an example of a computer. The safety diagnostic apparatus 2 is a computer as shown in FIG. 3, for example. A computer 1000 shown in FIG. 3 includes a CPU (Central Processing Unit) 1001, a main storage device 1002, and an auxiliary storage device (external storage device).
1003, a communication IF (Interface) 1004, an input / output IF (Interface) 1005, a drive device 1006, and a communication bus 1007. The CPU 1001 performs processing and the like according to this embodiment by executing a program. The main storage device 1002 caches programs and data read by the CPU 1001 and secures a work area for the CPU. Specifically, the main storage device is a RAM (Random Access Memory), a ROM (Read Only Memory), or the like. The auxiliary storage device 1003 stores programs executed by the CPU 1001, setting information used in the present embodiment, and the like. Specifically, the auxiliary storage device 1003 is an HDD (Hard-disk Drive), an SSD (Solid State Drive), an eMMC (embedded Multi-Media Card), a flash memory, or the like. The main storage device 1002 and the auxiliary storage device 1003 function as the storage unit 202. The communication IF 1004 transmits / receives data to / from other computers. The communication IF 1004 is specifically a wired or wireless network card or the like. The input / output IF 1005 is connected to the input / output device and accepts input from the user or outputs information to the user. Specifically, the input / output device is a keyboard, a mouse, a display, a touch panel, the acceleration sensor 1, and the like. The drive device 1006 reads data recorded on a storage medium such as a magnetic disk, a magneto-optical disk, and an optical disk, and writes data to the storage medium. The above components are connected by a communication bus 1007. A plurality of these components may be provided, or some of the components (for example, the communication IF 1004 and the drive device 1006) may not be provided. Further, the input / output device may be integrated with the computer. In addition, the program executed in this embodiment is provided via a portable storage medium readable by the drive device 1006, a portable auxiliary storage device 1003 such as a flash memory, a communication IF 1004, and the like. It may be. Then, when the CPU 1001 executes the program, the above-described computer is operated as the safety diagnostic apparatus 2 shown in FIG. 2, for example. Note that at least a part of the configuration of the exemplified computer may exist on the network. For example, a so-called cloud-like mode in which a service is provided by one or more servers that bear at least a part of the functional units illustrated in FIG.

<Safety diagnosis processing>
FIG. 4 is a process flowchart showing an example of the safety diagnosis process. This process is executed, for example, when the acceleration sensor 1 detects an acceleration exceeding a predetermined threshold.

  First, the measurement value acquisition unit 201 of the safety diagnostic apparatus 2 acquires the acceleration value measured by the acceleration sensor 1 and stores it in the storage unit 202 (FIG. 4: S1). As described above, the acceleration sensor 1 according to the present embodiment measures a response value measured on the second or higher floors, not the input value measured on the ground (first floor).

FIG. 5 is a graph showing an example of acceleration values measured by the acceleration sensor 1. In the graph of FIG. 4, the horizontal axis represents time (seconds), and the vertical axis represents acceleration values (cm / s ² ). The measurement value acquisition unit 201 continuously acquires, for example, an acceleration value as shown in FIG. 4 until the magnitude of the acceleration becomes a predetermined threshold value or less, and stores it in the storage unit 202.

Next, the natural frequency estimation unit 203 performs frequency analysis on the acceleration value acquired by the measurement value acquisition unit 201 to estimate the natural frequency (FIG. 4: S2). Specifically, the natural frequency estimation unit 203 performs a Fourier transform on the acceleration values shown in FIG. 5 to obtain a Fourier amplitude spectrum as shown in FIG. In the graph of FIG. 6, the horizontal axis represents the frequency (Hz) and the vertical axis represents the Fourier amplitude ((cm / sec ² ) · sec). The natural frequency appears as a local maximum value in the Fourier amplitude spectrum. Note that in order from the lowest frequency, the primary natural frequency (FIG. 6 (1)), the secondary natural frequency (FIG. 6 (2)), and the third natural frequency (FIG. 6 (3)). ) ... The natural frequency estimating unit 203 obtains a Fourier amplitude spectrum, and searches for a local maximum value of the Fourier amplitude, for example, from a lower frequency to a higher frequency, and specifies a predetermined number of natural frequencies. For example, the natural frequencies from the first to the third order are obtained. For example, an upper limit of the frequency may be determined in advance such as 7 Hz, and the natural frequency appearing from 0 Hz to the upper limit may be searched.

  Further, the eigenvalue analysis unit 204 sets the mass and rigidity, performs eigenvalue analysis, and estimates the eigenmode and the stimulation function (FIG. 4: S3). The eigenvalue analysis is performed using a mass system model with degrees of freedom. That is, in this embodiment, the analysis model is not used, and instead, the mass and rigidity of each mass point (each rank) are set. In this step, the mass is set as a predetermined value, and an appropriate value is searched for while changing the stiffness value.

As for the mass, a mass matrix (M matrix) based on the mass ratio of each floor is determined with a predetermined floor as a reference. FIG. 7 is a diagram illustrating an example of setting the mass. FIG. 7 is a graph in which the horizontal axis represents the mass ratio between layers (floors) and the vertical axis represents layers (floors). In the example of FIG. 7, assuming that the first floor is 1 as a reference value, a ratio that is reduced by 1% is set every time the layer rises by one floor. This can be said to reflect the tendency of the cross-sectional area of building columns and beams to decrease toward the upper floor. In the example of FIG. 7, the rooftop mass is set to 1.5 times the mass of the top floor. This takes into account the load of equipment installed on the roof. In this embodiment, by using such simple values, safety diagnosis can be performed even when values based on the design model or the like cannot be set.

  As for the rigidity, a rigidity matrix (K matrix) based on the rigidity ratio of each floor is obtained based on a predetermined floor. Specifically, eigenvalue analysis is performed while changing the stiffness ratio, and an appropriate value that approximates the measurement value acquired in S1 is searched using a predetermined evaluation function. FIG. 8 is a diagram illustrating an example of the initial value of the stiffness ratio used for the search. FIG. 8 is a graph in which the horizontal axis represents the rigidity ratio between layers (floors) and the vertical axis represents layers (floors). The initial values in FIG. 8 are values set based on the seismic layer shear force distribution coefficient (that is, “Ai distribution”). Specifically, the shear force of each layer based on the Ai distribution obtained from the mass of each mass point shown in FIG. 7 is normalized by the value of the first floor.

ここで、Ｃは評価値、ｆａ_ｉは推定されたｉ次の固有振動数（Ｈｚ）、ｆｍ_ｉはＳ２において求められたｉ次の固有振動数（Ｈｚ）、βａ_ｉはＳ３において推定されたｉ次の刺激係数、ｎは予め定められたモード次数の上限値である。評価関数は、１次の固有振動数に対する２次以上の固有振動数の比に刺激係数による重みを乗じ、予め定められた上限次数までの二乗和を表すものである。刺激関数は、質量及び各次数のモード（変位）を用いて求めることができる。一般的に、モード次数が低いほど刺激係数（刺激関数の最上層での値に相当）の絶対値が大きくなる性質があるため、各次数の刺激係数の絶対値で重みづけしている。なお、以下の式（２）に基づいて固有値ω^２が求められる。
（［Ｋ］−ω^２［Ｍ］）｛φ｝＝０ ・・・（２）
［Ｋ］は剛性マトリクス、［Ｍ］は質量マトリクス、｛φ｝は固有ベクトルである。ωは固有円振動数であり、以下の式（３）により上述の振動数ｆが求められる。
ｆ＝ω／２π ・・・（３） The search procedure is as follows. First, the mass of each layer and the rigidity of the lowermost layer (one layer) are fixed values, and the rigidity of each layer excluding the uppermost layer (rooftop) is increased or decreased within a range of ± 10% of the initial value, for example. Eigenvalue analysis is performed on the combinations. Since the boundary conditions of the top floor are different from those of other layers, the fluctuation range is set to ± 20% with respect to the initial value, for example. Then, the evaluation function defined by the following formula (1) is repeatedly calculated to search for a combination of mass and rigidity that minimizes the evaluation value C.

Here, C is the evaluation value, fa _i is the estimated i-th natural frequency (Hz), fm _i is the _i-th natural frequency (Hz) obtained in S2, and βa _i is estimated in S3. An i-th order stimulation coefficient, n is an upper limit value of a predetermined mode order. The evaluation function represents a sum of squares up to a predetermined upper limit order by multiplying the ratio of the secondary or higher natural frequency with respect to the primary natural frequency by the weight of the stimulation coefficient. The stimulation function can be obtained using the mass and the mode (displacement) of each order. In general, the lower the mode order, the larger the absolute value of the stimulation coefficient (corresponding to the value at the top layer of the stimulation function), so weighting is performed with the absolute value of the stimulation coefficient of each order. The eigenvalue ω ² is obtained based on the following equation (2).
([K] −ω ² [M]) {φ} = 0 (2)
[K] is a stiffness matrix, [M] is a mass matrix, and {φ} is an eigenvector. ω is a natural circular frequency, and the above-described frequency f is obtained by the following equation (3).
f = ω / 2π (3)

  As described above, the rigidity of each layer can be specified by the least square method so that the error from the natural frequency ratio obtained from the actual measurement value in S2 is minimized. Further, the eigenmode shape as shown in FIG. 9 can be estimated based on the mass ratio and the rigidity ratio of each layer, and the stimulation function can be obtained. The absolute values of the rigidity and mass of each layer set here do not match the actual building values. However, since the ratio of the natural frequency is adopted as the evaluation value, if the ratio in the height direction of the set rigidity and mass is reasonable, the eigenvalue obtained from the eigenvalue analysis result using those values will be used. Mode shapes and stimulus functions are also reasonable values.

Further, the transfer function setting unit 205 multiplies the transfer function of each order obtained from the natural frequency and damping constant of each order obtained in S2 by the stimulation function obtained in S3 and develops it with multiple degrees of freedom. A transfer function in the complex number domain is derived (FIG. 4: S4). The attenuation constant of the first-order mode may be a general value such as 2% in the case of steel structure (S), 3% in the case of reinforced concrete structure (RC), or the increase in amplitude of measurement data. For example, a value that increases linearly may be adopted. For the attenuation constant of the second order mode or higher, the same concept as that used in the design such as the same as the first mode (constant mode) or the rigidity proportional type may be adopted. Furthermore, for the Fourier amplitude of the observation record, it is assumed that only the dominant component of each mode is flat in the Fourier amplitude of the input wave. A method for obtaining an attenuation constant may be used.

  FIG. 10 is a graph showing an example of a transfer function (amplitude) obtained for one degree of freedom. In the graph of FIG. 10, the horizontal axis represents the frequency (Hz) and the vertical axis represents the transfer function (amplitude). FIG. 11 is a graph showing an example of a transfer function (phase difference) obtained for one degree of freedom. In the graph of FIG. 11, the horizontal axis represents the frequency (Hz) and the vertical axis represents the transfer function (phase difference). The transfer function (amplitude) is the ratio of the response value at each floor to the input value at the lowest floor. In other words, it can be said that the transfer function (amplitude) represents the magnitude of fluctuation with reference to the first floor. . The transfer function (phase difference) is the directionality of the response at each floor with respect to the input at the lowest floor. When the phase difference is 0 degree, the response vibrates in the same direction with respect to the input, similarly 180 degrees. In the case of, the response to the input is oscillating in the opposite direction.

  FIG. 12 is a graph showing an example of a transfer function (amplitude) developed (combined) with multiple degrees of freedom. In the graph of FIG. 12, the horizontal axis represents the frequency (Hz) and the vertical axis represents the transfer function (amplitude). FIG. 13 is a graph showing an example of a transfer function (phase difference) developed with multiple degrees of freedom. In the graph of FIG. 13, the horizontal axis represents the frequency (Hz) and the vertical axis represents the transfer function (phase difference). By combining the transfer functions for one degree of freedom shown in FIGS. 10 and 11 in the complex number domain, the transfer functions as shown in FIGS. 12 and 13 can be obtained.

  Further, the input value estimation unit 206 calculates an input value on the first floor using the acceleration value obtained in S1 and the transfer function obtained in S4 (FIG. 4: S5). The input wave in the lowest layer (first floor) is obtained by dividing the frequency analysis result of the acceleration value obtained in S1 by the transfer function obtained in S4. Specifically, an input wave to the first floor is obtained by multiplying the Fourier amplitude spectrum subjected to frequency analysis in S2 by the inverse of the multi-degree-of-freedom transfer function in the complex domain and performing inverse Fourier transform.

FIG. 14 is a diagram illustrating an example of an estimated value of an input wave on the first floor. In the graph of FIG. 14, the horizontal axis represents time (seconds) and the vertical axis represents acceleration (cm / s ² ). In S5, for example, a value as shown in FIG. 14 is obtained.

  After that, the mode decomposing unit 207 multiplies the frequency analysis result obtained in S2 by the ratio of the absolute value of the transfer function of each mode to the multi-degree-of-freedom transfer function obtained in S4, and separates the waveform of each mode. (FIG. 4: S6). In this step, the waveform of each mode is obtained by multiplying the Fourier amplitude spectrum subjected to frequency analysis in S2 by the ratio of the absolute value of the transfer function of each mode to the transfer function of multiple degrees of freedom and performing inverse Fourier transform. FIGS. 15 (1) to (3) are examples of graphs obtained by separating the waveform measured by the acceleration sensor from the primary mode to the tertiary mode.

  Moreover, the response value estimation part 208 calculates the response value in each floor (FIG. 4: S7). First, the relative acceleration for each mode on each floor is obtained by multiplying the difference (that is, relative acceleration) between the waveform for each mode obtained in S6 and the first floor input value obtained in S5 by the eigenmode shape obtained in S3. It is done. And the estimated value of the absolute acceleration response in each floor is obtained by adding together the relative acceleration for every mode and the first floor input value calculated | required by S5 about each floor.

FIG. 16 is a diagram illustrating an example of an estimated value of an absolute acceleration response in a certain floor. In the graph of FIG. 16, the horizontal axis represents time (seconds) and the vertical axis represents acceleration (cm / s ² ). In S4, for example, values as shown in FIG. 16 are obtained.

  Then, the interlayer deformation angle estimation unit 209 calculates an interlayer deformation angle in each floor (FIG. 4: S8). In this step, the absolute displacement is obtained by second-order integration of the absolute acceleration response of each floor. And the difference of an absolute displacement is calculated | required about each between continuous upper and lower floors, and the amount of interlayer deformation is obtained. Further, the interlayer deformation angle is obtained by dividing the interlayer deformation amount by the floor height of each floor.

  FIG. 17 is a diagram illustrating an example of the maximum interlayer deformation angle calculated for each floor. In the graph of FIG. 17, the horizontal axis represents the maximum interlayer deformation angle (rad.), And the vertical axis represents the rank. In S8, for example, values as shown in FIG. 17 are obtained.

  Further, the safety diagnosis unit 210 determines the safety of each floor based on the interlayer deformation angle in each floor, and outputs the result (FIG. 4: S9). It is assumed that information indicating safety is stored in advance in the storage unit 202 in association with the range of the magnitude of the interlayer deformation angle. The safety diagnosis unit 210 reads out information indicating safety corresponding to the interlayer deformation angle calculated in S8 from the storage unit 202, for example, for each floor, and outputs the information.

  FIG. 18 is a diagram illustrating an output example of the safety diagnosis result. In the output example of FIG. 18, for each floor of the structure 3, the estimated seismic intensity, the diagnosis result regarding the structure, the diagnosis result regarding the ceiling / furniture, and the diagnosis result regarding the facility / equipment are displayed as precautions regarding the use of the structure 3. Has been. The diagnosis result is a display of the safety diagnosis apparatus 2 as shown in FIG. 13 in a plurality of stages S9 representing the degree of safety such as “safety”, “caution”, “danger”, etc. It may be displayed on the screen, or may be transmitted to a predetermined destination such as an owner of the structure 3, a user, a management company, a construction company, or other related parties by e-mail or the like. In addition, printing from a printer or the like may be used in preparation for a power failure.

<Effect>
According to the safety diagnostic device 2 according to the present embodiment, based on the measurement data from one acceleration sensor 1 and simple settings held in advance, the safety of the structure can be obtained without using an analysis model. Sex can be diagnosed. That is, by using the mass of each layer of the structure that is simply set and searching for the rigidity of each layer that approximates the natural frequency ratio based on the actual measurement value, the estimation of the first-order input wave and the response value of each layer Estimation is possible. The response value is a waveform that includes the characteristics of the actual shaking of the structure. According to the present embodiment, the estimated value can be obtained in such a manner that the response value measured on a certain floor is developed on other floors. In this respect, the validity of the estimated value according to the present embodiment is higher than that in the case where all the shaking of the structure is analytically obtained using the input value that is the waveform of the ground motion that does not reflect the shaking of the structure. I can say that.

<Modification>
Instead of the acceleration sensor 1, a speed sensor may be used. Also in this case, for example, by performing the first order integration, the displacement can be obtained from the velocity of each layer, and the interlayer deformation angle can be calculated.

  A plurality of acceleration sensors 1 may be used. Also in this case, a plurality of response values on the upper floor are acquired instead of the input values on the first floor. By using a plurality of measured values, a mode amplitude ratio between measurement floors can be added to an evaluation value that is only the natural frequency ratio, and the accuracy of approximation performed in S3 of FIG. 4 can be improved.

  The mass ratio based on the predetermined layer shown in FIG. Moreover, when the design book of the structure 3 exists, the ratio based on the design book may be set in advance.

The initial values of the stiffness search shown in FIG. However, the calculation amount for obtaining the approximate value can be reduced by using the Ai distribution ratio as the initial value as in the embodiment.

  In addition, microtremor measurement is performed on all floors in advance, and the natural frequency, mode shape, damping constant, and stimulation function of each order are obtained in advance from the analysis of the obtained measurement results. You may substitute suitably as a setting numerical value in a process.

<Others>
Further, the present invention is not limited to the above-described embodiment, and can be modified within the scope not departing from the gist of the present invention. In addition, the above-described embodiments and the modifications mentioned as appropriate can be implemented in combination as much as possible.

  The present invention also includes a computer program that executes the above-described processing and a computer-readable recording medium that records the program. The recording medium on which the program is recorded can perform the above-described processing by causing the computer to execute the program. The above-described processing may be performed by a seismometer, or may be performed by a smartphone or the like.

  Here, the computer-readable recording medium refers to a recording medium in which information such as data and programs is accumulated by electrical, magnetic, optical, mechanical, or chemical action and can be read from a computer. Examples of such a recording medium that can be removed from the computer include a flexible disk, a magneto-optical disk, an optical disk, a magnetic tape, and a memory card. In addition, examples of the recording medium fixed to the computer include a hard disk drive and a ROM.

DESCRIPTION OF SYMBOLS 1 Acceleration sensor 2 Safety diagnostic apparatus 201 Measurement value acquisition part 202 Storage part 203 Natural frequency estimation part 204 Eigenvalue analysis part 205 Transfer function setting part 206 Input value estimation part 207 Mode decomposition part 208 Response value estimation part 209 Interlayer deformation angle estimation Part 210 Safety diagnosis part 3 Structure

Claims (5)

A safety diagnostic device for diagnosing the safety of structures after an earthquake,
An eigenvalue analysis is performed using a response value acquired from a sensor for measuring a response value generated in the structure and a predetermined mass of the structure to obtain a stiffness value of the structure, and the mass and An eigenvalue analysis unit for obtaining an eigenmode and a stimulus function based on a combination with the obtained stiffness value;
Input value estimation for obtaining an input value of an earthquake wave for the structure using the natural frequency estimated from the response value, a transfer function obtained using the stimulus function obtained by the eigenvalue analysis unit, and the response value And
A relative value obtained by multiplying the difference between the waveform for each eigenmode of the response value obtained using the transfer function and the input value obtained by the input value estimation unit by the eigenmode shape obtained by the eigenvalue analysis. A response value estimation unit that obtains an absolute acceleration response or an absolute speed response by adding the acceleration or relative speed and the input value;
Integrating the absolute acceleration response or absolute velocity response obtained by the response value estimation unit, based on the absolute displacement obtained for each floor of the structure and the interlayer deformation angle obtained using the floor height of the structure , A safety diagnostic unit for judging the safety of the structure,
A safety diagnostic device comprising: The safety diagnostic apparatus according to claim 1, wherein the eigenvalue analysis unit obtains a value of a stiffness matrix such that a ratio of a secondary or higher natural frequency to a primary natural frequency approaches an actual measurement value. The transfer function is obtained by multiplying the stimulus function by a transfer function of one degree of freedom obtained based on each natural frequency estimated from the response value and a predetermined damping constant, and synthesizing it with multiple degrees of freedom. The safety diagnosis device according to claim 1, further comprising a transfer function setting unit. A safety diagnostic method for diagnosing the safety of a structure after an earthquake,
An eigenvalue analysis is performed using a response value acquired from a sensor for measuring a response value generated in the structure and a predetermined mass of the structure to obtain a stiffness value of the structure, and the mass and An eigenvalue analysis step for obtaining an eigenmode and a stimulation function based on a combination with the obtained stiffness value;
An input value for obtaining an input value of the seismic wave for the structure using the natural frequency estimated from the response value, a transfer function obtained using the stimulation coefficient obtained in the eigenvalue analysis step, and the response value An estimation step;
Obtained by multiplying the difference between the waveform for each eigenmode of the response value obtained using the transfer function and the input value obtained in the input value estimation step by the eigenmode shape obtained by the eigenvalue analysis. A response value estimation step of obtaining an absolute acceleration response or an absolute speed response by adding the relative acceleration or the relative speed and the input value;
Integrating the absolute acceleration response or absolute velocity response obtained in the response value estimation step, based on the interlayer displacement angle obtained using the absolute displacement obtained for each floor of the structure and the floor height of the structure A safety diagnosis step for determining safety of the structure;
A safety diagnostic method executed by a computer. A safety diagnostic program for diagnosing the safety of structures after an earthquake,
An eigenvalue analysis is performed using a response value acquired from a sensor for measuring a response value generated in the structure and a predetermined mass of the structure to obtain a stiffness value of the structure, and the mass and An eigenvalue analysis step for obtaining an eigenmode and a stimulation function based on a combination with the obtained stiffness value;
An input value for obtaining an input value of the seismic wave for the structure using the natural frequency estimated from the response value, a transfer function obtained using the stimulation function obtained in the eigenvalue analysis step, and the response value An estimation step;
Obtained by multiplying the difference between the waveform for each eigenmode of the response value obtained using the transfer function and the input value obtained in the input value estimation step by the eigenmode shape obtained by the eigenvalue analysis. A response value estimation step of obtaining an absolute acceleration response or an absolute speed response by adding the relative acceleration or the relative speed and the input value;
Integrating the absolute acceleration response or absolute velocity response obtained in the response value estimation step, based on the interlayer displacement angle obtained using the absolute displacement obtained for each floor of the structure and the floor height of the structure A safety diagnosis step for determining safety of the structure;
Is a safety diagnostic program that runs a computer.

Images