Method, device, equipment and storage medium for constructing finite element model of steel structure

Abstract

The present application is applicable to the technical field of dynamic response analysis of steel structures, and provides a method, device, equipment and storage medium for constructing a finite element model of a steel structure. The method includes: in the process of finite element modeling for a high-rise steel structure building, a node deformation ring is constructed at the intersection of a column unit and a beam unit through multiple peripheral nodes and multiple connection units, and two node deformation rings are connected with the rest of the column unit and the rest of the beam unit through multiple column connection nodes and multiple beam connection nodes to obtain a three-dimensional node area model at the intersection of the beam unit and the column unit. By constructing a three-dimensional node area model that can analyze the force and deformation of the beam-column node area in the process of finite element modeling for a high-rise steel structure building, the nonlinear response results of the steel structure can be accurately calculated, and at the same time, the problem of low calculation efficiency can be avoided.

Description

Technical Field

[0001] The present application belongs to the technical field of dynamic response analysis of steel structures, and in particular, relates to a method, device, equipment and storage medium for constructing a finite element model of a steel structure. Background Art

[0002] The basic dynamic characteristics (modes, periods, etc.) of the finite element model of steel structures and the nonlinear deformation of the structure under strong loads are important indicators for measuring the accuracy and reliability of the finite element model.

[0003] At present, in the process of finite element modeling for steel structure high-rise buildings, the calculation efficiency will be significantly reduced due to the excessive complexity of the steel structure or the excessive number of structural components, and the beam-column node area has little effect on the linear elastic response performance of the entire steel structure. Therefore, only the main load-bearing components such as beams and columns are modeled, and the influence of the beam-column node area on the dynamic response performance of the structure is ignored to improve the calculation efficiency. However, the beam-column node area has a great influence on the nonlinear response performance of the entire steel structure, which will cause the response calculation results to be significantly different from the actual situation. Summary of the invention

[0004] The embodiments of the present application provide a method, device, equipment and storage medium for constructing a finite element model of a steel structure, which can solve the problem that in the process of finite element modeling of a steel structure high-rise building, usually only the main load-bearing components such as beams and columns are modeled, while the influence of the node area where the beams and columns intersect on the nonlinear response performance of the entire steel structure is ignored, resulting in a significant difference between the simulated structural response results and the actual situation.

[0005] In a first aspect, an embodiment of the present application provides a method for constructing a finite element model of a steel structure, the method comprising:

[0006] In the process of modeling a finite element model for a steel structure high-rise building, a node deformation ring is constructed at the intersection of a column unit and a beam unit through a plurality of peripheral nodes and a plurality of connection units, wherein the shape of the node deformation ring is a parallelogram, the connection unit is used to represent the edge of the intersection of the column unit and the beam unit, and the peripheral node is used to represent the vertex or midpoint of the intersection of the column unit and the beam unit;

[0007] The two node deformation rings are connected with the rest of the column unit and the rest of the beam unit through a plurality of column connection nodes and a plurality of beam connection nodes, so as to obtain a three-dimensional node area model at the intersection of the beam unit and the column unit; wherein the rest is the portion of the column unit or the beam unit excluding the intersection of the column unit and the beam unit.

[0008] In a possible implementation of the first aspect, the node deformation ring includes an X-coordinate node deformation ring and a Y-coordinate node deformation ring, the beam unit includes an X-coordinate beam unit and a Y-coordinate beam unit; the peripheral node includes an upper node, a lower node, a left node, a right node, and a plurality of corner nodes;

[0009] The method connects the two node deformation rings with the rest of the column unit and the rest of the beam unit through a plurality of column connection nodes and a plurality of beam connection nodes to obtain a three-dimensional node area model at the intersection of the beam unit and the column unit, including:

[0010] Connecting a first X-beam connection node to a left-side node of the X-coordinate node deformation ring, and connecting a second X-beam connection node to a right-side node of the X-coordinate node deformation ring, so that the X-coordinate node deformation ring is connected to the X-coordinate beam unit;

[0011] Connecting a first Y-beam connection node to a Y-ring left node of the Y-coordinate node deformation ring, and connecting a second Y-beam connection node to a Y-ring right node of the Y-coordinate node deformation ring, so that the Y-coordinate node deformation ring is connected to the Y-coordinate beam unit;

[0012] Connecting the first column connection node to the upper X-ring node of the X-coordinate node deformation ring and the upper Y-ring node of the Y-coordinate node deformation ring, and connecting the second column connection node to the lower X-ring node of the X-coordinate node deformation ring and the lower Y-ring node of the Y-coordinate node deformation ring, so as to connect the X-coordinate node deformation ring, the Y-coordinate node deformation ring and the rest of the column unit;

[0013] The translational degrees of freedom of the upper X-ring node of the X-coordinate node deformation ring, the upper Y-ring node of the Y-coordinate node deformation ring, the lower X-ring node of the X-coordinate node deformation ring, and the lower Y-ring node of the Y-coordinate node deformation ring are constrained to obtain the three-dimensional node area model.

[0014] In a possible implementation manner of the first aspect, the X-coordinate node deformation ring includes four corner point positions, namely, an X-upper-left corner point position, an X-upper-right corner point position, an X-lower-right corner point position, and an X-lower-left corner point position, and each corner point position includes two corner nodes whose positions overlap;

[0015] The method of constructing a node deformation ring at the intersection of a column unit and two beam units by using a plurality of peripheral nodes and a plurality of connection units comprises:

[0016] Constraining the translational degrees of freedom of the two corner nodes at each corner point position so that the displacements corresponding to the two corner nodes are the same;

[0017] The rotational freedom of the corner nodes of the X upper left corner point position, the X lower right corner point position and the X lower left corner point position is not constrained, so that the X upper left corner point position, the X lower right corner point position and the X lower left corner point position have a zero stiffness hinge function;

[0018] The two corner nodes at the upper right corner point are connected by an X zero-length unit in the connection unit, and a force-deformation constitutive model is applied to the X zero-length unit so that the X upper right corner point has a rotation spring function;

[0019] Through the multiple X-bar units in the connection unit, each corner node at the X-upper left corner point position, the X-upper right corner point position, the X-lower right corner point position and the X-lower left corner point position is connected in sequence with the X-ring upper node, the X-ring right node, the X-ring lower node and the X-ring left node to obtain the X-coordinate node deformation ring.

[0020] In a possible implementation manner of the first aspect, the Y coordinate node deformation ring includes four corner point positions, namely, a Y upper left corner point position, a Y upper right corner point position, a Y lower right corner point position, and a Y lower left corner point position, and each of the corner point positions includes two corner nodes whose positions overlap;

[0021] The method of constructing a node deformation ring at the intersection of a column unit and two beam units by using a plurality of peripheral nodes and a plurality of connection units comprises:

[0022] Constraining the translational degrees of freedom of the two corner nodes at each corner point position so that the displacements corresponding to the two corner nodes are the same;

[0023] The rotational freedom of the corner nodes of the Y upper left corner point position, the Y lower right corner point position and the Y lower left corner point position is not constrained, so that the Y upper left corner point position, the Y lower right corner point position and the Y lower left corner point position have a zero stiffness hinge function;

[0024] The two corner nodes at the upper right corner point are connected by a Y zero-length unit in the connection unit, and a force-deformation constitutive model is applied to the Y zero-length unit so that the Y upper right corner point has a rotation spring function;

[0025] Through the multiple Y rod units in the connecting unit, each corner node at the Y upper left corner point position, the Y upper right corner point position, the Y lower right corner point position and the Y lower left corner point position is connected with the upper node of the Y ring, the right node of the Y ring, the lower node of the Y ring and the left node of the Y ring in sequence to obtain the Y coordinate node deformation ring.

[0026] In a possible implementation manner of the first aspect, after constructing a node deformation ring on each beam unit by using a plurality of peripheral nodes and a plurality of connection units, the method further includes:

[0027] The length values ​​of the plurality of connection units are determined according to the height value of the beam unit, the length value of the column unit in the X direction, and the length value of the column unit in the Y direction.

[0028] In a possible implementation manner of the first aspect, the method further includes:

[0029] Using an elastic unit model, simulating a linear elastic section in the middle of the beam unit and the column unit;

[0030] The plastic hinge sections at both ends of the beam element and the column element are simulated by using a non-elastic unit model.

[0031] In a possible implementation of the first aspect, the stiffness of the connection unit is a first preset multiple of the stiffness of the beam unit and the stiffness of the column unit; the strength of the connection unit is a second preset multiple of the strength of the beam unit and the strength of the column unit.

[0032] In a second aspect, an embodiment of the present application provides a device for constructing a finite element model of a steel structure, comprising:

[0033] A node deformation ring construction module is used to construct a node deformation ring at the intersection of a column unit and a beam unit through a plurality of peripheral nodes and a plurality of connection units during finite element modeling of a steel structure high-rise building, wherein the shape of the node deformation ring is a parallelogram, the connection unit is used to represent the intersection of the column unit and the beam unit, and the peripheral node is used to represent the vertex or midpoint of the intersection of the column unit and the beam unit;

[0034] The beam-column node area construction module is used to connect the two node deformation rings with the rest of the column unit and the rest of the beam unit through multiple column connection nodes and multiple beam connection nodes to obtain a three-dimensional node area model at the intersection of the beam unit and the column unit; wherein the rest is the portion of the column unit or the beam unit excluding the intersection of the column unit and the beam unit.

[0035] In a third aspect, an embodiment of the present application provides a terminal device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the method for constructing a finite element model of a steel structure as described in any one of the above items when executing the computer program.

[0036] In a fourth aspect, an embodiment of the present application provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the method for constructing a finite element model of a steel structure described in any one of the above items is implemented.

[0037] In a fifth aspect, an embodiment of the present application provides a computer program product. When the computer program product is run on a terminal device, the terminal device executes the method for constructing a finite element model of a steel structure as described in any one of the above-mentioned first aspects.

[0038] Compared with the prior art, the embodiments of the present invention have the following beneficial effects:

[0039] The embodiment of the present application provides a method for constructing a finite element model of a steel structure. In the process of finite element modeling for a high-rise steel structure building, a node deformation ring is constructed at the intersection of a column unit and a beam unit through multiple peripheral nodes and multiple connection units, wherein the shape of the node deformation ring is a parallelogram, the connection unit is used to characterize the edge of the intersection of the column unit and the beam unit, and the peripheral node is used to characterize the vertex or midpoint of the intersection of the column unit and the beam unit; through multiple column connection nodes and multiple beam connection nodes, the two node deformation rings are connected with the rest of the column unit and the rest of the beam unit to obtain a three-dimensional node area model at the intersection of the beam unit and the column unit; wherein the rest is the part of the column unit or the beam unit except the intersection of the column unit and the beam unit. By constructing a three-dimensional node area model that can analyze the stress deformation of the beam-column node area in the process of finite element modeling for a high-rise steel structure building, the nonlinear response results of the steel structure can be accurately calculated, and at the same time, the problem of low calculation efficiency can be avoided. BRIEF DESCRIPTION OF THE DRAWINGS

[0040] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0041] Figure 1 It is a flow chart of a method for constructing a finite element model of a steel structure provided in one embodiment of the present application;

[0042] Figure 2 It is a structural schematic diagram of a beam-column node area provided by an embodiment of the present application;

[0043] Figure 3 is a schematic diagram of a three-dimensional node area model provided by an embodiment of the present application;

[0044] Figure 4 It is a decomposition schematic diagram of a three-dimensional node area model provided by an embodiment of the present application;

[0045] Figure 5 is a schematic diagram of a node deformation ring provided by an embodiment of the present application;

[0046] Figure 6 is a schematic diagram of a force-deformation constitutive model provided in one embodiment of the present application;

[0047] Figure 7 This is a schematic diagram of modeling a beam-column component provided by an embodiment of the present application;

[0048] Figure 8 is a schematic diagram of a mIMK model provided in an embodiment of the present application;

[0049] Fig. 9 It is a structural schematic diagram of a device for constructing a finite element model of a steel structure provided in one embodiment of the present application;

[0050] Fig.10 It is a structural schematic diagram of a comparison model provided in one embodiment of the present application;

[0051] Fig.11 It is a schematic diagram of the first two vibration modes of a structure under different model schemes provided by an embodiment of the present application;

[0052] Fig.12 It is a schematic diagram of the maximum inter-story displacement angle response of different model schemes under the action of preliminary nonlinear seismic loads provided by an embodiment of the present application;

[0053] Fig.13 It is a schematic diagram of the calculation time consumption of different model schemes under the action of preliminary nonlinear seismic loads provided by an embodiment of the present application;

[0054] Fig.14 It is a schematic diagram of the maximum inter-story displacement angle response of different model schemes under the action of a completely elastic-plastic earthquake load provided by an embodiment of the present application;

[0055] Fig.15 It is a structural diagram of a terminal device provided in one embodiment of the present application. DETAILED DESCRIPTION

[0056] In the following description, specific details such as specific system structures, technologies, etc. are provided for the purpose of illustration rather than limitation, so as to provide a thorough understanding of the embodiments of the present application. However, it should be clear to those skilled in the art that the present application may also be implemented in other embodiments without these specific details. In other cases, detailed descriptions of well-known systems, devices, circuits, and methods are omitted to prevent unnecessary details from obstructing the description of the present application.

[0057] It should be understood that when used in the present specification and the appended claims, the term "comprising" indicates the presence of described features, wholes, steps, operations, elements and / or components, but does not exclude the presence or addition of one or more other features, wholes, steps, operations, elements, components and / or combinations thereof.

[0058] It should also be understood that the term “and / or” used in the specification and appended claims refers to and includes any and all possible combinations of one or more of the associated listed items.

[0059] As used in the specification and appended claims of this application, the term "if" can be interpreted as "when" or "uponce" or "in response to determining" or "in response to detecting", depending on the context. Similarly, the phrase "if it is determined" or "if [described condition or event] is detected" can be interpreted as meaning "uponce it is determined" or "in response to determining" or "uponce [described condition or event] is detected" or "in response to detecting [described condition or event]", depending on the context.

[0060] In addition, in the description of the present application specification and the appended claims, the terms "first", "second", "third", etc. are only used to distinguish the descriptions and cannot be understood as indicating or implying relative importance.

[0061] References to "one embodiment" or "some embodiments" etc. described in the specification of this application mean that one or more embodiments of the present application include specific features, structures or characteristics described in conjunction with the embodiment. Therefore, the statements "in one embodiment", "in some embodiments", "in some other embodiments", "in some other embodiments", etc. that appear in different places in this specification do not necessarily refer to the same embodiment, but mean "one or more but not all embodiments", unless otherwise specifically emphasized in other ways. The terms "including", "comprising", "having" and their variations all mean "including but not limited to", unless otherwise specifically emphasized in other ways.

[0062] See also Figure 1 , Figure 11 is a flow chart of a method for constructing a finite element model of a steel structure provided in one embodiment of the present application. The method comprises:

[0063] S11. In the process of modeling a finite element model for a steel structure high-rise building, a node deformation ring is constructed at the intersection of a column unit and a beam unit through a plurality of peripheral nodes and a plurality of connection units, wherein the node deformation ring is in the shape of a parallelogram, the connection unit is used to represent the edge of the intersection of the column unit and the beam unit, and the peripheral node is used to represent the vertex or midpoint of the intersection of the column unit and the beam unit;

[0064] S12. Connect the two node deformation rings with the rest of the column unit and the rest of the beam unit through multiple column connection nodes and multiple beam connection nodes to obtain a three-dimensional node area model at the intersection of the beam unit and the column unit; wherein the rest is the portion of the column unit or the beam unit excluding the intersection of the column unit and the beam unit.

[0065] It should be noted that this method is mainly used in the finite element model of steel structure high-rise buildings, mainly for modeling the three-dimensional node area where beam components and column components intersect.

[0066] It should be noted that in the process of analyzing the structural mechanical behavior of steel high-rise buildings, it is usually simulated and analyzed by establishing a finite element model for the steel structure. The column element is the part that represents the column in the steel high-rise building in the finite element model, that is, the column component in the Z coordinate direction. The beam element is the part that represents the beam in the steel high-rise building in the finite element model, that is, the beam component in the X coordinate direction or the Y coordinate direction.

[0067] like Figure 2 As shown, Figure 2 This is a structural schematic diagram of a beam-column node area provided by an embodiment of the present application. The intersection of beam-column components in a steel structure high-rise building is the beam-column node area. However, in the process of modeling the steel structure high-rise building through the finite element model, the influence of the beam-column node area on the structural dynamic response is usually ignored. In this embodiment, a method for constructing a finite element model of the beam-column node area is proposed, such as Figure 3 As shown, Figure 3 is a schematic diagram of a three-dimensional node area model provided by an embodiment of the present application. The node deformation ring is a three-dimensional structure constructed at the intersection of the column unit and the beam unit in the finite element model, such as Figure 3 The parallelogram in the figure is the node deformation ring. The connection element is used to connect the peripheral nodes to form the edge of the node deformation ring, which represents the boundary of the intersection of the column element and the beam element. Figure 3 The peripheral nodes are the vertices or midpoints of the node deformation loop that represent the intersection of the column element and the beam element, such as Figure 3 The four vertices and four midpoints of the parallelogram in the 3D node area are the intersection of the beam element and the column element. In the modeling process, the node deformation ring is used as the intersection of the beam element and the column element, which is equivalent to splitting the beam element and the column element, splitting the intersection of the two, and then connecting the intersection with the rest of the intersection through the column connection node and the beam connection node. Among them, the column connection node and the beam connection node represent the nodes used to connect the two node deformation rings with the rest of the column element and the rest of the beam element.

[0068] Specifically, simulation models for three-dimensional node areas usually include octahedral models, scissor-shaped models, and parallelogram models. After many simulation tests, it was found that it is difficult to arrange zero-stiffness hinges and rotation springs on the octahedral node area model, and the expected deformation target cannot be achieved; in the scissor-shaped node area model, the node area boundary and the connection position of the beam-column unit cannot maintain a right angle, which is different from the actual node area deformation; the parallelogram node area model can better simulate the deformation of the actual three-dimensional node area. Therefore, in this embodiment, in order to more accurately simulate the deformation behavior of the intersection area of ​​the beam unit and the column unit and improve the accuracy of the model, the parallelogram model is used as the modeling model of the three-dimensional node area model, and the shape of the node deformation ring is set to a parallelogram. Figure 3 As shown in the figure, in the three-dimensional node area model, the beam unit and the column unit are respectively connected to the peripheral nodes of the node area model to transmit the axial force and bending moment. Through the three zero-stiffness hinges and one rotation spring arranged at the corner points of the parallelogram, the model can produce a certain shear deformation under the action of external force, which can better simulate the deformation behavior of the intersection area of ​​the beam unit and the column unit.

[0069] It should be noted that, at present, general finite element software is usually used for structural calculation and verification. If the structure is too complex or the number of structural components is too large, the calculation efficiency will be significantly reduced. Compared with the main components of the structure (beam components and column components), the influence of the beam-column node area on the overall structure is small under the linear elastic condition. Therefore, the node area is usually ignored to improve the calculation efficiency. In this embodiment, in the process of modeling the finite element model of the steel structure high-rise building, the overall influence of the beam-column node area on the structure will be considered. Therefore, the open source finite element software OpenSEES is used for simulation and analysis. Compared with general finite element software, the simulation analysis speed is faster and the degree of freedom is higher.

[0070] It can be understood that the method for constructing a finite element model of a steel structure provided in an embodiment of the present application, in the process of finite element modeling for a steel structure high-rise building, a node deformation ring is constructed at the intersection of a column unit and a beam unit through multiple peripheral nodes and multiple connection units, wherein the shape of the node deformation ring is a parallelogram, the connection unit is used to characterize the edge of the intersection of the column unit and the beam unit, and the peripheral node is used to characterize the vertex or midpoint of the intersection of the column unit and the beam unit; through multiple column connection nodes and multiple beam connection nodes, the two node deformation rings are connected with the rest of the column unit and the rest of the beam unit to obtain a three-dimensional node area model at the intersection of the beam unit and the column unit; wherein the rest is the part of the column unit or the beam unit except the intersection of the column unit and the beam unit. By constructing a three-dimensional node area model that can analyze the stress deformation of the beam-column node area in the process of finite element modeling for a steel structure high-rise building, the nonlinear response results of the steel structure can be accurately calculated, and at the same time, the problem of low calculation efficiency can be avoided.

[0071] In a possible implementation, the node deformation ring includes an X-coordinate node deformation ring and a Y-coordinate node deformation ring, the beam unit includes an X-coordinate beam unit and a Y-coordinate beam unit; the peripheral node includes an upper node, a lower node, a left node, a right node and a plurality of corner nodes;

[0072] Through multiple column connection nodes and multiple beam connection nodes, the two node deformation rings are connected with the rest of the column element and the rest of the beam element to obtain a three-dimensional node area model at the intersection of the beam element and the column element, including:

[0073] Connecting the first X-beam connection node to the X-ring left node of the X-coordinate node deformation ring, and connecting the second X-beam connection node to the X-ring right node of the X-coordinate node deformation ring, so that the X-coordinate node deformation ring is connected to the X-coordinate beam unit;

[0074] Connecting the first Y beam connection node to the Y ring left node of the Y coordinate node deformation ring, and connecting the second Y beam connection node to the Y ring right node of the Y coordinate node deformation ring, so that the Y coordinate node deformation ring is connected to the Y coordinate beam unit;

[0075] Connecting the first column connection node to the upper X-ring node of the X-coordinate node deformation ring and the upper Y-ring node of the Y-coordinate node deformation ring, and connecting the second column connection node to the lower X-ring node of the X-coordinate node deformation ring and the lower Y-ring node of the Y-coordinate node deformation ring, so as to connect the X-coordinate node deformation ring, the Y-coordinate node deformation ring and the rest of the column unit;

[0076] The translational degrees of freedom of the upper X-ring node of the X-coordinate node deformation ring, the upper Y-ring node of the Y-coordinate node deformation ring, the lower X-ring node of the X-coordinate node deformation ring, and the lower Y-ring node of the Y-coordinate node deformation ring are constrained to obtain a three-dimensional node area model.

[0077] It should be noted that if Figure 4 As shown, Figure 4 This is a decomposition schematic diagram of a three-dimensional node area model provided by an embodiment of the present application. The connection between the three-dimensional node area model and the beam unit and the column unit can be decomposed into three parts, namely the X-coordinate node deformation ring, the Y-coordinate node deformation ring and the column unit, wherein the X-coordinate node deformation ring and the Y-coordinate node deformation ring are the main deformation parts of the three-dimensional node area model, which simulate the shear deformation of the beam-column node area in two orthogonal directions through two intersecting parallelograms. The X-coordinate node deformation ring and the Y-coordinate node deformation ring and the column unit are connected through two nodes at the top and bottom of the deformation ring, and the X-coordinate beam unit and the Y-coordinate beam unit are connected through two nodes on both sides of the deformation ring.

[0078] Specifically, Figure 4 In the three-dimensional node zone model, there are two node deformation rings, namely, the X-coordinate node deformation ring and the Y-coordinate node deformation ring. The X-coordinate node deformation ring and the Y-coordinate node deformation ring represent the node deformation rings constructed in the X direction and the Y direction in the standard coordinate system, respectively. The X-coordinate node deformation ring and the Y-coordinate node deformation ring are used to simulate the deformation behavior of the intersection of the beam element and the column element in their respective directions. Accordingly, there are two beam elements, namely, the X-coordinate beam element and the Y-coordinate beam element. The X-coordinate beam element and the Y-coordinate beam element represent the beam components constructed in the X direction and the Y direction, respectively. The column element represents the column component constructed in the Z direction.

[0079] like Figure 5 As shown, Figure 5 This is a schematic diagram of a node deformation ring provided in an embodiment of the present application. Figure 5 (5-a) in the figure shows the schematic diagram of the X-coordinate node deformation ring. Figure 5 (5-b) in FIG. 5 is a schematic diagram of the Y-coordinate node deformation ring. The peripheral nodes include the upper node, the lower node, the left node, the right node and multiple corner nodes, which are used to represent the nodes at different positions in the node deformation ring. The upper node is the node located at the upper midpoint of the node deformation ring, such as the upper node of the X ring ( Figure 5 The lower node is the node located at the lower midpoint of the node deformation ring, such as the lower node of the X ring ( Figure 5 Node 09 in (5-a) and the lower node of the Y ring ( Figure 5 The left node is the node located at the midpoint of the left side of the node deformation ring, such as the left node of the X ring ( Figure 5 The right node is the node located at the midpoint of the right side of the node deformation ring, such as the right node of the X ring ( Figure 5 The node 06 in (5-a) and the node on the right side of the Y ring (node ​​26 in (5-b)); the corner nodes are the nodes located at the four vertices of the node deformation ring, such as Figure 5 In (5-a), nodes 01, 02, 04, 05, 07, 08, 10, 11, and Figure 5 In (5-b), nodes 21, 22, 24, 25, 27, 28, 30, and 31.

[0080] The column connection node is a node that connects the node deformation ring to the rest of the column element, including the first column connection node and the second column connection node, such as Figure 4 Midpoint 41, node 42. Beam connection nodes are nodes that connect the node deformation ring to the rest of the beam unit, including the first X-beam connection node and the second X-beam connection node (such as Figure 4 mid-node 43, node 44), and the first Y-beam connection node and the second Y-beam connection node (such as Figure 4 Middle nodes 45, 46).

[0081] Specifically, first, the first Y-beam connection node is connected to the left node of the Y-ring of the Y-coordinate node deformation ring, and the second Y-beam connection node is connected to the right node of the Y-ring of the Y-coordinate node deformation ring, so that the Y-coordinate node deformation ring is connected to the Y-coordinate beam unit; secondly, the first column connection node is connected to the upper node of the X-ring of the X-coordinate node deformation ring and the upper node of the Y-ring of the Y-coordinate node deformation ring, and the second column connection node is connected to the lower node of the X-ring of the X-coordinate node deformation ring and the lower node of the Y-ring of the Y-coordinate node deformation ring, so that the X-coordinate node deformation ring, the Y-coordinate node deformation ring and the rest of the column unit are connected; in this way, the two node deformation rings are connected to the rest of the column unit and the rest of the beam unit. Then, the EqualDOF command in OpenSEES is used to constrain the translational degrees of freedom of the upper nodes of the X-coordinate deformation ring of the node, the upper nodes of the Y-coordinate deformation ring of the node, the lower nodes of the X-coordinate deformation ring of the node, and the lower nodes of the Y-coordinate deformation ring of the node to ensure the integrity of the three-dimensional node after being subjected to force. At this time, the rotational degrees of freedom of the nodes are ignored to avoid the abnormal internal force of the structure caused by the rotation of the node causing the rotation of the connection unit, thereby obtaining a three-dimensional node area model. Among them, the translational degrees of freedom refer to the number of directions in which an object can move freely in space. By restricting the movement of these peripheral nodes in specific directions, the stability and accuracy of the three-dimensional node area model are ensured.

[0082] In a possible implementation, the X-coordinate node deformation ring includes four corner point positions, namely, an X-upper left corner point position, an X-upper right corner point position, an X-lower right corner point position, and an X-lower left corner point position, and each corner point position includes two corner nodes that overlap;

[0083] A node deformation ring is constructed at the intersection of the column element and two beam elements through multiple peripheral nodes and multiple connection elements, including:

[0084] The translational degrees of freedom of the two corner nodes at each corner point are constrained so that the displacements corresponding to the two corner nodes are the same;

[0085] The rotational freedom of the corner nodes at the X-upper left corner point, the X-lower right corner point and the X-lower left corner point is not constrained, so that the X-upper left corner point, the X-lower right corner point and the X-lower left corner point have zero stiffness hinge functions;

[0086] The two corner nodes at the upper right corner position are connected by the X zero-length unit in the connection unit, and the force-deformation constitutive model is applied to the X zero-length unit to make the X upper right corner position have the function of a rotation spring;

[0087] By connecting multiple X-bar units in the unit, each corner node at the X-upper left corner point position, the X-upper right corner point position, the X-lower right corner point position and the X-lower left corner point position is connected with the upper node of the X ring, the right node of the X ring, the lower node of the X ring and the left node of the X ring in turn to obtain the X-coordinate node deformation ring.

[0088] It should be noted that the corner point position refers to the positions of the four vertices of the node deformation ring, and each corner point position is composed of two corner nodes that overlap. Figure 5 In (5-a), in the X-coordinate node deformation ring, the corner point positions include the X-upper left corner point position (including nodes 01 and 02), the X-upper right corner point position (including nodes 04 and 05), the X-lower right corner point position (including nodes 07 and 08), and the X-lower left corner point position (including nodes 10 and 11). Among them, the translational degree of freedom refers to the number of directions in which an object can move freely in space; the rotational degree of freedom refers to the number of directions in which an object can rotate freely in space. By constraining the translational degree of freedom or rotational degree of freedom of the corner node at the corner point position, its movement or rotation in a specific direction can be restricted. A zero stiffness hinge refers to a hinge that allows the node to rotate freely in a specific direction while restricting the translation of the node, and the node has zero stiffness in the direction of rotation. The connection unit includes 1 X zero length unit (such as Figure 5 The unit 109 in (5-a) and 8 X-bar units (such as Figure 5 The zero-length element is a special connection element with zero length, on which a force-deformation constitutive model can be applied to simulate the function of a rotational spring. The bar element is a connection element that connects each peripheral node and is used to simulate the rigidity and deformation behavior of the node deformation ring.

[0089] Specifically, first, the rotational freedom of the corner nodes at the X upper left corner point position, the X lower right corner point position and the X lower left corner point position is not constrained, so that the X upper left corner point position, the X lower right corner point position and the X lower left corner point position have the function of a zero stiffness hinge; secondly, the two corner nodes at the upper right corner point position are connected through the X zero length unit in the connection unit, and the X upper right corner point position is made to have a rotation spring function by applying a force-deformation constitutive model on the X zero length unit; then, each corner node at the X upper left corner point position, the X upper right corner point position, the X lower right corner point position and the X lower left corner point position is connected in sequence with the upper node of the X ring, the right node of the X ring, the lower node of the X ring and the left node of the X ring through multiple X rod units in the connection unit to obtain an X coordinate node deformation ring.

[0090] like Figure 5In (5-a), the X-coordinate node deformation ring consists of 12 peripheral nodes and 9 connection units. There are 2 nodes at each corner point of the X-coordinate node deformation ring. Nodes 01 and 02, nodes 04 and 05, nodes 07 and 08, and nodes 10 and 11 are located at the four corner points of the X-coordinate node deformation ring, respectively. The positions of the two nodes at the same corner point completely overlap. Among them, nodes 01 and 02, nodes 04 and 05, nodes 07 and 08, nodes 10 and 11 are respectively restricted by the EqualDOF command in OpenSEES to make the displacements of the two corresponding nodes completely consistent, ensuring that the three-dimensional node area model maintains its original shape under external loads; the constraints on the rotational freedom of the nodes are slightly different. The rotational freedom of nodes 01 and 02, nodes 07 and 08, nodes 10 and 11 are not constrained, and the two nodes can produce relative rotation, thereby achieving the three corner positions of the X coordinate to the node deformation ring (i.e., the X upper left corner position, the X lower right corner position and the X lower left corner position) Figure 5 The zero stiffness hinge function of (5-a) in the figure is used, while nodes 04 and 05 are connected through unit 109 (zero length unit) in OpenSEES, and a force-deformation constitutive model is assigned to the zero length unit to simulate the shear deformation of the node under the action of external force. Nodes 03, 06, 09, and 12 are connecting nodes for beam units, column units, and X-coordinate node deformation rings. Their main function is to realize the force transmission and displacement coordination between the three-dimensional node area model and the beam-column unit.

[0091] It should be noted that the force-deformation constitutive model is a mathematical model that describes the deformation behavior of a material or structure when subjected to force. In this embodiment, the force-deformation constitutive model of the rotation spring is intended to adopt a trilinear model, such as Figure 6 As shown, Figure 6 is a schematic diagram of a force-deformation constitutive model provided by an embodiment of the present application. Figure 6 In the figure, the V axis represents shear force, the γ axis represents torsion, and V y is the yield shear force, γ y K is the yield torsion e is the elastic shear stiffness, V p is the fully plastic shear force at the node, γ p is the torsion corresponding to the fully plastic shear force at the node, is the ratio of the shear stiffness to the elastic shear stiffness after the node enters the fully plastic stage. y , yield torsion γ y , elastic shear stiffness K e 、Node full plastic shear force V p and its corresponding torsion γ p Determined according to the following formula:

[0092] ;

[0093] ;

[0094] ;

[0095] ;

[0096] ;

[0097] wherein, A eff is the effective shear area, d c is the height of the column element, t p is the web thickness of the column element, G is the shear modulus of the column element material, b c is the width of the flange of the column element, t cf is the thickness of the flange of the column element, d b is the height of the beam element, is the yield stress, is the shear stiffness after the node yields.

[0098] In a possible implementation, the Y - coordinate - direction node deformation ring includes four corner - point positions: the Y upper - left - corner point position, the Y upper - right - corner point position, the Y lower - right - corner point position, and the Y lower - left - corner point position. Each corner - point position includes two coincident corner nodes;

[0099] Construct a node deformation ring at the intersection of the column element and two beam elements through multiple peripheral nodes and multiple connecting elements, including:

[0100] Constrain the translational degrees of freedom of the two corner nodes at each corner - point position so that the displacements corresponding to the two corner nodes are the same;

[0101] Do not constrain the rotational degrees of freedom of the corner nodes at the Y upper - left - corner point position, the Y lower - right - corner point position, and the Y lower - left - corner point position so that the Y upper - left - corner point position, the Y lower - right - corner point position, and the Y lower - left - corner point position have the function of a zero - stiffness hinge;

[0102] Connect the two corner nodes at the upper - right - corner point position through the Y zero - length element in the connecting element, and apply a force - deformation constitutive model on the Y zero - length element so that the Y upper - right - corner point position has the function of a rotational spring;

[0103] Connect each corner node at the Y upper - left - corner point position, the Y upper - right - corner point position, the Y lower - right - corner point position, and the Y lower - left - corner point position to the Y - ring upper node, the Y - ring right node, the Y - ring lower node, and the Y - ring left node in sequence through multiple Y rod elements in the connecting element to obtain the Y - coordinate - direction node deformation ring.

[0104] It should be noted that, as in Figure 5 in (5-b), in the node deformation ring in the Y coordinate direction, the corner positions include the Y upper left corner position (including node 21 and node 22), the Y upper right corner position (including node 24 and node 25), the Y lower right corner position (including node 27 and node 28), and the Y lower left corner position (including node 30 and node 31). Among them, the translational degree of freedom refers to the number of directions in which an object can freely move in space; the rotational degree of freedom refers to the number of directions in which an object can freely rotate in space. By constraining the translational degree of freedom or rotational degree of freedom of the corner nodes at the corner positions, its movement or rotation in a specific direction can be restricted. A zero stiffness hinge means allowing the node to freely move in a specific direction, that is, having zero stiffness. The connecting unit includes 1 Y zero-length unit (such as Figure 5 unit 209 in (5-b)) and 8 Y rod units (such as Figure 5 units 201 to 208 in (5-b)). Among them, the zero-length unit is a special connecting unit with a length of zero, and a force-deformation constitutive model can be applied to it to simulate the function of a rotational spring. The rod unit is a connecting unit that connects each peripheral node and is used to simulate the rigidity and deformation behavior of the node deformation ring.

[0105] Specifically, first, the rotational degrees of freedom of the corner nodes at the Y upper left corner position, the Y lower right corner position, and the Y lower left corner position are not constrained, so that the Y upper left corner position, the Y lower right corner position, and the Y lower left corner position have the function of a zero stiffness hinge; second, through the Y zero-length unit in the connecting unit, the two corner nodes at the upper right corner position are connected, and by applying a force-deformation constitutive model to the Y zero-length unit, the Y upper right corner position has the function of a rotational spring; then, through multiple Y rod units in the connecting unit, each corner node at the Y upper left corner position, the Y upper right corner position, the Y lower right corner position, and the Y lower left corner position is sequentially connected to the Y ring upper node, the Y ring right node, the Y ring lower node, and the Y ring left node to obtain the node deformation ring in the X coordinate direction.

[0106] As in Figure 5In (5-b), the Y-coordinate node deformation ring is composed of 12 peripheral nodes and 9 connection units. There are 2 nodes at each corner point of the Y-coordinate node deformation ring. Nodes 21 and 22, nodes 24 and 25, nodes 27 and 28, and nodes 30 and 31 are located at the four corner points of the Y-coordinate node deformation ring, respectively. The positions of the two nodes at the same corner point completely overlap. Among them, nodes 21 and 22, nodes 24 and 25, nodes 27 and 28, nodes 30 and 31 are respectively restricted by the EqualDOF command in OpenSEES to make the displacements of the two corresponding nodes completely consistent, ensuring that the three-dimensional node area model maintains its original shape under external loads; the constraints on the rotational freedom of the nodes are slightly different. The rotational freedom of nodes 21 and 22, nodes 27 and 28, nodes 30 and 31 are not constrained, and the two nodes can produce relative rotation, thereby achieving the three corner positions of the Y coordinate to the node deformation ring (i.e., the Y upper left corner position, the Y lower right corner position and the Y lower left corner position) Figure 5 The zero stiffness hinge function of (5-b) in the figure is used, while nodes 24 and 25 are connected through unit 209 (zero length unit) in OpenSEES, and a force-deformation constitutive model is assigned to the zero length unit to simulate the shear deformation of the node under the action of external force. Nodes 23, 26, 29, and 32 are connecting nodes for beam units, column units, and Y-coordinate node deformation rings. Their main function is to realize the force transmission and displacement coordination between the three-dimensional node area model and the beam-column unit.

[0107] In a possible implementation, after constructing a node deformation ring on each beam element through a plurality of peripheral nodes and a plurality of connection units, the method for constructing the steel structure finite element model further includes:

[0108] The length values ​​of the plurality of connection units are determined according to the height value of the beam unit, the length value of the column unit in the X direction, and the length value of the column unit in the Y direction.

[0109] It should be noted that the height value of the beam unit refers to the size of the beam unit perpendicular to the XY plane in the finite element model, that is, the size of the beam unit in the thickness direction; the length value of the column unit in the X direction refers to the size of the column unit in the X direction in the finite element model; the length value of the column unit in the Y direction refers to the size of the column unit in the Y direction in the finite element model. These dimensions are related to the geometric shape and size of the beam unit and the column unit. The length value of the connection unit is determined according to the size of the beam unit and the column unit to ensure the accuracy and reliability of the finite element model.

[0110] Specifically, in this embodiment, Figure 5As shown in the figure, in the node deformation ring in the X-coordinate direction and the node deformation ring in the Y-coordinate direction, the length of each connection unit is related to the size of the three-dimensional node area. Assume that the size of the three-dimensional node area is d cx ×d cy ×d b , where d b is the height of the beam element, d cx is the length of the column element in the X direction, d cy is the length of the column unit in the Y direction, then the length of unit 101, unit 102, unit 105, and unit 106 can be 0.5d cx , the length of unit 103, unit 104, unit 107, and unit 108 can be 0.5d b ; The length of unit 201, unit 202, unit 205, and unit 206 can be 0.5d cy , the length of unit 203, unit 204, unit 207, and unit 208 can be 0.5d b .

[0111] In a possible implementation, the method for constructing the steel structure finite element model further includes:

[0112] The elastic unit model is used to simulate the linear elastic section in the middle of the beam unit and the column unit;

[0113] The inelastic unit model is used to simulate the plastic hinge sections at both ends of the beam unit and column unit.

[0114] It should be noted that, based on the mechanical response mechanism of beam-column components in actual working conditions, in the simulation process of the finite element model of the steel structure, all beam units and column units are divided into two sections, namely the linear elastic section located in the middle of the beam unit and the column unit, and the plastic hinge section located at the two ends of the beam unit or the column unit. This is because in actual working conditions, the two ends of the beam-column components are often the sections with the largest stress concentration and deformation. Therefore, it is necessary to set plastic hinges at the two ends of the beam-column components in order to accurately capture the deformation when entering the nonlinear stage. Figure 7 As shown, Figure 7 This is a schematic diagram of modeling a beam-column component provided by an embodiment of the present application. Figure 7 In the structure, the middle part of the beam-column component is an elastic rod, i.e., a linear elastic section; the two ends of the beam-column component are finite length plastic hinges, i.e., plastic hinge sections. The linear elastic section refers to the area where the structure is always in the elastic deformation stage when subjected to force. The plastic hinge section refers to the area where the structure enters the plastic deformation stage after the force reaches a certain degree of yield.

[0115] The elastic unit model is used in finite element analysis to simulate the elastic deformation of a structure when it is subjected to stress. In this embodiment, the linear elastic section in the middle of the beam unit and the column unit is simulated by the elastic unit model. The inelastic unit model is used in finite element analysis to simulate the plastic deformation, yielding or other inelastic behavior of a structure when it is subjected to stress. The restoring force model is used to simulate the relationship between the bending moment and the rotation angle of a structural member after being subjected to stress. In this embodiment, the plastic hinge sections at the two end parts of the beam unit and the column unit are simulated by the modified IMK restoring force model (modified Ibaara-Medina-Krawinkler, mIMK model) to give it nonlinear deformation capability.

[0116] Specifically, for the concentrated plastic model, the linear elastic section between the beam element and the column element is modeled using the elastic element model. bc Modeling the linear elastic section will reduce the stiffness of the beam and column elements, resulting in changes in the dynamic response characteristics of the structure. It is necessary to use the stiffness correction method to correct it. The moment of inertia of the beam-column structure after correction is I mod = I bc ×(n+1) / n, where n can be 10. The plastic hinge sections at both ends of the beam element and column element can be modeled using zero-length elements, using the modified IMK restoring force model (modified Ibaara-Medina-Krawinkler, mIMK model). Figure 8 As shown, Figure 8 It is a schematic diagram of a mIMK model provided in one embodiment of the present application. The mIMK model can effectively simulate the stiffness and strength degradation of components under repeated loading. It has three different loading degradation rules: bilinear, peak-pointing, and pinching. Among them, the bilinear mIMK model is widely used to simulate the hysteresis performance of steel structure components. The bilinear mIMK model can be used to simulate the plastic hinge of steel components. The parameters of the mIMK model can be determined according to the following formula, namely:

[0117] ;

[0118] ;

[0119] in: is the modified plastic hinge stiffness; is the modified plastic hinge strength; K 0 is the stiffness of the beam-column element; is the plastic development coefficient. For the strong axis and weak axis of the H-section, The values ​​of are 1.05 and 1.2 respectively. For hollow sections, The value is 1.05; W is the elastic section modulus; F y is the yield stress; P g is the axial load caused by gravity; P ye is the axial yield load of the column.

[0120] For the mIMK model of H-section members, the plastic rotation angle θ p 、The post-peak rotation angle θ pc And the cumulative plastic deformation capacity Λ can be determined according to the following formula:

[0121] ;

[0122] ;

[0123] ;

[0124] Where h is the web height of the H-section member, t w is the web thickness of the H-section member, b f is the flange width of the H-section member, t f is the flange thickness of the H-section member, L is the length of the H-section member, d is the height of the H-section member, F y is the yield stress.

[0125] For the mIMK model of box-section members, the plastic rotation angle θ p 、The post-peak rotation angle θ pc And the cumulative plastic deformation capacity Λ can be determined according to the following formula:

[0126] ;

[0127] ;

[0128] ;

[0129] Where D is the height of the box section member, t is the thickness of the box section member, P g is the axial load caused by gravity, P ye is the axial yield load of the column, F y is the yield stress.

[0130] Specifically, the peak moment to yield moment ratio of the mIMK model is c / M y It is determined according to the type of beam unit and column unit. In this embodiment, after multiple tests, it is found that the M of the beam unit and the column unit is c / M yThe values ​​can be 1.17 and 1.10 respectively; the limiting rotation angle is 0.055rad; and the residual strength ratio is 0.4.

[0131] In a possible implementation, the stiffness of the connection unit is a first preset multiple of the stiffness of the beam unit and the stiffness of the column unit; the strength of the connection unit is a second preset multiple of the strength of the beam unit and the strength of the column unit.

[0132] It should be noted that in the finite element model, the connection unit is an element used to connect different structural units (such as beam units, column units, etc.). The connection unit simulates the mechanical behavior of the connection part in the actual steel structure, mainly including strength, stiffness, etc. The first preset multiple is a preset value for calculating the stiffness of the connection unit; the second preset multiple is a preset value for calculating the strength of the connection unit. The first preset multiple and the second preset multiple are determined after multiple tests. Since the deformation of the three-dimensional node area at the intersection of the beam and column of the steel structure is mainly shear deformation, and its bending deformation and axial deformation are very small, the strength and stiffness of all connection units can be taken as larger values ​​to avoid axial deformation and bending deformation in the three-dimensional node area, but the stiffness of the connection unit will affect the dynamic characteristics of the steel structure. If the stiffness of the connection unit in the node is not appropriate, it may cause the dynamic characteristics (modal, period, etc.) of the finite element model of the rigid structure to deviate from the dynamic characteristics of the real steel structure, resulting in inaccurate calculation results of the finite element model. After many tests, the strength and stiffness of all connection units in the three-dimensional node area model can be set to 20 times the strength and stiffness of the beam unit and the column unit, that is, in this embodiment, the first preset multiple and the second preset multiple can be 20. At this time, the strength and stiffness of the three-dimensional node area are large enough to ensure that the axial deformation and bending deformation of the three-dimensional node area are small and will not affect the nonlinear deformation of the structure, and can also ensure that the dynamic characteristics of the finite element model of the steel structure are basically consistent with the dynamic characteristics of the actual steel structure.

[0133] It should be understood that the size of the serial numbers of the steps in the above embodiments does not mean the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of the present application.

[0134] Corresponding to the method for constructing a finite element model of a steel structure in the above embodiment, Fig. 9 A schematic structural diagram of a device for constructing a finite element model of a steel structure provided in one embodiment of the present application is shown. For ease of explanation, only the parts related to the embodiment of the present application are shown.

[0135] Reference Fig. 9 , the construction device 3 of the steel structure finite element model of this embodiment includes:

[0136] A node deformation ring construction module 31 is used to construct a node deformation ring at the intersection of a column unit and a beam unit through a plurality of peripheral nodes and a plurality of connection units during the finite element modeling process for a steel structure high-rise building, wherein the node deformation ring is in the shape of a parallelogram, the connection unit is used to represent the intersection of the column unit and the beam unit, and the peripheral node is used to represent the vertex or midpoint of the intersection of the column unit and the beam unit;

[0137] The beam-column node area construction module 32 is used to connect the two node deformation rings with the rest of the column unit and the rest of the beam unit through multiple column connection nodes and multiple beam connection nodes to obtain a three-dimensional node area model at the intersection of the beam unit and the column unit; wherein the rest is the part of the column unit or the beam unit excluding the intersection of the column unit and the beam unit.

[0138] It can be understood that in the construction device 3 of a steel structure finite element model provided in the embodiment of the present application, in the process of finite element modeling for a steel structure high-rise building, the node deformation ring construction module 31 constructs a node deformation ring at the intersection of the column unit and the beam unit through multiple peripheral nodes and multiple connection units, wherein the shape of the node deformation ring is a parallelogram, the connection unit is used to characterize the edge of the intersection of the column unit and the beam unit, and the peripheral node is used to characterize the vertex or midpoint of the intersection of the column unit and the beam unit; the beam-column node area construction module 32 connects the two node deformation rings with the rest of the column unit and the rest of the beam unit through multiple column connection nodes and multiple beam connection nodes to obtain a three-dimensional node area model at the intersection of the beam unit and the column unit; wherein the rest is the part of the column unit or the beam unit except the intersection of the column unit and the beam unit. By constructing a three-dimensional node area model that can analyze the stress deformation of the beam-column node area in the process of finite element modeling for a steel structure high-rise building, the nonlinear response results of the steel structure can be accurately calculated, and at the same time, the problem of low calculation efficiency can be avoided.

[0139] It should be noted that the information interaction, execution process, etc. between the modules in the above-mentioned steel structure finite element model construction device 3 are based on the same concept as the method embodiment of the present application. Their specific functions and technical effects can be found in the method embodiment part and will not be repeated here.

[0140] Specifically, in order to verify the feasibility of the three-dimensional node area model in this application, the model and test results published on the website are selected for comparison, which is called the comparison model. The comparison model is a series of tests to evaluate the seismic performance of steel structure high-rise buildings under long-period ground motion and near-fault earthquakes, including large-scale shaking table tests of the comparison model. The plane layout, overall structure and test vibration direction of the comparison model are as follows Fig.10 As shown, Fig.10It is a structural schematic diagram of a comparison model provided in one embodiment of the present application. Fig.10 (10-a) in the figure is the plane layout of the comparison model. Fig.10 (10-b) is the overall three-dimensional schematic diagram of the comparison model. Fig.10 (10-b) is a schematic diagram of the test of the comparison model. Fig.10 In (10-a), C1 represents a column, B1 represents a longitudinal beam, and G1 represents a transverse beam. The specific parameters of the structural components of the comparison model can be seen in Table 1. The mass of each layer is 219 kN. Table 1 is a table of the size parameters of the structural components of the comparison model.

[0141] Table 1 Structural component size parameters of the comparison model (unit: mm)

[0142]

[0143] In order to comprehensively compare the influence of the three-dimensional node area model of the beam-column unit on the structural response simulation results under different structural modeling schemes, the beam-column modeling scheme of the finite length plastic hinge model is adopted. There are three node area modeling schemes for the beam-column unit, namely, the node area-free scheme, the scissor-shaped node area scheme, and the parallelogram node area scheme. The above beam-column modeling scheme and the node area modeling scheme of the beam-column unit are combined, and there are three modeling schemes in total. The corresponding numbers of the modeling schemes are shown in Table 2. Table 2 shows the corresponding numbers of the three model schemes.

[0144] Table 2 Corresponding numbers of the three model schemes

[0145]

[0146] The mechanical parameters (such as yield strength) of the beam-column unit need to be determined according to the specific steel type used in the component. Material property tests were conducted on steels of different sizes and grades, and the relevant parameters in the finite element model were taken as the average value of the material property tests. The C1 columns, all B1 beams, and G1 beams of the 1-6F floors are made of SM490A steel, and the C1 columns of the 7-18F floors are made of BCR295 steel. The elastic modulus of all steels is E = 210 GPa, and the Poisson's ratio is ν = 0.3, the shear modulus is G = 79.3 GPa, and the yield strength of each component is averaged according to the material property test results as shown in Table 3. Table 3 shows the yield strength of steel of different materials and thicknesses.

[0147] Table 3 Yield strength of steel of different materials and thicknesses

[0148]

[0149] As shown in Table 4, Table 4 compares the period of the first two-order modes of the structure under different model schemes with the period of the test model. Fig.11 As shown, Fig.11 It is a schematic diagram of the first two vibration modes of a structure under different model schemes provided in an embodiment of the present application. Fig.11 In the figure, C-0 represents a finite-length plastic hinge model without nodes, CS represents a finite-length plastic hinge model with three-dimensional scissor-type nodes, and CP represents a finite-length plastic hinge model with three-dimensional parallelogram nodes.

[0150] Table 4 Comparison of the period of the first two modes of the structure under different model schemes and the period of the test model

[0151]

[0152] The finite element model was established by applying four synthetic seismic waves respectively, and the difference between the inter-layer displacement angle of the structure under different schemes and the test results was compared. The seismic intensity was selected from the pseudo velocity spectrum as pSV 110cm / s and pSV 340cm / s. Under these two seismic intensities, the structure corresponds to the initial nonlinear and fully elastic-plastic response conditions respectively.

[0153] (1) Preliminary comparative analysis of nonlinear seismic responses

[0154] like Fig.12 As shown, Fig.12 It is a schematic diagram of the maximum inter-story displacement angle response of different model schemes under the action of preliminary nonlinear seismic loads provided in an embodiment of the present application. Fig.12 (12-a) in the figure represents the maximum inter-story displacement angular response of the C-0 model under the preliminary nonlinear earthquake load; Fig.12 (12-b) in the figure represents the maximum inter-story displacement angular response of the CS model under the preliminary nonlinear earthquake load; Fig.12 (12-c) in the figure represents the maximum inter-story displacement angle response of the CP model under the preliminary nonlinear seismic load. In order to quantitatively compare the calculation accuracy of different model schemes, the average error and maximum inter-story displacement angle error of all floors of the three model schemes under different seismic waves are given, as shown in Tables 5 and 6. Table 5 shows the average error of the average inter-story displacement angle of all floors of different model schemes under different seismic waves under the preliminary nonlinear seismic load, and Table 6 shows the maximum inter-story displacement angle error of different model schemes under different seismic waves under the preliminary nonlinear seismic load.

[0155] Table 5 Average errors of inter-story displacement angles of all floors under different model schemes under different earthquake waves under preliminary nonlinear earthquake loads

[0156]

[0157] Table 6 Maximum inter-story displacement angle errors of different model schemes under different seismic waves under preliminary nonlinear seismic loads

[0158]

[0159] Overall, the calculated structural responses are consistent with the test results, with or without the three-dimensional node zone model, and the maximum inter-story displacement angle error does not exceed 7%.

[0160] In order to compare the computational efficiency of different model schemes, the time consumption of calculating 1000 time steps in three model schemes is given, such as Fig.13 As shown, Fig.13 It is a schematic diagram of the calculation time consumption of different model schemes under the action of preliminary nonlinear seismic loads provided by an embodiment of the present application. Fig.13 In the figure, C-0 represents a finite-length plastic hinge model without nodes, CS represents a finite-length plastic hinge model with three-dimensional scissor-type nodes, and CP represents a finite-length plastic hinge model with three-dimensional parallelogram nodes. Fig.13 In the present study, as the complexity of the modeling scheme increases, the computational time also increases significantly, and the model with three-dimensional node area is more time-consuming than the model without three-dimensional node area.

[0161] (2) Comparative analysis of nonlinear seismic response

[0162] like Fig.14 As shown, Fig.14 It is a schematic diagram of the maximum inter-story displacement angle response of different model schemes under the action of completely elastic-plastic seismic loads provided in one embodiment of the present application. Fig.14 (14-a) in the figure represents the maximum inter-story displacement angle response of the C-0 model under the fully elastic-plastic earthquake load; Fig.14 (14-b) in the figure represents the maximum inter-story displacement angular response of the CS model under fully elastic-plastic earthquake loads; Fig.14 (14-c) in the figure represents the maximum inter-story displacement angle response of the CP model under the fully elastic-plastic earthquake load. In order to quantitatively compare the calculation accuracy of different modeling schemes, the average error and maximum inter-story displacement angle error of all floors of the three model schemes under the fully elastic-plastic earthquake load are given, as shown in Tables 7 and 8. Table 7 shows the average error of the inter-story displacement angle of all floors of different model schemes under different earthquake waves under the fully elastic-plastic earthquake load, and Table 8 shows the maximum inter-story displacement angle error of different model schemes under different earthquake waves under the fully elastic-plastic earthquake load.

[0163] Table 7 Average error of inter-story displacement angles of all floors under different model schemes under completely elastic-plastic earthquake loads and different earthquake waves

[0164]

[0165] Table 8 Maximum inter-story displacement angle errors of different model schemes under completely elastic-plastic earthquake loads under different seismic waves

[0166]

[0167] During the scaled model test, the structure collapsed under an earthquake intensity of pSV 340cm / s; although no collapse failure occurred in the three model schemes in the finite element analysis, the maximum inter-story displacement angle of the structure exceeded 0.04rad, and the structure had entered a critical failure state.

[0168] It can be seen that the simulation accuracy of the finite-length plastic hinge model is related to the modeling method of the three-dimensional node area. The average displacement angle of all floors and the maximum inter-story displacement angle errors of the model without three-dimensional node area are the largest. Compared with the scissor-shaped node and the parallelogram node, although the maximum displacement angle error of the parallelogram node is greater than that of the scissor-shaped node, the average displacement angle error of all floors of the parallelogram node model is less than 5%, while that of the scissor-shaped node model reaches 32%.

[0169] In summary, both the scissor-shaped node and parallelogram node modeling schemes have higher accuracy than the model without three-dimensional node area, among which the parallelogram node modeling scheme has the highest calculation accuracy.

[0170] The present application also provides a terminal device, such as Fig.15 As shown, Fig.15 This is a schematic diagram of the structure of a terminal device provided in one embodiment of the present application. Fig.15 The terminal device 4 of this embodiment includes: a memory 41, a processor 42, and a computer program stored in the memory 41 and executable on the processor 42. When the processor 42 executes the computer program, the steps in any one of the above-mentioned methods for constructing a finite element model of a steel structure are implemented.

[0171] The embodiment of the present application further provides a computer-readable storage medium, which stores a computer program. When the computer program is executed by a processor, the steps in the above-mentioned method embodiments can be implemented.

[0172] An embodiment of the present application provides a computer program product. When the computer program product runs on a mobile terminal, the mobile terminal can implement the steps in the above-mentioned method embodiments when executing the computer program product.

[0173] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the present application implements all or part of the processes in the above-mentioned embodiment method, which can be completed by instructing the relevant hardware through a computer program. The computer program can be stored in a computer-readable storage medium. When the computer program is executed by the processor, the steps of the above-mentioned method embodiments can be implemented. Among them, the computer program includes computer program code, and the computer program code can be in source code form, object code form, executable file or some intermediate form. The computer-readable medium may at least include: any entity or device that can carry the computer program code to the camera / terminal device, recording medium, computer memory, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), electric carrier signal, telecommunication signal and software distribution medium. For example, a USB flash drive, a mobile hard disk, a magnetic disk or an optical disk. In some jurisdictions, according to legislation and patent practice, computer-readable media cannot be electric carrier signals and telecommunication signals.

[0174] In the above embodiments, the description of each embodiment has its own emphasis. For parts that are not described or recorded in detail in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

[0175] Those of ordinary skill in the art will appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of this application.

[0176] In the embodiments provided in the present application, it should be understood that the disclosed devices / network equipment and methods can be implemented in other ways. For example, the device / network equipment embodiments described above are only schematic. For example, the division of modules or units is only a logical function division. There may be other division methods in actual implementation, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some interfaces, indirect coupling or communication connection of devices or units, which can be electrical, mechanical or other forms.

[0177] The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0178] The above embodiments are only used to illustrate the technical solutions of the present application, rather than to limit them. Although the present application has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. These modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present application, and should all be included in the protection scope of the present application.

Claims

1. A method for constructing a finite element model of a steel structure, characterized in that: include: In the process of modeling a finite element model for a steel structure high-rise building, a node deformation ring is constructed at the intersection of a column unit and a beam unit through a plurality of peripheral nodes and a plurality of connection units, wherein the shape of the node deformation ring is a parallelogram, the connection unit is used to represent the edge of the intersection of the column unit and the beam unit, and the peripheral node is used to represent the vertex or midpoint of the intersection of the column unit and the beam unit; The two node deformation rings are connected with the rest of the column unit and the rest of the beam unit through a plurality of column connection nodes and a plurality of beam connection nodes to obtain a three-dimensional node area model at the intersection of the beam unit and the column unit; wherein the rest is the portion of the column unit or the beam unit excluding the intersection of the column unit and the beam unit; wherein, The node deformation ring includes an X-coordinate node deformation ring and a Y-coordinate node deformation ring, the beam unit includes an X-coordinate beam unit and a Y-coordinate beam unit; the peripheral nodes include an upper node, a lower node, a left node, a right node and a plurality of corner nodes; The method connects the two node deformation rings with the rest of the column unit and the rest of the beam unit through a plurality of column connection nodes and a plurality of beam connection nodes to obtain a three-dimensional node area model at the intersection of the beam unit and the column unit, including: Connecting a first X-beam connection node to a left-side node of the X-coordinate node deformation ring, and connecting a second X-beam connection node to a right-side node of the X-coordinate node deformation ring, so that the X-coordinate node deformation ring is connected to the X-coordinate beam unit; Connecting a first Y-beam connection node to a Y-ring left node of the Y-coordinate node deformation ring, and connecting a second Y-beam connection node to a Y-ring right node of the Y-coordinate node deformation ring, so that the Y-coordinate node deformation ring is connected to the Y-coordinate beam unit; Connecting the first column connection node to the upper X-ring node of the X-coordinate node deformation ring and the upper Y-ring node of the Y-coordinate node deformation ring, and connecting the second column connection node to the lower X-ring node of the X-coordinate node deformation ring and the lower Y-ring node of the Y-coordinate node deformation ring, so as to connect the X-coordinate node deformation ring, the Y-coordinate node deformation ring and the rest of the column unit; Constraining the translational degrees of freedom of the upper X-ring node of the X-coordinate node deformation ring, the upper Y-ring node of the Y-coordinate node deformation ring, the lower X-ring node of the X-coordinate node deformation ring, and the lower Y-ring node of the Y-coordinate node deformation ring to obtain the three-dimensional node area model; and, The X-coordinate node deformation ring includes four corner point positions, namely, the X-upper left corner point position, the X-upper right corner point position, the X-lower right corner point position, and the X-lower left corner point position, and each corner point position includes two corner nodes that overlap; The method of constructing a node deformation ring at the intersection of a column unit and two beam units by using a plurality of peripheral nodes and a plurality of connection units comprises: Constraining the translational degrees of freedom of the two corner nodes at each corner point position so that the displacements corresponding to the two corner nodes are the same; The rotational freedom of the corner nodes of the X upper left corner point position, the X lower right corner point position and the X lower left corner point position is not constrained, so that the X upper left corner point position, the X lower right corner point position and the X lower left corner point position have a zero stiffness hinge function; The two corner nodes at the upper right corner point are connected by an X zero-length unit in the connection unit, and a force-deformation constitutive model is applied to the X zero-length unit so that the X upper right corner point has a rotation spring function; Through the multiple X-bar units in the connection unit, each corner node at the X-upper left corner point position, the X-upper right corner point position, the X-lower right corner point position and the X-lower left corner point position is connected in sequence with the X-ring upper node, the X-ring right node, the X-ring lower node and the X-ring left node to obtain the X-coordinate node deformation ring.

2. The method for constructing a finite element model of a steel structure according to claim 1, characterized in that: The Y coordinate node deformation ring includes four corner point positions, namely, a Y upper left corner point position, a Y upper right corner point position, a Y lower right corner point position, and a Y lower left corner point position, and each of the corner point positions includes two corner nodes that overlap; The method of constructing a node deformation ring at the intersection of a column unit and two beam units by using a plurality of peripheral nodes and a plurality of connection units comprises: Constraining the translational degrees of freedom of the two corner nodes at each corner point position so that the displacements corresponding to the two corner nodes are the same; The rotational freedom of the corner nodes of the Y upper left corner point position, the Y lower right corner point position and the Y lower left corner point position is not constrained, so that the Y upper left corner point position, the Y lower right corner point position and the Y lower left corner point position have a zero stiffness hinge function; The two corner nodes at the upper right corner point are connected by a Y zero-length unit in the connection unit, and a force-deformation constitutive model is applied to the Y zero-length unit so that the Y upper right corner point has a rotation spring function; Through the multiple Y rod units in the connecting unit, each corner node at the Y upper left corner point position, the Y upper right corner point position, the Y lower right corner point position and the Y lower left corner point position is connected with the upper node of the Y ring, the right node of the Y ring, the lower node of the Y ring and the left node of the Y ring in sequence to obtain the Y coordinate node deformation ring.

3. The method for constructing a finite element model of a steel structure according to claim 1, characterized in that: After constructing a node deformation ring on each beam element through a plurality of peripheral nodes and a plurality of connection elements, the method further comprises: The length values ​​of the plurality of connection units are determined according to the height value of the beam unit, the length value of the column unit in the X direction, and the length value of the column unit in the Y direction.

4. The method for constructing a finite element model of a steel structure according to claim 1, wherein: The method further comprises: Using an elastic unit model, simulating a linear elastic section in the middle of the beam unit and the column unit; The plastic hinge sections at both ends of the beam element and the column element are simulated by using a non-elastic unit model.

5. The method for constructing a finite element model of a steel structure according to any one of claims 1 to 4, characterized in that: The stiffness of the connection unit is a first preset multiple of the stiffness of the beam unit and the stiffness of the column unit; the strength of the connection unit is a second preset multiple of the strength of the beam unit and the strength of the column unit.

6. A device for constructing a finite element model of a steel structure, characterized in that: include: A node deformation ring construction module is used to construct a node deformation ring at the intersection of a column unit and a beam unit through a plurality of peripheral nodes and a plurality of connection units during finite element modeling of a steel structure high-rise building, wherein the shape of the node deformation ring is a parallelogram, the connection unit is used to represent the intersection of the column unit and the beam unit, and the peripheral node is used to represent the vertex or midpoint of the intersection of the column unit and the beam unit; The beam-column node area construction module is used to connect the two node deformation rings with the rest of the column unit and the rest of the beam unit through multiple column connection nodes and multiple beam connection nodes to obtain a three-dimensional node area model at the intersection of the beam unit and the column unit; wherein the rest is the part of the column unit or the beam unit excluding the intersection of the column unit and the beam unit; wherein, The node deformation ring includes an X-coordinate node deformation ring and a Y-coordinate node deformation ring, the beam unit includes an X-coordinate beam unit and a Y-coordinate beam unit; the peripheral nodes include an upper node, a lower node, a left node, a right node and a plurality of corner nodes; The beam-column node area construction module is specifically used for: Connecting a first X-beam connection node to a left-side node of the X-coordinate node deformation ring, and connecting a second X-beam connection node to a right-side node of the X-coordinate node deformation ring, so that the X-coordinate node deformation ring is connected to the X-coordinate beam unit; Connecting a first Y-beam connection node to a Y-ring left node of the Y-coordinate node deformation ring, and connecting a second Y-beam connection node to a Y-ring right node of the Y-coordinate node deformation ring, so that the Y-coordinate node deformation ring is connected to the Y-coordinate beam unit; Connecting the first column connection node to the upper X-ring node of the X-coordinate node deformation ring and the upper Y-ring node of the Y-coordinate node deformation ring, and connecting the second column connection node to the lower X-ring node of the X-coordinate node deformation ring and the lower Y-ring node of the Y-coordinate node deformation ring, so as to connect the X-coordinate node deformation ring, the Y-coordinate node deformation ring and the rest of the column unit; Constraining the translational degrees of freedom of the upper X-ring node of the X-coordinate node deformation ring, the upper Y-ring node of the Y-coordinate node deformation ring, the lower X-ring node of the X-coordinate node deformation ring, and the lower Y-ring node of the Y-coordinate node deformation ring to obtain the three-dimensional node area model; and, The X-coordinate node deformation ring includes four corner point positions, namely, the X-upper left corner point position, the X-upper right corner point position, the X-lower right corner point position, and the X-lower left corner point position, and each corner point position includes two corner nodes that overlap; The beam-column node area construction module is also specifically used for: Constraining the translational degrees of freedom of the two corner nodes at each corner point position so that the displacements corresponding to the two corner nodes are the same; The rotational freedom of the corner nodes of the X upper left corner point position, the X lower right corner point position and the X lower left corner point position is not constrained, so that the X upper left corner point position, the X lower right corner point position and the X lower left corner point position have a zero stiffness hinge function; The two corner nodes at the upper right corner point are connected by an X zero-length unit in the connection unit, and a force-deformation constitutive model is applied to the X zero-length unit so that the X upper right corner point has a rotation spring function; Through the multiple X-bar units in the connection unit, each corner node at the X-upper left corner point position, the X-upper right corner point position, the X-lower right corner point position and the X-lower left corner point position is connected in sequence with the X-ring upper node, the X-ring right node, the X-ring lower node and the X-ring left node to obtain the X-coordinate node deformation ring.

7. A terminal device, characterized in that: The method comprises a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the method according to any one of claims 1 to 5 when executing the computer program.

8. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the method according to any one of claims 1 to 5 is implemented.

Images