JP6916682B2 - How to estimate the seismic response of any part of the building - Google Patents
How to estimate the seismic response of any part of the building Download PDFInfo
- Publication number
- JP6916682B2 JP6916682B2 JP2017139146A JP2017139146A JP6916682B2 JP 6916682 B2 JP6916682 B2 JP 6916682B2 JP 2017139146 A JP2017139146 A JP 2017139146A JP 2017139146 A JP2017139146 A JP 2017139146A JP 6916682 B2 JP6916682 B2 JP 6916682B2
- Authority
- JP
- Japan
- Prior art keywords
- building
- floor
- mode
- measurement
- vibration
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Description
本発明は、地震により建物が振動したときの応答を推定する方法である。 The present invention is a method of estimating the response when a building vibrates due to an earthquake.
従来の地震観測システムでは、各階のセンサー設置位置など代表点のみについて、応答評価や応答予測を行い、被災度を評価している。同じ階でも偏心がある建物や、ねじれが生じやすい建物においては、平面上の位置による応答に差が大きくなるが、より応答の大きいと考えられる個所を評価することができなかった。 In the conventional seismic observation system, the degree of damage is evaluated by performing response evaluation and response prediction only for representative points such as sensor installation positions on each floor. In a building with eccentricity or a building that is prone to twisting even on the same floor, the difference in response depending on the position on the plane becomes large, but it was not possible to evaluate the part that is considered to have a larger response.
特許文献1(特開2014−211397号公報)には、建物の最下階と最上階に設置したセンサー2つの計測結果である加速度波形を記録し、地震記録応答推定装置により建物の最下階と最上階の加速度波形の差分を求めるとともに最上階の基本応答波を求め、建物固有の振動数及び振動モード形の係数に基づき、基本応答波を建物の複数次の固有振動数毎の応答波形に分離し、また、各階の対応する次数のモード係数を掛け合せてバンドパス波形を求め、各階毎のバンドパス波形を足し合わせて合成応答波形を求め、さらに各階における相対加速度を示す合成応答波形に最下階センサー3の計測結果である加速度波形を足し合わせ絶対加速度波形を求め、この絶対加速度波形を数値積分し、地震時の速度、変位の波形に変換することによって、限られた階に設置したセンサーで得られた建物の地震時応答情報に基づき、より精度よく建物各階の応答を推定することを可能にする建物の地震時応答/健全性確認方法が開示されている。
この先行発明は、各階の代表点について、応答推定を行うシステムであり、あらかじめ設定する振動モード情報は、設計パラメータとしての建物固有の振動数および振動モード形としている。
この先行発明では、(1)各階の代表点についてのみしか評価できないため、偏心のある建物、平面形状が細長い建物など、平面上の位置による応答の差が大きい建物は、被害状況を正確に伝えられず、(2)設計パラメータを利用しているため、設計データを得られない建物については推定ができない、また、竣工後、年数が経過している建物や、地震を経験した後の建物は、振動特性が変化することが知られており、築年数の経った建物では、設計時のモデルによる振動モード情報は、実状とは異なる可能性が高い。
In Patent Document 1 (Japanese Unexamined Patent Publication No. 2014-21137), the acceleration waveform which is the measurement result of the two sensors installed on the bottom floor and the top floor of the building is recorded, and the bottom floor of the building is recorded by the seismic record response estimation device. The difference between the acceleration waveform of the top floor and the acceleration waveform of the top floor is obtained, and the basic response waveform of the top floor is obtained. The band path waveform is obtained by multiplying the mode coefficients of the corresponding orders of each floor, and the band pass waveforms of each floor are added to obtain the composite response waveform. Installed on a limited floor by adding the acceleration waveforms that are the measurement results of the
This prior invention is a system that estimates the response for a representative point of each floor, and the vibration mode information set in advance is a building-specific frequency and a vibration mode type as design parameters.
In this prior invention, (1) only the representative points of each floor can be evaluated, so that a building with a large difference in response depending on the position on the plane, such as an eccentric building or a building with an elongated plane shape, accurately conveys the damage situation. It is not possible to estimate for buildings for which design data cannot be obtained because (2) design parameters are used, and for buildings that have been completed for many years or have experienced an earthquake. , It is known that the vibration characteristics change, and in an old building, the vibration mode information by the model at the time of design is likely to be different from the actual situation.
本発明は、解析モデルを持たない建物や設計時から振動性状が変化した振動モードが不明な建物について、実情に基づく地震応答を推定できる地震応答システムを開発することを目的とする。 An object of the present invention is to develop an earthquake response system capable of estimating an earthquake response based on actual conditions for a building having no analysis model or a building whose vibration properties have changed from the time of design and whose vibration mode is unknown.
本発明は、常時微振動や加振機を用いた人工振動などを測定して、建物の上下の限られた高さ方向のモード情報と水平方向のねじれモード情報に基づいて、建物の任意箇所の地震応答を推定する方法であり、この推定に基づいて地震応答や地震の被災状況を判定するシステムである。
本発明の主な構成は次のとおりである。
1.建物に印加する常時微振動又は加振動を測定した上下の限られた高さ方向のモード情報と水平方向のねじれモード情報に基づいて、建物の任意箇所の地震応答を推定する方法であって、高さ方向のモード情報は、上下2点間を一層ずつずらして測定し、層毎に伝達関数を掛け合わせて建物上下全体の変形を推定して得られるモード情報であることを特徴とする建物の任意箇所の地震応答を推定する方法。
2.水平方向のねじれモード情報は、各フロアについてそれぞれのフロアの2点間を測定し、各フロアの2点間の伝達関数と回転中心を求め、フロアのねじれ変形を推定して得られるモード情報であることを特徴とする1.記載の建物の任意箇所の地震応答を推定する方法。
3.1.または2.に記載の建物の任意箇所の地震応答を推定する方法により得られた、建物の任意箇所の地震応答を推定するシステム。
4.1.または2.に記載の建物の任意箇所の地震応答を推定する方法により得られた、建物の任意箇所の地震応答に基づいて、構造物の被災度を判定するシステム。
The present invention constantly measures microearthquakes, artificial vibrations using a vibration exciter, etc., and based on the limited height direction mode information and horizontal twist mode information of the top and bottom of the building, any part of the building. It is a method of estimating the seismic response of the above, and is a system for judging the seismic response and the damage situation of the earthquake based on this estimation.
The main configurations of the present invention are as follows.
1. 1. It is a method of estimating the seismic response of any part of the building based on the limited vertical mode information in the vertical direction and the horizontal twist mode information that measured the constant micro-vibration or excitation applied to the building. mode information in the height direction is between the upper and lower two points were measured by shifting one by one, you characterized in that the mode information obtained by multiplying the transfer function for each layer to estimate the deformation of the entire building vertical how to estimate the seismic response of any part of the building.
2 . The horizontal twist mode information is mode information obtained by measuring between two points on each floor for each floor, obtaining the transfer function and rotation center between the two points on each floor, and estimating the twist deformation of the floor. It is characterized by being 1. How to estimate the seismic response of any part of the serial placement of the building.
3 . 1. 1. Or 2. A system for estimating the seismic response of an arbitrary part of a building obtained by the method of estimating the seismic response of an arbitrary part of a building described in 1.
4 . 1. 1. Or 2. A system for determining the degree of damage to a structure based on the seismic response of an arbitrary part of a building obtained by the method of estimating the seismic response of an arbitrary part of a building described in 1.
1.本発明は、解析モデルを持たない建物や設計時から振動性状が変化した振動モードが不明な建物について、常時微振動などを測定して、建物の上下の限られた高さ方向のモード情報と水平方向のねじれモード情報に基づいて、建物の任意箇所の地震応答を実情に基づいた応答として推定できる。
2.特に、隣接する上下階の2点を測定セットとして1層ずつずらしながら測定して得られた層毎の振動情報に伝達関数を掛け合わせることで1階などの基準階に対する各階の伝達関数を利用して基準階に対する振幅情報と位相情報から高さ方向のモード情報を設定する。
3.また、各階の平面的なねじれは、任意の2点間の振動から回転中心を求めて、高さ方向に積み重ねることにより建物のねじれモードを設定する。
4.本発明は、定常時の微振動に着目して、上下方向も水平方向も2点間の測定情報に基づいて、建物の高さ方向の振動モードも、ねじれモードも推定して、任意箇所の地震応答に利用できるシステムである。
5.各層に基づく高さ方向のモード情報とねじれモード情報が得られるので、建物の任意の箇所の地震応答が推定でき、その推定された地震応答に基づく建物の各部位、部材や設備などに対する地震の被災度を判定することができる。
6.本発明は、把握された建物任意箇所や設備に対する地震応答性を利用して、建物に地震モデル情報を適用することにより、それぞれの箇所の安全度や危険度を予測することができ、改修などに利用することができる。また、実際に遭遇した地震の情報に基づいて、建物のダメージを瞬時に推定することができ、地震直後の対処情報として有用である。例えば、避難の必要の有無、設備の安全稼働などを判断することができる。
1. 1. The present invention constantly measures micro-vibration and the like for a building that does not have an analysis model or a building whose vibration properties have changed from the time of design and whose vibration mode is unknown, and obtains mode information in the limited height direction above and below the building. Based on the twist mode information in the horizontal direction, the seismic response at any part of the building can be estimated as a response based on the actual situation.
2. In particular, the transfer function of each floor with respect to the reference floor such as the first floor is used by multiplying the vibration information for each layer obtained by measuring the two points on the adjacent upper and lower floors as a measurement set while shifting them one layer at a time. Then, the mode information in the height direction is set from the amplitude information and the phase information with respect to the reference floor.
3. 3. Further, for the planar twist of each floor, the twist mode of the building is set by finding the center of rotation from the vibration between arbitrary two points and stacking them in the height direction.
4. The present invention pays attention to the micro-vibration in the steady state, estimates the vibration mode and the twist mode in the height direction of the building based on the measurement information between two points in the vertical direction and the horizontal direction, and estimates the vibration mode and the twist mode at an arbitrary position. It is a system that can be used for earthquake response.
5. Since the height direction mode information and the twist mode information based on each layer can be obtained, the seismic response of any part of the building can be estimated, and the earthquake response to each part, member, equipment, etc. of the building based on the estimated seismic response can be estimated. The degree of damage can be determined.
6. According to the present invention, it is possible to predict the safety level and the degree of danger of each part by applying the seismic model information to the building by utilizing the seismic response to the grasped arbitrary part of the building and the equipment, such as repair. Can be used for. In addition, the damage to the building can be estimated instantly based on the information of the earthquake actually encountered, which is useful as the coping information immediately after the earthquake. For example, it is possible to determine whether or not evacuation is necessary and the safe operation of equipment.
本発明は、建物の地震観測において、地震時の応答や被災度を評価するシステム、限られた観測点から観測していない箇所の応答を推定するシステムに関する。
常時微振動や加振動などを測定して、建物の上下の限られた高さ方向のモード情報と水平方向のねじれモード情報に基づいて、建物の任意箇所の地震応答を推定する方法及び応答推定システムである。なお、測定対象の振動は、風や交通など建物に自然に印加されている振動情報あるいは加振機をもちいて意識的に発生した振動である。
The present invention relates to a system for evaluating the response at the time of an earthquake and the degree of damage in earthquake observation of a building, and a system for estimating the response of an unobserved portion from a limited number of observation points.
A method and response estimation that constantly measures micro-vibration and vibration, and estimates the seismic response of any part of the building based on the limited height direction information and horizontal twist mode information above and below the building. It is a system. The vibration to be measured is vibration that is consciously generated by using vibration information or a vibration exciter that is naturally applied to a building such as wind or traffic.
建物の振動特性を常時微振動測定により把握する際には、建物全体にセンサーを配置し、同時計測を行うのが一般的であるが、測定機器や配線作業が大がかりとなりコストと時間がかかる。センサー数削減のため基準測点を固定したうえで他の測点を移動させながら測定を行う方法もあるが、本発明では、より簡易に建物全体挙動を把握するために、高さ方向の地震応答特性を、上下階をセットとして1層毎にセンサーを移動させながら部分測定を行ったデータを用いる方法を提案する。
また、水平方向では、各階毎の回転中心を把握するために、任意の2点を測定して回転中心をもとめることを提案する。
そして、1階などの基準階を基準にして、建物全体の高さ方向の特性とねじれ特性を高さ方向のモード情報と水平方向のねじれモード情報として、建物の地震応答特性を把握する。
この建物の地震応答特性を利用して、地震被災時に建物に設置したセンサーが感知した情報に基づいて、建物の部材や設備の受けるダメージを推定することができる。建物ダメージを現場確認する前に緊急対応ができることとなり、被災対策を速やかにとることができる。
また、本発明は、建物が受けるダメージを予測することができ、既存建物の地震改修などに利用することができる。
本発明は、地震応答に関する建物のデータが無い建物や設計当初のデータが利用できない建物などにも適用ができる。また、設計データがあっても、竣工後設計との整合性を確認する手法としても活用できる。
したがって、本発明は、常時微振動や加振機を用いた人工振動などを測定して、建物の上下の限られた高さ方向のモード情報と水平方向のねじれモード情報に基づいて、建物の任意箇所の地震応答を推定して、この推定に基づいて地震応答や地震の被災状況を判定するシステムである。
When grasping the vibration characteristics of a building by constant micro-vibration measurement, it is common to place sensors in the entire building and perform simultaneous measurement, but the measuring equipment and wiring work are large and costly and time consuming. In order to reduce the number of sensors, there is a method of measuring while moving other stations after fixing the reference station, but in the present invention, in order to more easily grasp the behavior of the entire building, an earthquake in the height direction We propose a method using data obtained by partially measuring the response characteristics while moving the sensor for each layer with the upper and lower floors as a set.
Further, in the horizontal direction, it is proposed to measure any two points to find the center of rotation in order to grasp the center of rotation for each floor.
Then, with reference to the reference floor such as the first floor, the seismic response characteristics of the building are grasped by using the characteristics in the height direction and the twist characteristics of the entire building as the mode information in the height direction and the twist mode information in the horizontal direction.
Using the seismic response characteristics of this building, it is possible to estimate the damage to the building members and equipment based on the information detected by the sensors installed in the building at the time of the earthquake. Emergency response can be taken before the building damage is confirmed on site, and disaster countermeasures can be taken promptly.
Further, the present invention can predict the damage to a building and can be used for earthquake repair of an existing building.
The present invention can be applied to a building for which there is no building data regarding the seismic response or a building for which the initial design data cannot be used. Moreover, even if there is design data, it can be used as a method for confirming consistency with the design after completion.
Therefore, the present invention constantly measures microearthquakes, artificial vibrations using a vibration exciter, and the like, and based on the limited height direction mode information and horizontal twist mode information of the upper and lower parts of the building, the building It is a system that estimates the seismic response at an arbitrary location and determines the seismic response and the damage status of the earthquake based on this estimation.
<高さ方向のモード情報について>
隣接する上下階2点の測定をセットとして1層ずつずらしながら測定し、得られた層毎の伝達関数を掛け合わせることで、1階などの基準階に対する各階の伝達関数を求める。
得られた伝達関数から、卓越する振動数について基準階に対する振幅情報と位相情報から、各次のモードの形状を求める。
<About mode information in the height direction>
The transfer function of each floor with respect to the reference floor such as the first floor is obtained by measuring the measurements of two adjacent upper and lower floors as a set while shifting each layer one by one, and multiplying the obtained transfer functions of each layer.
From the obtained transfer function, the shape of each next mode is obtained from the amplitude information and phase information with respect to the reference floor for the predominant frequency.
(高さ方向の振動モード把握)
本検討で想定する簡易な部分移動測定を図1に示す。隣接する上下階2点の測定をセットとして1層ずつずらしながら測定する。図示の例では、左から屋上と9階をセットで測定し、以降順次9階と8階のセット・・・・2階と1階のセットとして部分移動して測定することを示している。
得られた層毎の伝達関数を掛け合わせることで、図2に示す1階などの基準階に対する各階の伝達関数を求める。
(Understanding the vibration mode in the height direction)
Figure 1 shows a simple partial movement measurement assumed in this study. The measurement is performed by shifting the measurements of two adjacent upper and lower floors one layer at a time as a set. In the illustrated example, the rooftop and the 9th floor are measured as a set from the left, and then the 9th and 8th floors are set ... The 2nd and 1st floors are partially moved and measured.
By multiplying the obtained transfer functions for each layer, the transfer function of each floor with respect to the reference floor such as the first floor shown in FIG. 2 is obtained.
<平面的なモード情報>
ねじれ振動モード形状把握には、平面上の2点について、卓越振動数と振動方向が同時に得られる等高線図を利用する。
<Plane mode information>
In order to grasp the torsional vibration mode shape, a contour diagram is used in which the dominant frequency and the vibration direction can be obtained at the same time for two points on the plane.
ねじれを伴う平面的な振動モード把握の方法について概念を図3に示す。ねじれ振動モード形状把握には、平面上の2点について、卓越振動数と振動方向が同時に得られる等高線図を利用する。ここで利用する等高線図は次の文献(a)(b)等で、回転スペクトルとよばれており、本明細書でも同様である。
(a)日向仁,肥田剛典,高田毅士:地震観測記録を用いた偏心建物の固有モード形の同定,日本建築学会大会学術講演梗概集,構造II,pp.263〜264,2015.9
(b)江藤公信,林正司,太田勤,田子茂:常時微動測定による建物のねじれ軸の検討-公会堂建築を例にした動的解析-,日本建築学会大会学術講演梗概集,pp.973〜974,1998.9
回転スペクトルは、同時に記録した直交する2方向の時刻歴データ(図3(a)(b))を、ベクトル合成し(図3(c))、さらに0度から180度方向に一定の角度刻みで回転させた座標軸(図3(d))に対するベクトルの正射影成分の時刻歴データをフーリエ変換して求めたフーリエ振幅を、横軸を振動数、縦軸を角度(振動方向)として並べ、等高線グラフで表すものである。卓越する振動数において、どの方向で大きく振動しているか把握できる。
剛床を仮定すると、平面上の2点の、振動卓越方向と直交するそれぞれの線が交わる点が、そのモードにおける振動の回転中心となる。
2点で振動方向に差がでるよう、通常は端部の測定データを用いるが、本発明では、階段室まわりのなどの測定しやすい場所を選定する。比較的近い2点の測定箇所から、回転中心を求めることができる。テナントビルなどでは、セキュリティ上問題の少ない共用部である階段室まわりなど利用しやすいところを測定点として、柔軟に選定する。
FIG. 3 shows a concept of a method for grasping a planar vibration mode accompanied by twisting. Twist vibration mode To grasp the shape, use a contour map that can obtain the dominant frequency and vibration direction at the same time for two points on the plane. The contour maps used here are referred to as rotation spectra in the following documents (a) and (b), and the same applies to the present specification.
(A) Hitoshi Hinata, Takenori Hida, Takeshi Takada: Identification of eccentric building eigenmodes using seismic observation records, Architectural Institute of Japan Conference Academic Lecture Abstract, Structure II, pp.263-264, September 2015
(B) Kiminobu Eto, Masashi Hayashi, Tsutomu Ota, Shigeru Tako: Examination of the twist axis of a building by constant microtremor measurement-Dynamic analysis using public hall architecture as an example-, Architectural Institute of Japan Conference Academic Lecture Abstracts, pp.973- 974, 1998.9
For the rotation spectrum, time history data (FIGS. 3 (a) and (b)) recorded at the same time in two orthogonal directions are vector-combined (FIG. 3 (c)), and further, in a constant angle step from 0 degree to 180 degree. The Fourier amplitude obtained by Fourier transforming the time history data of the normal projection component of the vector with respect to the coordinate axes rotated in (Fig. 3 (d)) is arranged with the horizontal axis as the frequency and the vertical axis as the angle (vibration direction). It is represented by a contour graph. It is possible to grasp in which direction the vibration is large at an outstanding frequency.
Assuming a rigid floor, the intersection of two lines on the plane that are orthogonal to the vibration predominant direction is the center of rotation of vibration in that mode.
Normally, the measurement data at the end is used so that there is a difference in the vibration direction at the two points, but in the present invention, a place that is easy to measure, such as around the staircase, is selected. The center of rotation can be obtained from two relatively close measurement points. For tenant buildings, etc., flexibly select a place that is easy to use, such as around the staircase, which is a common area with few security problems, as a measurement point.
<測定概要>
常時微動の部分移動測定のイメージ
本発明でイメージする簡易な部分移動測定を図1に示す。隣接する上下階をセットとして1層ずつずらしながら測定し、層毎の伝達関数を掛け合わせることで、1階などの基準階に対する各階の伝達関数を求める。また、同じ階で基準点と端部の2点を測定し水平の伝達関数を算出して同様に掛け合わせを行うことで建物全体の立体挙動を把握する。無線で同期とデータ伝送をするセンサーの利用を念頭にしており、例えば階段室で1層ずつの測定を行えば、無線通信が途切れにくいため配線不要となり、共用部の測定のためテナントビルでもセキュリティ上比較的許容されやすいなど、負担の少ない測定が可能である。
ここでは検証のため全層同時計測を行い、違う時間帯のデータを用いることで本手法の有効性を検証する。
<Measurement outline>
Image of partial movement measurement of constant fine movement Fig. 1 shows a simple partial movement measurement imaged in the present invention. The transfer function of each floor with respect to the reference floor such as the first floor is obtained by measuring while shifting the adjacent upper and lower floors one layer at a time and multiplying them by the transfer function of each layer. In addition, the three-dimensional behavior of the entire building is grasped by measuring two points, the reference point and the end, on the same floor, calculating the horizontal transfer function, and performing the same multiplication. We are considering the use of sensors that synchronize and transmit data wirelessly. For example, if we measure one layer at a time in the staircase, wiring is not required because wireless communication is not interrupted easily, and security is required even in tenant buildings for measurement of common areas. In addition, it is possible to perform measurements with less burden, such as being relatively easy to tolerate.
Here, for verification, simultaneous measurement of all layers is performed, and the effectiveness of this method is verified by using data in different time zones.
建物概要
対象建物は、昭和50年以前の建築で地上9階、地下2階、軒高31m のSRC造の事務所建物である。基準階の平面形状を測定点と併せて図7に示す。対象建物では、耐震補強工事が実施されている。また、建物北側(吹き抜け部分)には立体駐車場が併設されている。
Building Overview The target building is an SRC office building with 9 floors above ground, 2 floors below ground, and an eaves height of 31 m, which was built before 1975. The planar shape of the reference floor is shown in FIG. 7 together with the measurement points. Seismic retrofitting work is being carried out at the target building. In addition, a multi-storey car park is located on the north side of the building (atrium).
常時微動測定
全層のコア部分のXY方向(CX,CY)で計測するケースと、代表階(3、7、9階)において水平展開するケースに分けて行った。センサーを配置した各ケースの測定点配置を図7に示す。図7(a)では全層コアに一カ所、水平配置した図7(b)では5カ所である。
1回の計測は30分間とした。センサーは東京測振製のサーボ型速度計(VSE-12)を用いた。
Constant fine movement measurement The case was divided into the case where the core part of all layers was measured in the XY direction (CX, CY) and the case where it was horizontally deployed on the representative floors (3rd, 7th, and 9th floors). FIG. 7 shows the arrangement of measurement points in each case in which the sensor is arranged. In FIG. 7 (a), there is one location on the core of all layers, and in FIG. 7 (b), which is arranged horizontally, there are five locations.
One measurement was for 30 minutes. The sensor used was a servo-type speedometer (VSE-12) manufactured by Tokyo Seisakusho.
<建物の振動特性の分析>
伝達関数の掛け合わせ(全層コア部分)
測定した30分間のデータを9等分し、3分20秒ごとにそれぞれ層毎の伝達関数を求め(図10)、違う時間帯の伝達関数の掛け合わせ(図2の式)で、1階に対する各階の伝達関数を求めた。Y(短辺)方向の結果について示す。掛け合わせの組み合わせを変えた9通りの伝達関数と、同時計測データにより直接1階に対する各階の伝達関数を求めたものを重ねて図11に示す。掛け合わせによる伝達関数のばらつきは若干であり、同時計測の場合とほぼ同じ形状をしている。ピーク付近の形状もよく合っており、固有振動数は1次モード1.5Hz、2次モード4.5Hzであった。
<Analysis of vibration characteristics of buildings>
Multiplication of transfer functions (all-layer core part)
Divide the measured 30-minute data into 9 equal parts, obtain the transfer function for each layer every 3 minutes and 20 seconds (Fig. 10), and multiply the transfer functions in different time zones (Equation in Fig. 2) on the first floor. The transfer function of each floor was obtained. The result in the Y (short side) direction is shown. FIG. 11
伝達関数の掛け合わせ(水平展開)
センサーを水平展開させた3、7、9階において、建物端部の測点の、同じ階の基準測点に対する水平伝達関数(9階の場合9NY/9CY,9SY/9CY)(図4)と、前節で求めた1階に対する上下伝達関数(9CY/1CY)とを掛け合わせ、建物端部の1階に対する伝達関数(9NY/1CY, 9SY/1CY)を求めた。9階について、掛け合わせによる伝達関数を同時計測での伝達関数と重ねて図5に示す。同時計測の場合と同様の形状が得られ、3.0Hz付近でねじれのモードが確認できる。
Multiplication of transfer functions (horizontal expansion)
On the 3rd, 7th, and 9th floors where the sensor is horizontally deployed, the horizontal transfer function of the station at the end of the building with respect to the reference station on the same floor (9NY / 9CY, 9SY / 9CY for the 9th floor) (Fig. 4). , The transfer function (9NY / 1CY, 9SY / 1CY) for the first floor at the end of the building was obtained by multiplying it with the vertical transfer function (9CY / 1CY) for the first floor obtained in the previous section. For the 9th floor, the transfer function by multiplication is shown in FIG. 5 overlaid with the transfer function in simultaneous measurement. The same shape as in the case of simultaneous measurement can be obtained, and the twist mode can be confirmed around 3.0 Hz.
モード形状
求めた掛け合わせパターンを変えた9通りの伝達関数について、1次固有振動数(1.5Hz)と2次固有振動数(4.5Hz)におけるピーク値を拾い、振幅が最大の階を1と基準化して振動モードとして求めた。同時計測の場合と比較して図6に示す。階による揺れの大きさの大小関係は正確に把握できることがわかる。
Mode shape For 9 transfer functions with different multiplication patterns, pick up the peak values at the primary natural frequency (1.5Hz) and the secondary natural frequency (4.5Hz), and set the order with the maximum amplitude to 1. It was standardized and obtained as the vibration mode. FIG. 6 shows a comparison with the case of simultaneous measurement. It can be seen that the magnitude relationship of the magnitude of the shaking depending on the floor can be accurately grasped.
常時微動の部分移動測定を想定して、部分毎の伝達関数の掛け合わせにより、1階に対する建物全体の伝達関数を算出した。1回あたり3分余りの短いデータを用いているが、掛け合わせによる伝達関数は同時計測によるものとほぼ同じ形状となり、部分移動測定で精度よく全体挙動が把握できる。 Assuming the partial movement measurement of constant fine movement, the transfer function of the entire building for the first floor was calculated by multiplying the transfer functions for each part. Although short data of about 3 minutes is used for each time, the transfer function by multiplication has almost the same shape as that by simultaneous measurement, and the overall behavior can be grasped accurately by partial movement measurement.
<既存建物の常時微動測定>
1.1 対象建物
対象建物は、地上9階、地下2階、搭屋3階、軒高31mのSRC造の事務所建物である。基準階の平面形状を測定点と併せて図7に示す。基準階平面形状は、右側に多少凹凸があるが、概形としては34.7m×24.5mの長方形である。昭和50年より前に建築され、柱補強などの耐震補強工事が実施されている。
<Constant tremor measurement of existing buildings>
1.1 Target building The target building is an SRC office building with 9 floors above ground, 2 floors below ground, 3 floors of the building, and an eaves height of 31 m. The planar shape of the reference floor is shown in FIG. 7 together with the measurement points. The standard floor plan shape is a rectangle of 34.7m x 24.5m, although there are some irregularities on the right side. It was built before 1975 and is undergoing seismic retrofitting work such as pillar reinforcement.
1.2 測定条件
検証に用いるため、長時間の同時測定を行い、4.1(1)に後述する方法で、2点の部分移動測定を想定し検証を行った。測定ケースは、階段室において1階から搭屋1階(R階)の全階のXY方向を測定するケース(全層コア測定)と、代表階(7、9階)において、端部を含めた複数個所を測定するケース(水平展開測定)に分けて実施した。各ケースの測定点配置を図7に示す。
前面道路は片道2車線で交通量は比較的多いが、測定当日の風は強くなく、常時微動が定常に近いと考えられる条件で実施している。
1回の計測は30分間とし、サンプリング振動数は100Hzとした。センサーは東京測振製のサーボ型速度計(VSE-12)を用いた。
1.2 Measurement conditions In order to use for verification, simultaneous measurement was performed for a long time, and verification was performed assuming partial movement measurement of two points by the method described later in 4.1 (1). The measurement cases include the case of measuring the XY direction of all floors from the 1st floor to the 1st floor (R floor) of the staircase (all-layer core measurement) and the representative floors (7th and 9th floors) including the ends. It was divided into cases (horizontal expansion measurement) in which multiple points were measured. The arrangement of measurement points in each case is shown in FIG.
The front road has two lanes each way and the traffic volume is relatively heavy, but the wind on the day of measurement is not strong, and it is carried out under conditions where tremors are considered to be close to steady at all times.
One measurement was for 30 minutes, and the sampling frequency was 100 Hz. The sensor used was a servo-type speedometer (VSE-12) manufactured by Tokyo Seisakusho.
3.同時測定データによる振動特性の分析
得られた各測定点の時刻歴データをフーリエ変換し、1階に対する各階のフーリエスペクトルの振幅比(以降、伝達関数とよぶ)を求めた。フーリエスペクトルの算出では、時刻歴データを40.96秒ずつに分割して平均化処理をし、0.1HzのParzenウィンドウを施した。
端部を測定している9階と、中間階である5階の伝達関数を図8に示す。
固有振動数は、X方向は、1次モード2.2Hz、2次モード7.6Hz、Y方向は、1次モード1.4Hz、2次モード4.5Hzであった。また、ねじれ(θ)1次モードは、端部で振幅の大きい3.2Hzであった。
なお、参考として、1次モードの減衰定数を、伝達関数から1/√2法で求めたところ、X1次モードで7.9%、Y1次モードで4.4%であった。
3. 3. Analysis of vibration characteristics using simultaneous measurement data The time history data of each measurement point obtained was Fourier transformed to obtain the amplitude ratio of the Fourier spectrum of each floor to the first floor (hereinafter referred to as the transfer function). In the calculation of the Fourier spectrum, the time history data was divided into 40.96 seconds and averaged, and a 0.1 Hz Parzen window was applied.
The transfer function of the 9th floor measuring the end and the 5th floor, which is the middle floor, is shown in FIG.
The natural frequency was 2.2 Hz in the primary mode in the X direction, 7.6 Hz in the secondary mode, and 1.4 Hz in the primary mode and 4.5 Hz in the secondary mode in the Y direction. The twist (θ) primary mode was 3.2 Hz, which had a large amplitude at the end.
As a reference, when the attenuation constant of the first-order mode was obtained from the transfer function by the 1 / √2 method, it was 7.9% in the X-first-order mode and 4.4% in the Y-first-order mode.
<既存建物の高さ方向のモード情報の把握例>
4.1 高さ方向の振動モード把握の検証
(1)手法検証のためのデータの取り扱い
本測定は、30分間の同時測定を行っているが、本発明では部分移動測定は、各2点の測定を、すべて異なる時間に実施しても良い。同時測定のデータを利用して部分移動測定の検証を行うため、高さ方向の振動モード把握の検証では、同時測定のデータを時刻歴上で分割してデータを取り扱う。30分間(1800秒)のデータを時刻歴上で9等分して3分20秒(200秒)ずつのデータとし、それぞれ層毎の伝達関数を求め、伝達関数の掛け合わせを行う場合は、すべてが異なる時間の掛け合わせとする。
掛け合わせのパターン(どの時間帯にどの層を測定したかの想定)を、図9に示すように階を順番に降りていく形となる9通りとした。実際に部分移動測定を実施する場合と比較すると、30分間という短い時間内での測定ではあるが、これら9通りの掛け合わせを同時測定の場合と比較する。
<Example of grasping mode information in the height direction of an existing building>
4.1 Verification of grasping vibration mode in the height direction (1) Handling of data for method verification Although this measurement is performed simultaneously for 30 minutes, in the present invention, partial movement measurement is performed at two points each. The measurements may all be performed at different times. In order to verify the partial movement measurement using the data of the simultaneous measurement, in the verification of grasping the vibration mode in the height direction, the data of the simultaneous measurement is divided in the time history and the data is handled. When the data of 30 minutes (1800 seconds) is divided into 9 equal parts in the time history to obtain the data of 3 minutes and 20 seconds (200 seconds) each, the transfer function for each layer is obtained, and the transfer functions are multiplied. Everything is a cross of different times.
As shown in FIG. 9, the crossing pattern (assuming which layer was measured at which time zone) was set to 9 patterns in which the floors were descended in order. Compared with the case where the partial movement measurement is actually performed, the measurement is performed within a short time of 30 minutes, but these nine combinations are compared with the case of the simultaneous measurement.
(2)伝達関数と固有振動数
全層コア測定を9等分した時刻歴データを、それぞれフーリエ変換し、層毎の伝達関数を求めた。フーリエスぺクトルの算出では、40.96秒ごとに分割して平均化処理し、0.1HzのParzenウィンドウを施している。
部分移動測定による層毎の伝達関数を、X方向を図10(a)に、Y方向を図11(a)にそれぞれ9本重ねて示す。
層毎の伝達関数はローカルな特性を表し、建物全体の卓越振動数はあらわれていない。
これら層毎の伝達関数を、図9に示すようにして9パターンを想定して、図2で表す式により掛け合わせ、1階に対する各階の伝達関数を求めた。
このようにして部分移動測定の伝達関数の掛け合わせにより求めた1階に対する各階の伝達関数と、同時測定データにより求めた対応する伝達関数を重ねて図10(b)と図11(b) に示す。部分移動測定による伝達関数のばらつきは若干であり、同時測定の場合とほぼ同じ形状をしている。ピーク付近の形状もよく合っており、固有振動数は、X1次モードで2.2Hz、X2次モードで7.6Hz、Y1次モードで1.4Hz、Y2次モードで4.5Hzと同じ値が得られている。参考として、部分移動測定による伝達関数から1/√2法で1次モードの減衰定数を求めたところ、X1次モードで7.3〜8.4%、Y1次モードで4.0〜4.9%であった。
(2) Transfer function and natural frequency The time history data obtained by dividing the full-layer core measurement into nine equal parts was Fourier-transformed to obtain the transfer function for each layer. In the calculation of the Fourier spectrum, the Parzen window of 0.1 Hz is applied by dividing and averaging every 40.96 seconds.
The transfer functions for each layer by the partial movement measurement are shown in FIG. 10 (a) in the X direction and in FIG. 11 (a) in the Y direction.
The transfer function for each layer represents a local characteristic, and the predominant frequency of the entire building does not appear.
As shown in FIG. 9, the transfer functions for each layer were multiplied by the equation shown in FIG. 2, assuming 9 patterns, and the transfer function for each floor with respect to the 1st floor was obtained.
In this way, the transfer function of each floor for the first floor obtained by multiplying the transfer functions of the partial movement measurement and the corresponding transfer function obtained by the simultaneous measurement data are superimposed and shown in FIGS. 10 (b) and 11 (b). show. The variation of the transfer function due to the partial movement measurement is slight, and the shape is almost the same as that of the simultaneous measurement. The shape near the peak also matches well, and the natural frequency is the same as 2.2Hz in X1st mode, 7.6Hz in X2nd mode, 1.4Hz in Y1st mode, and 4.5Hz in Y2rd mode. .. As a reference, when the attenuation constant of the first-order mode was obtained by the 1 / √2 method from the transfer function by partial movement measurement, it was 7.3 to 8.4% in the X-first-order mode and 4.0 to 4.9% in the Y-first-order mode.
(3)振動モード
部分移動測定による掛け合わせで求めた9通りの1階に対する各階の伝達関数において、1次固有振動数(X方向2.2Hz、Y方向1.4Hz)と2次固有振動数(X方向7.6Hz,Y方向4.5Hz)におけるピーク値を拾い、振動モード形状を求めた。同時測定の伝達関数からピーク値を拾ったものと比較して図12に示す。階による揺れの大きさの大小関係は正確に把握できた。
(3) Vibration mode The first-order natural frequency (2.2Hz in the X direction, 1.4Hz in the Y direction) and the second-order natural frequency (X) in the transmission functions of each floor for the nine first floors obtained by multiplying by partial movement measurement. The peak value in the direction (7.6Hz in the direction and 4.5Hz in the Y direction) was picked up to determine the vibration mode shape. It is shown in FIG. 12 in comparison with the peak value picked up from the transfer function of the simultaneous measurement. I was able to accurately grasp the magnitude relationship of the magnitude of the shaking depending on the floor.
<既存建物の平面的な振動モードの把握>
4.2 平面的な振動モード把握の検証
(1)回転スペクトル
水平展開測定を実施した7階と9階で、同時に測定した2点のX方向とY方向の記録を用いて、回転スペクトルを求める。回転スペクトルを求める評価点の位置を図13に●と■で示す。測定点CXとCYのデータから中央部評価点CC、測定点CXとC’Yのデータから中央部評価点CC’の位置の回転スペクトルを求める。また、測定点WXとSY のデータから評価点WSの位置の回転スペクトルを求めるなど、端部の測定点(WX,EX,NY,SY)のデータより隅部評価点(WS,WN,ES,EN)の回転スペクトルを求める。
回転スペクトルの算定では、まず、X方向の波形とY方向の波形の合成により、1度毎の角度に射影させ、0度から180度までの180本の波形を求めた。それらをフーリエ変換し、横軸を振動数、縦軸を振動方向(角度)として、振幅をコンター図で表現した回転スペクトルを図14に示す。
0度および180度はX方向、90度はY方向を示す。1.4Hzでみられるピークは、どの評価点においても、90度にピークがあることから、Y方向の並進モードであることがわかる。一方、2.2Hzと3.2Hzのピークは、測点により、ピークとなる振動方向が異なり、ねじれを伴う振動をしていることがわかる。
図14は、コンター図の性格上、ピークとなる振動方向と振幅を詳細に読み取ることが難しいため、X1次(2.2Hz)、Y1次(1.4Hz)、θ1次(3.2Hz)の各モードの振動数での断面をとり、図15に横軸を振動方向、縦軸をフーリエ振幅として示す。図15(a)(b)より、1.4HzのY1次モードでは、どの評価点においても90度付近にピークがあり、同程度の振幅で振動している。
一方、ねじれを伴う、例えば図15(e)に示す9階の3.2Hzのモードでは、評価点CCでは20度方向で最も振動が大きいのに対して、評価点CC’では150度方向で最も振動が大きい。9階のこのモードにおける回転中心は、図16に示すように、評価点を通る最も振動が大きい方向に直交する線の交点として、幾何的に求めることができる。
<Understanding the planar vibration mode of an existing building>
4.2 Verification of grasping the plane vibration mode (1) Rotational spectrum The rotation spectrum is obtained by using the records of two points measured at the same time in the X and Y directions on the 7th and 9th floors where the horizontal expansion measurement was performed. .. The positions of the evaluation points for obtaining the rotation spectrum are shown by ● and ■ in FIG. The rotation spectrum of the position of the central evaluation point CC is obtained from the data of the measurement points CX and CY, and the position of the central evaluation point CC'is obtained from the data of the measurement points CX and C'Y. In addition, the rotation spectrum of the position of the evaluation point WS is obtained from the data of the measurement points WX and SY, and the corner evaluation points (WS, WN, ES, Obtain the rotation spectrum of EN).
In the calculation of the rotation spectrum, first, by synthesizing the waveform in the X direction and the waveform in the Y direction, the waveform was projected at an angle of 1 degree, and 180 waveforms from 0 degree to 180 degrees were obtained. FIG. 14 shows a rotation spectrum obtained by Fourier transforming them, with the horizontal axis representing the frequency and the vertical axis representing the vibration direction (angle), and expressing the amplitude in a contour diagram.
0 degrees and 180 degrees indicate the X direction, and 90 degrees indicate the Y direction. Since the peak seen at 1.4 Hz has a peak at 90 degrees at any evaluation point, it can be seen that it is a translation mode in the Y direction. On the other hand, it can be seen that the peaks of 2.2 Hz and 3.2 Hz vibrate with twisting because the peak vibration direction differs depending on the station.
Since it is difficult to read the peak vibration direction and amplitude in detail in FIG. 14 due to the nature of the contour diagram, the X1 (2.2 Hz), Y1 (1.4 Hz), and θ1 (3.2 Hz) modes are shown. A cross section in terms of frequency is taken, and FIG. 15 shows the horizontal axis as the vibration direction and the vertical axis as the Fourier amplitude. From FIGS. 15 (a) and 15 (b), in the 1.4 Hz Y primary mode, there is a peak near 90 degrees at every evaluation point, and the vibration is performed with the same amplitude.
On the other hand, in the 3.2 Hz mode on the 9th floor shown in FIG. 15 (e), which involves twisting, the vibration is the largest in the 20-degree direction at the evaluation point CC, whereas it is the largest in the 150-degree direction at the evaluation point CC'. The vibration is large. As shown in FIG. 16, the center of rotation of the ninth floor in this mode can be geometrically obtained as the intersection of the lines passing through the evaluation point and orthogonal to the direction of the largest vibration.
(2)回転中心とねじれ振動モード
4.2(1)で述べたようにして、ねじれを伴うモード(2.2Hz、3.2Hz)について、回転中心を求めた。7階と9階における中央部評価点(CC,CC’)と隅部評価点(WS,WN,ES,EN)より、それぞれ求めた回転中心を、図17と図18に●と○で示す。
回転スペクトルを評価する2点と回転中心との位置関係において、2点と回転中心を結ぶ線が、平行に近い場合(角度が小さい場合)、交点に誤差が生じやすい。剛床を仮定しているので、対象とするモードのフーリエ振幅の比(振幅比)は、回転中心と2点との距離の比(回転半径比)に比例すると考えられることから、振幅比/回転半径比を求め、これが1に近いものを確からしい回転中心として求めた。図17には、振幅比/回転半径比が0.8〜1.2のものをプロットしている。中央部評価点より求めた回転中心では、振幅比/回転半径比の値が、7階および9階のそれぞれ2つのモードについて、0.9〜1.1であり確からしいと考えられる。中央部評価点および隅部評価点より求めた回転中心の位置は、図17に示す、回転中心が建物平面の外にあるX1次モード(2.2Hz)では若干ばらつきがあるが、図18に示す回転中心が建物平面内にあるθ1次モード(3.2Hz)ではよく合っている。
中央部評価点より求めた回転中心によるねじれモード形状と、同時測定による端部の測定点の伝達関数より求めたねじれモード形状を比較して図18に実線と点線で示している。回転中心によるねじれモード形状は、同時測定の伝達関数から求めたねじれモードとほぼ同じ形状となっている。
(2) Rotation center and torsional vibration mode As described in 4.2 (1), the rotation center was determined for the modes with twist (2.2Hz, 3.2Hz). The centers of rotation obtained from the central evaluation points (CC, CC') and corner evaluation points (WS, WN, ES, EN) on the 7th and 9th floors are shown by ● and ○ in FIGS. 17 and 18, respectively. ..
In the positional relationship between the two points for evaluating the rotation spectrum and the center of rotation, when the lines connecting the two points and the center of rotation are close to parallel (when the angle is small), an error is likely to occur at the intersection. Since a rigid bed is assumed, the ratio of Fourier amplitudes in the target mode (amplitude ratio) is considered to be proportional to the ratio of the distance between the center of rotation and the distance between the two points (radial rotation ratio). The turning radius ratio was found, and the one close to 1 was found as the probable center of rotation. FIG. 17 plots the amplitude ratio / turning radius ratio of 0.8 to 1.2. At the center of rotation obtained from the central evaluation point, the value of the amplitude ratio / radius of gyration ratio is 0.9 to 1.1 for each of the two modes of the 7th and 9th floors, which is considered to be probable. The positions of the rotation centers obtained from the central evaluation points and the corner evaluation points are slightly different in the X-first mode (2.2Hz) in which the rotation center is outside the building plane, as shown in FIG. 17, but are shown in FIG. It fits well in the θ1 primary mode (3.2Hz) where the center of rotation is in the building plane.
The twist mode shape obtained from the evaluation point at the center and the twist mode shape obtained from the transfer function of the measurement points at the ends by simultaneous measurement are compared and shown by a solid line and a dotted line in FIG. The twist mode shape due to the center of rotation is almost the same as the twist mode obtained from the transfer function of simultaneous measurement.
測定点の設置に制約を受けることの多い、既存建物の常時微動測定において、本実施例では簡易な測定で建物全体挙動を把握するために、全体の多点同時測定に代わる方法として、2点ずつの測定を移動させながら繰り返す部分移動測定方法につい有効性が確認された。
9階建て既存建物の測定のデータをもとに、1回あたりの計測時間が3分余りのデータで、各モードの固有振動数と振動モード形状を求め、同時測定と同様の振動性状が得られることを検証し、次の知見が得られた。
(1)層毎の伝達関数を掛け合わせることにより各階の1階に対する伝達関数と、1次および2次の高さ方向のモード形状を、精度よく把握できた。
(2)平面上2点の回転スペクトルを利用して、ねじれの回転中心とモード形状を求めることができ、同時測定による伝達関数から求めたモード形状とよく整合している。
(3)対象建物においては、階段室付近の同じ階の1スパン(7m程度)のみ離れた、2点3成分からねじれモードの回転中心が求まり、階段室まわりのみの測定で、ねじれモード形状を把握できる。
In the constant fine movement measurement of an existing building, which is often restricted by the installation of measurement points, in this embodiment, in order to grasp the behavior of the entire building by simple measurement, two points are used as an alternative method to simultaneous multipoint measurement of the entire building. The effectiveness of the partial movement measurement method, in which each measurement is repeated while moving, was confirmed.
Based on the measurement data of the existing 9-story building, the natural frequency and vibration mode shape of each mode can be obtained from the data of about 3 minutes per measurement time, and the same vibration properties as the simultaneous measurement can be obtained. The following findings were obtained by verifying that.
(1) By multiplying the transfer function for each layer, the transfer function for the first floor of each floor and the mode shapes in the primary and secondary height directions could be grasped accurately.
(2) The rotation center of the twist and the mode shape can be obtained by using the rotation spectra of two points on the plane, and the mode shape is well matched with the mode shape obtained from the transfer function by simultaneous measurement.
(3) In the target building, the rotation center of the twist mode can be obtained from the two points and three components separated by one span (about 7 m) on the same floor near the staircase, and the twist mode shape can be obtained by measuring only around the staircase. I can grasp it.
<地震の影響度判定方法>
建物について地震応答性能が決定され、さらに、この建物の限られた階にセンサーを配置して実際に暴露した地震データを観測データとして地震時に全層に渡る地震応答(影響度、ダメージ、被災度)を推定することができる。被災時に建物の各箇所を実際に確認することなく、直後地震対策を取ることができる手法であり、システムを実現する。
<Earthquake impact determination method>
Seismic response performance is determined for the building, and sensors are placed on the limited floors of this building, and the actual exposed seismic data is used as observation data to cover all layers of seismic response (impact, damage, damage). ) Can be estimated. It is a method that can take earthquake countermeasures immediately after the disaster without actually checking each part of the building, and realizes the system.
(限られた観測階による地震時の全層応答推定手法の開発)
N階建て、応答観測階がS箇所(地動観測点1点、応答観測点S点)の建物について、全N階の応答を推定する方法を示す。
1〜S次の固有モードΦは、振動測定、設計モデル等により、あらかじめ設定しておくものとする。ここで、対象とするモード次数はSとし、応答観測点数と一致させる。
観測データに基づく全層応答推定算出法を図19に示す。
(Development of full-thickness response estimation method at the time of earthquake with limited observation floor)
The method of estimating the response of all N floors for a building with N floors and response observation floors at S points (1 ground motion observation point, response observation point S point) is shown.
The 1st to Sth specific modes Φ shall be set in advance by vibration measurement, design model, etc. Here, the target mode order is S, which matches the number of response observation points.
FIG. 19 shows a full-layer response estimation calculation method based on the observation data.
本手法は理論上、全観測階について観測値と推定値が一致するため、防災システムとして利用者の混乱を招かず実用性が高いと考えられる。 Theoretically, this method has the same observed and estimated values for all observation floors, so it is considered to be highly practical as a disaster prevention system without causing confusion for users.
建物に応用した地震応答判定システム例を図20に示す。
建物の4カ所に地震センサーが配置されている。被災時にこれらのセンサーから得られる情報を演算処置装置に入力して建物の任意箇所の被災度を推定して表示部に出力する。
演算処理装置では、この建物に関してあらかじめ常時微振動などを利用して高さ方向と水平方向のモード情報が取得されている。このモード情報に地震センサーからの観測データを適用して、モード重畳法による応答演算を行い、任意箇所の被災度を推定する。
FIG. 20 shows an example of an earthquake response determination system applied to a building.
Seismic sensors are placed at four locations in the building. The information obtained from these sensors at the time of a disaster is input to the arithmetic treatment device, the degree of damage at any part of the building is estimated, and the information is output to the display unit.
In the arithmetic processing unit, mode information in the height direction and the horizontal direction is acquired in advance for this building by constantly using micro-vibration or the like. The observation data from the seismic sensor is applied to this mode information, and the response calculation by the mode superimposition method is performed to estimate the degree of damage at any location.
本実施例に示すように、常時微振動を利用して、上下階毎の測定と各層の任意の2箇所を測定をすることにより、全層に渡る高さ方向のモード情報とねじれモードによる平面的なモード情報が正確に得られることとなる。したがって、設計情報が無い建物や、竣工当初から改修などを経て地震応答が変化している建物についても、全層に渡って同じ箇所にセンサーを配置して、一斉に測定せずとも、高さ方向とねじれ情報が分かることとなる。
高さ方向のモードとねじれモードが決定できることにより、これに、想定される地震情報を入力することにより、建物の任意地点の地震応答を推定することができる。
そして、この任意地点に存在する柱、梁、壁、階段などの建築部材、あるいは、設備などに与える影響も推定できることとなり、それぞれの固有の強度などの物性を反映して、それぞれの建築部材や設備などの建物に関する構造物の被災度を判定することができる。
すなわち、得られた建物の地震応答推定にモード重畳法を適用して判定することができる。
地震応答に基づいて、構造物の被災度を判定するシステムは特開2016−109607号公報、特開2016−197013号公報、特許第6001740号公報に開示されるような公知のシステムを利用することができる。
As shown in this embodiment, by constantly using micro-vibration to measure each floor above and below and at any two points of each layer, the mode information in the height direction over all layers and the plane by the twist mode Mode information can be obtained accurately. Therefore, even for buildings for which there is no design information or for which the seismic response has changed due to repairs from the beginning of completion, sensors are placed at the same location over all floors, and the height does not need to be measured all at once. The direction and twist information will be known.
Since the mode in the height direction and the twist mode can be determined, the seismic response at an arbitrary point of the building can be estimated by inputting the assumed seismic information into the mode.
Then, it is possible to estimate the influence on building members such as columns, beams, walls, stairs, etc. existing at this arbitrary point, or equipment, etc., reflecting the physical properties such as the unique strength of each building member and It is possible to determine the degree of damage to structures related to buildings such as equipment.
That is, it can be determined by applying the mode superimposition method to the seismic response estimation of the obtained building.
As a system for determining the degree of damage to a structure based on an earthquake response, use a known system as disclosed in JP-A-2016-109607, JP-A-2016-97013, and Patent No. 6001740. Can be done.
Claims (4)
高さ方向のモード情報は、上下2点間を一層ずつずらして測定し、層毎に伝達関数を掛け合わせて建物上下全体の変形を推定して得られるモード情報であることを特徴とする建物の任意箇所の地震応答を推定する方法。 It is a method of estimating the seismic response of any part of the building based on the limited vertical mode information in the vertical direction and the horizontal twist mode information that measured the constant micro-vibration or vibration applied to the building.
Mode information in the height direction is between the upper and lower two points were measured by shifting one by one, you characterized in that the mode information obtained by multiplying the transfer function for each layer to estimate the deformation of the entire building vertical how to estimate the seismic response of any part of the building.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2017139146A JP6916682B2 (en) | 2017-07-18 | 2017-07-18 | How to estimate the seismic response of any part of the building |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2017139146A JP6916682B2 (en) | 2017-07-18 | 2017-07-18 | How to estimate the seismic response of any part of the building |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2019020261A JP2019020261A (en) | 2019-02-07 |
| JP6916682B2 true JP6916682B2 (en) | 2021-08-11 |
Family
ID=65355668
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2017139146A Active JP6916682B2 (en) | 2017-07-18 | 2017-07-18 | How to estimate the seismic response of any part of the building |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP6916682B2 (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP7276611B2 (en) * | 2020-06-19 | 2023-05-18 | 三菱電機株式会社 | shaking estimation system |
| JP7482001B2 (en) * | 2020-11-09 | 2024-05-13 | 株式会社東芝 | Plant evaluation system, method, and program |
| JP7623235B2 (en) * | 2021-06-25 | 2025-01-28 | 清水建設株式会社 | Building response estimation system and building response estimation method |
| JP2023159679A (en) * | 2022-04-20 | 2023-11-01 | 清水建設株式会社 | Soundness evaluation system and soundness evaluation method |
| JP7716698B1 (en) | 2024-04-26 | 2025-08-01 | 国立大学法人広島大学 | Earthquake risk assessment device, earthquake risk assessment method and program |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3876247B2 (en) * | 2000-12-28 | 2007-01-31 | 構造品質保証研究所株式会社 | Diagnosis method and diagnosis system of structure by microtremor observation |
| JP2016109607A (en) * | 2014-12-09 | 2016-06-20 | 凸版印刷株式会社 | Strong motion seismograph, measuring system, and damage state determination method |
| US10648881B2 (en) * | 2015-03-23 | 2020-05-12 | 9420410 Canada Inc. | Seismic response assessment of man-made structures |
-
2017
- 2017-07-18 JP JP2017139146A patent/JP6916682B2/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| JP2019020261A (en) | 2019-02-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6916682B2 (en) | How to estimate the seismic response of any part of the building | |
| JP3952851B2 (en) | Seismic performance evaluation method and apparatus for buildings | |
| JP5569900B2 (en) | Seismic performance evaluation method, seismic performance evaluation device, and seismic performance evaluation system | |
| JP3105450B2 (en) | Method and apparatus for determining risk of earthquake damage to structures | |
| CN105445756B (en) | Safety Diagnosis System for Structures | |
| WO2014109324A1 (en) | Building safety verification system and building safety verification method | |
| JP6001740B1 (en) | High precision evaluation of structure transfer functions, earthquake response prediction, deterioration diagnosis system and method | |
| JP2011095237A5 (en) | ||
| CN105928675B (en) | System and method for performing modal analysis on a structure | |
| JP2014211397A (en) | Method for confirming earthquake response and soundness of building | |
| JP2016018389A (en) | Comprehensive monitoring device, and comprehensive monitoring program | |
| JP7343380B2 (en) | Building health monitoring system | |
| JP3956302B2 (en) | Damaged part detection device for structure and method for detecting damaged part of structure | |
| JP5799183B2 (en) | Building safety verification system, building safety verification method and program | |
| Tuhta et al. | System identification of rc building using n4sid | |
| Zhang et al. | Ambient vibration testing & modal identification of an office building | |
| JP4822337B2 (en) | Diagnosis method for building structures | |
| Mahmoudabadi et al. | Identification of modal parameters of non‐classically damped linear structures under multi‐component earthquake loading | |
| JP6363539B2 (en) | Estimation method of damaged part of building | |
| JP6800049B2 (en) | Building vibration measurement method | |
| JP6438745B2 (en) | Estimation method for inter-building displacement | |
| RU2292433C1 (en) | Method to determine change in deflected mode of building structure or structure having complicated three-dimensional shape | |
| Balafas et al. | Extension of the rotation algorithm for earthquake damage estimation of complex structures | |
| JP3285787B2 (en) | Countermeasure effect judgment method and countermeasure effect judgment device | |
| Rendon | Determination Of The Modal Parameters Of A Five Story Reinforced Concrete Structure Using Ultra-Low Level Excitation And Computational Analysis |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A80 | Written request to apply exceptions to lack of novelty of invention |
Free format text: JAPANESE INTERMEDIATE CODE: A80 Effective date: 20170815 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20200623 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20210310 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20210316 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20210428 Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20210428 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20210629 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20210716 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6916682 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |