JP7211286B2 - Buckling stress estimation device, buckling stress estimation method, and buckling stress estimation program - Google Patents

Description

The present invention relates to a buckling stress estimation device, a buckling stress estimation method, and a buckling stress estimation program.

Conventionally, the buckling stress (local buckling stress) of a plate element to which a compressive force acts is estimated. In this buckling stress estimation method, it is common to use the Fourier series represented by formula (1), and when estimating the out-of-plane displacement _Ww of the plate element, formula (1) is changed to -1) The degree of buckling stress is estimated after changing the expression of the equation (see, for example, Non-Patent Document 1).
However, _bw is the length in the direction orthogonal to the thickness direction of the plate element.

Then, by applying the energy method to this Fourier series, the degree of buckling stress is estimated. Since the basis of the Fourier series has orthogonality, many curved surfaces can be approximated with high accuracy by the Fourier series.

Stephen P. Timoshenko and James M. Gere, "Theory of Elastic Stability" Second Edition

However, Fourier series suffer from the slow convergence of approximations. For example, when trying to approximate the buckling deformation that occurs when a compressive force acts on a building member or a plate element that constitutes a member, the partial sum of the 15th term of the Fourier series, that is, the formula (1-1) It is necessary to set N to about 15. In this case, the approximation function contains 30 trigonometric terms other than the undetermined coefficient _a0 .
As described above, the Fourier series requires a large number of terms for approximation convergence, so analysis using the Fourier series cannot be performed by manual calculation, and numerical calculation using a computer is required. Furthermore, the conventional method for estimating the degree of buckling stress has a technical problem that the estimation results cannot be formulated in a manner that is physically easy to understand.

The present invention has been made in view of the above problems. It is an object of the present invention to provide a buckling stress estimating device, a buckling stress estimating method, and a buckling stress estimating program that can formulate an estimation result in a physically easy-to-understand manner.

In order to solve the above problems, the present invention proposes the following means.
The apparatus for estimating the degree of buckling stress of the present invention spreads along the x-axis extending in the direction of the material axis and the y-axis extending in the width direction of the plate, and the plate element buckles when a compressive force acts in the direction along the x-axis. A buckling stress estimating device for estimating the buckling stress of the plate element, wherein the axis extending in the plate thickness direction of the plate element is the z-axis, and the direction along the x-axis of the plate element is directed toward the first end According to the above, when the half wavelength in the direction along the x-axis of the plate element that alternately displaces in the positive direction of the z-axis and the negative direction of the z-axis in a wave shape is a, along the z-axis a displacement estimating unit for estimating the out-of- _plane displacement _Ww of the plate element directed in the direction using equation (2); and a stress degree calculation unit that obtains based on the energy method.
However, _N is a natural number of 2 or more, a ₀ , an and _bn are undetermined coefficients, and _bw is the length of the plate element in the direction along the y-axis.

In addition, the buckling stress estimation method of the present invention spreads along the x-axis extending in the material axial direction and the y-axis extending in the plate width direction, and the compressive force acts in the direction along the x-axis to cause buckling. A buckling stress estimation method for estimating the buckling stress of a plate element, wherein the axis extending in the plate thickness direction of the plate element is the z-axis, and the first end of the plate element in the direction along the x-axis. Assuming that the half wavelength in the direction along the x-axis of the plate element displaced alternately in the positive direction of the z-axis and the negative direction of the z-axis in a wave-like manner toward the side is a, the z a displacement estimating step of estimating the out-of- _plane displacement _Ww of the plate element in the direction along the axis using equation (3); and a stress calculation step of obtaining the flexural stress based on the energy method.
However, _N is a natural number of 2 or more, a ₀ , an and _bn are undetermined coefficients, and _bw is the length of the plate element in the direction along the y-axis.

In addition, the buckling stress estimation program of the present invention spreads along the x-axis extending in the material axial direction and the y-axis extending in the plate width direction, and the compressive force acts in the direction along the x-axis to cause buckling. A buckling stress estimation program for an estimation device for estimating the buckling stress of a plate element, wherein the axis extending in the plate thickness direction of the plate element is the z-axis, and the direction along the x-axis of the plate element The half wavelength in the direction along the x-axis of the plate element alternately displaced in the positive direction of the z-axis and the negative direction of the z-axis toward the first end of the plate element is set to a. When the estimating device is oriented in the direction along the z-axis, a displacement estimating unit for estimating the out-of-plane displacement _Ww of the plate element by equation (4), and the out-of-plane displacement based on equation (4) The buckling stress degree of the plate element that produces _Ww is characterized by functioning as a stress degree calculation unit that obtains the degree of stress based on the energy method.
However, _N is a natural number of 2 or more, a ₀ , an and _bn are undetermined coefficients, and _bw is the length of the plate element in the direction along the y-axis.

The inventors have studied a plurality of functions that can estimate the out-of-plane displacement of the plate element on which the compressive force acts, with a smaller number of terms than the Fourier series while using trigonometric functions. As a result, it was found that the out-of-plane displacement _Ww of the plate element can be estimated using equation (2) with fewer terms than the Fourier series. The formulas (3) and (4) are the same as the formula (2).
Equation (2) can estimate the out-of-plane displacement _Ww of the plate element with a smaller number of terms than the Fourier series, so the calculations required for estimation can be performed more easily than when the Fourier series is used. Furthermore, when the number of terms required to estimate the out-of-plane displacement W _w with a certain accuracy is estimated using formula (2) than when estimating using the Fourier series of formula (1-1) Since there are fewer, the estimation result can be formulated physically easier to understand than the case of using the Fourier series.

In the buckling stress estimation device, the buckling stress estimation method, and the buckling stress estimation program, the N may be two.
According to these inventions, in equation (2), undetermined coefficients a ₁ , a ₂ , b ₁ , and b ₂ are obtained with the half wavelength a treated as a constant, and then the half wavelength a is treated as a variable. , the half wavelength a in the equation (2) is obtained. Since the number of undetermined coefficients to be obtained at one time is four or less when obtaining the undetermined coefficient a ₁ etc. in the equation (2), the undetermined coefficient a ₁ , a ₂ , b ₁ , b ₂ , etc. can be easily obtained.

Further, in the buckling stress estimation device, the stress calculation unit uses the equations (5) and (6) to calculate the minimum positive value of the buckling stress σ _cr according to the equation (7). The buckling stress magnitude σ _cr may be determined based on the a _n , b _n and the half wavelength a, which are real numbers giving values.
where E is the Young's modulus of the plate element, _v is the Poisson's ratio of the plate element, and tw is the thickness of the plate element.

Further, in the method for estimating the degree of buckling stress, in the stress degree calculation step, the minimum positive value of the buckling stress degree σ _cr according to the expression (10) is calculated using the expressions (8) and (9). The buckling stress magnitude σ _cr may be determined based on the a _n , b _n and the half wavelength a, which are real numbers giving values.
where E is the Young's modulus of the plate element, _v is the Poisson's ratio of the plate element, and tw is the thickness of the plate element.

Further, in the buckling stress estimation program, the stress calculation unit uses the equations (11) and (12) to calculate the minimum positive value of the buckling stress σ _cr according to the equation (13). The buckling stress magnitude σ _cr may be determined based on the a _n , b _n and the half wavelength a, which are real numbers giving values.
where E is the Young's modulus of the plate element, _v is the Poisson's ratio of the plate element, and tw is the thickness of the plate element.

According to these inventions, when obtaining the buckling stress degree σ _cr based on the energy method, it is possible to obtain the buckling stress degree σ _cr accurately using a mathematical formula.

According to the apparatus for estimating buckling stress, the method for estimating buckling stress, and the program for estimating buckling stress of the present invention, the buckling stress of the plate element where the compressive force acts more than when using the Fourier series Calculations necessary for estimating degrees can be easily performed, and estimation results can be formulated in a physically easy-to-understand manner.

1 is a perspective view of a building having a web to which a buckling stress estimation device according to an embodiment of the present invention is applied; FIG. It is a figure which shows the outline|summary of the estimation apparatus of the same buckling stress degree. FIG. 10 is a perspective view schematically showing a state in which the web is buckled under compressive force when the web is sufficiently long in the direction along the x-axis. 4 is an enlarged perspective view of the web of FIG. 3 for half a wavelength along the x-axis; FIG. FIG. 5 is a diagram showing an example of estimation results of the out-of-plane displacement of the web using the y-axis coordinates when a compressive force acts on the web. FIG. 4 is a diagram showing an example of the change in buckling coefficient with respect to N when a compressive force acts on the web;

An embodiment of an apparatus for estimating buckling stress according to the present invention will be described below with reference to FIGS. 1 to 6. FIG.
This buckling stress estimation device (estimation device, hereinafter simply referred to as an estimation device) is installed in a building 1 shown in FIG. Used to estimate flexural strength. The H-section steel 10 includes an upper flange (first flange) 11, a lower flange (second flange) 12, and a web 13 connecting the upper flange 11 and the lower flange 12 to each other. In FIG. 1, a floor slab 20, which will be described later, is indicated by a chain double-dashed line.
The upper flange 11 and lower flange 12 of the H-section steel 10 are also plate elements. A plate element is an elastic element. An elastic element is an element that does not consider material nonlinearity.

The H-section steel 10 extends, for example, along a horizontal plane. The upper flange 11 is formed in a flat plate shape, and arranged such that the thickness direction of the upper flange 11 extends along the vertical direction. The lower flange 12 is formed in a flat plate shape and arranged below the upper flange 11 . The lower flange 12 is arranged such that the thickness direction of the lower flange 12 extends along the vertical direction.
The web 13 is formed in a rectangular flat plate shape when viewed in the thickness direction of the web 13 . The web 13 is arranged so that the thickness direction of the web 13 is along the horizontal plane. The web 13 connects the widthwise center of the lower surface of the upper flange 11 and the widthwise center of the upper surface of the lower flange 12 .

Longitudinal ends of the H-section steel 10 are fixed to columns 15 and the like. The H-section steel 10 supports the floor slab 20 from below. An upper flange 11 of the H-section steel 10 is provided with a shear connector 21 such as a headed stud. A shear connector 21 is embedded in the floor slab 20 .
The building 1 is used by installing equipment (not shown) on the floor slab 20 .

FIG. 2 shows the estimation device 50 of this embodiment. The estimation device 50 is a computer, and includes a CPU (Central Processing Unit) 51, a main storage device 55, an auxiliary storage device 60, an input/output interface (IO/I/F) 65, and a recording/reproducing device 70. I have. The CPU 51 , main storage device 55 , auxiliary storage device 60 , input/output interface 65 , and recording/reproducing device 70 are connected to each other via a bus 75 .
The main storage device 55 is a RAM (Random Access Memory) or the like that serves as a work area for the CPU 51 or the like.
The input/output interface 65 is connected to an input device 66 such as a keyboard and mouse, and a display device 67 .
A recording/reproducing device 70 records and reproduces data on a recording medium 71 such as a USB (Universal Serial Bus) memory.

The auxiliary storage device 60 is a hard disk drive or the like that stores various data, programs, and the like. The auxiliary storage device 60 stores a buckling stress degree estimation program (hereinafter referred to simply as an estimation program) 61 for causing the computer to function as the estimation device 50, various programs such as an OS program, and the like. Various programs including the estimation program 61 are loaded from the recording medium 71 to the auxiliary storage device 60 via the recording/reproducing device 70 . The estimation program 61 and the like are stored in the recording medium 71 .
Note that these programs may be loaded into the auxiliary storage device 60 from an external device via a disk-type recording medium such as a CD or DVD, or a communication device (not shown).

The CPU 51 executes various arithmetic processing. The CPU 51 functionally includes a displacement estimator 52 for estimating the out-of-plane displacement of the web 13 and a stress calculator 53 for obtaining the buckling stress of the web 13 that causes the out-of-plane displacement based on the energy method. I have. The displacement estimation unit 52 and the stress degree calculation unit 53, which are functional components of the CPU 51, function when the CPU 51 executes the estimation program 61 and the like stored in the auxiliary storage device 60. FIG. The estimation program 61 and the like are programs for the estimation device 50 . The estimation program 61 causes the estimation device 50 to function as the displacement estimation section 52 and the stress degree calculation section 53 .

As shown in FIG. 3, the displacement estimating unit 52 and the stress degree calculating unit 53 of the estimating device 50 of the present embodiment set the positional coordinates of the web 13 to right-handed orthogonal Recognize based on coordinate system. In addition, in FIG. 3 and FIG. 5 described later, the out-of-plane displacement of the web 13 is shown larger than the estimated value.
The web 13 spreads along the x-axis extending in the axial direction of the web 13 and the y-axis extending in the sheet width direction of the web 13 . Of the x-axis and the y-axis, for example, the x-axis extends along the horizontal plane and the y-axis extends along the vertical direction. Let the axis extending in the plate thickness direction of the web 13 be the z-axis. The x-, y-, and z-axes are orthogonal to each other. When viewed in the direction along the z-axis (hereinafter referred to as the z-axis direction), the web 13 has sides extending in the direction along the x-axis (hereinafter referred to as the x-axis direction) and sides extending in the direction along the y-axis (hereinafter referred to as the y-axis direction). each have sides extending in the axial direction). It is assumed that the web 13 is sufficiently long in the x-axis direction. The term "web 13 sufficiently long in the x-axis direction" as used herein means that the boundary conditions of surfaces 13a arranged at each end of the web 13 in the x-axis direction and extending in the y-axis direction (hereinafter referred to as end surfaces in the x-axis direction) are This means that the web 13 has such a length that its influence on buckling deformation can be ignored.
The out-of-plane displacement of the web 13 is the displacement of the web 13 in the z-axis direction.

When a compressive force F1 acts in the x-axis direction on the end faces 13a of the web 13 in the x-axis direction, the web 13 may buckle (locally buckle). In this example, the compressive force F1 acting on each end face of the web 13 is a uniformly distributed load acting on the end face of the web 13 in the x-axis direction over the entire length in the y-axis direction. The magnitude of the compressive force F1 acting on each end surface of the web 13 is equal to each other.

In this case, the web 13 alternately displaces in the positive direction of the z-axis and in the negative direction of the z-axis toward the first end of the web 13 in the x-axis direction (one of the end faces 13a in the x-axis direction). As a whole, it may be displaced in a wave shape corresponding to a plurality of wavelengths (hereinafter referred to as a wave shape in the x-axis direction).
In one wavelength portion of the web 13 displaced along the x-axis, the second end opposite to the first end in the x-axis direction is the origin of the x-axis, and from this second end to the first end in the x-axis direction Let the direction be the positive direction of the x-axis.
Let the origin of the y-axis be the lower end (first end) of the web 13 in the y-axis direction. The positive direction of the y-axis is assumed to be upward (direction from the first end to the second end).
Let the origin of the z-axis be the center of the web 13 in the z-axis direction (the center in the thickness direction). The positive direction of the z-axis is defined as a direction forming a right-handed orthogonal coordinate system with respect to the positive directions of the x-axis and the positive direction of the y-axis.

It should be noted that SI units such as "m" are preferably used for the units of length and the like to be described below. Let bw be the length of the web 13 in the y-axis direction, and _tw be the thickness of the web 13 (the length in the _z -axis direction). Let E be the Young's modulus of the web 13 and ν be the Poisson's ratio of the web 13 .
The estimating device 50 estimates the degree of buckling stress of the web 13 when the compressive force F1 acts on the web 13 to which the flanges 11 and 12 are not connected and the web 13 buckles. Note that the estimating device 50 can be used when the thickness of the flanges 11 and 12 in the H-section steel 10 is negligibly thin, the rigidity of the flanges 11 and 12 is negligibly low, and so on. The flexural strength may be estimated as the buckling strength of the web 13 .

The displacement estimator 52 and the stress calculator 53 of the estimation device 50 of the present embodiment make the following assumptions 1 to 6 when estimating the buckling stress of the web 13 .
1. The thickness of the web 13 is thin, and the thickness of the web 13 is short compared to the length of the web 13 in the x-axis direction and the length in the y-axis direction.
2. The deflection (out-of-plane displacement due to buckling) of the web 13 is small and smaller than the thickness of the web 13 .
3. The central plane in the thickness direction of the web 13 does not expand or contract due to the bending of the web 13 and maintains a neutral plane.
4. The cross-section of the web 13 makes the assumption of planarity with respect to bending.
5. The material of web 13 is homogeneous and isotropic.
6. The displacement when an external force acts on the web 13 follows Hooke's law.

As shown in FIG. 3, when the compressive force F1 acts on the web 13, the web 13 may be displaced in a wavy shape in the x-axis direction. Assuming that the out-of-plane displacement in the x-axis direction of the web 13 undulatingly displaced in the x-axis direction is represented by the equation sin(πx/a), x The wavelength in the axial direction becomes 2a. A half wavelength (half the length of the wavelength) in the x-axis direction is a.
FIG. 4 is a diagram showing the out-of-plane displacement of the portion whose length in the x-axis direction is half the wavelength a in the central portion in the x-axis direction of the web 13 .

When the compressive force F1 is acting on the web 13, conventionally, the Fourier series represented by the formula (1-1) is used to apply the The out-of-plane displacement _ww of the web 13 was estimated.
Although the inventors used trigonometric functions, the out-of-plane displacement of the web 13 toward the z-axis direction (first A plurality of functions capable of estimating the out-of-plane displacement) _w were examined. Note that the out-of-plane displacement _ww is a function of the y-axis coordinates, not of the x-axis coordinates (it is a function at some coordinate on the x-axis).
As a result, when a compressive force F1 acts on the web 13, the out-of-plane displacement _ww of the web 13 in the z-axis direction at a certain coordinate on the x-axis and a certain coordinate on the y-axis is expressed by Fourier We found that it can be estimated with fewer terms than the series. However, _N is a natural number of 2 or more, and a ₀ , an and b _n are undetermined coefficients.
Equation (16) includes trigonometric terms using the power function of the y-axis coordinates. cos(π(2y/b _w −1) ⁿ ) and cos(π/2×(2y/b _w −1) ⁿ ) are the basis. Equation (16) can be preferably used for estimating out-of-plane deformation of plate elements such as the web 13 .

On the other hand, the out-of-plane displacement (second out-of-plane displacement) Ww of the web 13 in the _z -axis direction at an arbitrary coordinate on the x-axis and a certain coordinate on the y-axis is estimated by equation (18). Equation (18) is based on the assumption that the out-of-plane displacement in the z-axis direction at a coordinate on the x-axis of the web 13 is represented by the equation sin(πx/a), as described above. The out-of-plane displacement _Ww is a function of the y-axis and x-axis coordinates and is used to estimate the buckling stress of the web 13 . The out-of-plane displacement _Ww in FIG. 3 is obtained by repeating the out-of-plane displacement _Ww in FIG. 4 in the x-axis direction. Since estimating the buckling stress of the web 13 in FIG. 3 is synonymous with estimating the _buckling stress of the web 13 in FIG. is estimated.
By substituting the formula (16) into the formula (18), the formula (19) is obtained. When N is 2 in equation (19), it can be transformed into equation (20).

For example, before the compressive force _{F1 acts} on the web 13 in FIG. The coordinates are 0. After the compressive force F1 acts on the web 13, the out-of-plane displacement Ww in this estimation target portion is estimated by equation (19-1), and the _z -axis coordinates at this time are 0 and (19-1) It is a value obtained by adding the out-of-plane displacement _Ww estimated by the formula. That is, after the compressive force F1 acts, the estimation target portion is estimated to be arranged at a position where the x-axis coordinate is x ₀ , the y-axis coordinate is y ₀ , and the _z -axis coordinate is Ww. be.

The displacement estimator 52 estimates the out-of-plane displacement _Ww of the web 13 according to the equation (19).

The web 13 is arranged in a closed section in which the _coordinates of the x-axis are between 0 and a, and in which the coordinates of the y-axis are between 0 and bw. At this time, the strain energy U generated within the web 13 due to buckling deformation is expressed by equation (21), and the external force potential V given by the compressive force F1 is expressed by equation (22).
The total potential energy Π is represented by the equation (23) as the sum of the strain energy U and the external force potential V.

Here, _Dw is the plate stiffness of the web 13 and is represented by the formula (25). _σw is the stress function of the web 13 on which the compressive force F1 acts, and is expressed by the equation (26). However, σ _cr is the degree of buckling stress of the web 13, and compression is positive. The function _δw is a function that expresses the deformation of the web 13 in the x-axis direction when buckling occurs, and is expressed by Equation (27). is the displacement of the web 13 in the x-axis direction.
For example, (∂ ² W _w /∂x ² ) means a function obtained by second-order partial differentiation of W _w with respect to x, and (∂ ² W _w /(∂ x ∂ y)) means that W _w Denotes the function with partial differentiation.

The stress calculation unit 53 obtains the buckling stress σ _cr of the web 13 that causes the out-of-plane displacement W _w based on the equation (19) based on the energy method. The stress calculation unit 53 calculates the buckling stress σ _cr of the web 13 using the equations (30) and (31) to give the minimum positive value to the buckling stress σ _cr obtained by the equation (32). A buckling stress degree σ _cr is obtained based on the real numbers a _n , b _n and the half wavelength a.
However, the buckling coefficient k is obtained by the formula (33).

Specifically, in a state where the half wavelength a is treated as a constant, simultaneous equations representing that the function obtained by partially differentiating the total potential energy Π with the undetermined coefficients a _n and b _n are equal to 0, are real numbers _Obtain an and _bn . When there are multiple pairs of an and _bn that are solutions to the simultaneous equations, among the multiple pairs of an and _bn , the _buckling stress degree σ _cr according to equation (32) is _set to the smallest positive value. Based on the _given set of an and _bn (by _substituting the set of an and _bn into the equation (32)), the degree of buckling stress σ _cr is obtained. Next, the half-wavelength a is treated as a variable, and the half-wavelength a is obtained from the equation representing that the function obtained by partially differentiating the total potential energy _Π from which the set of an and _bn is obtained is equal to 0 by the half-wavelength a. demand. The degree of buckling stress σ _cr obtained based on the set of an and _bn obtained as described above and the half wavelength a is the degree of _{buckling stress σ cr to} be obtained.
When there is only one set of undetermined coefficients a _n and b _n that is the solution of the simultaneous equations, the set of a _n and b _n gives the minimum positive value to the buckling stress σ _cr according to equation (32). If _given , the degree of buckling stress σ _cr is determined based on the set of an and b _n . Next, the half wavelength a is treated as a variable, and the degree of buckling stress σ _cr is obtained as described above.

In the method of estimating the degree of buckling stress of the present embodiment (hereinafter simply referred to as the estimation method), the displacement estimation step of estimating the out-of-plane displacement _Ww of the web 13 by equation (19) and the out-of-plane displacement _Ww A stress calculation step is performed to obtain the buckling stress σ _cr of the web 13 that causes
The compressive force F1 acting on the web 13 when the buckling stress σ _cr occurs in the web 13 can be determined from the known equations for plate elements. Buckling of the web 13 can be suppressed by making the compressive force F1 acting on the web 13 smaller than the compressive force F1.

[Evaluation of estimation accuracy of out-of-plane displacement and buckling stress]
(Comparison by out-of-plane displacement)
FIG. 5 shows the results of estimating the out-of-plane displacement _Ww of the web 13 according to the y-axis coordinates when the compressive force F1 acts on the web 13. As shown in FIG. A surface 13b arranged at each end of the web 13 in the y-axis direction and extending in the x-axis direction (hereinafter referred to as an end face in the y-axis direction) is fixedly supported (non-rotatably supported, rigidly joined). This out-of-plane displacement _Ww is the out-of-plane displacement in the cross section taken along the cutting line AA in FIG.
In FIG. 5 , the horizontal axis represents the y-axis coordinate against the length _bw of the web 13 and the vertical axis represents the estimated out-of-plane displacement _Ww of the web 13 . 0 on the horizontal axis represents the bottom edge of the web 13 and 1 on the horizontal axis represents the top edge of the web 13 . On the vertical axis, the value of the out-of-plane displacement _Ww is made dimensionless by the maximum value of the out-of-plane displacement _Ww .

In FIG. 5, a solid line L1 represents the analysis result by FEM (Finite Element Method). The dotted line L2 represents the result of estimation using equation (19) of the present embodiment (where N is 2), and the dashed-dotted line L3 represents the result of estimation using equation (1-1) (where N is 15). Since the lines L2 and L3 substantially coincide with the line L1, the lines L2 and L3 overlap the line L1.
N in equation (19) of this embodiment, which is necessary to approximate the analysis result by FEM with the same accuracy as when N is 15 in the Fourier series of equation (1-1), is 2 I found out.

(Comparison by buckling stress)
The buckling stress σ _cr of the web 13 was determined by the formula (1-1) and the formula (19) of the present embodiment, and the results of comparison by buckling coefficient are shown in Table 1 and FIG.

N is set to 2 when the formula (19) of this embodiment is used. On the other hand, N was changed to 4, 5, 7, .
In Table 1, the buckling coefficient when the end face 13b in the y-axis direction of the web 13 is simply supported converges to 4.0 as the out-of-plane displacement _Ww and the buckling stress degree _σcr increase in estimation accuracy. (see line L11 in FIG. 6). The buckling coefficient when the end surface 13b in the y-axis direction of the web 13 is fixedly supported converges to 6.98 as the out-of-plane displacement _{Ww and the buckling stress degree σcr are} estimated with higher accuracy (Fig. 6 See line L12 in the middle).

In FIG. 6, the horizontal axis represents N, and the vertical axis represents buckling coefficient k. Outlined circles represent estimation results in the case of simple support using equation (19) of this embodiment. Black circles represent estimation results in the case of simple support using the Fourier series of formula (1-1). Outlined square marks represent estimation results in the case of fixed support using equation (19) of the present embodiment. Black squares represent estimation results in the case of fixed support using the Fourier series of formula (1-1).
The accuracy when the buckling coefficient k is evaluated using the equation (19) of this embodiment is the accuracy when N is 15 in the Fourier series of the equation (1-1) for both simple support and fixed support. was found to be the same.

As described above, according to the estimation device 50, the estimation method, and the estimation program 61 of the present embodiment, the inventors can use a trigonometric function, but have fewer terms than the Fourier series of formula (1-1) A number of functions that can estimate the out-of-plane displacement of the web 13 on which the compressive force F1 acts were investigated. As a result, it was found that the out-of-plane displacement _Ww of the web 13 can be estimated using equation (19) with fewer terms than the Fourier series.
Equation (19) can estimate the out-of-plane displacement _Ww of the web 13 with a smaller number of terms than the Fourier series of Equation (1-1). It can be carried out. Furthermore, when the number of terms required to estimate the out-of-plane displacement W _w with a certain accuracy is estimated using Equation (19) than when estimating using the Fourier series of Equation (1-1) Since there are fewer, the estimation result can be formulated physically easier to understand than the case of using the Fourier series.

Using equation (19), the out-of-plane displacement _Ww of the web 13 can be quickly approximated. That is, it is possible to quickly approximate a buckling waveform when a member of a building or a plate element constituting a member undergoing in-plane bending undergoes local buckling.

Also, N is two. First, in equation (19), undetermined coefficients a ₁ , a ₂ , b ₁ , and b ₂ are obtained with the half wavelength a treated as a constant. Obtain the wavelength a. The number of undetermined coefficients to be determined at _one time is four or less when determining the undetermined coefficients a1 and the like in the equation (19). Therefore, the undetermined coefficients a ₁ , a ₂ , b ₁ , b ₂ , etc. can be easily obtained using equations of order 4 or less whose solution formulas are known.
Further, the stress degree calculation unit 53 (in the stress degree calculation step) uses the equations (30) and (31) to obtain the minimum positive value for the buckling stress degree σ _cr according to the equation (32). A buckling stress degree σ _cr is obtained based on certain a _n , b _n and half wavelength a. When obtaining the buckling stress degree σ _cr based on the energy method, the buckling stress degree σ _cr can be accurately obtained using the formulas (30) to (32).

As described above, one embodiment of the present invention has been described in detail with reference to the drawings, but the specific configuration is not limited to this embodiment, and the configuration can be changed, combined, or deleted without departing from the scope of the present invention. etc. are also included.
For example, in the above embodiment, N in formula (19) may be 3 or more.
When obtaining the buckling stress degree σ _cr based on the energy method, the equations (30) to (32) need not be used.
In this embodiment, the estimating device 50 estimates the buckling stress of the web 13, which is a plate element, but the plate element is not limited to the web 13, and may be the upper flange 11, the lower flange 12, or the like.

13 web (plate element)
50 Buckling stress estimation device 52 Displacement estimator 53 Stress calculator 61 Buckling stress estimation program F1 Compressive force

Claims (9)

buckling stress for estimating the buckling stress of a plate element that spreads along the x-axis extending in the material axial direction and the y-axis extending in the width direction of the plate and buckles when a compressive force acts in the direction along the x-axis An estimator for
Let the axis extending in the plate thickness direction of the plate element be the z-axis,
Along the x-axis, the plate element displaces alternately in a wave-like manner in the positive direction of the z-axis and the negative direction of the z-axis toward the first end of the plate element in the direction along the x-axis. When the half wavelength in the direction is a,
a displacement estimating unit for estimating the out-of-plane displacement Ww of the plate element in the direction along the _z -axis using equation (1);
(1) a stress calculation unit that obtains the buckling stress of the plate element that causes the out-of-plane displacement _Ww based on the formula based on the energy method;
A device for estimating buckling stress, comprising:
However, _N is a natural number of 2 or more, a ₀ , an and _bn are undetermined coefficients, and _bw is the length of the plate element in the direction along the y-axis.

2. The buckling stress estimation device according to claim 1, wherein said N is two. The _stress degree calculation unit uses the formulas (2) and (3) to calculate the a _n , b _n and 3. The buckling stress degree estimating device according to claim 1, wherein the buckling stress degree _σcr is obtained based on the half wavelength a.
where E is the Young's modulus of the plate element, _v is the Poisson's ratio of the plate element, and tw is the thickness of the plate element.

buckling stress for estimating the buckling stress of a plate element that spreads along the x-axis extending in the material axial direction and the y-axis extending in the width direction of the plate and buckles when a compressive force acts in the direction along the x-axis A method for estimating the
Let the axis extending in the plate thickness direction of the plate element be the z-axis,
Along the x-axis, the plate element alternately displaces in the positive direction of the z-axis and the negative direction of the z-axis in a wave-like manner toward the first end side of the plate element in the direction along the x-axis. When the half wavelength in the direction along is a,
a displacement estimating step of estimating the out-of-plane displacement Ww of the plate element in the direction along the _z -axis using equation (5);
(5) a stress degree calculation step of obtaining the buckling stress degree of the plate element that causes the out-of-plane displacement _Ww based on the equation based on the energy method;
method for estimating buckling stress.
However, _N is a natural number of 2 or more, a ₀ , an and _bn are undetermined coefficients, and _bw is the length of the plate element in the direction along the y-axis.

5. The buckling stress estimation method according to claim 4, wherein said N is two. In the _stress degree calculation step, using formulas (6) and (7), the a _n , b _n and 6. The buckling stress degree estimation method according to claim 4 or 5, wherein the buckling stress degree _σcr is obtained based on the half wavelength a.
where E is the Young's modulus of the plate element, _v is the Poisson's ratio of the plate element, and tw is the thickness of the plate element.

An estimating device for estimating the degree of buckling stress of a plate element that spreads along the x-axis extending in the material axial direction and the y-axis extending in the width direction of the plate and buckles when a compressive force acts in the direction along the x-axis A buckling stress estimation program comprising:
Let the axis extending in the plate thickness direction of the plate element be the z-axis,
Along the x-axis, the plate element displaces alternately in a wave-like manner in the positive direction of the z-axis and the negative direction of the z-axis toward the first end of the plate element in the direction along the x-axis. When the half wavelength in the direction is a,
the estimation device,
a displacement estimating unit for estimating an out-of-plane displacement Ww of the plate element along the _z -axis using equation (9);
(9) a stress degree calculation unit that obtains the buckling stress degree of the plate element that causes the out-of-plane displacement _Ww based on equation (9) based on the energy method;
buckling stress estimation program that functions as
However, _N is a natural number of 2 or more, a ₀ , an and _bn are undetermined coefficients, and _bw is the length of the plate element in the direction along the y-axis.

8. The buckling stress estimation program according to claim 7, wherein said N is two. The stress degree calculation unit uses the _formulas (10) and (11) to calculate the a _n , b _n and 9. The buckling stress estimation program according to claim 7 or 8, wherein the buckling stress _σcr is calculated based on the half wavelength a.
where E is the Young's modulus of the plate element, _v is the Poisson's ratio of the plate element, and tw is the thickness of the plate element.

Images