Structural response determination method and device based on nonlinear self-finding secant iteration

Abstract

The present application relates to the technical field of nonlinear solution of structural response, and specifically to a method and device for determining structural response based on nonlinear self-finding secant iteration. The method comprises the following steps: establishing a structural nonlinear solution model based on the structural dynamic equilibrium state and the structural response relationship; obtaining two corresponding preset structural response value errors based on the structural nonlinear solution model according to two preset structural response values; obtaining an approximate structural response value according to the preset structural response value and the corresponding preset structural response value error; if the approximate structural response error value does not meet the error tolerance requirement, determining the next preset structural response value according to the preset structural response value error and the approximate structural response error value, until the approximate structural response error value meets the error tolerance requirement. This method can solve the technical problem that the initial iterative value deviates far from the true solution, causing the iterative process to diverge or converge to a non-true solution, thereby affecting the reliability of the structural response result.

Description

Technical Field

[0001] The present application relates to the technical field of nonlinear solution of structural response, and in particular to a method and device for determining structural response based on nonlinear self-finding secant iteration. Background Art

[0002] The response of a structure to external forces, such as the structural response of a seismic isolation device, requires a nonlinear dynamic analysis of the structure.

[0003] In existing techniques, nonlinear dynamic analysis of structures often relies on iterations of trial values, as the mechanical equilibrium equations are nonlinear and difficult to solve directly. However, if the initial iteration value deviates significantly from the true solution, or if the function exhibits complex behavior such as multiple extreme points or non-monotonicity within the initial value neighborhood, the iteration process may diverge or converge to an undesirable local solution, thus affecting the reliability of the structural response results. Summary of the Invention

[0004] The present application provides a method and device for determining structural response based on nonlinear self-finding secant iteration, which can solve the technical problem in the prior art that the initial iteration value deviates far from the true solution, which may cause the iterative process to diverge or converge to an undesired local solution, thereby affecting the reliability of the structural response results.

[0005] In order to achieve the above purpose, the technical solution adopted by the present invention is:

[0006] In one aspect, the present invention provides a method for determining structural response based on nonlinear self-finding secant iteration, comprising the following steps:

[0007] Based on the structural dynamic equilibrium state and structural response relationship, a structural nonlinear solution model is established;

[0008] Substituting the structural design parameters and the equivalent external load into the structural nonlinear solution model, and obtaining two corresponding preset structural response value errors according to the two preset structural response values;

[0009] Based on the secant method, an approximate structural response value is obtained according to two preset structural response values ​​and their corresponding preset structural response value errors;

[0010] If the approximate structural response error value corresponding to the approximate structural response value meets the error tolerance requirement, the approximate structural response value is used as the output;

[0011] If the approximate structural response error value corresponding to the approximate structural response value does not meet the error tolerance requirement, the next two preset structural response values ​​are determined according to the two preset structural response value errors and the approximate structural response error value until the obtained approximate structural response error value meets the error tolerance requirement, and the approximate structural response value is used as the output.

[0012] In some optional solutions, establishing a structural nonlinear solution model based on the structural dynamic equilibrium state and the structural response relationship includes:

[0013] The structural dynamic equilibrium equation is established based on the structural dynamic equilibrium state, and the structural response equation is established based on the structural response relationship;

[0014] According to the structural dynamic equilibrium equation and structural response equation, a structural dynamic equilibrium model is established;

[0015] Based on the structural dynamics equilibrium model, a structural nonlinear solution model is established.

[0016] In some optional solutions, the structural dynamics equilibrium equation is ;

[0017] The structural response equation is ;

[0018] Equivalent structural stiffness ;

[0019] Equivalent external load ;

[0020] The structural dynamics equilibrium model is ;

[0021] The structural nonlinear solution model is: ;

[0022] in, u is the displacement of the structure, v is the displacement velocity of the structure, Δt is the time step, Q is the equivalent external load of the structure, F d is the internal force of the seismic isolation structure, K is the equivalent stiffness of the structure excluding the isolation bearing components, is the structural stiffness matrix, is the structural mass matrix, 、 and are the displacement, velocity and acceleration of the structure at the previous moment, is the external load on the structure at the current moment, is the structural response error.

[0023] In some optional solutions, the two corresponding preset structure response value errors are respectively:

[0024] as well as ;

[0025] in, v0 is the first preset structure response value, v 1 is the second preset structure response value, is the preset structural response value error corresponding to the first preset structural response value, is the preset structural response value error corresponding to the second preset structural response value.

[0026] In some alternative solutions, according to the formula , to obtain the approximate structural response value v 2.

[0027] In some optional solutions, the error tolerance requirement is ,in, is the preset structural response value error corresponding to the approximate structural response value, ε is the error tolerance.

[0028] In some optional solutions, the error tolerance ε =10 -6 .

[0029] In some optional solutions, two preset structural response values ​​for the next step are determined based on the same-sign or different-sign relationship between the two preset structural response value errors and the approximate structural response error values.

[0030] In some optional solutions, determining the next two preset structural response values ​​based on the same-sign or different-sign relationship between the two preset structural response value errors and the approximate structural response error values ​​includes:

[0031] like f ( v 1) and f ( v 2) Different signs, the next step is to preset the response values ​​of the two structures v 01 = v 1, v 11 = v 2;

[0032] like f ( v 1) and f ( v 2) Same sign, and f ( v 2) with f ( v 0) different signs, the next step is to preset the response values ​​of the two structures v 01 = v 0, v 11 = v 2;

[0033] like f ( v 0), f ( v 1) and f ( v 2) With the same sign, the next step is to set the response value of the two preset structures to v 01 = v 1, v 11 = v 1+0.5( v 2- v 1);

[0034] in, v 0 is the first preset structure response value, v 1 is the second preset structure response value, v 2 is the approximate structural response value, v 01 For the next step, the first preset structure response value is v 11 Next, the second preset structure response value, is the preset structural response value error corresponding to the first preset structural response value, is the preset structural response value error corresponding to the second preset structural response value, is the preset structural response value error corresponding to the approximate structural response value.

[0035] On the other hand, the present invention further provides a structural response determination device based on nonlinear self-finding secant iteration, comprising:

[0036] A model building module is used to build a structural nonlinear solution model based on the structural dynamic equilibrium state and structural response relationship;

[0037] An error acquisition module is used to substitute the structural design parameters and the equivalent external load into the structural nonlinear solution model, and obtain two corresponding preset structural response value errors according to two preset structural response values;

[0038] An approximate value acquisition module, which is used to obtain an approximate structural response value based on two preset structural response values ​​and their corresponding preset structural response value errors based on a secant method;

[0039] The judgment iteration module is used to output the approximate structural response value when the approximate structural response error value corresponding to the approximate structural response value meets the error tolerance requirement; when the approximate structural response error value corresponding to the approximate structural response value does not meet the error tolerance requirement, determine the next two preset structural response values ​​according to the two preset structural response value errors and the approximate structural response error value, until the obtained approximate structural response error value meets the error tolerance requirement, and use the approximate structural response value as the output.

[0040] Compared with the prior art, the advantages of the present invention are: the present invention obtains an approximate structural response value based on two preset structural response values ​​and their corresponding preset structural response value errors; determines whether to use the approximate structural response value as output by judging whether the approximate structural response error value corresponding to the approximate structural response value meets the error tolerance requirement; when the approximate structural response error value corresponding to the approximate structural response value does not meet the error tolerance requirement, determines the next two preset structural response values ​​based on the two preset structural response value errors and the approximate structural response error value, until the obtained approximate structural response error value meets the error tolerance requirement, and uses the approximate structural response value as output. This can solve the problem in the prior art that when the trial value iteration method is used for solving, if the initial iteration value deviates far from the true solution, or the function presents complex states such as multiple extreme points and non-monotonicity in the neighborhood of the initial value, it may cause the iterative process to diverge, thereby affecting the reliability of the structural response result. BRIEF DESCRIPTION OF THE DRAWINGS

[0041] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the following briefly introduces the drawings required for use in the description of the embodiments. 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 any creative work.

[0042] Figure 1 Flowchart of a method for determining structural response based on nonlinear self-finding secant iteration in an embodiment of the present invention;

[0043] Figure 2 Schematic diagram of the hardware structure of a device for determining a structural response based on nonlinear self-finding secant iteration in an embodiment of the present invention;

[0044] Figure 3 4 is a comparison diagram of the convergence curves of the embodiment of the present invention and the traditional secant method. DETAILED DESCRIPTION

[0045] In order to enable those skilled in the art to better understand the present invention, the following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0046] The present invention provides a method for determining structural response based on nonlinear self-finding secant iteration, comprising the following steps:

[0047] S1: Based on the structural dynamic equilibrium state and structural response relationship, a structural nonlinear solution model is established.

[0048] In some optional embodiments, step S1 includes the following steps:

[0049] S11: The structural dynamic equilibrium equation is established based on the structural dynamic equilibrium state, and the structural response equation is established based on the structural response relationship.

[0050] The structural dynamic equilibrium equation is .

[0051] The structural response equation is .

[0052] Equivalent structural stiffness

[0053] The equivalent external load

[0054] In this example, a single-degree-of-freedom viscous damper is used as an example. The internal force of the viscous damper is ,

[0055] in, u is the displacement of the structure, v is the displacement velocity of the structure, Δt is the time step, Q is the equivalent external load of the structure, F d is the internal force of the seismic isolation structure, K is the equivalent stiffness of the structure excluding the isolation bearing components, is the structural stiffness matrix, is the structural mass matrix, 、 and is the displacement, velocity and acceleration of the structure at the previous moment, is the external load on the structure at the current moment, is the structural response error. C is the damping coefficient of the viscous damper, α is the viscous damper velocity index.

[0056] S12: Establish a structural dynamic equilibrium model based on the structural dynamic equilibrium equation and the structural response equation.

[0057] In this example, the structural dynamic equilibrium equation and the structural response equation are combined to establish a structural dynamic equilibrium model.

[0058] The structural dynamic equilibrium model is , taking the single degree of freedom viscous damper as an example, the structural dynamic equilibrium model is .

[0059] S13: Based on the structural dynamic equilibrium model, a structural nonlinear solution model is established.

[0060] In this example, the structural dynamic equilibrium model is deformed to obtain a structural nonlinear solution model.

[0061] The structural nonlinear solution model is: Taking a single-degree-of-freedom viscous damper as an example, the structural nonlinear solution model is: ;

[0062] S2: Substitute the structural design parameters and equivalent external loads into the structural nonlinear solution model, and obtain two corresponding preset structural response value errors based on the two preset structural response values.

[0063] In this example, two initial preset structural response values ​​are selected, namely the first preset structural response value v 0, and the second preset structure response value v 1.

[0064] Structural design parameters include internal forces of the seismic isolation structure F d And the equivalent stiffness K of the structure except the isolation support components, the equivalent external load is Q .

[0065] The two corresponding preset structure response value errors are:

[0066] as well as .

[0067] Taking a single-degree-of-freedom viscous damper as an example, the internal force of the seismic isolation structure is obtained according to the mechanical formula of the viscous damper: , the two corresponding preset structure response value errors are:

[0068] as well as ;

[0069] in, v 0 is the first preset structure response value, v1 is the second preset structure response value, is the preset structural response value error corresponding to the first preset structural response value, is the preset structural response value error corresponding to the second preset structural response value.

[0070] Taking a single degree of freedom viscous damper as an example, the damping coefficient C of the viscous damper is 100kN / (m / s) 0.2 , viscous damper velocity index α is 0.2, the equivalent stiffness K of the structure except the single degree of freedom damper is 1000kN / m, the displacement, velocity and acceleration of the structure at the previous moment 、 and are all zero, the time step Δt is 0.02, the external load is 10kN, and the equivalent external load Q =10kN. The mechanical equilibrium equation is: , the structural nonlinear solution model is obtained as .

[0071] The first preset structure response value and the second preset structure response value are respectively v 0=0 and v 1=1.

[0072] We obtain: the first preset structure response value v 0=0 corresponding to the preset structural response value error , and the second preset structural response value v 1=1 corresponding to the preset structural response value error .

[0073] S3: Based on the secant method, an approximate structural response value is obtained according to two preset structural response values ​​and their corresponding preset structural response value errors.

[0074] Specifically, the secant equation is .

[0075] In this example, according to the formula , and obtain the approximate structural response value. v 2 is the approximate structural response value.

[0076] Based on the first preset structural response value v 0=0 and the second preset structure response value v 1=1, get the approximate structural response value .

[0077] S41: If the approximate structural response error value corresponding to the approximate structural response value meets the error tolerance requirement, the approximate structural response value is used as the output.

[0078] In this example, the error tolerance is set to ε =10 -6 , if satisfied , then it is believed that v i for The approximate solution of is output as the approximate structural response value.

[0079] For example, in this embodiment, , damping coefficient is 100kN / (m / s) 0.2 , the velocity index is 0.2, the equivalent stiffness of the structure except the single degree of freedom damper is 1000kN / m, the displacement, velocity and acceleration of the structure at the previous moment are zero, the time step is 0.02s, and the equivalent external load Q =10kN, the first preset structural response value and the second preset structural response value are respectively v 0=0 and v 1=1, the error tolerance requirement is .

[0080] Approximate structural response value v 2 Substitution ,get , it is obvious that the approximate structural response error value corresponding to the approximate structural response value does not meet the error tolerance requirement.

[0081] S42: If the approximate structural response error value corresponding to the approximate structural response value does not meet the error tolerance requirement, determine the next two preset structural response values ​​according to the two preset structural response value errors and the approximate structural response error value until the obtained approximate structural response error value meets the error tolerance requirement, and use the approximate structural response value as output.

[0082] In some optional embodiments, when the approximate structural response error value does not meet the error tolerance requirement, the next two preset structural response values ​​are determined according to the two preset structural response value errors and the approximate structural response error value, that is, according to f ( v 0), f ( v 1) and f ( v 2) The symbol state searches for the next iteration value, including:

[0083] A: If f ( v 1) and f ( v 2) Different signs, the next step is to preset the response values ​​of the two structures v 01 = v 1, v 11 =v 2.

[0084] B: If f ( v 1) and f ( v 2) Same sign, and f ( v 2) with f ( v 0) different signs, the next step is to preset the response values ​​of the two structures v 01 = v 0, v 11 = v 2.

[0085] C: If f ( v 0), f ( v 1) and f ( v 2) With the same sign, the next step is to set the response value of the two preset structures to v 01 = v 1, v 11 = v 1+0.5( v 2- v 1).

[0086] in, v 0 is the first preset structure response value, v 1 is the second preset structure response value, v 2 is the approximate structural response value, v 01 For the next step, the first preset structure response value is v 11 Next, the second preset structure response value, is the preset structural response value error corresponding to the first preset structural response value, is the preset structural response value error corresponding to the second preset structural response value, is the preset structural response value error corresponding to the approximate structural response value.

[0087] In this embodiment, , , , meeting the condition B, enter the second step of iteration, that is, enter the first step to select the next preset structural response value according to the preset structural response value error: the next preset structural response value v 01 =0, v 11 =0.0833.

[0088] Calculate the next approximate structural response v 3.

[0089]

[0090] Determine the convergence conditions and check If the conditions are met, v 3 is the approximate solution of the structural dynamic equilibrium model, that is, the approximate structural response, output v 3. If the conditions are not met, select the preset structural response value before the next iteration.

[0091] Proceed to the next step and select the preset structural response value before iteration.

[0092] according to f ( v 01 ), f ( v 11 )and f ( v 3) The symbol state searches for the next preset structural response value and repeats steps S2 to S3.

[0093] , , , meeting the condition B, enter the third step of iteration, that is, enter the second step to select the next preset structural response value according to the preset structural response value error: the next iteration value v 02 =0, v 12 =0.0133.

[0094] Calculate approximate structural response v 4.

[0095] ;

[0096] Determine the convergence conditions and check If the conditions are met, v 4 is the approximate solution of the structural dynamic equilibrium model, that is, the approximate structural response, output v 4. If the conditions are not met, proceed to the next step and select the preset structural response value before iteration.

[0097] Proceed to the next step and select the preset structural response value before iteration.

[0098] according to f ( v 02 ), f (v 12 )and f ( v 4) The symbol state searches for the next iteration value and repeats steps S2 to S3.

[0099] And so on, it converges at the 76th iteration. , iterative convergence. The damper result speed

[0100] By using the structural nonlinear self-finding secant iterative solution method embodiment of the present invention, the structural displacement of the viscous damper with the degree of freedom of the embodiment is obtained. m, the speed is m / s, the convergence accuracy is , the total number of iteration steps is 76. The traditional secant method is used in the initial value of the iteration v 0=0 and v When 1=1, the iterative value oscillates in the interval (-150, 150) and cannot converge, such as Figure 3 As shown, it can be seen that the structural nonlinear self-finding secant iterative solution method of the present invention has no dependence on the initial iteration value of the structure and can better solve the nonlinear structural response.

[0101] In summary, the present invention obtains two corresponding preset structural response value errors based on two preset structural response values; obtains an approximate structural response value based on the two preset structural response values ​​and their corresponding preset structural response value errors; determines whether to use the approximate structural response value as output by judging whether the approximate structural response error value corresponding to the approximate structural response value meets the error tolerance requirement; when the approximate structural response error value corresponding to the approximate structural response value does not meet the error tolerance requirement, determines the next two preset structural response values ​​based on the two preset structural response value errors and the approximate structural response error value, until the obtained approximate structural response error value meets the error tolerance requirement, and uses the approximate structural response value as output, so as to solve the problem that the prior art adopts a trial value iteration method for solving, and if the initial iteration value deviates far from the true solution, or the function presents complex states such as multiple extreme points and non-monotonicity in the neighborhood of the initial value, it may cause the iterative process to diverge, thereby affecting the reliability of the structural response result.

[0102] In a second aspect, an embodiment of the present application further provides a structural response determination device based on nonlinear self-finding secant iteration, comprising: a model building module, an error acquisition module, an approximate value acquisition module and a judgment iteration module.

[0103] The model establishment module is used to establish a structural nonlinear solution model based on the structural dynamic equilibrium state and the structural response relationship; the error acquisition module is used to substitute the structural design parameters and the equivalent external load into the structural nonlinear solution model, and obtain two corresponding preset structural response value errors according to the two preset structural response values; the approximate value acquisition module is used to obtain the approximate structural response value based on the secant method according to the two preset structural response values ​​and their corresponding preset structural response value errors; the judgment iteration module is used to output the approximate structural response value when the approximate structural response error value corresponding to the approximate structural response value meets the error tolerance requirement; when the approximate structural response error value corresponding to the approximate structural response value does not meet the error tolerance requirement, the next two preset structural response values ​​are determined according to the two preset structural response value errors and the approximate structural response error value, until the obtained approximate structural response error value meets the error tolerance requirement, and the approximate structural response value is used as the output.

[0104] Among them, the functional implementation of each module in the above-mentioned structural response determination device based on nonlinear self-finding secant iteration corresponds to the various steps in the above-mentioned structural response determination method embodiment based on nonlinear self-finding secant iteration, and their functions and implementation processes will not be repeated here one by one.

[0105] In a third aspect, an embodiment of the present application provides a structural response determination device based on nonlinear self-finding secant iteration. The structural response determination device based on nonlinear self-finding secant iteration can be a personal computer (PC), a laptop computer, a server, or other device with data processing capabilities.

[0106] Reference Figure 2 , Figure 2 Schematic diagram of the hardware structure of the device for determining the structural response based on nonlinear self-finding secant iteration involved in the embodiment of the present application. In the embodiment of the present application, the device for determining the structural response based on nonlinear self-finding secant iteration may include a processor, a memory, a communication interface, and a communication bus.

[0107] The communication bus may be of any type and is used to interconnect the processor, memory, and communication interface.

[0108] Communication interfaces include input / output (I / O), physical, and logical interfaces, which interconnect components within the nonlinear self-secant iteration-based structural response determination device and other devices (such as other computing devices or user devices). Physical interfaces can include Ethernet, fiber, and ATM interfaces; user devices can include displays and keyboards.

[0109] The memory can be various types of storage media, such as random access memory (RAM), read-only memory (ROM), non-volatile RAM (NVRAM), flash memory, optical storage, hard disk, programmable ROM (PROM), erasable PROM (EPROM), electrically erasable PROM (EEPROM), etc.

[0110] The processor may be a general-purpose processor that can invoke a nonlinear self-secant iteration-based structural response determination program stored in a memory and execute the nonlinear self-secant iteration-based structural response determination method provided in the embodiments of the present application. For example, the general-purpose processor may be a central processing unit (CPU). The method executed when the nonlinear self-secant iteration-based structural response determination program is invoked can be referenced from the various embodiments of the nonlinear self-secant iteration-based structural response determination method of the present application and will not be further described here.

[0111] Those skilled in the art will understand that Figure 2 The hardware structure shown in the figure does not constitute a limitation to the present application and may include more or fewer components than shown in the figure, or a combination of certain components, or a different arrangement of components.

[0112] In a fourth aspect, an embodiment of the present application also provides a computer-readable storage medium.

[0113] The computer-readable storage medium of the present application stores a structural response determination program based on nonlinear self-finding secant iteration, wherein when the structural response determination program based on nonlinear self-finding secant iteration is executed by a processor, the steps of the structural response determination method based on nonlinear self-finding secant iteration as described above are implemented.

[0114] Among them, the method implemented when the structural response determination program based on nonlinear self-finding secant iteration is executed can refer to the various embodiments of the structural response determination method based on nonlinear self-finding secant iteration in this application, and will not be repeated here.

[0115] It should be noted that the serial numbers of the above-mentioned embodiments of the present application are for description only and do not represent the advantages or disadvantages of the embodiments.

[0116] The terms "including" and "having" and any variations thereof in the specification and claims of this application and the above-mentioned drawings are intended to cover non-exclusive inclusion. For example, a process, method, system, product or device that includes a series of steps or units is not limited to the listed steps or units, but optionally includes steps or units that are not listed, or optionally includes other steps or units inherent to these processes, methods, products or devices. The terms "first", "second" and "third" are used to distinguish different objects, etc., and do not represent a sequence, nor do they limit the "first", "second" and "third" to different types.

[0117] In the description of the embodiments of this application, the words "exemplary," "for example," or "for example" are used to indicate examples, illustrations, or descriptions. Any embodiment or design described as "exemplary," "for example," or "for example" in the embodiments of this application should not be construed as being preferred or advantageous over other embodiments or designs. Rather, the use of words such as "exemplary," "for example," or "for example" is intended to present the relevant concepts in a concrete manner.

[0118] In the description of the embodiments of the present application, unless otherwise specified, “ / ” means or, for example, A / B can mean A or B; “and / or” in the text is merely a description of the association relationship of associated objects, indicating that three relationships may exist, for example, A and / or B can mean: A exists alone, A and B exist at the same time, and B exists alone. In addition, in the description of the embodiments of the present application, “multiple” refers to two or more than two.

[0119] In some processes described in the embodiments of the present application, multiple operations or steps are included that appear in a specific order. However, it should be understood that these operations or steps may not be performed in the order in which they appear in the embodiments of the present application or may be performed in parallel. The sequence numbers of the operations are only used to distinguish between different operations, and the sequence numbers themselves do not represent any order of execution. In addition, these processes may include more or fewer operations, and these operations or steps may be performed in sequence or in parallel, and these operations or steps may be combined.

[0120] Through the description of the above embodiments, those skilled in the art can clearly understand that the above-mentioned embodiment methods can be implemented by means of software plus the necessary general hardware platform. Of course, they can also be implemented by hardware, but in many cases the former is a better implementation method. Based on this understanding, the technical solution of this application, or the part that contributes to the existing technology, can be embodied in the form of a software product. This computer software product is stored in a storage medium (such as ROM / RAM, magnetic disk, optical disk) as described above and includes a number of instructions for enabling a terminal device to execute the methods described in each embodiment of this application.

[0121] The above are only preferred embodiments of the present application and do not limit the patent scope of the present application. Any equivalent structure or equivalent process transformation made using the contents of the description and drawings of this application, or directly or indirectly applied in other related technical fields, are also included in the patent protection scope of the present application.

Claims

1. A method for determining structural response based on nonlinear self-finding secant iteration, characterized in that: The following steps are involved: Based on the structural dynamic equilibrium state and structural response relationship, a structural nonlinear solution model is established, including: Establish a structural dynamic equilibrium equation based on the structural dynamic equilibrium state, and establish a structural response equation based on the structural response relationship; establish a structural dynamic equilibrium model based on the structural dynamic equilibrium equation and the structural response equation; establish a structural nonlinear solution model based on the structural dynamic equilibrium model; The structural dynamics equilibrium equation is ; The structural response equation is ; Equivalent structural stiffness ; Equivalent external load ; The structural dynamics equilibrium model is ; The structural nonlinear solution model is ; in, u is the displacement of the structure, v is the displacement velocity of the structure, Δt is the time step, Q is the equivalent external load for structural dynamic analysis, F d is the internal force of the seismic isolation structure, K is the equivalent stiffness of the structure excluding the isolation bearing components, is the structural stiffness matrix, is the structural mass matrix, 、 and are the displacement, velocity and acceleration of the structure at the previous moment, is the external load on the structure at the current moment, is the structural response error; Substituting the structural design parameters and the equivalent external loads into the structural nonlinear solution model, and obtaining two corresponding preset structural response value errors according to the two preset structural response values; Based on the secant method, an approximate structural response value is obtained according to two preset structural response values ​​and their corresponding preset structural response value errors; If the approximate structural response error value corresponding to the approximate structural response value meets the error tolerance requirement, the approximate structural response value is used as the output; If the approximate structural response error value corresponding to the approximate structural response value does not meet the error tolerance requirement, the next two preset structural response values ​​are determined according to the two preset structural response value errors and the approximate structural response error value until the obtained approximate structural response error value meets the error tolerance requirement, and the approximate structural response value is used as the output.

2. The structural response determination method based on nonlinear self-finding secant iteration according to claim 1, characterized in that: The two corresponding preset structure response value errors are: as well as ; in, v 0 is the first preset structure response value, v 1 is the second preset structure response value, is the preset structural response value error corresponding to the first preset structural response value, is the preset structural response value error corresponding to the second preset structural response value.

3. The method for determining structural response based on nonlinear self-finding secant iteration according to claim 2, characterized in that: According to the formula , to obtain the approximate structural response value v 2.

4. The structural response determination method based on nonlinear self-finding secant iteration according to claim 1, characterized in that: The error tolerance requirement is ,in, is the preset structural response value error corresponding to the approximate structural response value, ε is the error tolerance.

5. The structural response determination method based on nonlinear self-finding secant iteration according to claim 4, characterized in that: The error margin ε =10 -6 .

6. The structural response determination method based on nonlinear self-finding secant iteration according to claim 1, characterized in that: According to the same-sign or different-sign relationship between the two preset structural response value errors and the approximate structural response error values, the two preset structural response values ​​for the next step are determined.

7. The structural response determination method based on nonlinear self-finding secant iteration according to claim 6, characterized in that: The step of determining the next two preset structural response values ​​based on the same-sign or different-sign relationship between the two preset structural response value errors and the approximate structural response error values ​​includes: like f ( v 1) and f ( v 2) Different signs, the next step is to preset the response values ​​of the two structures v 01 = v 1, v 11 = v 2; like f ( v 1) and f ( v 2) Same sign, and f ( v 2) with f ( v 0) different signs, the next step is to preset the response values ​​of the two structures v 01 = v 0, v 11 = v 2; like f ( v 0), f ( v 1) and f ( v 2) With the same sign, the next step is to set the response value of the two preset structures to v 01 = v 1, v 11 = v 1+0.5( v 2- v 1); in, v 0 is the first preset structure response value, v 1 is the second preset structure response value, v 2 is the approximate structural response value, v 01 For the next step, the first preset structure response value is v 11 Next, the second preset structure response value, is the preset structural response value error corresponding to the first preset structural response value, is the preset structural response value error corresponding to the second preset structural response value, is the preset structural response value error corresponding to the approximate structural response value.

8. A structural response determination device based on nonlinear self-finding secant iteration, characterized in that: include: The model building module is used to establish a structural nonlinear solution model based on the structural dynamic equilibrium state and structural response relationship, including: Establish a structural dynamic equilibrium equation based on the structural dynamic equilibrium state, and establish a structural response equation based on the structural response relationship; establish a structural dynamic equilibrium model based on the structural dynamic equilibrium equation and the structural response equation; establish a structural nonlinear solution model based on the structural dynamic equilibrium model; The structural dynamics equilibrium equation is ; The structural response equation is ; Equivalent structural stiffness ; Equivalent external load ; The structural dynamics equilibrium model is ; The structural nonlinear solution model is ; in, u is the displacement of the structure, v is the displacement velocity of the structure, Δt is the time step, Q is the equivalent external load for structural dynamic analysis, F d is the internal force of the seismic isolation structure, K is the equivalent stiffness of the structure excluding the isolation bearing components, is the structural stiffness matrix, is the structural mass matrix, 、 and are the displacement, velocity and acceleration of the structure at the previous moment, is the external load on the structure at the current moment, is the structural response error; An error acquisition module is used to substitute the structural design parameters and the equivalent external load into the structural nonlinear solution model, and obtain two corresponding preset structural response value errors according to two preset structural response values; An approximate value acquisition module, which is used to obtain an approximate structural response value based on two preset structural response values ​​and their corresponding preset structural response value errors based on a secant method; The judgment iteration module is used to output the approximate structural response value when the approximate structural response error value corresponding to the approximate structural response value meets the error tolerance requirement; when the approximate structural response error value corresponding to the approximate structural response value does not meet the error tolerance requirement, determine the next two preset structural response values ​​according to the two preset structural response value errors and the approximate structural response error value, until the obtained approximate structural response error value meets the error tolerance requirement, and use the approximate structural response value as the output.

Images