Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP7590908B2 - Floor vibration analysis system - Google Patents
[go: Go Back, main page]

JP7590908B2 - Floor vibration analysis system - Google Patents

Floor vibration analysis system Download PDF

Info

Publication number
JP7590908B2
JP7590908B2 JP2021063145A JP2021063145A JP7590908B2 JP 7590908 B2 JP7590908 B2 JP 7590908B2 JP 2021063145 A JP2021063145 A JP 2021063145A JP 2021063145 A JP2021063145 A JP 2021063145A JP 7590908 B2 JP7590908 B2 JP 7590908B2
Authority
JP
Japan
Prior art keywords
vibration
floor
partition wall
acceleration
evaluation point
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
Application number
JP2021063145A
Other languages
Japanese (ja)
Other versions
JP2022158338A (en
Inventor
安紀子 更谷
剛 平松
充教 足立
伸介 上大門
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tokkyokiki Corp
Daiwa House Industry Co Ltd
Original Assignee
Tokkyokiki Corp
Daiwa House Industry Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tokkyokiki Corp, Daiwa House Industry Co Ltd filed Critical Tokkyokiki Corp
Priority to JP2021063145A priority Critical patent/JP7590908B2/en
Publication of JP2022158338A publication Critical patent/JP2022158338A/en
Application granted granted Critical
Publication of JP7590908B2 publication Critical patent/JP7590908B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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 relates to a floor vibration analysis system that analyzes floor vibrations.

従来から、床に外力が加わった際の振動を予測する床振動解析システムが提案されている。例えば、特許文献1では、建物概要データ、通りデータ、梁データ、壁データ、床データ、柱データなどの各種データなどの入力データから、レイリーリッツ法などの振動解析の手法を用いて、床振動を予測するシステムが提案されている。 Conventionally, floor vibration analysis systems have been proposed that predict vibrations when an external force is applied to a floor. For example, Patent Document 1 proposes a system that predicts floor vibrations using vibration analysis techniques such as the Rayleigh-Ritz method from input data such as various types of data, including building outline data, street data, beam data, wall data, floor data, and column data.

特開2017-182408号公報JP 2017-182408 A

しかしながら、特許文献1のシステムで、床振動を予測する際に、各種データとして、各部材の寸法以外にも、これを構成する材料の物性値等、様々なデータを入力しなければならず、解析に膨大な時間を要する。さらに、床振動解析は動的解析であるので、実測した結果と一致しないことも多い。これに加えて、建物の躯体を構成しない非構造部位(二次部材)としての間仕切り壁は、床振動の解析には反映されていないが、実際のところ、間仕切り壁の有無、間仕切り壁により、実測した振動の結果は、大きく変わることがある。 However, when predicting floor vibrations with the system of Patent Document 1, various data must be input, including not only the dimensions of each component, but also the physical properties of the materials that make them up, and the analysis takes an enormous amount of time. Furthermore, because floor vibration analysis is a dynamic analysis, it often does not match the results of actual measurements. In addition, partition walls, which are non-structural parts (secondary components) that do not form part of the building's framework, are not reflected in the floor vibration analysis, but in reality, the results of actually measured vibrations can vary greatly depending on whether or not there are partition walls, and the presence or absence of partition walls.

本発明はこのような点を鑑みてなされたものであり、その目的とするところは、より簡単な入力条件で、間仕切り壁の影響を考慮した、建物の床振動を解析することができる床振動解析システムを提供することにある。 The present invention was made in consideration of these points, and its purpose is to provide a floor vibration analysis system that can analyze the floor vibration of a building, taking into account the effects of partition walls, with simpler input conditions.

本発明に係る床振動解析システムは、鉄骨造の建物の大梁に囲まれた空間に小梁が配置された梁構造と、梁構造に支持された矩形状の床スラブとを含む、平面視が矩形状の評価領域において、前記評価領域に存在する振動評価点の振動を解析する床振動解析システムであって、前記床振動解析システムは、間仕切り壁が配置された前記評価領域の前記振動評価点において、時間変化を伴う加振力が外力として作用したときに、前記時間変化を伴う前記振動評価点の上下方向の加速度を算出するものであり、前記梁構造に作用する前記評価領域の前記床スラブの重量から、前記振動評価点の振動における質量を算出する質量算出部と、前記梁構造の情報と前記床スラブの情報とから、前記評価領域における床の固有振動数を算出し、前記固有振動数に基づいて、前記振動評価点の振動におけるばね定数を算出するばね定数算出部と、少なくとも前記間仕切り壁の構の情報から、前記振動評価点の振動における減衰定数を算出する減衰定数算出部と、前記時間変化を伴う加振力、前記ばね定数、前記質量、および前記減衰定数に基づいて、前記床スラブおよび前記梁構造の1自由度系の振動として、時間変化を伴う前記振動評価点の加速度を算出する加速度算出部と、を備えることを特徴とする。 The floor vibration analysis system according to the present invention is a floor vibration analysis system that analyzes the vibration of a vibration evaluation point present in an evaluation area having a rectangular shape in plan view, the evaluation area including a beam structure in which small beams are arranged in a space surrounded by the main beams of a steel-framed building, and a rectangular floor slab supported by the beam structure. The floor vibration analysis system calculates the vertical acceleration of the time-varying vibration evaluation point when a time-varying excitation force acts as an external force at the vibration evaluation point of the evaluation area where a partition wall is arranged, and calculates the vertical acceleration of the time-varying vibration evaluation point from the weight of the floor slab of the evaluation area acting on the beam structure. The system is characterized by having a mass calculation unit that calculates the mass at the floor slab and the beam structure, a spring constant calculation unit that calculates the natural frequency of the floor in the evaluation area from the beam structure information and the floor slab information, and calculates the spring constant at the vibration of the vibration evaluation point based on the natural frequency, a damping constant calculation unit that calculates the damping constant at the vibration of the vibration evaluation point from at least the partition wall structure information, and an acceleration calculation unit that calculates the acceleration of the vibration evaluation point with time change as the vibration of a one-degree-of-freedom system of the floor slab and the beam structure based on the time-varying excitation force, the spring constant, the mass, and the damping constant.

発明者らは、床振動を解析する際に、間仕切り壁の有無により、床振動の解析結果がことなることが分かった。さらに、調査および検討を進める中で、間仕切り壁は、床振動の減衰への寄与率が高いことがわかった。これまでは、柱、外壁等の構造も、床振動解析において、考慮されてきたが、したがって、本発明によれば、床スラブと梁構造をバネとし、床スラブの重量を質点とし、間仕切り壁等をダンパとしたバネ・マス・ダンパ系の簡易的なモデルの振動として、前記床スラブおよび前記梁構造の1自由度系の振動における振動評価点の振動解析を行うことができる。これにより、間仕切り壁による床振動の影響を簡単に把握することができる。 The inventors found that when analyzing floor vibration, the analysis results differ depending on whether or not a partition wall is present. Furthermore, in the course of their investigation and study, they found that partition walls have a high contribution rate to attenuation of floor vibration. Until now, structures such as columns and exterior walls have also been taken into consideration in floor vibration analysis, but according to the present invention, it is possible to perform vibration analysis of vibration evaluation points in the vibration of a one-degree-of-freedom system of the floor slab and beam structure as a simple model of a spring-mass-damper system in which the floor slab and beam structure are springs, the weight of the floor slab is a mass point, and partition walls, etc. are dampers. This makes it easy to grasp the effect of floor vibration caused by partition walls.

より好ましい態様としては、前記間仕切り壁の情報は、間仕切り壁の長さと、前記振動評価点から前記間仕切り壁までの距離とであり、前記減衰定数算出部は、少なくとも前記間仕切り壁の長さおよび前記振動評価点から前記間仕切り壁までの距離に基づいて、前記減衰定数を算出する。 In a more preferred embodiment, the information about the partition wall is the length of the partition wall and the distance from the vibration evaluation point to the partition wall, and the damping constant calculation unit calculates the damping constant based on at least the length of the partition wall and the distance from the vibration evaluation point to the partition wall.

この態様によれば、振動評価点の振動の減衰特性は、間仕切り壁の長さと振動評価点までの距離とに大きく依存することから、これらの値の基づいた重回帰式や対応表等から、減衰定数を算出することで、より正確に床振動の解析を行うことができる。 According to this aspect, since the vibration damping characteristics of the vibration evaluation point are highly dependent on the length of the partition wall and the distance to the vibration evaluation point, it is possible to more accurately analyze floor vibration by calculating the damping constant from a multiple regression equation or correspondence table based on these values.

より好ましい態様としては、前記間仕切り壁は、複数あり、前記間仕切り壁の情報は、前記各間仕切り壁の長さと、前記振動評価点から前記各間仕切り壁までの距離と、前記複数の間仕切り壁の配置状態とであり、前記減衰定数算出部は、少なくとも前記各間仕切り壁の長さ、および前記振動評価点から前記各間仕切り壁までの距離に基づいて、前記減衰定数を算出し、前記複数の間仕切り壁の配置状態に対応した数値で、前記減衰定数を補正する。 In a more preferred embodiment, there are multiple partition walls, the partition wall information includes the length of each partition wall, the distance from the vibration evaluation point to each partition wall, and the arrangement of the multiple partition walls, and the damping constant calculation unit calculates the damping constant based on at least the length of each partition wall and the distance from the vibration evaluation point to each partition wall, and corrects the damping constant with a numerical value corresponding to the arrangement of the multiple partition walls.

ここで、振動評価点の振動の減衰特性は、間仕切り壁の長さと振動評価点までの距離とに大きく依存することから、これらの値の基づいた重回帰式や、間仕切り壁のパラメータ等と減衰定数とを対応付けた対応テーブル等を含むデータから、減衰特性を算出することで、より正確に床振動の解析を行うことができる。特に、間仕切り壁が複数存在する場合には、それらの間仕切り壁の配置状態が、減衰特性に大きく寄与することから、これらの配置状態に対応付けられた補正係数等で補正を行うことにより、より正確に床振動の解析を行うことができる。間仕切り壁の長さおよび配置状態(レイアウト)の最適化を図ることができる。 Here, since the vibration damping characteristics of the vibration evaluation point depend greatly on the length of the partition wall and the distance to the vibration evaluation point, floor vibration can be analyzed more accurately by calculating the damping characteristics from data including a multiple regression equation based on these values and a correspondence table that associates the parameters of the partition wall with damping constants. In particular, when there are multiple partition walls, the layout of the partition walls greatly contributes to the damping characteristics, so floor vibration can be analyzed more accurately by making corrections using correction coefficients that correspond to these layouts. The length and layout of the partition walls can be optimized.

より好ましい態様としては、前記床振動解析システムは、前記加速度算出部で算出した前記振動評価点の加速度の波形を周波数応答解析することにより、振動周波数ごとの加速度を算出する解析部を備える。この態様によれば、解析部により、加速度の波形を解析することにより、振動の知覚確率を算出することができる。 In a more preferred embodiment, the floor vibration analysis system includes an analysis unit that calculates the acceleration for each vibration frequency by performing frequency response analysis on the acceleration waveform of the vibration evaluation point calculated by the acceleration calculation unit. According to this embodiment, the analysis unit can calculate the perception probability of vibration by analyzing the acceleration waveform.

より好ましい態様としては、前記建物は、前記床スラブの上に、制振装置が配置された構造であり、前記評価領域には、前記制振装置を含み、前記加速度算出部は、前記制振装置と、前記床スラブおよび前記梁構造との2自由度系の振動における前記加速度を算出する。この態様によれば、制振装置が無い場合の、床スラブおよび梁構造の1自由度系の振動解析と、制振装置が有る場合の、制振装置と、前記床スラブおよび前記梁構造との2自由度系の振動解析とを行うことにより、制振装置の有無による床振動の影響を確認することができる。 In a more preferred embodiment, the building has a structure in which a vibration control device is placed on the floor slab, the evaluation area includes the vibration control device, and the acceleration calculation unit calculates the acceleration in vibration of a two-degree-of-freedom system consisting of the vibration control device, the floor slab, and the beam structure. According to this embodiment, the effect of floor vibration due to the presence or absence of a vibration control device can be confirmed by performing a vibration analysis of a one-degree-of-freedom system of the floor slab and beam structure when there is no vibration control device, and a vibration analysis of a two-degree-of-freedom system of the vibration control device, the floor slab, and the beam structure when there is a vibration control device.

本発明によれば、より簡単な入力条件で、間仕切り壁の影響を考慮した、対象建物の床振動を予測することができる。 According to the present invention, it is possible to predict floor vibrations of a target building, taking into account the effects of partition walls, with simpler input conditions.

本発明の床振動解析システムの装置構成図である。1 is a diagram showing the configuration of a floor vibration analysis system according to the present invention; 図1に示す演算装置のブロック図である。FIG. 2 is a block diagram of the arithmetic unit shown in FIG. 1 . 図1に示す床振動解析システムが行う解析モデルの一例である。2 is an example of an analysis model performed by the floor vibration analysis system shown in FIG. 1 . (a)は、図2に示すばね定数算出部の床スラブの有効幅を説明するための模式図であり、(b)は、評価領域を構成する床全体の断面二次モーメントの算出を説明するための模式図である。3A is a schematic diagram for explaining the effective width of the floor slab of the spring constant calculation section shown in FIG. 2, and FIG. 3B is a schematic diagram for explaining the calculation of the second moment of area of the entire floor constituting the evaluation area. (a)および(b)は、複数の間仕切り壁の配置状態を説明するための解析モデルの一例である。13(a) and (b) are examples of analytical models for explaining the arrangement of multiple partition walls. 図2に示す加速度算出部による加速度の算出を説明するための振動方程式の1自由度系の振動モデルである。3 is a vibration model of a one-degree-of-freedom system of a vibration equation for explaining the calculation of acceleration by the acceleration calculation unit shown in FIG. 2 . (a)は、評価点における加振力の波形であり、(b)は、評価点における加速度の波形である。4A is the waveform of the excitation force at the evaluation point, and FIG. 4B is the waveform of the acceleration at the evaluation point. 図2に示す解析部による振動周波数ごとの加速度を算出した結果を示したグラフである。3 is a graph showing the results of calculation of acceleration for each vibration frequency by the analysis unit shown in FIG. 2 . 変形例に係る図2に示す加速度算出部による加速度の算出を説明するための振動方程式の2自由度系の振動モデルである。3 is a vibration model of a two-degree-of-freedom system of a vibration equation for explaining calculation of acceleration by the acceleration calculation unit shown in FIG. 2 according to a modified example.

以下、本発明の本実施形態に係る床振動解析システム1を、図面を参照しながら説明する。 The floor vibration analysis system 1 according to this embodiment of the present invention will be described below with reference to the drawings.

〔第1実施形態〕
1.床振動解析システム1の装置構成について
本実施形態に係る床振動解析システム1は、図1に示すように、対象建物の床の振動評価点における床振動を解析するものである。以下に、床振動解析システム1を説明する。なお、振動評価点は、床振動を評価したい、床スラブ上の1点のことであり、床スラブの上に後述する制振装置が設けられている場合も、同様に床スラブ上の1点である。本実施形態で、床振動解析を行う建物の構造形式は、鉄骨造である。
First Embodiment
1. Device configuration of floor vibration analysis system 1 As shown in Fig. 1, the floor vibration analysis system 1 according to this embodiment analyzes floor vibration at a vibration evaluation point on a floor of a target building. The floor vibration analysis system 1 will be described below. Note that the vibration evaluation point is a point on a floor slab where floor vibration is to be evaluated, and even if a vibration control device (described later) is provided on the floor slab, it is also a point on the floor slab. In this embodiment, the structural type of the building where floor vibration analysis is performed is a steel frame structure.

本実施形態で、解析するモデルは、図3に示すような鉄骨造の建物5の大梁52に囲まれた空間に小梁54が配置された梁構造55と、梁構造55に支持された矩形状の床スラブ51とを含む、平面視が矩形状の評価領域56において、その評価領域56に存在する振動評価点Cの振動を床解析するためのモデルである。 In this embodiment, the model to be analyzed is a model for floor analysis of vibrations at a vibration evaluation point C that exists in an evaluation area 56 that is rectangular in plan view, and includes a beam structure 55 in which minor beams 54 are arranged in a space surrounded by major beams 52 of a steel-framed building 5 as shown in Figure 3, and a rectangular floor slab 51 supported by the beam structure 55.

ここで、振動評価点Cは、平面視において、評価領域56の中央である。大梁52は、第1方向L1に沿った大梁52と、第1方向A1に直交する第2方向A2に沿った大梁52とがある。矩形状の評価領域56の周縁は、第1方向A1と第2方向A2に沿った縁部で構成される。本実施形態では、小梁54は、3本設けられており、各小梁54は、第2方向A2に沿って、大梁52.52の間にわたされており、3本の小梁54は、第1方向A1に等間隔に設けられている。 Here, the vibration evaluation point C is the center of the evaluation area 56 in a plan view. The main beams 52 include a main beam 52 aligned along the first direction L1 and a main beam 52 aligned along the second direction A2 perpendicular to the first direction A1. The periphery of the rectangular evaluation area 56 is composed of edges along the first direction A1 and the second direction A2. In this embodiment, three small beams 54 are provided, and each small beam 54 is disposed between the main beams 52, 52 along the second direction A2, and the three small beams 54 are provided at equal intervals in the first direction A1.

評価領域56は、平面視において、長さL1の長辺と、長さL2の短辺とで構成される矩形状の評価領域であり、評価領域56の四隅には、柱57が設けられているが、本実施形態では、柱57を考慮せずに解析を行う。これは、なお、床スラブ51の上に、制振装置を設けた場合には、この評価領域56は、制振装置も含むものである。本実施形態では、床スラブ51の上に、間仕切り壁60Aが配置されている。 In plan view, the evaluation area 56 is a rectangular evaluation area consisting of a long side of length L1 and a short side of length L2. Columns 57 are provided at the four corners of the evaluation area 56, but in this embodiment, the analysis is performed without taking the columns 57 into consideration. Note that if a vibration control device is provided on the floor slab 51, this evaluation area 56 also includes the vibration control device. In this embodiment, a partition wall 60A is provided on the floor slab 51.

床振動解析システム1は、入力装置2、演算装置3、および表示装置4を少なくとも備えている。入力装置2は、演算装置3に解析条件などを入力するキーボードなどの装置である。演算装置3は、床振動解析プログラムが記録されたRAMまたはROMなどのメモリーと、床振動解析プログラムを実行するCPUと、を備えている。表示装置4は、演算装置3で演算された結果を出力するディスプレイである。 The floor vibration analysis system 1 comprises at least an input device 2, a calculation device 3, and a display device 4. The input device 2 is a device such as a keyboard for inputting analysis conditions and the like to the calculation device 3. The calculation device 3 comprises a memory such as a RAM or ROM in which a floor vibration analysis program is recorded, and a CPU that executes the floor vibration analysis program. The display device 4 is a display that outputs the results calculated by the calculation device 3.

なお、本実施形態では、入力装置2と表示装置4とを個別に設けたが、例えば、これらの機能が一体化したタッチパネルディスプレイを備えているもよい。また、入力装置2および表示装置4と、演算装置3とが、ネットワークを介して接続されていてもよい。 In this embodiment, the input device 2 and the display device 4 are provided separately, but for example, a touch panel display that integrates these functions may be provided. Also, the input device 2 and the display device 4 may be connected to the calculation device 3 via a network.

床振動解析システム1の演算装置3は、間仕切り壁60Aが配置された評価領域56の振動評価点Cにおいて、時間変化を伴う加振力F(t)が外力として作用したときに、時間変化を伴う振動評価点Cの上下方向の加速度A(t)を算出するものである。 The calculation device 3 of the floor vibration analysis system 1 calculates the vertical acceleration A(t) of the time-varying vibration evaluation point C when a time-varying excitation force F(t) acts as an external force at the vibration evaluation point C of the evaluation area 56 where the partition wall 60A is located.

具体的には、このような算出を行うべく、図2に示すように、演算装置3は、ソフトウエアとして、質量算出部31、ばね定数算出部32、減衰定数算出部33、加速度算出部34、および解析部35を主に備えている。 Specifically, to perform such calculations, as shown in FIG. 2, the calculation device 3 mainly comprises, as software, a mass calculation unit 31, a spring constant calculation unit 32, a damping constant calculation unit 33, an acceleration calculation unit 34, and an analysis unit 35.

2.質量算出部31について
梁構造55に作用する評価領域56の床スラブ51の重量から、振動評価点Cの振動における質量を算出する。たとえば、図4に示すように、小梁54に作用する床スラブ51の質量mは、床スラブ51の全体の質量に対して、0.25~0.3倍程度であるので、この範囲における値を、床スラブ51の全体の質量に乗じたものを質量mとして算出する。質量算出部31で算出した質量mは、加速度算出部34に入力される。
2. Mass calculation unit 31 The mass in the vibration of vibration evaluation point C is calculated from the weight of floor slab 51 in evaluation area 56 acting on beam structure 55. For example, as shown in FIG. 4, the mass m of floor slab 51 acting on small beam 54 is about 0.25 to 0.3 times the total mass of floor slab 51, so a value in this range is multiplied by the total mass of floor slab 51 to calculate mass m. The mass m calculated by mass calculation unit 31 is input to acceleration calculation unit 34.

3.ばね定数算出部32について
ばね定数算出部32は、梁構造55の情報と床スラブ51の情報とから、大梁52に囲まれた評価領域56内の固有振動数ωを算出し、固有振動数ωに基づいて、振動評価点Cの振動におけるばね定数kを算出する。具体的には、以下に示す式(1)を用いて、以下の質量算出部31で算出した有効質量mと、固有振動数(一次固有振動数)ωに基づいて、ばね定数kを算出する。
3. Regarding the spring constant calculation unit 32 The spring constant calculation unit 32 calculates the natural frequency ω in the evaluation area 56 surrounded by the girder 52 from the information on the beam structure 55 and the information on the floor slab 51, and calculates the spring constant k in the vibration of the vibration evaluation point C based on the natural frequency ω. Specifically, the spring constant k is calculated based on the effective mass m calculated by the mass calculation unit 31 below and the natural frequency (first-order natural frequency) ω, using the following formula (1).

Figure 0007590908000001
Figure 0007590908000001

ここで、式(1)に示す固有振動数ωは、以下に示す式(2)を用いて算出することができる。式(2)に示す、αは、算出対象となる大梁52または小梁54(以下「梁」という)の支持条件であり、単純支持梁の場合には、α=1であり、固定支持梁の場合には、α=1.5を用いる。 The natural frequency ω shown in formula (1) can be calculated using formula (2) below. In formula (2), α is the support condition of the main beam 52 or sub-beam 54 (hereinafter referred to as "beam") being calculated. In the case of a simply supported beam, α = 1, and in the case of a fixed supported beam, α = 1.5 is used.

Lは、算出対象となる梁の長さであり、Eは、その梁のヤング率である。Iは、その梁とこれが支持する床スラブ51を含む断面二次モーメント(完全合成梁の有効等価断面二次モーメント)である。ρは、梁の単位長さあたりの質量である。したがって、大梁52および小梁54の固有振動数ωを算出し、この中の最小の固有振動数ωを、大梁52に囲まれた評価領域56内の床の固有振動数(スパン全体の一次固有振動数)とみなす。 L is the length of the beam to be calculated, and E is the Young's modulus of that beam. I is the second moment of area including the beam and the floor slab 51 it supports (effective equivalent second moment of area of a completely composite beam). ρ is the mass per unit length of the beam. Therefore, the natural frequencies ω of the main girder 52 and the secondary beam 54 are calculated, and the smallest natural frequency ω among these is regarded as the natural frequency of the floor within the evaluation area 56 surrounded by the main girder 52 (the primary natural frequency of the entire span).

Figure 0007590908000002
Figure 0007590908000002

ここで、式(2)に示す完全合成梁の有効等価断面二次モーメントIは、以下に示す式(3)を用いて算出することができる。式(3)に示す、Bは、床スラブ51の有効幅であり(図4(a)参照)、tは、床スラブ51の厚みであり、Xは、床スラブ51の表面から中立軸までの長さであり、Iは、算出対象となる梁(大梁52または小梁54)の二次モーメントであり、Aは、その梁(大梁52または小梁54)の断面である。なお、図4(a)では、説明の便宜上、H形鋼である大梁52と小梁54の大きさを同じ大きさで描いているが、実際はこれらの大きさは異なる。 Here, the effective equivalent moment of inertia I of the complete composite beam shown in formula (2) can be calculated using formula (3) shown below. In formula (3), B is the effective width of the floor slab 51 (see FIG. 4(a)), t is the thickness of the floor slab 51, Xn is the length from the surface of the floor slab 51 to the neutral axis, I s is the moment of inertia of the beam (girder 52 or sub-girder 54) to be calculated, and A s is the cross section of the beam (girder 52 or sub-girder 54). In FIG. 4(a), the girder 52 and sub-girder 54, which are H-shaped steel, are drawn to have the same size for the sake of convenience, but in reality, their sizes are different.

さらに、式(3)に示す、dは、床スラブ51の表面から、梁(大梁52または小梁54)の重心までに長さである。なお、デッキプレートを用いた床スラブでは、デッキプレート山上の平板状のコンクリートの厚さを床スラブの有効厚さとし、その長さは、たとえば、80mm以上とする。 Furthermore, ds in formula (3) is the length from the surface of the floor slab 51 to the center of gravity of the beam (main beam 52 or sub-beam 54). In floor slabs using deck plates, the thickness of the flat concrete on the top of the deck plate is regarded as the effective thickness of the floor slab, and this length is, for example, 80 mm or more.

Figure 0007590908000003
Figure 0007590908000003

ここで、合成梁の弾性剛性算出に用いる、床スラブ51の有効幅Bは、図4(a)に示すように、梁の上フランジの幅bに、上フランジの両側または片側の床スラブ協力幅baを加えたものである。 Here, the effective width B of the floor slab 51 used in calculating the elastic rigidity of the composite beam is the width b of the upper flange of the beam plus the effective width b a of the floor slab on both sides or one side of the upper flange, as shown in Figure 4 (a).

床スラブ協力幅baは、以下の式(4-1)~式(4-4)のいずれか1つを用いて算出することができる。式(4-1)と式(4-2)は、固定支持梁(すなわちラーメン材および連続梁)の場合に用いる式であり、式(4-3)と式(4-4)は、単純支持梁の場合に用いる式である。 The floor slab joint width b a can be calculated using one of the following formulas (4-1) to (4-4). Formulas (4-1) and (4-2) are used in the case of fixed support beams (i.e. rigid frame members and continuous beams), while formulas (4-3) and (4-4) are used in the case of simply supported beams.

ここで、aは、図4(a)に示すように、並列梁では鉄骨上フランジ端部から相隣る梁の上フランジ端部までの距離であり、単独T形梁または片持ちスラブよりなる梁では、その片側の上フランジ端面から床スラブ端までの距離の2倍である。lは、固定支持梁のスパンの長さであり、lは、単純支持梁のスパンの長さである。 Here, a is the distance from the end of the steel frame top flange to the end of the adjacent top flange in the case of parallel beams, as shown in Figure 4(a), and in the case of a single T-beam or a beam consisting of a cantilever slab, it is twice the distance from the end of the top flange on one side to the end of the floor slab. l is the span length of the fixed support beam, and l 0 is the span length of the simply supported beam.

Figure 0007590908000004
Figure 0007590908000004

ここで、xおよびd-xは、以下の式(5)および式(6)によって算出することができる。なお、中立軸が床スラブ外にある。ここで、式(5)に示す、tは、t/dであり、Pは、A/(B・d)であり、nは、ヤング係数比であり、15である。 Here, xn and ds - xn can be calculated by the following formulas (5) and (6). Note that the neutral axis is outside the floor slab. Here, in formula (5), t1 is t/ ds , Pt is As / (B· ds ), and n is the Young's modulus ratio, which is 15.

Figure 0007590908000005
Figure 0007590908000005

なお、上の式(2)~式(6)を用いて固有振動数ωを算出したが、例えば、同じ方向に並設される小梁の本数が2本以下である場合には、以下に示す式(7)を用いてもよい。ここで、Lは、算出対象となる梁と直交する方向の評価領域56(直交スラブ)の長さである。 Although the natural frequency ω was calculated using the above formulas (2) to (6), for example, when the number of small beams arranged in parallel in the same direction is two or less, the following formula (7) may be used, where Ln is the length of the evaluation area 56 (orthogonal slab) in the direction perpendicular to the beam to be calculated.

Figure 0007590908000006
Figure 0007590908000006

このようにして、ばね定数算出部32では、大梁52のサイズ(具体的には、長さ、梁幅、梁成など)、小梁54のサイズ(具体的には、長さ、梁幅、梁成など)、小梁54の本数を含む、上の式で算出に要する梁構造55の情報と、床スラブ51の質量、厚さ、および必要に応じてその仕様等の含む床スラブ51の情報を入力し、ばね定数kを算出することができる。なお、小梁54の本数は、図4(a)に示すように、並列梁では鉄骨上フランジ端部から相隣る梁の上フランジ端部までの距離aを算出する際に、簡易的に用いることができる。 In this way, the spring constant calculation unit 32 can calculate the spring constant k by inputting information on the beam structure 55 required for calculation using the above formula, including the size of the main girder 52 (specifically, length, beam width, beam depth, etc.), the size of the secondary beams 54 (specifically, length, beam width, beam depth, etc.), and the number of secondary beams 54, as well as information on the floor slab 51, including the mass, thickness, and, if necessary, specifications of the floor slab 51. Note that the number of secondary beams 54 can be used simply when calculating the distance a from the end of the steel frame upper flange to the end of the upper flange of the adjacent beam in parallel beams, as shown in Figure 4(a).

4.減衰定数算出部33について
減衰定数算出部33は、少なくとも間仕切り壁60Aの情報から、振動評価点Cの振動における減衰定数hを算出する。減衰定数hを算出する際には、大梁52の長さの情報(すなわち、スパン)をさらに用いてもよい。本実施形態では、減衰係数hは、床スラブ51および梁構造55に、間仕切り壁60Aの情報をさらに加えて算出されることが好ましい。すなわち、算出される減衰係数hは、さらに間仕切り壁60Aによる減衰特性と、床スラブ51および梁構造55の減衰特性との影響より算出したものであることが好ましい。なお、1自由度系の振動の場合には、後述する加速度算出部34において、算出した減衰定数hから、後述する式(8)に示す減衰係数cを算出する。
4. Regarding the damping constant calculation unit 33 The damping constant calculation unit 33 calculates the damping constant h in the vibration of the vibration evaluation point C from at least the information of the partition wall 60A. When calculating the damping constant h, information on the length of the main girder 52 (i.e., the span) may be further used. In this embodiment, it is preferable that the damping coefficient h is calculated by further adding information on the partition wall 60A to the floor slab 51 and the beam structure 55. That is, it is preferable that the calculated damping coefficient h is calculated from the influence of the damping characteristics of the partition wall 60A and the damping characteristics of the floor slab 51 and the beam structure 55. In the case of vibration of a one-degree-of-freedom system, the acceleration calculation unit 34 described later calculates the damping coefficient c shown in Equation (8) described later from the calculated damping constant h.

ここで、間仕切り壁60Aの情報は、間仕切り壁60Aの長さG1と、振動評価点Cから間仕切り壁60A(たとえば中心g1)までの距離LG1とである。具体的には、減衰定数算出部33は、間仕切り壁60Aの長さG1および振動評価点Cから間仕切り壁60Aまでの距離LG1に基づいて、減衰定数hを算出する。 Here, the information on the partition wall 60A is the length G1 of the partition wall 60A and the distance LG1 from the vibration evaluation point C to the partition wall 60A (for example, the center g1). Specifically, the damping constant calculation unit 33 calculates the damping constant h based on the length G1 of the partition wall 60A and the distance LG1 from the vibration evaluation point C to the partition wall 60A.

振動評価点Cの振動の減衰特性は、間仕切り壁60Aの長さG1と振動評価点Cまでの距離とに大きく依存すことが、発明者の調査によりわかった。実際の建物において測定した減衰定数hと、長さG1および距離LG1とから、重回帰分析を行い、これらを変数とした回帰式や、長さG1および距離LG1と、減衰定数hとを対応付けた対応表から、減衰定数hを算出することができる。これにより、正確に床振動の解析を行うことができる。なお、この回帰式または対応表に、上述したスパンを変数として加えてもよい。さらに、これらの回帰式または対応表に、床スラブ51の物性値、短辺または長辺の長さ、梁構造55における大梁の長さ、または、スパン等のパラメータを変数として加えてもよい。 The inventor's research has revealed that the vibration damping characteristics of vibration evaluation point C are highly dependent on the length G1 of the partition wall 60A and the distance to the vibration evaluation point C. A multiple regression analysis is performed from the damping constant h measured in an actual building, the length G1, and the distance LG1, and the damping constant h can be calculated from a regression equation using these as variables or a correspondence table that associates the length G1 and the distance LG1 with the damping constant h. This allows floor vibration to be analyzed accurately. The above-mentioned span may be added as a variable to this regression equation or correspondence table. Furthermore, parameters such as the physical properties of the floor slab 51, the length of the short side or long side, the length of the main girder in the beam structure 55, or the span may be added as variables to these regression equations or correspondence tables.

さらに、図5(a)および(b)に示すように、間仕切り壁60A、60Bなど、間仕切り壁が複数ある場合、間仕切り壁の情報は、各間仕切り壁60A、60Bの長さG1、G2と、振動評価点Cから各間仕切り壁60A、60Bまでの距離LG1、LG2と、複数の間仕切り壁60A、60Bの配置状態とである。図5(a)では、2つの間仕切り壁60A、60Bが平行となるように配置されており、図5(b)では、2つの間仕切り壁60A、60Bが直交となるように配置されている。 Furthermore, as shown in Figures 5(a) and (b), when there are multiple partition walls such as partition walls 60A and 60B, the partition wall information includes the lengths G1 and G2 of each partition wall 60A and 60B, the distances LG1 and LG2 from the vibration evaluation point C to each partition wall 60A and 60B, and the arrangement of the multiple partition walls 60A and 60B. In Figure 5(a), the two partition walls 60A and 60B are arranged so as to be parallel, and in Figure 5(b), the two partition walls 60A and 60B are arranged so as to be perpendicular to each other.

減衰定数算出部33は、各間仕切り壁60A、60Bの長さG1、G2、および振動評価点Cから前記各間仕切り壁60A、60Bまでの距離LG1、LG2に基づいて、減衰定数hを算出する。ここで、実際の建物において測定した減衰定数hと、長さG1、G2および距離LG1、LG2とから、重回帰分析を行い、これらを変数とした重回帰式や、長さG1、G2および距離LG1、LG2と、減衰定数hとを対応付けた対応表から、減衰定数hを算出することができる。なお、重回帰式、対応表では、上に示したように、長さG1、G2をそれぞれの変数として取り扱わず、これらの総長さ(総和の値)を変数としてもよい。また、距離LG1、LG2も同様に総距離(総和の値)を変数としてもよい。さらに、上述した場合と同様に、これらの回帰式または対応表に、床スラブ51の物性値、短辺または長辺の長さ、梁構造55における大梁の長さ、または、スパン等のパラメータを変数として加えてもよい。 The damping constant calculation unit 33 calculates the damping constant h based on the lengths G1, G2 of the partition walls 60A, 60B and the distances LG1, LG2 from the vibration evaluation point C to the partition walls 60A, 60B. Here, a multiple regression analysis is performed from the damping constant h measured in an actual building, the lengths G1, G2, and the distances LG1, LG2, and the damping constant h can be calculated from a multiple regression equation using these as variables or a correspondence table that associates the lengths G1, G2, the distances LG1, LG2, and the damping constant h. Note that in the multiple regression equation and the correspondence table, as shown above, the lengths G1, G2 may not be treated as individual variables, but the total length (sum value) of these may be used as a variable. Similarly, the total distance (sum value) of the distances LG1, LG2 may be used as a variable. Furthermore, as in the above case, parameters such as the physical properties of the floor slab 51, the length of the short or long sides, the length of the girder in the beam structure 55, or the span may be added as variables to these regression equations or correspondence tables.

図5(a)(b)に示すように、間仕切り壁60A、60Bの配置状態が変われば、減衰定数hも変わる。したがって、これらの配置状態に対応した数値を補正係数として設定し、算出した減衰定数hを補正係数で乗じることにより、減衰定数hを補正してもよい。なお、この補正は、間仕切り壁が2個の場合に限られず、間仕切り壁の個数とこれらの配置状態に応じて、設定すればよい。 As shown in Figures 5(a) and (b), if the arrangement of the partition walls 60A and 60B changes, the damping constant h also changes. Therefore, a numerical value corresponding to these arrangements may be set as a correction coefficient, and the damping constant h may be corrected by multiplying the calculated damping constant h by the correction coefficient. Note that this correction is not limited to the case where there are two partition walls, and may be set according to the number of partition walls and their arrangement.

このように、振動評価点Cの振動の減衰特性は、間仕切り壁60A、60Bの長さと振動評価点Cまでの距離とに大きく依存することから、これらの値の基づいた重回帰式や、間仕切り壁のパラメータ等と減衰定数とを対応付けた対応テーブル等を含むデータから、減衰定数hを算出することで、より正確に床振動の解析を行うことができる。 As such, since the vibration damping characteristics of vibration evaluation point C are highly dependent on the length of partition walls 60A, 60B and the distance to vibration evaluation point C, floor vibration can be analyzed more accurately by calculating the damping constant h from data including a multiple regression equation based on these values and a correspondence table that associates the parameters of the partition walls with damping constants.

特に、間仕切り壁60A、60Bが複数存在する場合には、それら間仕切り壁60A、60Bの配置状態が、減衰特性に大きく寄与することから、これらの配置状態に対応付けられた補正係数等で補正を行うことにより、より正確に床振動の解析を行うことができる。これにより、間仕切り壁60A、60Bの長さG1、G2および配置状態(レイアウト)の最適化を図ることができる。 In particular, when there are multiple partition walls 60A, 60B, the arrangement of these partition walls 60A, 60B contributes greatly to the damping characteristics, so by making corrections using correction coefficients and the like associated with these arrangements, floor vibrations can be analyzed more accurately. This allows optimization of the lengths G1, G2 and arrangement (layout) of the partition walls 60A, 60B.

5.加速度算出部34について
次に、加速度算出部は、時間変化を伴う加振力F(T)、ばね定数k、質量m、および減衰係数c(減衰定数hから算出された値)に基づいて、床スラブ51および梁構造55の1自由度系の振動として、時間変化を伴う振動評価点の加速度A(T)を算出する。なお、加速度算出部34において、減衰係数cは、減衰定数算出部33で算出された、減衰定数hから算出される値である。具体的には、図6に示すような、1自由度系の振動モデルを想定し、式(8)に示す振動方程式が得られる。ここで、カッコ内のTは時間であり、(T)は、時間Tの関数であること意味する。
5. Acceleration Calculation Unit 34 Next, the acceleration calculation unit calculates the acceleration A(T) of the vibration evaluation point with time change as the vibration of the one-degree-of-freedom system of the floor slab 51 and the beam structure 55 based on the excitation force F(T) with time change, the spring constant k, the mass m, and the damping coefficient c (value calculated from the damping constant h). In the acceleration calculation unit 34, the damping coefficient c is a value calculated from the damping constant h calculated by the damping constant calculation unit 33. Specifically, a vibration model of the one-degree-of-freedom system as shown in FIG. 6 is assumed, and the vibration equation shown in Equation (8) is obtained. Here, T in parentheses is time, and (T) means that it is a function of time T.

時間変化を伴う加振力F(T)は、床スラブ51の評価領域56に付与される外力の条件である。たとえば評価領域56において、特定の人数の人が、飛び跳ね、歩行、小走り、かかと衝撃動作、またはエアロビクス屈伸運動などを行った際に、これらの条件に対応する数値である。 The time-varying excitation force F(T) is a condition of the external force applied to the evaluation area 56 of the floor slab 51. For example, it is a numerical value that corresponds to the conditions when a certain number of people jump, walk, jog, perform heel strikes, aerobic bending and stretching, etc. in the evaluation area 56.

A(T)は、上述した如く時間変化を伴う振動評価点の加速度であり、V(T)は、時間変化を伴う振動評価点Cの速度であり、X(T)は、時間変化を伴う振動評価点Cの変位である。なお、変位X(T)を時間で微分すれば、速度V(T)となり、速度V(T)を時間で微分すれば、加速度A(T)となる。 As described above, A(T) is the acceleration of the vibration evaluation point that changes with time, V(T) is the velocity of the vibration evaluation point C that changes with time, and X(T) is the displacement of the vibration evaluation point C that changes with time. Note that if the displacement X(T) is differentiated with respect to time, it becomes the velocity V(T), and if the velocity V(T) is differentiated with respect to time, it becomes the acceleration A(T).

Figure 0007590908000007
Figure 0007590908000007

式(8)の微分方程式に基づいて、数値計算により、時間変化を伴う振動評価点Cの加速度A(T)を算出することができる。たとえば、このような手法としては、平均加速度法、線形加速度法、ニューマークβ法、ウィルソンθ法、ルンゲクッタ法、など一般的に知られている方法により、加速度A(T)を算出することができる。 Based on the differential equation of formula (8), the acceleration A(T) of the vibration evaluation point C, which changes over time, can be calculated by numerical calculation. For example, the acceleration A(T) can be calculated by commonly known methods such as the average acceleration method, the linear acceleration method, the Newmark β method, the Wilson θ method, and the Runge-Kutta method.

ここで、ニューマークβ法を簡単に説明すると、解析したい加振力F(T)を、時間間隔ΔTごとの数値データで構成される波形として準備し、m回分の時間間隔ΔTの経過した加振力を、加振力Fとし、その時の加速度をAとし、速度をVとし、変位をXとすると、m+1回分の時間間隔ΔTの経過した、速度をVm+1とし、変位をXm+1を以下の式(9)、式(10)で仮定することができる。なお、式(10)のβが、1/4のときは、平均加速度法であり、1/6のときは、線形加速度法であり、加速度の実測値と対比し、より精度の高い加速度の解析値となるようにβを設定すればよい。 Here, to briefly explain the Newmark β method, the excitation force F(T) to be analyzed is prepared as a waveform composed of numerical data for each time interval ΔT, and if the excitation force after m time intervals ΔT is defined as excitation force Fm , the acceleration at that time is defined as Am , the velocity is defined as Vm , and the displacement is defined as Xm , then the velocity after m+1 time intervals ΔT is defined as Vm +1 and the displacement is defined as Xm +1 can be assumed by the following formulas (9) and (10). Note that when β in formula (10) is 1/4, it is the average acceleration method, and when it is 1/6, it is the linear acceleration method. β can be set so as to obtain an analysis value of acceleration with higher accuracy by comparing it with the actual measured value of acceleration.

Figure 0007590908000008
Figure 0007590908000008

この式(9)および式(10)を、式(8)をm+1回分の時間間隔ΔTだけ経過した振動方程式である式(11)に代入する。 These equations (9) and (10) are substituted into equation (11), which is the vibration equation after a time interval ΔT equivalent to m+1 times of equation (8).

Figure 0007590908000009
Figure 0007590908000009

m+1で整理すると、m+1回分の時間間隔ΔT経過した加速度をAm+1を、先に求めたばね定数k、質量m、1つ前のM回分経過した、加振力F、加速度A、速度V、変位Xを用いて算出することができる。このようにして、ΔTごとの、加速度を算出することができる。 When rearranged in terms of V m+1 , the acceleration A m+1 after the time interval ΔT for m+1 times has elapsed can be calculated using the previously calculated spring constant k, mass m, and the excitation force F m , acceleration A m , velocity V m , and displacement X m for the previous M times. In this way, the acceleration for each ΔT can be calculated.

このようにして、図7(a)に示す時間変化を伴う加振力の波形から、図7(b)に示す、時間変化をともなる加速度の波形を演算することができる。 In this way, the waveform of the acceleration, which also changes with time, shown in Figure 7(b) can be calculated from the waveform of the excitation force, which also changes with time, shown in Figure 7(a).

6.解析部35について
解析部35は、加速度算出部34で算出した振動評価点Cの加速度A(T)の波形を周波数応答解析することにより(FFTにより)、振動周波数ごとの加速度を算出する。これにより、図8に示すように、振動周波数ごとの加速度を算出した結果を示した振動の知覚確率のグラフ(具体的には、制振装置なしのグラフ)を得ることができる。
6. Analysis Unit 35 The analysis unit 35 calculates the acceleration for each vibration frequency by performing frequency response analysis (by FFT) on the waveform of the acceleration A(T) of the vibration evaluation point C calculated by the acceleration calculation unit 34. This makes it possible to obtain a graph of the vibration perception probability (specifically, a graph without a vibration control device) showing the results of calculating the acceleration for each vibration frequency, as shown in FIG.

本実施形態によれば、床スラブ51と梁構造55をバネとし、床スラブ51の重量を質点とし、間仕切り壁60A(60A、60B)をダンパ等としたバネ・マス・ダンパ系の簡易的なモデルの振動として、振動評価点Cの振動解析を行うことができる。これにより、間仕切り壁60A(60A、60B)による床振動の影響を簡単に把握することができる。 According to this embodiment, vibration analysis of vibration evaluation point C can be performed as the vibration of a simple model of a spring-mass-damper system in which the floor slab 51 and beam structure 55 are springs, the weight of the floor slab 51 is a mass point, and the partition wall 60A (60A, 60B) is a damper. This makes it easy to grasp the effect of floor vibration caused by the partition wall 60A (60A, 60B).

たとえば、変形例として、建物は、床スラブ51の上に、制振装置80が配置された構造であり、評価領域56には、制振装置を含む。加速度算出部34は、制振装置80と、床スラブ51および梁構造55との2自由度系の振動における加速度A(T)を算出する。 For example, as a modified example, the building has a structure in which a vibration control device 80 is placed on the floor slab 51, and the evaluation area 56 includes the vibration control device. The acceleration calculation unit 34 calculates the acceleration A(T) in the vibration of the two-degree-of-freedom system consisting of the vibration control device 80, the floor slab 51, and the beam structure 55.

この場合には、図9に示すモデルが想定でき、床スラブ51の質量をm1とし、ばね定数をk1とし、減衰定数をh1とし、これらの値は、上述した方法で梁構造55の情報と床スラブ51の情報とから算出することができる。さらに、制振装置80の質量をm2とし、ばね定数をk2とし、減衰定数をh2とし、これらの値は、装置固有の予め決定された値である。なお、制振装置80の有りによる床振動解析を行うことから、上に示す制振装置80が無い場合と同様に、床スラブ51に加振力F(t)を作用させ、床スラブ51の振動評価点Cにおける加速度A(T)から、制振装置80の有りによる床振動の評価を行う。 In this case, the model shown in FIG. 9 can be assumed, where the mass of the floor slab 51 is m1, the spring constant is k1, and the damping constant is h1. These values can be calculated from the information of the beam structure 55 and the information of the floor slab 51 using the method described above. Furthermore, the mass of the vibration control device 80 is m2, the spring constant is k2, and the damping constant is h2, and these values are predetermined values specific to the device. Note that, since floor vibration analysis is performed with the vibration control device 80, an excitation force F(t) is applied to the floor slab 51, and floor vibration with the vibration control device 80 is evaluated from the acceleration A(T) at the vibration evaluation point C of the floor slab 51, just as in the case without the vibration control device 80 shown above.

具体的には、図9に示す2自由度系の振動モデルから、これらの質量m1、m2、減衰定数h1、h2、ばね定数k1、k2で構成される、質量行列m、減衰行列c、ばね定数行列k(剛性行列)を用いて、式(8)と同様の振動方程式を作成し、上に示した方法と同様の方法で、加速度算出部34で、振動評価点Cの加速度A(T)を算出し、解析部35で、振動周波数ごとの加速度を算出してもよい。たとえば、減衰行列cは、床スラブ51および梁構造55の振動系、制振装置80の振動系の減衰定数h1、h2から算出される。これにより、制振装置80が無い場合の、床スラブ51および梁構造55の1自由度系の振動解析と、制振装置80が有る場合の、制振装置80と、床スラブ51および梁構造55との2自由度系の振動解析とを行うことで、図8に示すように、制振装置80の有無による床振動の影響を確認することができる。 Specifically, from the vibration model of the two-degree-of-freedom system shown in FIG. 9, a vibration equation similar to that of Equation (8) is created using the mass matrix m, damping matrix c, and spring constant matrix k (stiffness matrix) consisting of the masses m1 and m2, damping constants h1 and h2, and spring constants k1 and k2, and the acceleration A(T) of the vibration evaluation point C is calculated in the acceleration calculation unit 34 in the same manner as described above, and the acceleration for each vibration frequency is calculated in the analysis unit 35. For example, the damping matrix c is calculated from the damping constants h1 and h2 of the vibration system of the floor slab 51 and the beam structure 55 and the vibration system of the vibration damping device 80. As a result, by performing a vibration analysis of the one-degree-of-freedom system of the floor slab 51 and the beam structure 55 when the vibration damping device 80 is not present, and a vibration analysis of the two-degree-of-freedom system of the vibration damping device 80, the floor slab 51, and the beam structure 55 when the vibration damping device 80 is present, the influence of floor vibration due to the presence or absence of the vibration damping device 80 can be confirmed as shown in FIG. 8.

以上、本発明の実施形態について詳述したが、本発明は、前記の実施形態に限定されるものではなく、特許請求の範囲に記載された本発明の精神を逸脱しない範囲で、種々の設計変更を行うことができるものである。 Although the embodiments of the present invention have been described in detail above, the present invention is not limited to the above-mentioned embodiments, and various design changes can be made without departing from the spirit of the present invention as described in the claims.

1:床振動解析システム、31:質量算出部、32:ばね定数算出部、33:減衰定数算出部、34:加速度算出部、35:解析部、51:床スラブ、52:大梁、54:小梁、55:梁構造、56:評価領域、60A、60B:間仕切り壁、C:振動評価点

1: Floor vibration analysis system, 31: Mass calculation section, 32: Spring constant calculation section, 33: Damping constant calculation section, 34: Acceleration calculation section, 35: Analysis section, 51: Floor slab, 52: Main beam, 54: Small beam, 55: Beam structure, 56: Evaluation area, 60A, 60B: Partition wall, C: Vibration evaluation point

Claims (5)

鉄骨造の建物の大梁に囲まれた空間に小梁が配置された梁構造と、梁構造に支持された矩形状の床スラブとを含む、平面視が矩形状の評価領域において、前記評価領域に存在する振動評価点の振動を解析する床振動解析システムであって、
前記床振動解析システムは、間仕切り壁が配置された前記評価領域の前記振動評価点において、時間変化を伴う加振力が外力として作用したときに、前記時間変化を伴う前記振動評価点の上下方向の加速度を算出するものであり、
前記梁構造に作用する前記評価領域の前記床スラブの重量から、前記振動評価点の振動における質量を算出する質量算出部と、
前記梁構造の情報と前記床スラブの情報とから、前記評価領域における床の固有振動数を算出し、前記固有振動数に基づいて、前記振動評価点の振動におけるばね定数を算出するばね定数算出部と、
少なくとも前記間仕切り壁の情報から、前記振動評価点の振動における減衰定数を算出する減衰定数算出部と、
前記時間変化を伴う加振力、前記ばね定数、前記質量、および前記減衰定数に基づいて、前記床スラブおよび前記梁構造の1自由度系の振動として、時間変化を伴う前記振動評価点の加速度を算出する加速度算出部と、を備えることを特徴とする床振動解析システム。
A floor vibration analysis system for analyzing vibrations of vibration evaluation points in an evaluation area having a rectangular shape in a plan view, the evaluation area including a beam structure in which a small beam is arranged in a space surrounded by a large beam of a steel-framed building and a rectangular floor slab supported by the beam structure,
the floor vibration analysis system calculates an upward and downward acceleration of the time-varying vibration evaluation point when a time-varying excitation force acts as an external force at the vibration evaluation point of the evaluation area where a partition wall is arranged,
A mass calculation unit that calculates a mass in vibration at the vibration evaluation point from the weight of the floor slab in the evaluation area acting on the beam structure;
A spring constant calculation unit that calculates a natural frequency of the floor in the evaluation area from the beam structure information and the floor slab information, and calculates a spring constant in vibration of the vibration evaluation point based on the natural frequency;
a damping constant calculation unit that calculates a damping constant for vibration at the vibration evaluation point from at least information on the partition wall;
and an acceleration calculation unit that calculates the acceleration of the vibration evaluation point, which changes over time, as a vibration of a one-degree-of-freedom system of the floor slab and the beam structure based on the time-varying excitation force, the spring constant, the mass, and the damping constant.
前記間仕切り壁の情報は、間仕切り壁の長さと、前記振動評価点から前記間仕切り壁までの距離とであり、
前記減衰定数算出部は、少なくとも前記間仕切り壁の長さおよび前記振動評価点から前記間仕切り壁までの距離に基づいて、前記減衰定数を算出することを特徴とする請求項1に記載の床振動解析システム。
The information on the partition wall includes a length of the partition wall and a distance from the vibration evaluation point to the partition wall,
2. The floor vibration analysis system according to claim 1, wherein the damping constant calculation unit calculates the damping constant based on at least a length of the partition wall and a distance from the vibration evaluation point to the partition wall.
前記間仕切り壁は、複数あり、
前記間仕切り壁の情報は、前記各間仕切り壁の長さと、前記振動評価点から前記各間仕切り壁までの距離と、前記複数の間仕切り壁の配置状態とであり、
前記減衰定数算出部は、少なくとも、前記各間仕切り壁の長さ、および前記振動評価点から前記各間仕切り壁までの距離に基づいて、前記減衰定数を算出し、前記複数の間仕切り壁の配置状態に対応した数値で、前記減衰定数を補正することを特徴とする請求項1に記載の床振動解析システム。
The partition wall is provided in a plurality of locations.
The partition wall information includes a length of each of the partition walls, a distance from the vibration evaluation point to each of the partition walls, and an arrangement state of the plurality of partition walls;
The floor vibration analysis system of claim 1, characterized in that the damping constant calculation unit calculates the damping constant based on at least the length of each partition wall and the distance from the vibration evaluation point to each partition wall, and corrects the damping constant with a numerical value corresponding to the arrangement state of the multiple partition walls.
前記床振動解析システムは、前記加速度算出部で算出した前記振動評価点の加速度の波形を周波数応答解析することにより、振動周波数ごとの加速度を算出する解析部を備えることを特徴とする請求項1~3のいずれか一項に記載の床振動解析システム。 The floor vibration analysis system according to any one of claims 1 to 3, characterized in that the floor vibration analysis system includes an analysis unit that calculates the acceleration for each vibration frequency by performing frequency response analysis on the acceleration waveform of the vibration evaluation point calculated by the acceleration calculation unit. 前記建物は、前記床スラブの上に、制振装置が配置された構造であり、前記評価領域には、前記制振装置を含み、
前記加速度算出部は、前記制振装置と、前記床スラブおよび前記梁構造との2自由度系の振動における前記加速度を算出することを特徴とする請求項1~4のいずれか一項に記載の床振動解析システム。
The building has a structure in which a vibration control device is arranged on the floor slab, and the evaluation area includes the vibration control device,
A floor vibration analysis system according to any one of claims 1 to 4, characterized in that the acceleration calculation unit calculates the acceleration in vibration of a two-degree-of-freedom system consisting of the vibration control device, the floor slab and the beam structure.
JP2021063145A 2021-04-01 2021-04-01 Floor vibration analysis system Active JP7590908B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2021063145A JP7590908B2 (en) 2021-04-01 2021-04-01 Floor vibration analysis system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2021063145A JP7590908B2 (en) 2021-04-01 2021-04-01 Floor vibration analysis system

Publications (2)

Publication Number Publication Date
JP2022158338A JP2022158338A (en) 2022-10-17
JP7590908B2 true JP7590908B2 (en) 2024-11-27

Family

ID=83638638

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2021063145A Active JP7590908B2 (en) 2021-04-01 2021-04-01 Floor vibration analysis system

Country Status (1)

Country Link
JP (1) JP7590908B2 (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116090051B (en) * 2022-12-27 2025-08-26 清华大学 Method and device for establishing rail transit-soil-structure coupling model
CN119167480A (en) * 2024-08-29 2024-12-20 宁波大学 A method for calculating the natural frequency of partial vibration modes of point-supported building curtain wall panels
KR102755079B1 (en) * 2024-11-05 2025-01-21 서울대학교산학협력단 Finite Element Model Updating Method and System using Improving System Identification through Sensor Fusion Techniques

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009174120A (en) 2008-01-21 2009-08-06 Takenaka Komuten Co Ltd Ceiling structure, method of constructing ceiling structure, and building
JP2016128640A (en) 2015-01-09 2016-07-14 株式会社竹中工務店 Building structure and vibration control device
JP2017182408A (en) 2016-03-30 2017-10-05 株式会社フジタ Floor vibration analysis method, floor vibration analysis program, and floor vibration analysis device
US20180106696A1 (en) 2015-03-23 2018-04-19 Ghyslaine MCCLURE Seismic response assessment of man-made structures

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009174120A (en) 2008-01-21 2009-08-06 Takenaka Komuten Co Ltd Ceiling structure, method of constructing ceiling structure, and building
JP2016128640A (en) 2015-01-09 2016-07-14 株式会社竹中工務店 Building structure and vibration control device
US20180106696A1 (en) 2015-03-23 2018-04-19 Ghyslaine MCCLURE Seismic response assessment of man-made structures
JP2017182408A (en) 2016-03-30 2017-10-05 株式会社フジタ Floor vibration analysis method, floor vibration analysis program, and floor vibration analysis device

Also Published As

Publication number Publication date
JP2022158338A (en) 2022-10-17

Similar Documents

Publication Publication Date Title
JP7590908B2 (en) Floor vibration analysis system
Zou et al. An optimal resizing technique for seismic drift design of concrete buildings subjected to response spectrum and time history loadings
El-Dardiry et al. Modelling of the dynamic behaviour of profiled composite floors
JP7368805B2 (en) Structure data generation device, structure data generation method, and structure data generation program
Kaveh et al. Cost optimization of steel-concrete composite floor systems with castellated steel beams
Yenidogan et al. Shake table test of a full-scale four-story reinforced concrete structure and numerical representation of overall response with modified IMK model
KR101655058B1 (en) Modal participation based flexural stiffness estimation method for building structure, and structural health assessment method using of it
Abdelaziz et al. Seismic evaluation of reinforced concrete structures infilled with masonry infill walls
JP5016468B2 (en) Design support system for damping buildings
Vaziri et al. The wall–frame interaction effect in corrugated steel plate shear walls systems
JP7336945B2 (en) Earthquake damage prediction system
Bhatt et al. Assessing the seismic response of existing RC buildings using the extended N2 method
JP2017182408A (en) Floor vibration analysis method, floor vibration analysis program, and floor vibration analysis device
Jeyaraj Vibro-acoustic behavior of an isotropic plate with arbitrarily varying thickness
Satheesh et al. Effect of in-plan eccentricity in vertically mass irregular RC framed buildings under seismic loads
JP6815547B2 (en) Floor vibration analysis program and floor vibration analysis device
Badri Rational decisions for seismic collapse-prevention design of low-rise asymmetric RC buildings considering deterioration properties
Latifi et al. A comparison of structural analyses procedures for earthquake-resistant design of building
JP7448405B2 (en) Floor vibration prediction system
JP2001349776A (en) Building vibration level prediction method and building design method
JP2025130163A (en) Method for predicting vibration of dry-jointed floors and method for modifying design based on the results of vibration prediction of dry-jointed floors
Badaoui et al. Estimation of accidental eccentricities for multi-storey buildings using artificial neural networks
JP2009059221A (en) Design simulation system, design simulation method, design simulation program, and recording medium
JP5155047B2 (en) Method for predicting vertical vibration of buildings
Lin et al. Seismic performance of high strength reinforced concrete buildings evaluated by nonlinear pushover and dynamic analyses

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20240319

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20241021

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: 20241029

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20241115

R150 Certificate of patent or registration of utility model

Ref document number: 7590908

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150