JP2976353B2 - Dynamic response property analysis system by wind force - Google Patents

Description

DETAILED DESCRIPTION OF THE INVENTION

[0001]

BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a dynamic model of an analytical model similar to a building model based on a wind tunnel test of a building model, generating directional wind pressure coefficient data from the observation results, and using the data as a basis. The present invention relates to a system that enables analysis of response characteristics.

[0002]

2. Description of the Related Art When a geometrically similar object, for example, a model that is geometrically similar to a full-scale building, is placed in a fluid having a similar flow, the similarity rule that the wind pressure around the object is similar. According to, a model of about one hundredth of a building is installed in the wind tunnel, and the wind is sent to the model, and the wind speed V considering the environment around the model
And the wind direction blown to the model,
The wind pressure P in the normal direction is measured. Then, the speed pressure q is calculated by the following equation (1) q = 1/2 (ρV ² ) (ρ is the air density), and the normal direction wind pressure coefficient C _q is calculated by the following equation (2) C _q = P / It is obtained by calculating q. Next, so that (the A pressure receiving area) the following equation _{(3) P W = qC q} A calculates Kazegairyoku P _w is used to grasp the dynamic response characteristics due to wind force of the wind direction.

[0003]

As described above, the design guidelines for grasping the dynamic response characteristics due to the external wind force in the wind direction have already been specified, but the guidelines for analyzing the dynamic response characteristics in the orthogonal direction have not been defined. The Architectural Institute of Japan does not specify anything, so there is a problem that it is not possible to analyze the response properties in all directions, and even if it is similar to a building model that was subjected to a wind tunnel experiment, it will be designed Previously, it is common practice to perform a wind tunnel experiment each time to obtain a wind pressure coefficient and an external wind force, so that there is a problem that it takes time and effort.

[0004] Furthermore, the dynamic response of a dome-shaped building due to an external wind force has not been analyzed yet.

SUMMARY OF THE INVENTION The present invention has been made in view of the above-mentioned problems, and has as its object to dynamically analyze an analysis model similar to a building model used in a wind tunnel test by an external wind force applied in all directions. It is an object of the present invention to provide a device that enables analysis of a response situation and that analyzes dynamic response characteristics due to an external wind force without newly performing a wind tunnel experiment on an analysis model similar to a building model once subjected to an experiment.

Further, the present invention obtains a two-dimensional directional component wind pressure coefficient and a torsional wind pressure coefficient on all sides of a high-rise building model in order to calculate an external wind force applied to an analysis model similar to the high-rise building model. It is an object of the present invention to provide an analyzing apparatus provided with a directional wind pressure coefficient data storage unit in which is stored.

Another object of the present invention is to provide an analysis model similar to a high-rise building model and an apparatus for analyzing the dynamic response state of its members by the wind main direction, the orthogonal direction and the torsional wind force. .

Another object of the present invention is to provide an analysis apparatus including a directional wind pressure coefficient data storage means for calculating an external wind force applied to an analysis model similar to a dome-shaped building.

Another object of the present invention is to provide an analysis model similar to a dome-shaped building model, and an apparatus for analyzing a dynamic response state of a member of the analysis model due to an external wind force.

[0010]

According to the present invention, as shown in the block diagram of the invention in FIG. 1, the normal direction for each time history at each observation point level obtained by rotating a building model by a wind tunnel experiment is shown. From the wind pressure coefficient, means for obtaining a directional wind pressure coefficient at each observation point level, directional wind pressure coefficient data storage means for storing the directional wind pressure coefficient output from the means for calculating the directional wind pressure coefficient, and read from the directional wind pressure coefficient data storage means Means for correcting the directional wind pressure coefficient to the height level of the analytical model similar to the building model, and means for calculating the external wind force applied to the analytical model similar to the building model from the directional wind pressure coefficient corrected to the height level Means for calculating the response displacement, velocity, and acceleration of the analysis model by inputting the external wind force, and the product of the response displacement of the member end obtained by inputting the external wind force and the rigidity of the member. The axial force of the member of the analysis model, the shear force of the member, and means for calculating the moment of the member, the response displacement, velocity, and acceleration and the axial force of the member, the shearing force of the member,
And means for selecting the maximum value from the values of the respective time histories with the moment of the member.

In the present invention, as shown in the block diagram of the invention in FIG. 2, the building model is a high-rise building model,
A means for calculating a wind pressure coefficient of a two-dimensional directional component at each observation point level, a means for calculating a torsional wind pressure coefficient at each observation point level, and a two-dimensional direction Directional wind pressure coefficient data storage means for storing the component wind pressure coefficient and the torsional wind pressure coefficient.

Further, as shown in the block diagram of the invention in FIG. 2, the present invention is based on the wind pressure coefficient in the normal direction at each observation point level on all sides of a high-rise building model obtained by a wind tunnel experiment. Means for calculating the two-dimensional directional component wind pressure coefficient, means for calculating the torsional wind pressure coefficient, directional wind pressure coefficient data storage means for storing the two-dimensional direction wind pressure coefficient and torsional wind pressure coefficient, and reading from the directional wind pressure coefficient data storage means Means for correcting the two-dimensional directional component wind pressure coefficient to a floor level of an analysis model similar to a high-rise building model, and a two-dimensional direction from the two-dimensional directional component wind pressure coefficient corrected to the floor height level and the torsion wind pressure coefficient. Means for calculating the component and torsional wind force, and the member shear force output from the means for calculating the member axial force, member shear force, and member moment of the analysis model. Means for calculating a layer shear analysis model is a force, which is constituted to comprise a means for calculating the overturning moment of the analysis model is input layer shear.

Further, according to the present invention, as shown in the block diagram of the invention in FIG. 3, the building model is a dome-shaped building model, and the dome-shaped building model is provided at each observation point level on the entire surface of the dome. Means for converting the normal direction wind pressure coefficient into a three-dimensional direction component wind pressure coefficient at each observation point level; and a direction wind pressure coefficient data storage means for storing the normal direction wind pressure coefficient and the three-dimensional direction component wind pressure coefficient from the conversion means. Is provided.

Further, according to the present invention, as shown in the block diagram of the invention in FIG. 3, the normal direction wind pressure coefficient at each observation point level on the entire surface of the dome of the dome-shaped building model obtained by the wind tunnel test is calculated. Means for converting the wind pressure coefficient into the three-dimensional component wind pressure coefficient, the direction wind pressure coefficient data storage means for storing the normal direction wind pressure coefficient and the three-dimensional direction component wind pressure coefficient from the conversion means, and the direction wind pressure coefficient data storage means. Means for correcting the normal direction wind pressure coefficient to each node position level of the analysis model similar to the dome-shaped building model, and the three-dimensional component wind external force from the normal direction wind pressure coefficient corrected to the node position level And a calculating means.

[0015]

[Action] From the wind tunnel experiment, the directional wind pressure coefficient at each observation point level is obtained from the normal direction wind pressure coefficient at each observation point level while rotating the building model at each observation point level.
This directional wind pressure coefficient is stored in the directional wind pressure coefficient data storage means, the read directional wind pressure coefficient is corrected to the height level of an analysis model similar to the model, the wind external force is calculated from the directional wind pressure coefficient, and the wind external force is calculated. Calculate the response displacement, velocity, and acceleration of the analysis model as input, and calculate the axial force, shear force, and moment of the analysis model member from the product of the response displacement at the end of the member and the rigidity of the member obtained by inputting the external wind force. Is calculated,
Find the maximum value from that value.

Further, the wind pressure coefficient of the two-dimensional direction component of the surface at each observation point level and the wind pressure coefficient of torsion are calculated from the normal direction wind pressure coefficients at each observation point level on all sides of the high-rise building model. Is stored in the direction wind pressure coefficient data storage means.

Further, the two-dimensional directional component wind pressure coefficient read from the directional wind pressure coefficient data storage means of the high-rise building model is corrected to a floor height level of an analysis model similar to the high-rise building model, and the floor height level is corrected. A two-dimensional wind pressure component and a torsional wind force are calculated from the corrected two-dimensional wind component wind pressure coefficient, and the response displacement of the analytical model calculated by inputting the wind force and the axial force and shear force of the member are input. , The moment is calculated, the shear force is input, the layer shear force of the analysis model is calculated, the layer shear force is input, the overturning moment of the analysis model is calculated, and the maximum value is obtained from these values.

Also, the normal direction wind pressure coefficient of the entire dome of the dome-shaped building model is converted into a three-dimensional direction component wind pressure coefficient at each observation point level, and the normal direction wind pressure coefficient and the three-dimensional direction component wind pressure coefficient are converted into directions. It is stored in the wind pressure coefficient data storage means.

Further, the normal wind pressure coefficient at each observation point level read from the directional wind pressure coefficient data storage means is corrected to each nodal position level of the analysis model similar to the dome-shaped building model, and the corrected normal is corrected. 3 from directional wind pressure coefficient
Calculate the external wind force of the dimensional component.

[0020]

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS The details of the present invention will be described below with reference to the accompanying drawings.

4 to 15 show an embodiment of the present invention. As shown in FIG. 4, in a wind tunnel 1 partially indicated by a dotted line, a model 2 of a hundredth of a building having a high-rise shape, having a three-dimensional frame by columns and beams, and having, for example, a square cross section is provided. Are arranged. Then, the observation point levels H1... Hk... Hn (Hk
Only one is shown. ) Has many holes 3 ...
One end (not shown) of a pipe 5... Is connected to these holes 3... And this pipe 5... Is connected to the wind pressure detector (not shown) built in the measuring device 4 from the inside of the model 2 through the bottom surface of the wind tunnel 1. It is connected. Further, the measuring device 4 is provided with a memory for storing the air density ρ and the calculation result, and an arithmetic circuit for calculating the normal direction wind pressure coefficient.

Although not shown, an anemometer, which is blown from the direction of the arrow and measures the wind speed V in consideration of the surrounding environment of the model 2, is also provided in the wind tunnel 1.

In the figure, reference numeral 6 denotes a peripheral building model installed on a turntable (not shown) together with the high-rise model 2, and A1, A2, and A3 indicate arrangements of the wind holes 3 at the observation point level Hk. At the installation position, the dominant area on which the intensively calculated wind pressure coefficient acts is shown.

When the wind is sent to the high-rise model 2 from the direction of the arrow for an appropriate time, the wind speed V affected by the surrounding environment by the surrounding buildings 6 is detected, while the wind holes 3 in the model 2 are detected.
Is sent to the measuring device 4 through the pipes 5... And the fluctuating wind pressure is detected by the wind pressure detector. At this time, tube 5
The measurement error due to the difference in the length is corrected by correcting the frequency characteristic, and an accurate wind pressure P in the normal direction for each time history is detected.

In the arithmetic circuit of the measuring device 4, the wind speed V
And the air density ρ called out from the memory, q = １ ／
Substitute ρV ² to calculate the design speed pressure q. Further, the normal direction wind pressure coefficient for each time history is calculated by substituting the above-mentioned design speed pressure q and wind pressure P into an equation, C = P / q, and adding a code for each observation point level. As shown in FIG. 6, this is read into the RAM 12A connected to the central processing unit 16 of the apparatus for analyzing the dynamic response situation due to the external wind force of the assumed analysis model, as shown in FIG.

The RAM 12A includes, via an input device 14, the dimensions of the building model 2 and the control area A1 of the building model 2 shown in FIG. .. A6 and distances L1, L2, L3 from the center line in the main direction of building model 2 shown in FIG.
L4 and the design speed pressure q are stored.

The setting data used when calculating the external wind force includes floor heights h1, h2,... Of the analysis model 2 '.
hn, the dominant height widths h'1, h'2,... of each floor height, and the width W in the horizontal direction are stored. Further, the analysis model 2 'is used as setting data used in dynamic response analysis. , The stiffness matrix K, the damping matrix C, the mass matrix M, the stiffness matrix Km ′ of the member in the member coordinate system of the member 17 of the analysis model 2 ′, and the heights e1, e2. 15 is stored.

On the other hand, in the RAM 12B detachably connected to the central processing unit 16, two-dimensional wind pressure coefficients FX and FY and a twist wind pressure coefficient F # are stored as directional wind pressure coefficient data for analysis.

In order to convert the normal wind pressure coefficient calculated from the result of the wind tunnel experiment into the directional wind pressure coefficient for analysis, the ROM 13 stores the normal wind pressure coefficient for each of the observation point levels H 1. Dominant areas A1, A2 at each observation point level
..Calculation to convert to normal wind pressure coefficient per hit and calculation of directional wind pressure coefficients FX and FY per unit area of two-dimensional dominant area from normal wind pressure coefficient per control area;
Control program for creating directional wind pressure coefficient data for calculating the torsional wind pressure coefficient FΘ, the two-dimensional directional wind pressure coefficients FX and FY, and the torsional wind pressure coefficient F
In order to convert Θ into a wind external force applied to the floor level of the analysis model 2 ′, the two-dimensional directional component wind pressure coefficients FX and FY for each observation point level are calculated for each floor level h1 of the analysis model 2 ′.
h2 .., and the two-dimensional wind pressure coefficients FX and FY
Is converted into a wind pressure coefficient distribution that is continuous in the vertical direction of the main direction plane and the orthogonal direction plane, respectively, and further, the dominant height widths h′1, h′2,. The applied wind pressure coefficient distribution is multiplied by the lateral width W of the analysis model to calculate its volume Dj, which is then multiplied by the design speed pressure q to continue the analysis model 2 ′ in the X and Y planes in the vertical direction. A control program for calculating wind external force Fx, Fy and calculating a wind external force Fθ from the vector sum of the wind external forces Fx, Fy is stored.

Further, the wind external force calculated as described above,
The stiffness matrix K, mass matrix M, and damping matrix C of the analysis model are expressed by the following incremental vibration equation (4): M (x ″) + C (x ′) + K (Δx) + (fs) = F (t) ( Here, Δx is a response displacement vector, x ′ is a velocity vector, x ″ is an acceleration vector, fs is a static balance vector, and F (t) is an external wind force Fx, Fy, Fθ for each time history.) Then, the vibration equation is solved by a numerical integration method to obtain the two-dimensional response and the torsional response displacement of the analysis model. Further, a1 (incremental displacement at the current time) + a3 (velocity at the previous time) + a
4 (Acceleration at the previous time) The response speed is calculated by the following equation: a0 (incremental displacement at the current time) + a2 (speed at the previous time) + a
3 (acceleration at the previous time), a control program for obtaining response acceleration is stored (a0 to a4 are constants, if Δt is the analysis time and β is 0.25, a0 = 1 / β (△ t)
² , a1 = 0.5 / (β △ t), a2 = −1 / (β △)
t), a3 = 1-0.5 / β, a4 = (1-0.25 /
β) △ t. Also, the static balance vector is
(Fs) _i2 = (fs) _i1 + K _i (Δx) It is obtained from ⁱ² _i1 . Here, i1 indicates a step at an arbitrary time, and i2 indicates a step next to i1. ).

Further, as shown in FIG. 13, the stresses generated at the i-end and the j-end of the members 17 such as columns and beams in each layer in the member coordinate system used in the analysis model, that is, the axial forces Ni and Nj, , Y-axis and z-axis shear forces Qyi, Qzi and Qyj, Qzj, and torsion moment My in the x-axis direction
To calculate i, Mzi and Myj, Mzj, R
The member stiffness matrix Km 'of the member coordinate system stored in the AM 12A is called, and this is expressed by the following equation (5), Km = Tg
^T · Km ′ · Tg (where Tg indicates a coordinate conversion matrix, and Tg ^T indicates a transposed coordinate conversion matrix), is converted to a member rigidity matrix Km of the whole coordinate system, and ui, vi, wi, θxi, θyi,
θzi, and uj, vj, wj, θxj, θyj, θz
And the total number of nodes of the analysis model × 6 components (ui, vi, wi, θx at the end of the member)
(6 degrees of freedom of i, θyi, θzi). This overall rigidity matrix Kg, the mass matrix M of the member,
The damping matrix C of the member is substituted into the above-mentioned vibration equation M (x ″) + C (x ′) + Kg (Δx) = F (t), an external wind force F (t) is added, and the total stiffness Kg is calculated by the Newmark β method. Thus, the following iterative equation: velocity vector, (x ′) _i2 = (x ′) _i1 + (Δt)
[(X ″) _i1 + (x ″) _i2 ] / 2 displacement increment vector, (Δx) ⁱ¹ _i2 = Δt (x ′) _i1
+ (Δt) ² [(1 / 2−β) (x ″) _i1 + β
(X ″) _i2 ] and solve (assuming i1 is a step and i2 is the next step after i1), calculate displacement, velocity, and acceleration vector, and calculate the displacement, velocity, and acceleration vector in the global coordinate system shown in Expression (6). Obtain displacement vectors x for all nodes,

[0032]

(Equation 6)

From the displacement vector x, the displacements at the i and j ends of the arbitrary member 17 are extracted. Since this displacement is a displacement of the whole coordinate system, the displacement conversion formula of the member coordinate system of the following equation (7) is used.

[0034]

(Equation 7)

That is, the displacement is substituted into the displacement conversion formula of (displacement of the member coordinate system) = Tg ^T (displacement of the whole coordinate system), and as shown in FIG. 14, the displacement of the whole coordinate system is converted into the displacement of the member coordinate system. The coordinate system displacement and the member stiffness matrix are expressed by equation (8).

[0036]

(Equation 8)

As shown in FIG. 15, the control program for calculating the member stress by calculating the following formula, and the layer shear force, which is the sum of the shear forces of the columns existing in each layer of the analysis model 2 ', are calculated as shown in FIG. The calculated shear forces Qzi, Qz in the z-axis direction for each layer
1Qa,... nQa, 1Qb,.
A control program that calculates by substituting nQb into the following equation (9): Shearing force of the first layer = 1Qa + 1Qb.... Shearing force of the Nth layer = nQa + nQb and, as shown in FIG. Heights e1, e2,... En of each layer of the called analysis model 2 ′;
The shear force 1Qa of each layer calculated by the above equation (9),...
NQa and 1Qb,... 1Qnb are calculated by the following equation (1).
0) Falling moment of one layer = (1Qa + 1Qb) × e1 Falling moment of two layers = (1Qa + 1Qb) × (e1 +
e2) + (2Qa + 2Qb) × e2... n-layer fall-down = (1Qa + 1Qb) × (e1 +
e2 + ... en) + (2Qa + 2Qb) × (e2 + ·
..En) +... + (NQa + nQb) × en and a control program to be calculated, and the above-mentioned response displacement, its speed and acceleration in each time history, axial force of member, shear force of member, member , And a control program for obtaining the maximum value from the values of the layer shearing force and the overturning moment and inputting the maximum value to the output device 15 is stored.

The operation of the system configured as described above will be described with reference to FIGS. 4 to 15 and flowcharts shown in FIGS.

In FIG. 16, a high-rise building model 2 is set on a turntable in a wind tunnel 1 (step S1), and while the table is being rotated, wind is blown for a predetermined period from the direction of an arrow to blow all sides of the model 2. Is detected (step S2), and the pipe 5 from the hole 3 to the measuring device 4 is detected.
Correcting the error of the wind pressure P due to the difference in length by correcting the frequency characteristic (step S3). Further, a wind speed V affected by the surrounding environment is detected (step S4).

Next, after calculating the design speed pressure q by the above equation (1) (step S5), using the above equation (2), the corrected wind pressure P is divided by the design speed pressure q to obtain each time history. Is calculated (step S6), and the calculation result is stored in a memory built in the measuring device 4 (step S7).

Next, in FIG. 17, the contents stored in the above-mentioned memory are read into the RAM 12A of the arithmetic unit shown in FIG.
8) Design speed pressure q and observation point level H1 of experimental model 2
..The dominant areas A1..A6 in the main and orthogonal directions at Hn
And the center distance L1... Ln from the main direction are stored (step S9).

Next, the observation point levels H1... H shown in FIG.
.. A6 in the main direction and the orthogonal direction in n are called from the RAM 12A, the normal direction wind pressure coefficients detected through the air holes 3... included in the controlled area are called, and these are aggregated at the observation point level. After calculating and obtaining wind pressure coefficient values CX1, CX2, CX4, CX5 and CY3, CY6 (step S10), the wind pressure coefficient values and the dominant areas A1, A2, A4, A5 are substituted into the following equation (11). , FX = ([CX1 · A1 + CX2 · A2] − [CX4 ·
A4 + CX5 · A5]) / [A1 + A2 + A4 + A5] The main direction wind pressure coefficient FX per unit area of the dominant area at each observation point level is calculated, and then the above wind pressure coefficient values CY3, CY
6 and the dominant areas A3, A6 into the following equation (12),
Y = [CY3 · A3-CY6 · A6] / [A3 + A6] The orthogonal wind pressure coefficient FY per unit area of the dominant area is calculated (step S11). Next, as shown in FIG. 8, distances L1, L2, L3, and L4 from the center line in the main direction are called, and the following equation (13) FΘ = ([CX2 · A2 · L2 + CX5 · A5 · L3]
-[CX1, A1, L1 + CX4, A4, L4]) /
[A1 + A2 + A3 + A4 + A5] to calculate a torsional wind pressure coefficient FΘ (step S
12). After performing the above calculation processing for each observation point level and for each time history, these directional wind pressure coefficient data FX, F
Y, F} are stored in the RAM 12B (step S1).
3).

In FIG. 18, the RAM 12A further stores
Via the input device 14, floor levels h1, h of the analysis model
2., and the dominant height width h'1, h'2 .. of each floor height,
The width W in the horizontal direction in the main and orthogonal directions, the heights e1, e2,... Of each layer of the analysis model, the rigidity matrix K, the damping matrix C, and the mass matrix M of the analysis model
And the member rigidity matrix Km 'in the member coordinate system are stored (step S14).

Next, the two-dimensional direction and the torsional wind pressure coefficient F
In order to convert X, FY, FΘ into the external wind force of the analysis model 2 ′, two-dimensional and torsional directional wind pressure coefficient data are read from the RAM 12B (step S15), and the main direction wind pressure coefficient FX is assumed as shown in FIG. 1 are set at the floor height levels h1, h2,... On the main direction plane of the analysis model 2 '.
As shown in FIG. 0, a continuous vertical wind pressure coefficient distribution is formed between these floor height levels (step S16).
As shown in FIG. 11, the dominant height h′1 of the wind pressure coefficient distribution,
(step S17), the wind pressure coefficient distribution applied to the dominant height widths h'1, h'2... is calculated (step S18), and as shown in FIG.
The volume D1..Dj .. of the continuous wind pressure coefficient distribution is calculated by multiplying the wind pressure coefficient distribution applied to 1 ', h'2 .. and the lateral width W of the analysis model (step S19), and the volume D
.., Dj,... Are multiplied by the called design speed pressure q, the external wind force Fx applied thereto is calculated, and the external wind force for each time history is calculated (step S20).

Similarly, the wind external force Fy in the orthogonal direction
, And Fθ from the vector sum of the external wind forces Fx and Fy.
Are calculated and stored in the RAM 12A as wind external force data (step S21).

Next, in FIG. 19, the above-mentioned vibration equation (4) is added to the wind external force Fx, which is called from the RAM 12A.
Fy, Fθ, the stiffness matrix K, the damping matrix C, and the mass matrix M of the analysis model are substituted and solved by a numerical integration method, and each response displacement, its speed and acceleration for each time history of the two-dimensional component and torsion are obtained. It is calculated (step S22).

Next, the member 17 used for the analysis model 2 '
Is retrieved from the RAM 12A, and this Km 'and a member end response displacement obtained by inputting an external wind force as described above and converting the obtained member end response displacement into a member coordinate system are expressed by the above equation. Then, the moment of the member, the axial force of the member, and the shearing force of the member are calculated (Step S23).

Next, the above-mentioned shearing force 1Qa ·· nQa,
1Qb ··· nQb is calculated by the above equation (9) to calculate the layer shearing force of each layer (step S24), and as shown in FIG. 15, the shearing force of each layer and the heights e1 and e of each layer of the analysis model.
2 ·· en is substituted into the above equation (10) to calculate the overturning moment (step S25). And the response displacement,
The maximum value is determined from the time history values of the speed, acceleration, moment of each member, axial force, shearing force, layer shearing force, and overturning moment (step S26).

As described above, in addition to the maximum value of the dynamic response state due to the external wind force, the continuous instantaneous value is output to the CRT or the output device 15 of the printer for recording and display (step S2).
7).

Further, it is also possible to calculate and record and display a two-dimensional direction and torsional wind pressure coefficient, a Fourier spectrum and a power spectrum of an external wind force, a Fourier spectrum and a power spectrum such as a response displacement and a speed.

The cylindrical high-rise model 20, the high-rise model 21, and the high-rise model 22 having the spire shown in FIGS.
For wind tunnel experiments, create a wind pressure coefficient database based on the experimental results, store it in memory, load it into the analyzer, and apply it to an analytical model with the same shape and scale ratio as the above model. Calculate and record Fourier spectrum, power spectrum, response displacement, velocity, acceleration, etc.
It can also be displayed.

Next, the analysis of the dynamic response situation of the dome-shaped analysis model will be described. The same parts as those in the description of the high-rise building model will be briefly described.

The dome-shaped building model 8 shown in FIG. 5 is formed so as to have a three-dimensional frame by columns and beams.
Is installed on a turntable (not shown) in the wind tunnel 1 shown in FIG. 1 and one end of a pipe 5 is connected to wind holes 3... Provided on the entire surface of the dome of the dome-shaped building model 8. The inside of the building model 8 passes through the bottom surface of the wind tunnel 1 and is connected to a wind pressure detector that detects wind pressure in the normal direction of the measuring device 4.

Then, similarly to the above-described high-rise building model, the measuring device 4 calculates the design speed pressure q at each observation point level J using the air density ρ and the wind speed V.
Further, the actually measured wind pressure P is corrected to an accurate wind pressure P by correcting the frequency characteristic, and using this and the design speed pressure q, FIG.
As shown in FIG. 3, the wind pressure coefficient FJ in the normal direction at each observation point level J for each time history is obtained over the entire dome, a code is attached for each observation point level, and stored in the memory of the measuring device 4. Let it.

The dome-shaped building model 8 is stored in the RAM 12A.
25, and the nodal position levels k1 and k1 of the skeleton 9 of the dome-shaped analysis model 8 'as shown in FIG.
k2... and each node position level k as shown in FIG.
1, assumed area of dominant surface at k2
And a stiffness matrix K, a damping matrix C, and a mass matrix M of the analysis model 8 ', a stiffness matrix Km' of a member in a member coordinate system of the analysis model 8 ', and a design speed pressure q are stored.

The RAM 12B stores a normal direction wind pressure coefficient FJ at each observation point level for each time history of the dome-shaped building model 8 and a directional wind pressure of a three-dimensional direction component, which will be apparent from the following description. Coefficients FX, FY, and FZ are stored.

In the ROM 13, as shown in FIG.
The normal direction wind pressure coefficient FJ at the observation point level J of the dome-shaped building model 8 exemplarily showing a part of the dome surface is called, and this is converted into the XYZ plane of the three-dimensional absolute coordinate system, and X,
A control program for creating wind pressure coefficients FX, FY, and FZ in the Y and Z directions, and a normal direction wind pressure coefficient FJ as shown in FIG.
Is converted into a wind external force applied to the nodal position levels k1, k2,... Of the analysis model 8 'by a normal direction wind pressure coefficient FJ as shown in FIG. A directional wind pressure coefficient distribution is formed, and as shown in FIG. 27, each node position level k1, k2 of the analysis model 8 '
・ ・ Call dominant area of dome surface at 10 ・ ・
As shown in FIG. 28, the dominant area 10 is translated in the normal direction, and the points Pj1... Pjn at which the translation plane 11 intersects the normal direction wind pressure distribution and the nodes Mj1.
• Calculate the volume Dj of the wind pressure coefficient distribution composed of Mjn and multiply this by the design speed pressure q to calculate the wind external force Fj at each node position level, and as shown in FIG. Convert into three-dimensional direction components Fx, Fy, Fz,
A control program for generating an external wind force for obtaining the external wind force in the three-dimensional direction, and the rigidity, mass,
A control program for calculating the three-dimensional response displacement, response speed, and response acceleration of the analysis model from the vibration equation (4) described above, and inputting the external wind force as described above. Then, the member end displacement of the member coordinate system obtained by performing the coordinate conversion process and the member rigidity matrix Km 'of the member coordinate system of the analysis model are substituted into the above equation (8), and the axial force of the member of the analysis model is substituted. And a control program for calculating a shear force and a moment; a response displacement, a speed, an acceleration, an axial force,
A control program for obtaining the maximum value of the shearing force and the moment is stored.

Next, the operation of the device configured as described above will be described.
This will be described with reference to the flowcharts shown in FIGS.

In FIG. 30, a dome-shaped building model 8 is set on a turntable (not shown) in the wind tunnel 1 (step S31), and while rotating, a wind is blown for a predetermined period from the direction of the arrow to rotate the entire dome. The wind pressure P in the normal direction is detected (step S32), and the wind pressure P is corrected by correcting the frequency characteristic to correct the error due to the difference in the length of the pipes 5 from the holes 3 to the measuring device 4 (step S32). Step S
33). Further, the wind speed V affected by the surrounding environment is detected (step S34).

Next, after calculating the design speed pressure q (step S35), the corrected wind pressure P is divided by the design speed pressure q to obtain the observation points of the dome-shaped building model 8 as shown in FIG. Normal wind pressure coefficient F for each time history acting on the level
J is calculated (step S36), and an identification code is attached to each observation point level and stored in a memory built in the measuring device 4 (step S37).

Next, in FIG. 31, the normal direction wind pressure coefficient FJ of the memory is read into the RAM 12B of the apparatus shown in FIG. 6 via the input device 14 (step S3).
8) The setting data of the design speed pressure q is input to the RAM 12A (step S39).

The normal wind pressure coefficient FJ at each observation point level on the dome surface of the dome-shaped building model 8 is stored in the RAM 12.
B, and this is coordinate-converted into directional wind pressure coefficients FX, FY, and FZ on the three-dimensional absolute coordinate plane as shown in FIG. 24, that is, converted into the main direction, orthogonal direction, and vertical direction wind pressure coefficients on the dome surface. After performing this process for each observation point level (step S40), the three-dimensional directional wind pressure coefficients FX and F together with the previously stored normal direction wind pressure coefficient FJ are stored.
Y and FZ are stored in the RAM 12B as direction wind pressure coefficient data (step S41).

Next, in FIG. 32, the nodal position levels k1, k2,...
The RAM 12A reads the setting data of the matrixes K, C, and M of stiffness, attenuation, and mass of the dome-shaped analysis model and the stiffness matrix Km of the member (step S42). Next, in order to calculate the external wind force acting on the surface of the skeleton 9 of the dome-shaped analysis model 8 'shown in FIG. 25, the above normal direction wind pressure coefficient FJ is called from the RAM 12B (step S43), and shown in FIG. Thus, the normal direction wind pressure coefficient distribution is formed on the entire dome of the analysis model (step S44). Then, the dominant area 10 of the nodal position levels k1, k2,.
12A, and calculates the node coordinates Mj1... Mjm surrounding the node position level kj as shown in FIG.
8, the control area 10 is translated in the normal direction (step S45), and the point coordinates Pj1... Pjm of the surface 11 intersecting with the control area 10 and the wind pressure coefficient are calculated, and the node Mj is calculated.
The volume Dj of the wind pressure coefficient distribution of the polyhedron surrounded by 1... Mjm and the intersection Pj1... Pjm is calculated (step S46), and the volume Dj is multiplied by the design speed pressure q to calculate the external wind force Fj. As shown in FIG. 29, this wind external force Fj
Is the normal direction to the dominant area 10, and is converted into three-dimensional coordinates to convert the wind force Fx in the main direction, the orthogonal direction, and the vertical direction into
Fy and Fz are obtained (step S47), this processing is performed for each time history over the entire dome of the analysis model, and these are stored in the RAM 12A to create a wind external force database (step S48).

Next, in FIG. 33, the called wind external forces Fx, Fy, Fz and the stiffness matrix K, the damping matrix C, and the mass matrix M of the dome-shaped analysis model are substituted into the vibration equation (4). The response displacement, the speed, and the acceleration of each time history of the three-dimensional component are calculated by the numerical integration method (step S49).

The member stiffness matrix Km 'of the member coordinate system used for the dome-shaped analysis model is stored in the RAM 12A.
, And the external force is input, and the member end displacement converted into the member coordinate system is substituted into the above equation (8) to calculate the moment of the member, the axial force of the member, and the shearing force of the member. (Step S50). Next, a maximum value is obtained from the calculated response displacement, velocity, acceleration, and the values of the axial force, shearing force, and moment of the member (step S51).

In this manner, the maximum value of the dynamic response status value of the dome-shaped analysis model due to the external wind force is output to the CRT or the output device 15 of the printer for recording and displaying, and also displays the continuous instantaneous value. (Step S52).

Further, a Fourier spectrum and a power spectrum of a three-dimensional direction wind pressure coefficient and a normal direction wind pressure coefficient and an external wind force, a Fourier spectrum such as a response displacement and a velocity,
It is also possible to calculate, record and display a power spectrum.

The directional wind pressure coefficient database of the above-mentioned dome-shaped building model was loaded into the analysis device, and the analysis model having the same scale ratio as the dome-shaped building model tested was subjected to dynamic analysis without new wind tunnel experiments. It is possible to analyze the target response situation.

[0069]

As described above, according to the present invention, the directional wind pressure coefficient at each observation point level is obtained from the normal wind pressure coefficient at each observation point level obtained by rotating the building model by a wind tunnel experiment. Means, a directional wind pressure coefficient data storage means for storing the directional wind pressure coefficient output from the means for obtaining the directional wind pressure coefficient, and a directional wind pressure coefficient read from the directional wind pressure coefficient data storage means, the height of the analytical model similar to the building model. Means for calculating the external force to be applied to the analysis model from the directional wind pressure coefficient corrected for the height level, and calculating the response displacement, velocity and acceleration of the analysis model by inputting the external force. Means and the axial force, shear force, and moment of the member of the analysis model from the product of the member end response displacement and the member rigidity of the member coordinate system obtained by inputting the external wind force. And calculating the response displacement, speed, and acceleration, the axial force of the member, shear force, and the means for selecting the maximum value from the value of each time history of the moment of the member, so that it is configured to include It is possible to calculate the external wind force applied in all directions of an analysis model similar to a building model, and therefore analyze the dynamic response situation due to the external force applied in all directions and the dynamic response situation of the members of the analysis model. It becomes possible.

Since the directional wind pressure coefficient data storage means is provided, the directional wind pressure coefficient data stored in the directional wind pressure coefficient data storage means is applied to the analysis of the dynamic response situation of an analytical model similar to a building model. It is possible to
For this reason, for an analysis model similar to the building model once subjected to the experiment, the dynamic response situation can be analyzed without performing a new wind tunnel test, thereby saving the trouble of performing a new wind tunnel experiment. It becomes possible.

Further, the maximum value is selected and output from the response displacement, velocity, and acceleration of the analysis model similar to the building model and the time history of the axial force, shearing force, and moment of the member. Can be done.

Further, according to the present invention, the building model is a high-rise building model, and the normal direction wind pressure coefficients at each observation point level on all sides of the high-rise building model are used to calculate a value of 2 at each observation point level. Since it is configured to include means for calculating the dimensional component wind pressure coefficient, means for calculating the torsional wind pressure coefficient, and directional wind pressure coefficient data storage means for storing the two-dimensional directional component wind pressure coefficient and the torsional wind pressure coefficient, The coefficient data storage means can store the wind pressure coefficients in the main direction, the orthogonal direction, and the torsional direction of the high-rise building model.

According to the present invention, a means for calculating a two-dimensional wind pressure coefficient from a normal wind pressure coefficient at each observation point level on all sides of a building model obtained by a wind tunnel experiment, and a torsional wind pressure coefficient Calculating means, a two-dimensional wind pressure coefficient and a torsional wind pressure coefficient storing means, and a two-dimensional wind pressure coefficient read from the storage means are used to store the two-dimensional wind pressure coefficient in an analysis model similar to a high-rise building model. From the means for correcting to the level, the means for calculating the two-dimensional directional component corrected to the floor height level, and the torsional wind force, and the means for calculating the axial force, shear force, and moment of the member of the analysis model. Means for calculating the layer shear force of the analysis model by inputting the output shear force, and means for calculating the overturning moment of the analysis model by inputting the layer shear force Therefore, for an analysis model similar to a high-rise building model, it is possible to analyze not only the main direction of the wind, but also the dynamic response situation in the orthogonal direction and the torsional external force. The situation can be analyzed.

Further, according to the present invention, means for converting the normal direction wind pressure coefficient at each observation point level on the entire surface of the dome of a dome-shaped building model into a wind pressure coefficient of a three-dimensional component at each observation point level, A directional wind pressure coefficient data storage means for storing the linear wind pressure coefficient and the three-dimensional directional component wind pressure coefficient is provided, so that the wind added to the analysis model similar to the dome-shaped building model is stored in the directional wind pressure coefficient data storage means. Directional wind pressure coefficient data serving as a basis for calculating an external force can be stored.

Further, according to the present invention, the normal direction wind pressure coefficient at each observation point level on the entire surface of the dome of the dome-shaped building model obtained by the wind tunnel experiment is converted into a three-dimensional component wind pressure coefficient at each observation point level. Means for storing the normal direction wind pressure coefficient and the three-dimensional direction wind pressure coefficient from the conversion means, and a normal direction wind pressure coefficient read from the direction wind pressure coefficient data storage means. It is configured to have a means for correcting each node position level of the analysis model similar to the model and a means for calculating the external wind force of the three-dimensional direction, so that it is similar to the dome-shaped building model that has not been analyzed yet. With respect to the dome surface of the analysis model to be analyzed and its members, the dynamic response state can be analyzed.

[Brief description of the drawings]

FIG. 1 is a configuration diagram of a first invention of the present invention.

FIG. 2 is a configuration diagram of the second and third inventions of the present invention.

FIG. 3 is a configuration diagram of a fourth and a fifth invention of the present invention.

FIG. 4 is a perspective view showing a part of a wind tunnel in which a high-rise building model is arranged.

FIG. 5 is a perspective view of a dome-shaped building model.

FIG. 6 is a block diagram of an arithmetic unit according to the present invention.

FIG. 7 is an explanatory diagram for calculating main and orthogonal wind pressure coefficients;

FIG. 8 is an explanatory diagram for calculating a torsional wind pressure coefficient.

FIG. 9 is a diagram showing an assumed analysis model.

FIG. 10 is a diagram showing a wind pressure coefficient distribution formed on the main direction surface of the analysis model.

FIG. 11 is a diagram showing a wind pressure coefficient distribution set to a dominant width of each floor height.

FIG. 12 is a diagram showing a wind pressure coefficient distribution volume.

FIG. 13 is a diagram showing displacements and rotation angles acting on both ends of a member of the analysis model.

FIG. 14 is a diagram illustrating a coordinate relationship between members in the overall coordinate system and the member coordinate system.

FIG. 15 is an explanatory diagram for obtaining a falling moment.

FIG. 16 is a flowchart for calculating a normal direction wind pressure coefficient applied to a high-rise building model installed in a wind tunnel.

FIG. 17 is a flowchart for obtaining a two-dimensional directional wind pressure coefficient and a torsional wind pressure coefficient from directional wind pressure coefficient data obtained from a result of a wind tunnel experiment.

FIG. 18 is a flowchart for obtaining wind external force data.

FIG. 19 is a flowchart for obtaining a response displacement and the like for an analysis model, an axial force of a member, and the like.

FIG. 20 is a diagram showing a cylindrical high-rise building model used for a wind tunnel experiment.

FIG. 21 is a view showing a high-rise building model used for a wind tunnel test.

FIG. 22 is a view showing a high-rise building model with a spire used in a wind tunnel experiment.

FIG. 23 is a diagram showing a normal direction wind pressure coefficient acting on the entire dome surface of the dome-shaped building model.

FIG. 24 is a diagram showing wind pressure coefficients acting on a part of the dome surface of the dome-shaped building model in the normal direction and the three-dimensional direction.

FIG. 25 is a diagram showing an external wind force acting on a skeleton forming a dome-shaped analysis model.

FIG. 26 is a diagram showing a wind pressure coefficient distribution on the dome surface of the analysis model.

FIG. 27 is a diagram showing a dominant surface at a nodal position level of the dome surface of the dome-shaped analysis model.

FIG. 28 is a diagram showing a wind pressure coefficient distribution volume on the dome surface.

FIG. 29 is a diagram showing an external wind force in a normal direction of a dominant surface and decomposing the external force into a three-dimensional component.

FIG. 30 is a flowchart for obtaining a normal direction wind pressure coefficient of a dome-shaped building model.

FIG. 31 is a flowchart for obtaining directional wind pressure coefficient data.

FIG. 32 is a flowchart for obtaining an external wind force.

FIG. 33 is a flowchart for obtaining an axial force of a member such as a response displacement of an analysis model;

[Explanation of symbols]

1 wind tunnel, 2 and 8 building model, 3 model, 4 measuring device, 5 pipe from air hole to measuring device, 16 central processing unit, 14 input device, 12 and 12A storage device, 1
5 Output device

──────────────────────────────────────────────────続 き Continued on the front page (58) Field surveyed (Int.Cl. ⁶ , DB name) G01M 9/00 G01M 9/06

Claims (5)

(57) [Claims] A means for obtaining a directional wind pressure coefficient at each observation point level from a normal direction wind pressure coefficient for each time history at each observation point level obtained by rotating a building model by a wind tunnel experiment; Directional wind pressure coefficient data storage means for storing the directional wind pressure coefficient output from the means for obtaining the directional wind pressure coefficient, and correcting the directional wind pressure coefficient read from the directional wind pressure coefficient data storage means to a height level of an analysis model similar to the building model. Means for calculating the wind force to be applied to the analysis model from the directional wind pressure coefficient corrected to the height level, and calculating the response displacement, speed, and acceleration of the analysis model by inputting the wind force. Means,
Means for calculating the axial force of the member of the analysis model, the shear force of the member, and the moment of the member from the product of the response displacement of the member end and the rigidity of the member obtained by inputting the external wind force, and the response displacement, speed, Means for selecting a maximum value from values of time histories of acceleration, an axial force of a member, a shear force of a member, and a moment of a member, and a dynamic response characteristic analysis system based on an external wind force. . 2. The building model is a high-rise building model, and a wind pressure of a two-dimensional component at each observation point level is obtained from a normal direction wind pressure coefficient at each observation point level on all sides of the high-rise building model. The apparatus according to claim 1, further comprising: means for calculating a coefficient; means for calculating a torsional wind pressure coefficient; and directional wind pressure coefficient data storage means for storing the two-dimensional directional component wind pressure coefficient and the torsional wind pressure coefficient. Dynamic response analysis system by wind force. 3. A means for calculating a two-dimensional component wind pressure coefficient at each observation point level from a normal direction wind pressure coefficient at each observation point level on all sides of a high-rise building model obtained by a wind tunnel experiment, and a torsional wind pressure coefficient. , A directional wind pressure coefficient data storage means for storing the two-dimensional directional wind pressure coefficient and the torsional wind pressure coefficient, and a two-dimensional directional component wind pressure coefficient read from the directional wind pressure coefficient data storage means. Means for correcting the floor level of the analysis model similar to the above, and means for calculating the two-dimensional direction component and the torsional wind force from the two-dimensional direction component wind pressure coefficient and the torsion wind pressure coefficient corrected to the floor level, The member shear force output from the means for calculating the axial force of the member, the member shear force, and the member moment of the analysis model is input, and the layer shear of the analysis model is input. 2. The dynamic response characteristic analysis system according to claim 1, further comprising: means for calculating a force; and means for calculating the overturning moment of the analysis model by inputting the layer shearing force. 4. The building model is a dome-shaped building model, and converts a normal direction wind pressure coefficient over the entire dome of the dome-shaped building model into a wind pressure coefficient of a three-dimensional component at each observation point level. Means, the normal wind pressure coefficient and 3
2. The dynamic response characteristic analysis system according to claim 1, further comprising: a direction wind pressure coefficient data storage unit configured to store a wind pressure coefficient of a dimensional direction component. 5. A conversion means for converting a normal direction wind pressure coefficient at each observation point level on the entire surface of the dome of a dome-shaped building model obtained by a wind tunnel experiment into a wind pressure coefficient of a three-dimensional component at each observation point level. A directional wind pressure coefficient data storage means for storing a normal wind pressure coefficient and a three-dimensional directional component wind pressure coefficient from the conversion means; and a normal wind pressure coefficient read from the directional wind pressure coefficient data storage means, the dome-shaped building. Means for correcting each node position level of the analysis model similar to the model, and means for calculating a three-dimensional component wind force from the normal direction wind pressure coefficient corrected to the node position level. The dynamic response characteristic analysis system based on an external wind force according to claim 1.