Train dynamic load equivalent simulation calculation method, device, equipment and storage medium

Abstract

The present invention discloses a train dynamic load equivalent simulation calculation method, device, equipment, and storage medium, relating to the field of bridge engineering technology. The method comprises: dividing nodes according to the main beam structure to obtain the number of nodes and the corresponding vertical coordinates; classifying the nodes according to the distance between the target node and adjacent nodes; calculating the equivalent load-time function of different node types based on the train length, speed, train uniform load concentration, and the distance between the target node and adjacent nodes; determining the loading moment and equivalent load-time function of each node based on the position of each node and the equivalent load-time function of different node types; and obtaining a time history curve of the internal force of the member and the node displacement based on the equivalent load-time function of each node type and the loading moment and equivalent load-time function type number of each node. The method solves the problems of inaccurate equivalent simulation calculation results and low computational efficiency in the prior art.

Description

Technical Field

[0001] The present invention relates to the technical field of bridge engineering, and in particular to a train dynamic load equivalent simulation calculation method, device, equipment and storage medium. Background Art

[0002] When engineers deal with the impact of train loads on bridge structures, they usually perform a static analysis of the train's effects on the structure. The specific method is to first calculate the influence line of the structure's stress and deformation under the action of a unit concentrated load, and then determine the most unfavorable position for loading based on the influence line. The entire process is automatically completed by the calculation software, which is fast and efficient. However, the disadvantage of the above method is that it cannot take into account the impact dynamic effect of the train on the bridge structure when the train passes through the bridge, and it also ignores the geometric nonlinear effects caused by the train loading, resulting in a large difference between the theoretical calculation results and the actual structural response. At present, dynamic analysis software is relatively mature. The difficulty in calculating the dynamic response of the structure when the train is running lies in the equivalent simulation of the dynamic load of the train.

[0003] In the train-bridge coupling analysis software, due to the need to consider the coupling effect of train, track, and bridge, the beam elements are divided into very short sections, usually less than 1 meter in length. The train load is applied according to the actual axle load. The train load is concentrated in the load group. Due to the small node spacing, the axle load at each moment can be directly applied to the node below the axle.

[0004] In bridge structure mechanical performance analysis, train load diagrams must be applied according to the ZK or ZKH load diagrams in the "Code for Design of Railway Bridges and Culverts." These load diagrams primarily utilize uniformly distributed loads, differing from the axle-load-based loading method used in actual trains. Therefore, train load simulation methods used in train-bridge coupling analysis cannot be used. Furthermore, mechanical performance analysis models are typically used for multiple scheme comparisons, resulting in lengthy unit lengths. Therefore, a novel equivalent simulation method for dynamic train loads is needed to simplify the equivalent process and improve computational efficiency while ensuring accuracy. Summary of the Invention

[0005] In view of the defects existing in the prior art, the purpose of the present invention is to provide a method, device, equipment and storage medium for equivalent simulation calculation of train dynamic load, which can solve the problems in the prior art of insufficient accuracy of equivalent simulation calculation results of train dynamic load or low calculation efficiency when the results are accurate.

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

[0007] In a first aspect, this solution provides a train dynamic load equivalent simulation calculation method, comprising the following steps:

[0008] Divide the nodes according to the main beam structure and obtain the number of nodes and the corresponding vertical coordinates;

[0009] According to the node number and the corresponding vertical coordinate, the nodes are classified according to the distance between the target node and the adjacent nodes;

[0010] Based on the train length, speed, uniform load concentration of the train, and the distance between the target node and the adjacent nodes, the equivalent load-time function of different types of nodes is calculated;

[0011] Determine the loading time and equivalent load-time function of each node based on the location of each node and the equivalent load-time function of different types of nodes;

[0012] According to the equivalent load-time function of each type of node, the loading moment of each node and the type number of the equivalent load-time function, the time history curves of the internal force of the member and the node displacement are obtained.

[0013] In some optional solutions, dividing the nodes according to the main beam structure and obtaining the number of nodes and the corresponding vertical coordinates include:

[0014] Extract the node numbers and corresponding ordinates of each node on the beam that the train passes in sequence according to the direction of travel, and form a node number vector and a corresponding ordinate vector;

[0015] According to the formula L={l1 … l n} T , obtain the distance vector from the node on the beam to the node at the end of the beam in the direction of train travel;

[0016] Among them, L is the distance vector from the node on the beam to the node at the end of the beam in the direction of train travel, n is the number of nodes, l i =|x i -x1|,X={x1 … x n} T , l i is the distance vector from the i-th node to the beam end node in the train travel direction, X is the vertical coordinate vector of the node, x i is the vertical coordinate vector of the i-th node, i=1~n.

[0017] In some optional solutions, the nodes are classified according to the distance between the target node and the adjacent nodes, including:

[0018] When the distance between the first target node and the rear adjacent node is the same as the distance between the second target node and the rear adjacent node, and the distance between the first target node and the front adjacent node is the same as the distance between the second target node and the front adjacent node, the first target node and the second target node are classified into the same category.

[0019] In some optional solutions, when the distance between the first target node and the rear adjacent node is the same as the distance between the second target node and the rear adjacent node, and the distance between the first target node and the front adjacent node is the same as the distance between the second target node and the front adjacent node, classifying the first target node and the second target node into the same category includes:

[0020] According to the spacing matrix d i =l i+1 -l i , row number vector w = {1 ... n-2} T ;

[0021] The rows with the same elements in the spacing matrix S are classified by two-layer loops, including:

[0022] Outer loop: Select row j S(j,:) from rows 1 to (n-3-k) in ascending order, where j = 1 to (n-3-k);

[0023] Inner loop: Select the i-th row S(i,:) from (n-2-k) to (j+1) in descending order and compare it with S(j,:). If S(i,:) = S(j,:), update the first variable m to make the first variable m+1, and set the w(i)th element f(w(i)) of the identification vector to j. Then delete the i-th row in the matrix S and vector w, and continue to the next inner loop, where f is the identification vector.

[0024] After the inner loop ends, jump out to the outer loop, update the second variable k, make the second variable k+m, and set the w(j)th element f(w(j))=j, update j, make j+1 and enter the inner loop again until all outer loops end;

[0025] After the two-layer loop is completed, a new matrix S without repeated rows, a row number vector w and an identification vector f are obtained;

[0026] Among them, d i is the distance between the i-th node and the adjacent node, k is the second variable, its initial value k=0, m is the first variable, its initial value m=0.

[0027] In some optional solutions, the equivalent load-time function of different types of nodes is calculated based on the train length, speed, train uniform load concentration, and the distance between the target node and the adjacent nodes, including:

[0028] right The element of the matrix S corresponding to the rth row of the rth type equivalent load-time function is S(r,:)=(d w(r)+1 d w(r)+2), calculate the load-time function of the target node, that is, the w(r)+2 node, where r is an element in [1,n-2-k];

[0029] Each type of equivalent load-time function curve starts from time 0, and the uniformly distributed load effect between adjacent nodes is equivalent according to the lever method, and the reaction force of the target node is applied to the target node in the reverse direction as a concentrated load;

[0030] At t∈[t r-1,1 ,t r,1 ] moment, the reaction force of the target node

[0031] At t∈[t r,1 ,t r+1,1 ] moment, the reaction force of the target node

[0032] At t∈[t r+1,1 ,t r-1,2 ] moment, the reaction force of the target node

[0033] At t∈[t r-1,2 ,t r,2 ] moment, the reaction force of the target node

[0034] At t∈[t r,2 ,t r+1,2 ] moment, the reaction force of the target node

[0035] At t∈[t r+1,2 ,∞], the reaction force R of the target node r (t) = 0;

[0036] Where (n-2-k) is the number of rows in the spacing matrix S, that is, there are a total of (n-2-k) equivalent load-time functions, d w(r)+1 is the distance between the w(r)+2th node and its adjacent node, d w(r)+2 is the distance between the w(r)+2th node and the adjacent node in front, D is the train length, v is the speed, p is the uniform load concentration of the train, t r-1,1 =0,t r-1,1 is the time when the vehicle head reaches the adjacent node behind, t r,1 is the time when the train head arrives at the target node, t r+1,1 is the time when the train head arrives at the adjacent node in front, t r-1,2is the time when the rear end of the train arrives at the adjacent node behind, t r,2 is the time when the rear end of the train arrives at the destination node, t r+1,2 is the time when the rear end of the train arrives at the adjacent node in front.

[0037] In some optional solutions, determining the loading time and equivalent load-time function of each node based on the position of each node and the equivalent load-time function of different types of nodes includes:

[0038] Iterate over all elements in the node vector. Its node number is n(q), and the corresponding loading time is The corresponding equivalent load-time function category is f(q-1), that is, the f(q-1)th category equivalent load-time function is taken as the equivalent load-time function of the n(q)th node, where q is an element in [2,n-1].

[0039] In some optional solutions, obtaining the time history curves of the member internal force and the node displacement according to the equivalent load-time function of each type of node, the loading moment of each node, and the equivalent load-time function type number includes:

[0040] A finite element model of the entire bridge was established, and the equivalent load-time function of each type of node, the loading moment of each node, and the type number of the equivalent load-time function were input into the finite element program. Dynamic time history analysis was performed, and the time history curves of the internal force of the member and the node displacement were calculated.

[0041] In a second aspect, the present invention provides a train dynamic load equivalent simulation calculation device, comprising:

[0042] The node number and vertical coordinate acquisition module is used to divide the nodes according to the main beam structure and obtain the node number and corresponding vertical coordinate;

[0043] The node classification module is used to classify nodes according to the node number and the corresponding vertical coordinate, and the distance between the target node and the adjacent nodes;

[0044] The first function acquisition module is used to calculate the equivalent load-time function of different types of nodes based on the train length, speed, train uniform load concentration and the distance between the target node and the adjacent nodes;

[0045] The second function acquisition module is used to determine the loading time and equivalent load-time function of each node according to the position of each node and the equivalent load-time function of different types of nodes;

[0046] The internal force and node displacement time history curve acquisition module is used to obtain the time history curves of the internal force and node displacement of the member according to the equivalent load-time function of each type of node and the loading moment and equivalent load-time function type number of each node.

[0047] In a third aspect, the present solution provides a computer device comprising a processor, a memory, and a computer program stored in the memory and executable by the processor, wherein when the computer program is executed by the processor, the steps of the train dynamic load equivalent simulation calculation method as described in any one of the above items are implemented.

[0048] In a fourth aspect, the present invention provides a computer-readable storage medium having a computer program stored thereon, wherein when the computer program is executed by a processor, the steps of the train dynamic load equivalent simulation calculation method as described in any one of the above items are implemented.

[0049] Compared with the prior art, the advantages of the present invention are as follows: this solution divides nodes according to the main beam structure, obtains the number of nodes and the corresponding vertical coordinates; classifies nodes according to the node number and the corresponding vertical coordinates, and according to the distance between the target node and the adjacent nodes; calculates the equivalent load-time function of different types of nodes based on the train length, speed, train uniform load concentration, and the distance between the target node and the adjacent nodes; determines the loading time and equivalent load-time function of each node based on the position of each node and the equivalent load-time function of different types of nodes; and obtains the time history curve of the internal force of the member and the node displacement based on the equivalent load-time function of each type of node and the loading time and equivalent load-time function type number of each node. By classifying equivalent load-time functions of the same type and applying different loading times to different nodes, the number of equivalent load-time functions is effectively reduced, the workload of data processing is reduced, the calculation process is simplified, and work efficiency is improved. This solves the problem in the prior art that the equivalent simulation calculation results of train dynamic loads are not accurate enough, or the calculation efficiency is low when the results are accurate. BRIEF DESCRIPTION OF THE DRAWINGS

[0050] 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.

[0051] Figure 1 Schematic diagram of a train crossing a bridge in an embodiment of the present invention;

[0052] Figure 2Schematic diagram of a flow chart of a train dynamic load equivalent simulation calculation method according to an embodiment of the present invention;

[0053] Figure 3 Schematic diagram of equivalent load-time function calculated in an example in an embodiment of the present invention;

[0054] Figure 4 is a schematic block diagram of the structure of a computer device in an embodiment of the present invention;

[0055] In the picture: 1. Train; 2. Main beam. DETAILED DESCRIPTION

[0056] To make the purpose, technical solutions, and advantages of the embodiments of this application more clear, the technical solutions in the embodiments of this application will be clearly and completely described below in conjunction with the drawings in the embodiments of this application. Obviously, the described embodiments are part of the embodiments of this application, not all of the embodiments. Based on the embodiments in this application, all other embodiments obtained by ordinary technicians in this field without making creative efforts are within the scope of protection of this application.

[0057] The embodiments of the present invention are further described in detail below with reference to the accompanying drawings.

[0058] Figure 1 Schematic diagram of a train crossing a bridge in an embodiment of the present invention, wherein a train 1 travels from one end of a main beam 2 to the other end. For a target node, a node close to the starting end is a rear node, and a node close to the ending end is a front node.

[0059] like Figure 2 As shown, in a first aspect, the present invention provides a method for calculating equivalent simulation of train dynamic load, comprising the following steps:

[0060] S1: Divide the nodes according to the main beam structure and obtain the number of nodes and the corresponding vertical coordinates.

[0061] Extract the node numbers and corresponding ordinates of each node on the beam that the train passes in sequence according to the direction of travel, and form a node number vector and a corresponding ordinate vector;

[0062] According to the formula L={l1 … l n} T , obtain the distance vector from the node on the beam to the node at the end of the beam in the direction of train travel;

[0063] Among them, L is the distance vector from the node on the beam to the node at the end of the beam in the direction of train travel, n is the number of nodes, l i =|x i -x1|,X={x1 … x n} T , l iis the distance vector from the i-th node to the beam end node in the train travel direction, X is the vertical coordinate vector of the node, x i is the vertical coordinate vector of the i-th node, i=1~n.

[0064] S2: Classify the nodes according to the node number and the corresponding vertical coordinate, and according to the distance between the target node and the adjacent nodes.

[0065] The beam ends are constrained by vertical supports, so the vertical loads on the beam end nodes will not affect the internal forces and deformations of the structure. Therefore, the equivalent load-time function of the end nodes is not considered in the spacing matrix.

[0066] When the distance between the first target node and the rear adjacent node is the same as the distance between the second target node and the rear adjacent node, and the distance between the first target node and the front adjacent node is the same as the distance between the second target node and the front adjacent node, the first target node and the second target node are classified into the same category.

[0067] According to the spacing matrix d i =l i+1 -l i , row number vector w = {1 ... n-2} T ;

[0068] The rows with the same elements in the spacing matrix S are classified by two-layer loops, including:

[0069] Outer loop: Select row j S(j,:) from rows 1 to (n-3-k) in ascending order, where j = 1 to (n-3-k);

[0070] Inner loop: Select the i-th row S(i,:) from (n-2-k) to (j+1) in descending order and compare it with S(j,:). If S(i,:) = S(j,:), update the first variable m to make the first variable m+1, and set the w(i)th element f(w(i)) of the identification vector to j. Then delete the i-th row in the matrix S and vector w, and continue to the next inner loop, where f is the identification vector.

[0071] After the inner loop ends, jump out to the outer loop, update the second variable k, make the second variable k+m, and set the w(j)th element f(w(j))=j, update j, make j+1 and enter the inner loop again until all outer loops end;

[0072] After the two-layer loop is completed, a new matrix S without repeated rows, a row number vector w and an identification vector f are obtained;

[0073] Among them, d iis the distance between the i-th node and the adjacent node, k is the second variable, its initial value k=0, m is the first variable, its initial value m=0.

[0074] S3: Based on the train length, speed, uniformly distributed train load concentration, and the distance between the target node and adjacent nodes, the equivalent load-time function of different types of nodes is calculated.

[0075] right The element of the matrix S corresponding to the rth row of the rth type equivalent load-time function is S(r,:)=(d w(r)+1 d w(r)+2 ), calculate the load-time function of the target node, that is, the w(r)+2 node, where r is an element in [1,n-2-k];

[0076] Each type of equivalent load-time function curve starts from time 0, and the uniformly distributed load effect between adjacent nodes is equivalent according to the lever method, and the reaction force of the target node is applied to the target node in the reverse direction as a concentrated load;

[0077] At t∈[t r-1,1 ,t r,1 ], the train load partially enters the interval between the target node and the adjacent node behind it, but does not enter the interval between the target node and the adjacent node in front of it. The reaction force of the target node is

[0078] At t∈[t r,1 ,t r+1,1 ], the train load covers the interval between the target node and the adjacent node behind it, and partially enters the interval between the target node and the adjacent node in front of it. The reaction force of the target node is

[0079] At t∈[t r+1,1 ,t r-1,2 ], the train load completely covers the interval of the adjacent nodes in front and behind the target node, and the reaction force of the target node

[0080] At t∈[t r-1,2 ,t r,2 ], the train load covers the interval between the target node and the adjacent node in front, and partially leaves the interval between the target node and the adjacent node in the back. The reaction force of the target node

[0081] At t∈[t r,2 ,t r+1,2 ], the train load completely leaves the interval between the target node and the adjacent node behind, and partially leaves the interval between the target node and the adjacent node in front. The reaction force of the target node

[0082] At t∈[t r+1,2 ,∞], the train load has completely left the interval of the adjacent nodes in front and behind the target node, and the reaction force R r (t) = 0;

[0083] Where (n-2-k) is the number of rows in the spacing matrix S, that is, there are a total of (n-2-k) equivalent load-time functions, d w(r)+1 is the distance between the w(r)+2th node and its adjacent node, d w(r)+2 is the distance between the w(r)+2th node and the adjacent node in front, D is the train length, v is the speed, p is the uniform load concentration of the train, t r-1,1 =0,t r-1,1 is the time when the vehicle head reaches the adjacent node behind, t r,1 is the time when the train head arrives at the target node, t r+1,1 is the time when the train head arrives at the adjacent node in front, is the time when the rear end of the train arrives at the adjacent node behind, t r,2 is the time when the rear end of the train arrives at the destination node, t r+1,2 is the time when the rear end of the train arrives at the adjacent node in front.

[0084] S4: Determine the loading time and equivalent load-time function of each node based on the position of each node and the equivalent load-time function of different types of nodes.

[0085] Iterate over all elements in the node vector. Its node number is n(q), and the corresponding loading time is The corresponding equivalent load-time function category is f(q-1), that is, the f(q-1)th category equivalent load-time function is taken as the equivalent load-time function of the n(q)th node, where q is an element in [2,n-1].

[0086] S5: According to the equivalent load-time function of each type of node, the loading moment of each node and the type number of the equivalent load-time function, the time history curves of the internal force of the member and the node displacement are obtained.

[0087] A finite element model of the entire bridge was established, and the equivalent load-time function of each type of node, the loading moment of each node, and the type number of the equivalent load-time function were input into the finite element program. Dynamic time history analysis was performed, and the time history curves of the internal force of the member and the node displacement were calculated.

[0088] In this embodiment, nodes are divided according to the main beam structure, and the node number and corresponding ordinate are obtained. Nodes are classified according to the node number and corresponding ordinate, based on the distance between the target node and adjacent nodes. Equivalent load-time functions for different node types are calculated based on train length, speed, uniformly distributed train load concentration, and the distance between the target node and adjacent nodes. The loading moment and equivalent load-time function of each node are determined based on its position and the equivalent load-time function of each node type. Time history curves of the internal force and node displacement of the member are obtained based on the equivalent load-time function of each node type, the loading moment of each node, and the type number of the equivalent load-time function. By categorizing equivalent load-time functions of the same type and applying different loading moments to different nodes, the number of equivalent-time functions is effectively reduced, the data processing workload is reduced, the calculation process is simplified, and work efficiency is improved. The proposed equivalent load-time function fully considers the load distribution relationship between adjacent nodes, avoids the problem of duplicate loading of two adjacent nodes, and improves calculation accuracy. The invention solves the problem in the prior art that the equivalent simulation calculation results of the train dynamic load are not accurate enough or the calculation efficiency is low when the results are accurate.

[0089] The following is a specific calculation example to facilitate understanding of this solution.

[0090] In this example, the main beam is divided into five units, with a total of six nodes, or n = 6. The train load concentration is p = 64 kN / m, the train length is D = 550 m, and the speed is v = 250 km / h. The train enters from the left, and the time when the front of the train enters the left node is 0. The vector formed by the node numbers on the beam from left to right is n = {101 102 103 104 105 106} T The distance vector L from the left-to-right beam node to the beam end node in the train travel direction is {14 26 40 52 66} T , the original spacing matrix The number of rows is n-2=4, and the row number vector w={1 2 3 4} T .

[0091] The rows with the same elements in the distance matrix S are classified by two-layer loop. Initially, the second variable k=0.

[0092] Outer loop: (n-3-k)=3, select the j-th row S(j,:) from rows 1 to 3 in ascending order.

[0093] In the first outer loop, let j = 1, get S(1,:), initially take the first variable m = 0, and enter the inner loop:

[0094] Select the i-th row from (n-2-k) to (j+1) in descending order.

[0095] First let i=n-2-k=4, and compare S(4,:) with S(1,:), they are not equal.

[0096] Then let i-1, so that i=3, compare S(3,:) and S(1,:), they are equal. Let the first variable m+1, so that m=1, and let the w(3)th element f(3) of the identity vector = 1, then delete the 3rd row in the matrix S and vector w. The new matrix is

[0097] Then let i-1, so that i=2, and compare S(2,:) with S(1,:), they are not equal.

[0098] The first inner loop ends and jumps to the outer loop, setting k+m, making k=1, and setting the w(1)th element f(1)=1 of the identification vector.

[0099] Second outer loop: Let j+1 again, make j=2, get S(2,:), initial variable m=0, enter the inner loop:

[0100] Select the i-th row from (n-2-k) to (j+1) in descending order,

[0101] First, let i = n-2-k = 3, and compare S(3,:) with S(2,:), which are equal. Let the first variable m+1, so that m=1, and let the w(3)th element f(4) of the identity vector = 2, then delete the third row in the matrix S and the vector w. The new matrix is

[0102] Since (n-2-k)=(j+1)=3, the second inner loop ends and jumps to the outer loop, setting k+m, making k=2, and setting the w(2)th element f(2)=2 of the identification vector.

[0103] At this point, (n-3-k)=j=2, the outer loop is completed, and all loops end. After the two layers of loops are completed, a new matrix with no repeated rows is obtained. Row number vector w = {12} T and identity vector f = {1212} T .

[0104] The first type of equivalent load-time function, i.e. r = 1, corresponds to the first row in the new matrix S. The distance from the node to the adjacent node behind is 14m, and the distance to the adjacent node in front is 12m. Let the time when the vehicle head reaches the adjacent node behind be t r-1,1 = 0, the time when the train head arrives at the node of interest is t r,1 = 0.2016s, the time when the train head arrives at the adjacent node in front is t r+1,1 = 0.3744s. The time when the rear of the train arrives at the adjacent node behind is tr-1,2 =7.9200s, the time when the rear end of the train arrives at the node of interest is t r,2 = 8.1216s, the time when the rear end of the train arrives at the adjacent node in front is t r+1,2 =8.2944s. The corresponding equivalent load-time function is as follows Figure 3 The second type of equivalent load-time function can be obtained by the same method.

[0105] Iterate over all elements in the node vector. The node number is n(q), and the loading time at the corresponding node is The corresponding equivalent load-time function category is f(q-1), that is, take the f(q-1)th category of equivalent load-time function.

[0106] The matrix consisting of the node number, loading time t(q), and equivalent load-time function type of all nodes except the beam end nodes is:

[0107] In a second aspect, the present invention provides a train dynamic load equivalent simulation calculation device, comprising:

[0108] The node number and vertical coordinate acquisition module is used to divide the nodes according to the main beam structure and obtain the node number and corresponding vertical coordinate;

[0109] The node classification module is used to classify nodes according to the node number and the corresponding vertical coordinate, and the distance between the target node and the adjacent nodes;

[0110] The first function acquisition module is used to calculate the equivalent load-time function of different types of nodes based on the train length, speed, train uniform load concentration and the distance between the target node and the adjacent nodes;

[0111] The second function acquisition module is used to determine the loading time and equivalent load-time function of each node according to the position of each node and the equivalent load-time function of different types of nodes;

[0112] The internal force and node displacement time history curve acquisition module is used to obtain the time history curves of the internal force and node displacement of the member according to the equivalent load-time function of each type of node and the loading moment and equivalent load-time function type number of each node.

[0113] It should be noted that those skilled in the art will clearly understand that, for the convenience and brevity of description, the specific working processes of the above-described devices and modules and units can refer to the corresponding processes in the aforementioned embodiments and will not be repeated here.

[0114] The apparatus provided in the above embodiment can be implemented in the form of a computer program. The computer program can be used in Figure 4 Runs on the computer equipment shown.

[0115] See also Figure 4 , Figure 4 This is a schematic block diagram of the structure of a computer device provided in an embodiment of the present application. The computer device may be a terminal.

[0116] like Figure 4 As shown, the computer device includes a processor, a memory, and a network interface connected via a system bus, wherein the memory may include a non-volatile storage medium and an internal memory.

[0117] The non-volatile storage medium can store an operating system and a computer program. The computer program includes program instructions, which, when executed, can enable the processor to execute any train dynamic load equivalent simulation calculation method.

[0118] The processor is used to provide computing and control capabilities and support the operation of the entire computer equipment.

[0119] The internal memory provides an environment for the operation of the computer program in the non-volatile storage medium. When the computer program is executed by the processor, the processor can execute any train dynamic load equivalent simulation calculation method.

[0120] The network interface is used for network communication, such as sending assigned tasks, etc. Those skilled in the art will understand that Figure 4 The structure shown in the figure is only a block diagram of a part of the structure related to the solution of the present application, and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may include more or fewer components than shown in the figure, or combine certain components, or have a different component arrangement.

[0121] It should be understood that the processor may be a central processing unit (CPU), or other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. The general-purpose processor may be a microprocessor or any conventional processor, etc.

[0122] In one embodiment, the processor is used to run a computer program stored in a memory to implement the steps of the above-mentioned train dynamic load equivalent simulation calculation method.

[0123] An embodiment of the present application further provides a computer-readable storage medium, on which a computer program is stored. The computer program includes program instructions. The method implemented when the program instructions are executed can refer to the various embodiments of the present application.

[0124] The computer-readable storage medium may be an internal storage unit of the computer device described in the aforementioned embodiment, such as a hard disk or memory of the computer device. The computer-readable storage medium may also be an external storage device of the computer device, such as a plug-in hard disk, a SmartMedia Card (SMC), a Secure Digital (SD) card, a flash memory card, etc., equipped on the computer device.

[0125] In the description of this application, it should be noted that the terms "upper" and "lower" and the like indicate orientations or positional relationships based on the orientations or positional relationships shown in the accompanying drawings, and are only for the convenience of describing this application and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore cannot be understood as a limitation on this application. Unless otherwise clearly specified and limited, the terms "installed", "connected", and "connected" should be understood in a broad sense, for example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection or an electrical connection; it can be a direct connection, or an indirect connection through an intermediate medium, or it can be internal communication between two elements. For those of ordinary skill in the art, the specific meanings of the above terms in this application can be understood according to the specific circumstances.

[0126] It should be noted that, in this application, relational terms such as "first" and "second" are used only to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any actual relationship or order between these entities or operations. Moreover, the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device comprising a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, an element defined by the sentence "comprising a ..." does not exclude the presence of other identical elements in the process, method, article or device comprising the element.

[0127] The foregoing is merely a list of specific embodiments of the present application, intended to enable those skilled in the art to understand or implement the present application. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the general principles defined herein may be implemented in other embodiments without departing from the spirit or scope of the present application. Therefore, the present application is not limited to the embodiments shown herein, but is intended to conform to the broadest scope consistent with the principles and novel features of the present application.

Claims

1. A train dynamic load equivalent simulation calculation method, characterized in that: The following steps are involved: Divide the nodes according to the main beam structure and obtain the number of nodes and the corresponding vertical coordinates; According to the node number and the corresponding vertical coordinate, the nodes are classified according to the distance between the target node and the adjacent nodes; Based on the train length, speed, uniform load concentration of the train, and the distance between the target node and the adjacent nodes, the equivalent load-time function of different types of nodes is calculated; Determine the loading time and equivalent load-time function of each node based on the location of each node and the equivalent load-time function of different types of nodes; According to the equivalent load-time function of each type of node, the loading moment of each node and the type number of the equivalent load-time function, the time history curve of the internal force of the member and the node displacement is obtained; The node classification according to the distance between the target node and the adjacent nodes includes: When the distance between the first target node and the rear adjacent node is the same as the distance between the second target node and the rear adjacent node, and the distance between the first target node and the front adjacent node is the same as the distance between the second target node and the front adjacent node, the first target node and the second target node are classified into the same category.

2. The train dynamic load equivalent simulation calculation method according to claim 1, characterized in that: The method of dividing the nodes according to the main beam structure and obtaining the number of nodes and the corresponding vertical coordinates includes: Extract the node numbers and corresponding ordinates of each node on the beam that the train passes in sequence according to the direction of travel, and form a node number vector and a corresponding ordinate vector; According to the formula , obtain the distance vector from the node on the beam to the node at the end of the beam in the direction of train travel; in, is the distance vector from the node on the beam to the node at the end of the beam in the direction of train travel, is the number of nodes, , , For the The distance vector from the node No. to the beam end node in the direction of train travel, is the vertical coordinate vector of the node, For the The vertical coordinate vector of node number, .

3. The train dynamic load equivalent simulation calculation method according to claim 1, characterized in that: When the distance between the first target node and the rear adjacent node is the same as the distance between the second target node and the rear adjacent node, and the distance between the first target node and the front adjacent node is the same as the distance between the second target node and the front adjacent node, classifying the first target node and the second target node into the same category includes: According to the spacing matrix d i =l i+1 -l i , row number vector ; Cycle the spacing matrix in two layers The rows with the same elements are grouped into categories, including: Outer loop: from Select the row from smallest to largest OK ,in, ; Inner loop: from Select the first one from largest to smallest OK , and with For comparison, if , update the first variable , so that the first variable , and the first Elements , then in the matrix and vector Delete the Line, continue to the next inner loop, where is the identity vector; After the inner loop ends, jump out to the outer loop and update the second variable , so that the second variable , and the first Elements ,renew ,make Enter the inner loop again until all outer loops are completed; After the two layers of loops, a new matrix with no repeated rows is obtained. , row number vector and logo vector ; in, For the The distance between the node and the adjacent node, is the second variable, and its initial value , is the first variable, and its initial value .

4. The train dynamic load equivalent simulation calculation method according to claim 3, characterized in that: The equivalent load-time functions of different types of nodes are calculated based on train length, speed, train uniform load concentration, and the distance between the target node and the adjacent nodes, including: right , No. The equivalent load-time function corresponds to The elements of the row matrix S are , calculate the target node, which is The load-time function of the node is: for Elements in Each type of equivalent load-time function curve starts from time 0, and the uniformly distributed load effect between adjacent nodes is equivalent according to the lever method, and the reaction force of the target node is applied to the target node in the reverse direction as a concentrated load; At t∈[t r-1，1 , t r，1 ] moment, the reaction force of the target node ; At t∈[t r，1 , t r+1，1 ] moment, the reaction force of the target node ; exist The reaction force of the target node during the time ; exist The reaction force of the target node during the time ; exist The reaction force of the target node during the time ; exist The reaction force of the target node during the time ; in, is the spacing matrix The number of rows, that is, there are An equivalent load-time function, For the The distance between the node and the adjacent node behind it, For the The distance between the node and the adjacent node in front, is the train length, is the vehicle speed, is the uniformly distributed load concentration of the train, , is the time when the vehicle head reaches the adjacent node behind, , is the time when the train head arrives at the target node, , is the time when the train head arrives at the adjacent node in front, , is the time when the rear end of the train arrives at the adjacent node behind, , is the time when the rear end of the train arrives at the destination node, , is the time when the rear end of the train arrives at the adjacent node in front.

5. The train dynamic load equivalent simulation calculation method according to claim 4, characterized in that: Determining the loading time and equivalent load-time function of each node based on the position of each node and the equivalent load-time function of different types of nodes includes: Traverse all elements in the node vector and , whose node number is , the corresponding loading time is , the corresponding equivalent load-time function category is , that is, take the Equivalent load-time function as the first The equivalent load-time function of node No. is: for Elements in .

6. The train dynamic load equivalent simulation calculation method according to claim 1, characterized in that: The method of obtaining the time history curve of the internal force of the member and the node displacement according to the equivalent load-time function of each type of node, the loading moment of each node, and the type number of the equivalent load-time function includes: A finite element model of the entire bridge was established, and the equivalent load-time function of each type of node, the loading moment of each node, and the type number of the equivalent load-time function were input into the finite element program. Dynamic time history analysis was performed, and the time history curves of the internal force of the member and the node displacement were calculated.

7. A train dynamic load equivalent simulation calculation device, characterized in that: include: The node number and vertical coordinate acquisition module is used to divide the nodes according to the main beam structure and obtain the node number and corresponding vertical coordinate; The node classification module is used to classify nodes according to the node number and the corresponding vertical coordinate, and the distance between the target node and the adjacent nodes; The first function acquisition module is used to calculate the equivalent load-time function of different types of nodes based on the train length, speed, train uniform load concentration and the distance between the target node and the adjacent nodes; The second function acquisition module is used to determine the loading time and equivalent load-time function of each node according to the position of each node and the equivalent load-time function of different types of nodes; The module for obtaining the time history curve of internal force and node displacement is used to obtain the time history curve of member internal force and node displacement according to the equivalent load-time function of each type of node, the loading moment of each node and the type number of the equivalent load-time function; The node classification according to the distance between the target node and the adjacent nodes includes: When the distance between the first target node and the rear adjacent node is the same as the distance between the second target node and the rear adjacent node, and the distance between the first target node and the front adjacent node is the same as the distance between the second target node and the front adjacent node, the first target node and the second target node are classified into the same category.

8. A computer device, characterized in that: The computer device includes a processor, a memory, and a computer program stored in the memory and executable by the processor, wherein when the computer program is executed by the processor, the steps of the train dynamic load equivalent simulation calculation method as described in any one of claims 1 to 6 are implemented.

9. A computer-readable storage medium, characterized in that The computer-readable storage medium stores a computer program, wherein when the computer program is executed by a processor, the steps of the train dynamic load equivalent simulation calculation method as claimed in any one of claims 1 to 6 are implemented.

Images