A method, system and device for calculating prestress of cable dome structure

Abstract

The present invention provides a method, system and device for calculating the prestress of a cable dome structure, which relates to the technical field of prestress analysis. The method mainly includes: collecting the structural information of the cable dome and expressing it in a parameterized manner to obtain a node dictionary and an adjacency matrix; initializing a data structure; traversing the node dictionary to extract the first node of each ring and each layer; obtaining the connection node of the first node; traversing the adjacency matrix to find other nodes connected to the specified node; calculating the direction vector and component parameters of the first node and its connection nodes; calculating the prestress of the first node; based on the ring where the node is located and the node type, calculating the internal force of the component on each node and constructing a prestress dictionary; through finite element analysis, cyclically iterating and calculating the prestress of each node until the final design value of the cable prestress that meets the requirements is obtained. The scheme can automatically, quickly and accurately perform prestress analysis, thereby helping to speed up the progress of the project, reduce the cost of the project, and improve the quality and safety of the project.

Description

Technical Field

[0001] The present invention relates to the technical field of prestress analysis, and in particular to a method, system and device for calculating the prestress of a cable dome structure. Background Art

[0002] Cable dome structure is a common large-span building structure. In theory, it mainly includes edges and nodes; the edges are various components, mainly cables and rods; the nodes are the connection points between various components. In actual engineering, Figure 1 As shown, the cable dome structure is often divided into two layers: the upper structure mainly refers to the geometric structure similar to a single-layer spherical lattice shell, which is composed of struts (central struts 1 and circumferential struts 2) and ridge cables connected to their corresponding upper end points 6; the lower structure is mainly composed of the lower end points of the struts and other rods, nodes (cable support nodes 5) and cables (circular cables 4 and oblique cables 7) outside the upper structure. The upper structure or the lower structure also includes a concentric ring structure centered on the dome, and the concentric ring structure is divided into a central ring structure and a non-central ring structure according to different positions.

[0003] Replacing the ridge cables in the cable dome structure with the cable-through ridge rods 3 can help greatly improve the overall performance of the cable dome. However, when calculating the prestress value of the cable-through ridge rods 3 for such a cable dome structure, complex manual calculations are often required, involving a large number of parameters and formulas, which consumes a lot of time and manpower, is prone to errors, and brings quality and safety risks to the construction. Summary of the invention

[0004] The object of the present invention is to provide a method, system and device for calculating the prestress of a cable dome structure, so as to solve at least one of the above-mentioned technical problems existing in the prior art.

[0005] In a first aspect, in order to solve the above technical problems, the present invention provides a method for calculating the prestress of a cable dome structure, comprising the following steps:

[0006] Step 1: Receive basic structural data of the cable dome structure; group and number each node according to the structure in which it is located; generate the point information node_dict of each node and the adjacency matrix between each node based on the node number ID, and build a visual inspection model of the cable dome structure; the basic structural data includes: structural topology information, span , Rise-Span Ratio , Ring division number and radial division number etc.; the structural topology information refers to the type of cable dome structure; the annular division number refers to the number of equal parts into which the entire circumference is divided along the annular direction of the dome; the radial division number refers to the number of rings into which the entire dome is divided from the center of the dome outward;

[0007] Step 2, receiving known mechanical parameters of the cable dome structure; the known mechanical parameters include the prestress value of the central strut, the node load and the prestress value of the prestressed ridge pole, etc.;

[0008] Step 3: Extract the nodes that need to be calculated from the point information node_dict and generate a dictionary of nodes that need to be calculated; traverse each node in the point information node_dict, extract the first node ID of each ring in each layer structure, and store it in the first node dictionary first_nodes;

[0009] Step 4: Extract the first nodes one by one from the first node dictionary first_nodes, and check whether the node belongs to the dictionary of nodes to be calculated: if so, execute step 5; if not, execute step 4;

[0010] Step 5: From the adjacency matrix, find other nodes connected to the current first node as slave nodes, and add the slave node ID to the sublist of the connected node dictionary point_dic_j_ids;

[0011] Step 6, obtain the coordinates of the node and the ring it is in through the node ID; based on the connected node dictionary point_dic_j_ids, traverse and calculate the direction vector, component length, height difference and other information between the current head node and each slave node; store the direction vector in the direction vector dictionary direction_vectors;

[0012] Step 7: Based on the point information node_dict and the connected node dictionary point_dic_j_ids, obtain the layer structure and ring where the node is located and make a judgment: if the current first node is located in the upper structure, execute step 8; if the current first node is located in the lower structure, execute step 9;

[0013] Step 8, if the node is located in the central ring, search the corresponding direction vector from the direction vector dictionary direction_vectors, and calculate the prestressed ridge rod internal force based on the known node load; if the node is located in the non-central ring, search the corresponding direction vector from the direction vector dictionary direction_vectors, and calculate the prestressed ridge rod internal force and the strut internal force based on the known node load; execute step 10; the internal force represents the force generated inside each component of the cable dome structure after the cable dome structure is balanced;

[0014] Step 9: If the node is located in the central ring, search the corresponding direction vector from the direction vector dictionary direction_vectors, and calculate the internal force of the oblique cable based on the known node load; if the node is located in the non-central ring, search the corresponding direction vector from the direction vector dictionary direction_vectors, and calculate the internal force of the oblique cable and the ring cable based on the known node load; execute step 10;

[0015] Step 10: Update the prestress dictionary based on the internal force values ; Iterate step 4 until all the first nodes are traversed; the prestress dictionary The spine prestressing assembly including each ring 、The prestressing force of each ring , the prestressed cable of each ring And the prestressing set of the struts of each ring :

[0016] in, It can be expressed as:

[0017] ;

[0018] It can be expressed as:

[0019] ;

[0020] It can be expressed as:

[0021] ;

[0022] It can be expressed as:

[0023] ;

[0024] Step 11: Traverse the prestress dictionary For each key-value pair in the key-value pair, the key in the key-value pair represents the component type, and the value in the key-value pair represents the prestress value of the component at different rings, that is, the internal force of the tensioning cable required in the tensioning construction; based on the key-value pair, a prestress distribution diagram is drawn and a data label is annotated on each data point of the prestress distribution diagram for display, so as to achieve visualization;

[0025] Step 12: Calculate the initial value and design value of the cable prestress in equilibrium by using a Q network (Q-network or Deep Q-network, DQN); the cable prestress includes oblique cable prestress and ring cable prestress, which is specifically determined by the design document and the construction plan; construct a geometric model of the cable dome structure by using a finite element tool; mesh the geometric model to obtain a finite element model; calculate physical quantities such as stress and deformation of the cable dome structure by using finite element analysis and then display them; adjust the geometric model parameters according to external feedback until the design value of the cable prestress that meets the preset conditions is obtained as the final design value.

[0026] Through the above steps, the final design value of the prestress of the cable dome structure can be automatically calculated, which saves a lot of manpower and time, reduces calculation errors, and improves design efficiency and engineering quality.

[0027] In a feasible implementation manner, the step 1 specifically includes:

[0028] Step 101, in the Python environment, based on the input parameters of the cable dome structure, the design parameters of the cable dome structure such as the spherical radius and the unit inclination angle are obtained through the proportional relationship of similar triangles and the spherical coordinate system calculation; based on the rectangular coordinate system, the design parameters are converted into the position information dictionary node_dict of the nodes in the cable dome structure for subsequent programming indexing, and stored in the node list; based on the node list, with reference to the cable dome structure, the component information of the cable-through ridge pole, the strut, the ring cable and the oblique cable is added to construct a three-dimensional list;

[0029] Step 102: Generate the components such as the cable-through ridge pole, the strut pole, the ring rope and the oblique rope of each ring in sequence based on the three-dimensional list; then construct a two-dimensional adjacency matrix of the connection relationship of each node in each ring; and then splice the two-dimensional adjacency matrix into a three-dimensional adjacency matrix according to the order of the rings;

[0030] Step 103, determining the positions of components such as cable-through ridge poles, struts, ring cables, and oblique cables layer by layer and ring by ring; performing three-dimensional visualization processing on the cable dome structure through the matplotlib library to obtain a visualization inspection model of the cable dome structure, thereby facilitating homomorphic adjustment and preview;

[0031] Step 104: Establishing a corresponding relationship between node numbers, node coordinates, a three-dimensional adjacency matrix, and a cable dome structure visualization inspection model and storing the relationship in a database so that the node coordinates and the adjacency matrix can be queried and updated according to the node numbers later;

[0032] Step 105: Based on the design requirement data, manually determine whether the cable dome structure visualization inspection model meets the requirements: if so, output information about the specified nodes and their connection relationships; if not, output unqualified information and receive adjustment values ​​of input parameters, and update the cable dome structure visualization inspection model until the requirements are met.

[0033] Through the above steps, the input parameters of the cable dome structure can be directly converted into an editable visual model in the Python environment, eliminating the data conversion steps between the modeling environment and the architectural software, which can greatly improve the design efficiency and quality of the cable dome structure.

[0034] In a feasible implementation manner, the input parameters in step 101 include span, aspect ratio, radial division number, circumferential division number, angle between the inclined cable and the horizontal plane, etc.; the radial division number refers to the number of rings into which the entire dome is divided from the center of the dome outward; the circumferential division number refers to the number of equal parts into which the entire circumference is divided along the circumferential direction of the dome.

[0035] In a feasible implementation, the nodes are divided into upper structure nodes and lower structure nodes, the point type mark of the upper structure node is 0, and the point type mark of the lower structure node is 1; the location information dictionary node_dict includes node numbers and node coordinates.

[0036] In a feasible implementation manner, the three-dimensional list specifically includes:

[0037] The first dimension represents the type of rod, including: cable-through ridge rod, strut, ring cable and oblique cable;

[0038] The second dimension represents the number of radial divisions, including the ordinal number of the ring;

[0039] The third dimension represents the connection relationship between nodes. A value of 0 indicates no connection, and a value of 1 indicates connection.

[0040] In a feasible implementation manner, step 101 specifically includes the following steps:

[0041] Step 1011: Calculate the spherical radius of the cable dome structure based on the spherical coordinate system. and unit inclination , the specific formula can be:

[0042] ;

[0043] ;

[0044] in, It represents the coverage radius of the cable dome structure. The specific calculation formula can be:

[0045] ;

[0046] It represents the sagittal height of the cable dome structure. The specific calculation formula can be:

[0047] ;

[0048] Step 1012: Calculate the vertical heights of nodes in adjacent rings based on the proportional relationship of similar triangles. The specific formula may be:

[0049] ;

[0050] ;

[0051] in, Indicates The vertical height of the nodes in the ring; Indicates The vertical height of the nodes in the ring; Indicates The specific formula for the proportional coefficient of the ring is:

[0052] ;

[0053] Indicates The specific formula for the proportional coefficient of the ring is:

[0054] ;

[0055] Step 1013: Calculate the radial coordinates of the nodes in the adjacent rings. The specific formula may be:

[0056] ;

[0057] ;

[0058] in, Indicates The radial coordinates of the nodes in the ring; Indicates The radial coordinates of the nodes in the ring;

[0059] Calculate the radial difference and store it in the radial difference list The specific formula can be:

[0060] ;

[0061] in, Indicates The nodes in the ring and The radial difference between nodes in the ring;

[0062] Step 1014, calculate the height of the support pole, the specific formula can be:

[0063] ;

[0064] in, Indicates The height of the poles in the ring; Indicates The height of the brace in the ring in the superstructure can be expressed as:

[0065] ;

[0066] in, Indicates The angle between the inclined cable and the horizontal plane in each ring;

[0067] Indicates The height of the brace in the lower structure in a ring can be expressed as:

[0068] ;

[0069] Step 1015: Calculate the node coordinates based on the rectangular coordinate system ( , ) and the corresponding coordinates of the lower nodes of the struts ( , ), the specific formula can be:

[0070] when When it is an even number:

[0071] ;

[0072] in, Indicates the x-coordinate of the node; Indicates the y-coordinate of the node; Indicates the current ring equal parts;

[0073] And there are:

[0074] ;

[0075] when When it is an odd number:

[0076] ;

[0077] And there are:

[0078] ;

[0079] Define the mapping relationship between node numbers and rectangular coordinates. The specific expression can be:

[0080] ;

[0081] in, Indicates the node number. The topmost node of the cable dome structure can be taken as the first node, and the numbers can be arranged in sequence in a counterclockwise spiral outward direction to the last node.

[0082] In a feasible implementation manner, step 102 specifically includes the following steps:

[0083] Step 1021, respectively define the adjacency matrices of the rod types, namely: the through-cable ridge rod adjacency matrix adj_matrix_RB, the ring-cable adjacency matrix adj_matrix_CC, the oblique-cable adjacency matrix adj_matrix_DC and the strut adjacency matrix adj_matrix_S;

[0084] Each adjacency matrix is ​​a three-dimensional array:

[0085] The first dimension represents the ordinal number of the ring;

[0086] The second and third dimensions form a square matrix, which represents the connection relationship between all nodes in the current ring. The specific setting rules are: when the element value is 0, it means no connection; when the element value is 1, it means there is a connection;

[0087] Initialize the three-dimensional array and set all elements to 0;

[0088] In this way, complex structural data can be clearly organized and managed through the three-dimensional array, so that the connection status of each ring can be independently stored and called for subsequent analysis and operation;

[0089] Step 1022, traverse and select any two nodes (point1 and point2) to form a node pair, find the node numbers of the two nodes based on the position information dictionary node_dict and determine the connection relationship, if there is a connection relationship, generate a unique member number ID according to the sequence number of the ring, the node number and the member type, so as to distinguish the members in the subsequent sequence; set the connection relationship and member number ID of the two nodes in the adjacency matrix of the ring where the two nodes are located;

[0090] Step 1023: Store the adjacency matrix of the member type in the location information dictionary node_dict for subsequent access and calling.

[0091] In a feasible implementation manner, the specific method for determining the connection relationship in step 1022 includes:

[0092] Step 10221, according to the ring where the node is located and the point type mark, define the node classification rules: a superstructure node in the same ring is connected to the nearest substructure node to form a strut; in adjacent rings, the two nearest superstructure nodes are connected to form a cable-through ridge pole; in adjacent rings, a substructure node in the previous ring is connected to the nearest superstructure node in the next ring to form an oblique cable; in the same ring, two adjacent substructure nodes are connected to form a ring cable;

[0093] Step 10222: In the location information dictionary node_dict, traverse and construct node pairs; based on the point type mark of the node and the equal part in which it is located, determine the connection relationship of the node pair through the node classification rule; according to the classification, mark the element value of the connection relationship in the corresponding adjacency matrix.

[0094] In a feasible implementation, the specific method for determining the position of the cable-threaded spine pole in step 103 includes: in the upper structure, traversing each ring one by one to calculate the polar angles of all nodes at the same height in each ring, and when the polar angle difference of the corresponding nodes on adjacent rings is equal to a preset value, determining the position of the cable-threaded spine pole.

[0095] In a feasible implementation, the specific method for determining the position of the strut in step 103 includes: at the same height, by using the solid geometry theorem, calculating the foot of the perpendicular of the midpoint of the line connecting the upper nodes of adjacent struts in the z plane of the lower node, wherein the z plane is the horizontal plane where the z coordinate of the lower node is located; then taking the foot of the perpendicular and the line connecting the upper nodes of a certain adjacent strut as the positive strut; changing the polar diameter of the positive strut and increasing the offset strut distance to determine the polar diameter of the lower node of the strut; connecting the upper node of the strut and the lower node of the strut to form a strut.

[0096] In a feasible implementation manner, the specific method for determining the positions of the hoop and the oblique cable in step 103 includes: at the same height, using the solid geometry theorem to convert the line between the midpoints of the line connecting the lower nodes of adjacent struts as the hoop;

[0097] At the same height, traverse all the lower nodes of the struts ring by ring and take them as an endpoint; among all the nodes, select the node with the same polar angle as the endpoint as the other endpoint; connect these two endpoints to form an oblique cable.

[0098] In a feasible implementation manner, the design requirement data in step 105 include design requirement values ​​of span, rise-to-span ratio, number of radial divisions, number of circumferential divisions, angle between the inclined cable and the horizontal plane, etc.

[0099] In a feasible implementation manner, in step 8, the method for calculating the internal force of the prestressed spinal rod in the central ring is:

[0100] In the cable dome structure visualization inspection model, select each upper node in the central ring; search for several rods connected to the node; search for the other end number ID of the rod to obtain the corresponding rod number ID; search for the rectangular coordinates of the rod based on the rod number ID; calculate the direction vector (x, y and z) of each rod; construct the force balance equations in the x, y and z directions, and solve them simultaneously to obtain the internal force of the prestressed ridge rod. The specific formula can be:

[0101] Number of struts * strut internal force * x-direction component + number of vertebrae * vertebrae internal force * x-direction component = 0;

[0102] Number of struts * strut internal force * y-direction component + number of vertebrae * vertebrae internal force * y-direction component = 0;

[0103] Number of struts*strut internal force*z-direction component+number of ridgepole*ridgepole internal force*z-direction component+equivalent load of the node=0.

[0104] In a feasible implementation, in step 8, the method for calculating the internal force of the prestressed vertebral rod and the internal force of the strut in the non-center ring is:

[0105] In the cable dome structure visualization inspection model, select each upper node in the non-central ring; search for several rods connected to the node; search for the other end number ID of the rod to obtain the corresponding rod number ID; search for the rectangular coordinates of the rod based on the rod number ID; calculate the direction vector (x, y and z) of each rod; construct the force balance equations in the x, y and z directions, and solve them simultaneously to obtain the internal force of the prestressed ridge rod. The specific formula can be:

[0106] Number of vertebrae*internal force of vertebrae*x-direction component+number of oblique cables*internal force of oblique cables*x-direction component+number of struts*internal force of struts*x-direction component=0;

[0107] Number of vertebrae*internal force of vertebrae*component in y direction+number of oblique cables*internal force of oblique cables*component in y direction+number of struts*internal force of struts*component in y direction=0;

[0108] Number of ridge poles*ridge pole internal force*z-direction component + number of oblique cables*oblique cable internal force*z-direction component + number of struts*strut internal force*z-direction component + equivalent load at this node = 0.

[0109] In a feasible implementation, the method for calculating the internal force of the oblique cable in the central ring is:

[0110] In the cable dome structure visualization inspection model, select each lower node in the central ring; search for several rods connected to the node; search for the other end number ID of the rod to obtain the corresponding rod number ID; search for the rectangular coordinates of the rod based on the rod number ID; calculate the direction vector (x, y and z) of each rod; construct the force balance equations in the x, y and z directions, and solve them simultaneously to obtain the internal force of the inclined cable. The specific formula can be:

[0111] Number of struts * internal force of struts * x-direction component + number of inclined cables * internal force of inclined cables * x-direction component = 0;

[0112] Number of struts * internal force of struts * component in y direction + number of inclined cables * internal force of inclined cables * component in y direction = 0;

[0113] Number of struts*internal force of struts*component in z direction+number of inclined cables*internal force of inclined cables*component in z direction+equivalent load of the node=0.

[0114] In a feasible implementation, the method for calculating the internal forces of the oblique cables and the ring cables in the non-central ring is:

[0115] In the cable dome structure visualization inspection model, select each lower node in the non-central ring; search for several rods connected to the node; search for the other end number ID of the rod to obtain the corresponding rod number ID; search for the rectangular coordinates of the rod based on the rod number ID; calculate the direction vector (x, y and z) of each rod; construct the force balance equations in the x, y and z directions, and solve them simultaneously to obtain the internal force of the inclined cable. The specific formula can be:

[0116] Number of struts * internal force of struts * x-direction component + number of oblique cables * internal force of oblique cables * x-direction component + number of hoop cables * internal force of hoop cables * x-direction component = 0;

[0117] Number of struts * internal force of struts * component in y direction + number of oblique cables * internal force of oblique cables * component in y direction + number of hoop cables * internal force of hoop cables * component in y direction = 0;

[0118] Number of struts*internal force of struts*component in z direction+number of inclined cables*internal force of inclined cables*component in z direction+number of hoop cables*internal force of hoop cables*component in z direction+equivalent load of this node=0.

[0119] In a feasible implementation manner, the specific training process of the Q network in step 12 includes:

[0120] Step 121, setting the initial value of the cable prestress as the initial state;

[0121] Step 122: According to the current status , select an action using the ε-greedy strategy (ε represents the probability of randomly selecting an action) or the policy gradient method (Adjust cable prestress), the specific formula can be:

[0122] ;

[0123] in, Indicates action The vector representation of ; Indicates the Q value of the Q network; state The vector representation of is: ;

[0124] in, Indicates the current prestress value of the cable; Indicates the total number of cables;

[0125] And define Indicates the range of change of prestress value, used to limit the action space;

[0126] Step 123: Execute the action , get a new state and rewards ; The reward , which is used to measure the closeness of the current prestress value to the target prestress value after each prestress adjustment. The specific formula can be:

[0127] ;

[0128] in, Indicates the current prestress value; Indicates the target prestress value;

[0129] Step 124: Update the Q value through the Q learning algorithm. The specific formula can be:

[0130] ;

[0131] in, represents the learning rate; represents the discount factor; Indicates a new action;

[0132] Step 125: After taking the new state as the current state, iterate step 122 until the termination condition is reached;

[0133] Step 126: After the Q network is trained, the Q network is tested using a validation set to evaluate the computing performance of the Q network under different initial states. The parameters of the Q network may also be adjusted based on the test results.

[0134] In a feasible implementation manner, the termination condition of step 125 includes reaching a target prestress value or a maximum number of iterations.

[0135] In the second aspect, based on the same inventive concept, the present application also provides a cable dome structure prestress calculation system, including a data receiving module, a data processing module and a result generating module;

[0136] The data receiving module is used to receive basic structural data, known mechanical parameters and nodes that need to be calculated of the cable dome structure;

[0137] The data processing module performs calculations using the cable dome structure prestress calculation method to obtain a final design value of the cable prestress;

[0138] The result generation module is used to send out the final design value of the cable prestressing.

[0139] In the third aspect, based on the same inventive concept, the present application also provides a cable dome structure prestress calculation device, including a processor, a memory and a bus, the memory stores instructions and data read by the processor, the processor is used to call the instructions and data in the memory to execute the above-mentioned cable dome structure prestress calculation method, and the bus connects the functional components for transmitting information.

[0140] By adopting the above technical solution, the present invention has the following beneficial effects:

[0141] The present invention provides a cable dome structure prestress calculation method, system and device, which can automatically, quickly and accurately calculate the final design value of the prestress of the cable dome structure, saving a lot of manpower and time, reducing calculation errors, and improving design efficiency, thereby helping to speed up project progress, reduce project costs, and improve project quality and safety. BRIEF DESCRIPTION OF THE DRAWINGS

[0142] In order to more clearly illustrate the specific implementation methods of the present invention or the technical solutions in the prior art, the drawings required for use in the specific implementation methods or the description of the prior art will be briefly introduced below. Obviously, the drawings in the following description are some implementation methods of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0143] Figure 1 A schematic diagram of a cable dome structure provided in the background technology;

[0144] Figure 2 A flow chart of a method for calculating prestress of a cable dome structure provided by an embodiment of the present invention;

[0145] Figure 3A prestress distribution diagram provided by an embodiment of the present invention;

[0146] Figure 4 for Figure 2 Specific flow chart of step 1 in FIG.

[0147] Figure 5 A top view of a visual inspection model of a cable dome structure provided by an embodiment of the present invention;

[0148] Figure 6 A stereogram of a visual inspection model of a cable dome structure provided by an embodiment of the present invention;

[0149] Figure 7 for Figure 6 A top view of

[0150] Figure 8 A diagram illustrating calculation of cable dome structural parameters provided by an embodiment of the present invention;

[0151] Fig. 9 A diagram illustrating the definition of a cable dome structural component provided by an embodiment of the present invention;

[0152] Fig.10 A z-plane top view illustrating calculation of the offset distance provided in an embodiment of the present invention;

[0153] Fig.11 A diagram of a prestress calculation system for a cable dome structure provided by an embodiment of the present invention;

[0154] Reference numerals:

[0155] 1-Center support rod; 2-Annular support rod; 3-Cable-through ridge rod; 4-Annular cable; 5-Cable support node; 6-Upper end point; 7-Oblique cable. DETAILED DESCRIPTION

[0156] The technical solution of the present invention will be described clearly and completely below 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 the present invention.

[0157] In the description of the present invention, it should be noted that the terms "center", "upper", "lower", "left", "right", "vertical", "horizontal", "inner", "outer", etc., indicating the orientation or positional relationship, are based on the orientation or positional relationship shown in the drawings, and are only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying 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 limiting the present invention. In addition, the terms "first", "second", and "third" are used for descriptive purposes only, and cannot be understood as indicating or implying relative importance.

[0158] In the description of the present invention, it should be noted that, 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 it can be indirectly connected through an intermediate medium, or it can be the internal communication of two components. For ordinary technicians in this field, the specific meanings of the above terms in the present invention can be understood according to specific circumstances.

[0159] The present invention is further explained below in conjunction with specific implementation modes.

[0160] It should also be noted that the following specific embodiments or specific implementations are a series of optimized settings listed in the present invention to further explain the specific content of the invention, and these settings can be used in combination or in association with each other.

[0161] Embodiment 1:

[0162] like Figure 2 As shown, the present embodiment provides a method for calculating the prestress of a cable dome structure, comprising the following steps:

[0163] Step 1, receiving the basic structural data of the cable dome structure; grouping and numbering each node according to the structure in which it is located; based on the node number ID, generating the point information node_dict of each node and the adjacency matrix between each node, and constructing a visual inspection model of the cable dome structure; the basic structural data includes: structural topological information, span, arrow span ratio, annular division number and radial division number; the structural topological information refers to the type of the cable dome structure; the annular division number refers to the number of equal parts into which the entire circumference is divided along the annular direction of the dome; the radial division number refers to the number of rings into which the entire dome is divided from the center of the dome outward;

[0164] Step 2, receiving known mechanical parameters of the cable dome structure; the known mechanical parameters include the prestress value of the central strut, the node load and the prestress value of the prestressed ridge pole, etc.;

[0165] Step 3: Extract the nodes that need to be calculated from the point information node_dict and generate a dictionary of nodes that need to be calculated; traverse each node in the point information node_dict, extract the first node ID of each ring structure in each layer structure, and store it in the first node dictionary first_nodes;

[0166] Step 4: Extract the first nodes one by one from the first node dictionary first_nodes, and check whether the node belongs to the dictionary of nodes to be calculated: if so, execute step 5; if not, execute step 4;

[0167] Step 5: From the adjacency matrix, find other nodes connected to the current first node as slave nodes, and add the slave node ID to the sublist of the connected node dictionary point_dic_j_ids;

[0168] Step 6, obtain the coordinates of the node and the ring structure where it is located through the node ID; based on the connected node dictionary point_dic_j_ids, traverse and calculate the direction vector, component length and height difference between the current head node and each slave node; store the direction vector in the direction vector dictionary direction_vectors;

[0169] Step 7: Based on the point information node_dict and the connected node dictionary point_dic_j_ids, obtain the layer structure and ring structure of the node and make a judgment: if the current first node is in the upper structure, execute step 8; if the current first node is in the lower structure, execute step 9;

[0170] Step 8, if the node is located in the central ring structure, search the corresponding direction vector from the direction vector dictionary direction_vectors, and calculate the prestressed ridge rod internal force based on the known node load; if the node is located in the non-central ring structure, search the corresponding direction vector from the direction vector dictionary direction_vectors, and calculate the prestressed ridge rod internal force and the strut internal force based on the known node load; execute step 10; the internal force represents the force generated inside each component of the cable dome structure after the cable dome structure is balanced;

[0171] Step 9: If the node is located in the central ring structure, search the corresponding direction vector from the direction vector dictionary direction_vectors, and calculate the internal force of the inclined cable based on the known node load; if the node is located in the non-central ring structure, search the corresponding direction vector from the direction vector dictionary direction_vectors, and calculate the internal force of the inclined cable and the ring cable based on the known node load; execute step 10;

[0172] Step 10: Update the prestress dictionary based on the internal force values ; Iterate step 4 until all the first nodes are traversed;

[0173] Step 11: Traverse the prestress dictionary In each key-value pair, the key in the key-value pair represents the component type, and the value in the key-value pair represents the prestress value of the component at different rings, that is, the internal force of the tensioning cable required in the tensioning construction; based on the key-value pair, the prestress distribution diagram is drawn by the plt.plot() function, and the data label is marked on each data point of the prestress distribution diagram by the plt.text() function, and then displayed by the plt.show() function, such as Figure 3 As shown;

[0174] Step 12, using the Q network, calculate the initial value and design value of the cable prestress in the equilibrium state; the cable prestress includes inclined cable prestress and ring cable prestress, which is specifically determined by the design documents and the construction plan; using finite element tools, such as calling ANSYS finite element software through the pyansys library in Python, to construct a geometric model of the cable dome structure; meshing the geometric model to obtain a finite element model; calculating physical quantities such as stress and deformation of the cable dome structure through finite element analysis and displaying them; adjusting the geometric model parameters according to external feedback until the design value of the cable prestress that meets the preset conditions is obtained as the final design value.

[0175] Through the above steps, the final design value of the prestress of the cable dome structure can be automatically calculated, which saves a lot of manpower and time, reduces calculation errors, and improves design efficiency and engineering quality.

[0176] Furthermore, if Figure 4 As shown, the step 1 specifically includes:

[0177] Step 101, in the Python environment, based on the input parameters of the cable dome structure, the design parameters of the cable dome structure such as the spherical radius and the unit inclination angle are obtained through the proportional relationship of similar triangles and the spherical coordinate system calculation; based on the rectangular coordinate system, the design parameters are converted into the position information dictionary node_dict of the nodes in the cable dome structure for subsequent programming indexing, and stored in the node list; based on the node list, with reference to the cable dome structure, the component information of the cable-through ridge pole, the strut, the ring cable and the oblique cable is added to construct a three-dimensional list;

[0178] Step 102: Generate the components such as the cable-through ridge pole, the strut pole, the ring rope and the oblique rope of each ring in sequence based on the three-dimensional list; then construct a two-dimensional adjacency matrix of the connection relationship of each node in each ring; and then splice the two-dimensional adjacency matrix into a three-dimensional adjacency matrix according to the order of the rings;

[0179] Step 103, determining the positions of components such as cable-through ridge poles, struts, ring cables, and oblique cables layer by layer and ring by ring; performing three-dimensional visualization processing on the cable dome structure through the matplotlib library to obtain a visualization inspection model of the cable dome structure, thereby facilitating homomorphic adjustment and preview;

[0180] The cable-through vertebral rod is a component in which the vertebral cord passes through the vertebral rod along the axial direction. The component applies prestress to make the vertebral cord play a tensile role, while the vertebral rod mainly bears pressure. The component fully utilizes the respective advantages of steel cables and steel pipes to enhance the overall stability and bearing capacity of the structure.

[0181] The component also includes an outer sleeve and connecting parts:

[0182] The outer sleeve can be specifically made by welding a seamless steel pipe and a partition; a circular hole is provided in the middle of the partition for positioning and tensioning the ridge cable;

[0183] The connecting piece is used to connect the two sections of the cable-threading spinal rod end to end along the axial direction into one piece;

[0184] Step 104: Figure 5-7 As shown, the corresponding relationship between the node number, the node coordinate, the three-dimensional adjacency matrix and the cable dome structure visual inspection model is established and stored in the database, so that the node coordinate and the adjacency matrix can be queried and updated according to the node number later;

[0185] Step 105: Based on the design requirement data, manually determine whether the cable dome structure visualization inspection model meets the requirements: if so, output information about the specified nodes and their connection relationships; if not, output unqualified information and receive adjustment values ​​of input parameters, and update the cable dome structure visualization inspection model until the requirements are met.

[0186] Through the above steps, the input parameters of the cable dome structure can be directly converted into an editable visual model in the Python environment, eliminating the data conversion steps between the modeling environment and the architectural software, which can greatly improve the design efficiency and quality of the cable dome structure.

[0187] Furthermore, the input parameters in step 101 include span, rise-to-span ratio, radial division number, circumferential division number, angle between the inclined cable and the horizontal plane, etc.; the radial division number refers to the number of rings into which the entire dome is divided from the center of the dome outward; the circumferential division number refers to the number of equal parts into which the entire circumference is divided along the circumferential direction of the dome.

[0188] Furthermore, the nodes are divided into upper structure nodes and lower structure nodes, the point type mark of the upper structure node is 0, and the point type mark of the lower structure node is 1; the location information dictionary node_dict includes the node number and the node coordinates.

[0189] Furthermore, the three-dimensional list specifically includes:

[0190] The first dimension represents the type of rod, including: cable-through ridge rod, strut, ring cable and oblique cable;

[0191] The second dimension represents the number of radial divisions, including the ordinal number of the ring;

[0192] The third dimension represents the connection relationship between nodes. A value of 0 indicates no connection, and a value of 1 indicates connection.

[0193] Furthermore, the step 101 specifically includes the following steps:

[0194] Step 1011: Figure 8 As shown, based on the spherical coordinate system, the spherical radius of the cable dome structure is calculated and unit inclination , the specific formula can be:

[0195] ;

[0196] ;

[0197] in, It represents the coverage radius of the cable dome structure. The specific calculation formula can be:

[0198] ;

[0199] It represents the sagittal height of the cable dome structure. The specific calculation formula can be:

[0200] ;

[0201] Step 1012: Calculate the vertical heights of nodes in adjacent rings based on the proportional relationship of similar triangles. The specific formula may be:

[0202] ;

[0203] ;

[0204] in, Indicates The vertical height of the nodes in the ring; Indicates The vertical height of the nodes in the ring; Indicates The specific formula for the proportional coefficient of the ring is:

[0205] ;

[0206] Indicates The specific formula for the proportional coefficient of the ring is:

[0207] ;

[0208] Step 1013: Calculate the radial coordinates of the nodes in the adjacent rings. The specific formula may be:

[0209] ;

[0210] ;

[0211] in, Indicates The radial coordinates of the nodes in the ring; Indicates The radial coordinates of the nodes in the ring;

[0212] Calculate the radial difference and store it in the radial difference list The specific formula can be:

[0213] ;

[0214] in, Indicates The nodes in the ring and The radial difference between the nodes in the ring;

[0215] Step 1014, calculate the height of the support pole, the specific formula can be:

[0216] ;

[0217] in, Indicates The height of the poles in the ring; Indicates The height of the brace in the ring in the superstructure can be expressed as:

[0218] ;

[0219] in, Indicates The angle between the inclined cable and the horizontal plane in each ring;

[0220] Indicates The height of the brace in the lower structure in a ring can be expressed as:

[0221] ;

[0222] Step 1015: Calculate the node coordinates based on the rectangular coordinate system ( , ) and the corresponding coordinates of the lower nodes of the struts ( , ), the specific formula can be:

[0223] when When it is an even number:

[0224] ;

[0225] in, Indicates the x-coordinate of the node; Indicates the y-coordinate of the node; Indicates the current ring equal parts;

[0226] And there are:

[0227] ;

[0228] when When it is an odd number:

[0229] ;

[0230] And there are:

[0231] ;

[0232] Define the mapping relationship between node numbers and rectangular coordinates. The specific expression can be:

[0233] ;

[0234] in, Indicates the node number. The topmost node of the cable dome structure can be taken as the first node, and the numbers can be arranged in sequence in a counterclockwise spiral outward direction to the last node.

[0235] Furthermore, the step 102 specifically includes the following steps:

[0236] Step 1021, respectively define the adjacency matrices of the rod types, namely: the through-cable ridge rod adjacency matrix adj_matrix_RB, the ring-cable adjacency matrix adj_matrix_CC, the oblique-cable adjacency matrix adj_matrix_DC and the strut adjacency matrix adj_matrix_S;

[0237] Each adjacency matrix is ​​a three-dimensional array:

[0238] The first dimension represents the ordinal number of the ring;

[0239] The second and third dimensions form a square matrix, which represents the connection relationship between all nodes in the current ring. The specific setting rules are: when the element value is 0, it means no connection; when the element value is 1, it means there is a connection;

[0240] Initialize the three-dimensional array and set all elements to 0;

[0241] In this way, complex structural data can be clearly organized and managed through the three-dimensional array, so that the connection status of each ring can be independently stored and called for subsequent analysis and operation;

[0242] Step 1022, traverse and select any two nodes (point1 and point2) to form a node pair, find the node numbers of the two nodes based on the position information dictionary node_dict and determine the connection relationship, if there is a connection relationship, generate a unique member number ID according to the sequence number of the ring, the node number and the member type, so as to distinguish the members in the subsequent sequence; set the connection relationship and member number ID of the two nodes in the adjacency matrix of the ring where the two nodes are located;

[0243] Step 1023: Store the adjacency matrix of the member type in the location information dictionary node_dict for subsequent access and calling.

[0244] Furthermore, the specific method of determining the connection relationship in step 1022 includes:

[0245] Step 10221, according to the ring where the node is located and the point type mark, define the node classification rules: an upper structure node in the same ring is connected to the nearest lower structure node to form a strut; in adjacent rings, the two nearest upper structure nodes are connected to form a cable-through ridge; in adjacent rings, a lower structure node in the previous ring is connected to the nearest upper structure node in the next ring to form an oblique cable; in the same ring, two adjacent lower structure nodes are connected to form a ring cable, such as Fig. 9 As shown;

[0246] Step 10222: In the location information dictionary node_dict, traverse and construct node pairs; based on the point type mark of the node and the equal part in which it is located, determine the connection relationship of the node pair through the node classification rule; according to the classification, mark the element value of the connection relationship in the corresponding adjacency matrix.

[0247] Furthermore, the specific method for determining the position of the cable-threaded spine pole in step 103 includes: in the upper structure, traversing each ring one by one to calculate the polar angles of all nodes at the same height in each ring, and when the polar angle difference of corresponding nodes on adjacent rings is equal to a preset value, determining the position of the cable-threaded spine pole.

[0248] Furthermore, if Fig.10 As shown, the specific method for determining the position of the strut in step 103 includes: at the same height, by solid geometry theorem, calculating the foot A of the perpendicular of the midpoint of the line connecting the upper nodes of adjacent struts in the z plane of the lower node, wherein the z plane is the horizontal plane where the z coordinate of the lower node is located; then taking the line AE between the foot of the perpendicular and a certain upper node of an adjacent strut as the positive strut; changing the polar diameter of the positive strut and increasing the offset distance AF to determine the polar diameter OF of the lower node of the strut; connecting the upper node of the strut and the lower node of the strut to form a strut FE.

[0249] Furthermore, the specific method for determining the positions of the hoop and the oblique cable in step 103 includes: at the same height, by converting the solid geometry theorem, taking the line between the midpoints of the line connecting the lower nodes of adjacent struts as the hoop;

[0250] At the same height, traverse all the lower nodes of the struts ring by ring and take them as an endpoint; among all the nodes, select the node with the same polar angle as the endpoint as the other endpoint; connect these two endpoints to form an oblique cable.

[0251] Furthermore, the design requirement data in step 105 include design requirement values ​​of span, rise-to-span ratio, number of radial divisions, number of circumferential divisions, angle between the inclined cable and the horizontal plane, etc.

[0252] Furthermore, in step 8, the method for calculating the internal force of the prestressed vertebral rod in the central ring is:

[0253] In the cable dome structure visualization inspection model, select each upper node in the central ring; search for several rods connected to the node; search for the other end number ID of the rod to obtain the corresponding rod number ID; search for the rectangular coordinates of the rod based on the rod number ID; calculate the direction vector (x, y and z) of each rod; construct the force balance equations in the x, y and z directions, and solve them simultaneously to obtain the internal force of the prestressed ridge rod. The specific formula can be:

[0254] Number of struts * strut internal force * x-direction component + number of vertebrae * vertebrae internal force * x-direction component = 0;

[0255] Number of struts * strut internal force * y-direction component + number of vertebrae * vertebrae internal force * y-direction component = 0;

[0256] Number of struts*strut internal force*z-direction component+number of ridgepole*ridgepole internal force*z-direction component+equivalent load of the node=0.

[0257] Furthermore, in step 8, the method for calculating the internal force of the prestressed vertebral rod and the internal force of the strut in the non-center ring is:

[0258] In the cable dome structure visualization inspection model, select each upper node in the non-central ring; search for several rods connected to the node; search for the other end number ID of the rod to obtain the corresponding rod number ID; search for the rectangular coordinates of the rod based on the rod number ID; calculate the direction vector (x, y and z) of each rod; construct the force balance equations in the x, y and z directions, and solve them simultaneously to obtain the internal force of the prestressed ridge rod. The specific formula can be:

[0259] Number of vertebrae*internal force of vertebrae*x-direction component+number of oblique cables*internal force of oblique cables*x-direction component+number of struts*internal force of struts*x-direction component=0;

[0260] Number of vertebrae*internal force of vertebrae*component in y direction+number of oblique cables*internal force of oblique cables*component in y direction+number of struts*internal force of struts*component in y direction=0;

[0261] Number of ridge poles*ridge pole internal force*z-direction component + number of oblique cables*oblique cable internal force*z-direction component + number of struts*strut internal force*z-direction component + equivalent load at this node = 0.

[0262] Furthermore, in step 9, the method for calculating the internal force of the oblique cable in the central ring is:

[0263] In the cable dome structure visualization inspection model, select each lower node in the central ring; search for several rods connected to the node; search for the other end number ID of the rod to obtain the corresponding rod number ID; search for the rectangular coordinates of the rod based on the rod number ID; calculate the direction vector (x, y and z) of each rod; construct the force balance equations in the x, y and z directions, and solve them simultaneously to obtain the internal force of the inclined cable. The specific formula can be:

[0264] Number of struts * internal force of struts * x-direction component + number of inclined cables * internal force of inclined cables * x-direction component = 0;

[0265] Number of struts * internal force of struts * component in y direction + number of inclined cables * internal force of inclined cables * component in y direction = 0;

[0266] Number of struts*internal force of struts*component in z direction+number of inclined cables*internal force of inclined cables*component in z direction+equivalent load of the node=0.

[0267] Furthermore, in step 9, the method for calculating the internal force of the oblique cable and the internal force of the ring cable in the non-center ring is:

[0268] In the cable dome structure visualization inspection model, select each lower node in the non-central ring; search for several rods connected to the node; search for the other end number ID of the rod to obtain the corresponding rod number ID; search for the rectangular coordinates of the rod based on the rod number ID; calculate the direction vector (x, y and z) of each rod; construct the force balance equations in the x, y and z directions, and solve them simultaneously to obtain the internal force of the inclined cable. The specific formula can be:

[0269] Number of struts * internal force of struts * x-direction component + number of oblique cables * internal force of oblique cables * x-direction component + number of hoop cables * internal force of hoop cables * x-direction component = 0;

[0270] Number of struts * internal force of struts * component in y direction + number of oblique cables * internal force of oblique cables * component in y direction + number of hoop cables * internal force of hoop cables * component in y direction = 0;

[0271] Number of struts*internal force of struts*component in z direction+number of inclined cables*internal force of inclined cables*component in z direction+number of hoop cables*internal force of hoop cables*component in z direction+equivalent load of this node=0.

[0272] Furthermore, the specific training process of the Q network in step 12 includes:

[0273] Step 121, setting the initial value of the cable prestress as the initial state;

[0274] Step 122: According to the current status , select an action using the ε-greedy strategy (ε represents the probability of randomly selecting an action) or the policy gradient method (Adjust cable prestress), the specific formula can be:

[0275] ;

[0276] in, Indicates action The vector representation of ; Indicates the Q value of the Q network; state The vector representation of is: ;

[0277] in, Indicates the current prestress value of the cable; Indicates the total number of cables;

[0278] And define Indicates the range of change of prestress value, used to limit the action space;

[0279] Step 123: Execute the action , get a new state and rewards ; The reward , which is used to measure the closeness of the current prestress value to the target prestress value after each prestress adjustment. The specific formula can be:

[0280] ;

[0281] in, Indicates the current prestress value; Indicates the target prestress value;

[0282] Step 124: Update the Q value through the Q learning algorithm. The specific formula can be:

[0283] ;

[0284] in, represents the learning rate; represents the discount factor; Indicates a new action;

[0285] Step 125: After taking the new state as the current state, iterate step 122 until the termination condition is reached;

[0286] Step 126: After the Q network is trained, the Q network is tested using a validation set to evaluate the computing performance of the Q network under different initial states. The parameters of the Q network may also be adjusted based on the test results.

[0287] Furthermore, the termination condition of step 125 includes reaching a target prestress value or a maximum number of iterations.

[0288] Embodiment 2:

[0289] like Fig.11 As shown, this embodiment provides a cable dome structure prestress calculation system, including a data receiving module, a data processing module and a result generating module;

[0290] The data receiving module is used to receive basic structural data, known mechanical parameters and nodes that need to be calculated of the cable dome structure;

[0291] The data processing module performs calculations using the cable dome structure prestress calculation method to obtain a final design value of the cable prestress;

[0292] The result generation module is used to send out the final design value of the cable prestressing.

[0293] Embodiment three:

[0294] The present embodiment provides a cable dome structure prestress calculation device, comprising a processor, a memory and a bus, wherein the memory stores instructions and data read by the processor, the processor is used to call the instructions and data in the memory to execute the cable dome structure prestress calculation method, and the bus connects the functional components to transmit information.

[0295] In another embodiment, the present solution can be implemented by a device, which may include a corresponding module for performing each or several steps in the above-mentioned embodiments. The module may be one or more hardware modules specially configured to perform the corresponding steps, or implemented by a processor configured to perform the corresponding steps, or stored in a computer-readable medium for implementation by a processor, or implemented by some combination.

[0296] The processor performs the various methods and processes described above. For example, the method implementation in the present solution can be implemented as a software program, which is tangibly contained in a machine-readable medium, such as a memory. In some embodiments, part or all of the software program can be loaded and / or installed via a memory and / or a communication interface. When the software program is loaded into the memory and executed by the processor, one or more steps in the method described above can be executed. Alternatively, in other embodiments, the processor can be configured to perform one of the above methods by any other appropriate means (for example, by means of firmware).

[0297] The device can be implemented using a bus architecture. The bus architecture can include any number of interconnecting buses and bridges, depending on the specific application and overall design constraints of the hardware. The bus connects various circuits including one or more processors, memories, and / or hardware modules together. The bus can also connect various other circuits such as peripherals, voltage regulators, power management circuits, external antennas, etc.

[0298] The bus can be an Industry Standard Architecture (ISA) bus, a Peripheral Component Interconnect (PCI) bus, or an Extended Industry Standard Component (EISA) bus, etc. The bus can be divided into an address bus, a data bus, a control bus, etc.

[0299] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or replace some or all of the technical features therein with equivalents. However, these modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for calculating the prestress of a cable dome structure, characterized in that: include: Step 1, receiving the structural basic data of the cable dome structure; Group and number each node according to its structure; Based on the node number, the node position information node_dict of each node and the adjacency matrix between each node are generated to build a visual inspection model of the cable dome structure; Step 2, receiving known mechanical parameters of the cable dome structure; Step 3: Extract the nodes that need to be calculated from the point information node_dict and generate a dictionary of nodes that need to be calculated; traverse each node in the point information node_dict, extract the first node number of each ring structure in each layer structure, and store it in the first node dictionary first_nodes; Step 4: Extract the first nodes one by one from the first node dictionary first_nodes, and check whether the node belongs to the dictionary of nodes to be calculated: if so, execute step 5; if not, execute step 4; Step 5: From the adjacency matrix, find other nodes connected to the current first node as slave nodes, and add the slave node number to the sublist of the connected node dictionary point_dic_j_ids; Step 6: Obtain the coordinates of the node and the ring structure where it is located through the node number; based on the connected node dictionary point_dic_j_ids, traverse and calculate the direction vector, component length and height difference between the current head node and each slave node; Store the direction vector in the direction vector dictionary direction_vectors; Step 7: Based on the point information node_dict and the connected node dictionary point_dic_j_ids, obtain the layer structure and ring structure of the node and make a judgment: if the current first node is in the upper structure, execute step 8; if the current first node is in the lower structure, execute step 9; Step 8: If the node is located in the central ring structure, search the corresponding direction vector from the direction vector dictionary direction_vectors, and calculate the prestressed ridge rod internal force based on the known node load; if the node is located in the non-central ring structure, search the corresponding direction vector from the direction vector dictionary direction_vectors, and calculate the prestressed ridge rod internal force and the brace rod internal force based on the known node load; execute step 10; Step 9: If the node is located in the central ring structure, search the corresponding direction vector from the direction vector dictionary direction_vectors, and calculate the internal force of the inclined cable based on the known node load; If the node is located in a non-central ring structure, the corresponding direction vector is searched from the direction vector dictionary direction_vectors, and the internal force of the inclined cable and the internal force of the ring cable are calculated based on the known node load; step 10 is executed; Step 10: Update the prestress dictionary based on the internal force values ; Iterate step 4 until all the first nodes are traversed; Step 11: Traverse the prestress dictionary For each key-value pair in the graph, a prestress distribution diagram is drawn and a data label is annotated on each data point of the prestress distribution diagram for display; Step 12, calculating the initial value and design value of the cable prestress in the equilibrium state through the Q network; Use finite element tools to build a geometric model of the cable dome structure; Meshing the geometric model to obtain a finite element model; Through finite element analysis, the stress and deformation of the cable dome structure are calculated and displayed; According to external feedback, the geometric model parameters are adjusted until the design value of the cable prestress that meets the preset conditions is obtained as the final design value.

2. The method according to claim 1, characterized in that The basic structural data include: structural topological information, span, span-to-width ratio, annular division number and radial division number.

3. The method according to claim 1, characterized in that The known mechanical parameters include the prestress value of the central support rod, the node load and the prestress value of the prestressed ridge rod.

4. The method according to claim 1, characterized in that: In step 8, the method for calculating the internal force of the prestressed vertebral rod in the central ring is: In the cable dome structure visualization inspection model, select each upper node in the central ring; search for several rods connected to the node; search for the other end number of the rod to obtain the corresponding rod number; search for the rectangular coordinates of the rod based on the rod number; calculate the direction vector of each rod; construct the force balance equations in the x-, y- and z-directions, and solve them simultaneously to obtain the internal force of the prestressed ridge rod. The specific formula is: Number of struts * strut internal force * x-direction component + number of vertebrae * vertebrae internal force * x-direction component = 0; Number of struts * strut internal force * y-direction component + number of vertebrae * vertebrae internal force * y-direction component = 0; Number of struts*strut internal force*z-direction component+number of ridgepole*ridgepole internal force*z-direction component+equivalent load of the node=0.

5. The method according to claim 1, characterized in that In step 8, the method for calculating the internal force of the prestressed ridge rod and the internal force of the strut in the non-central ring is: In the cable dome structure visualization inspection model, select each upper node in the non-central ring; search for several rods connected to the node; search for the other end number of the rod to obtain the corresponding rod number; search for the rectangular coordinates of the rod based on the rod number; calculate the direction vector of each rod; construct the force balance equations in the x-, y- and z-directions, and solve them simultaneously to obtain the internal force of the prestressed ridge rod. The specific formula is: Number of vertebrae*internal force of vertebrae*x-direction component+number of oblique cables*internal force of oblique cables*x-direction component+number of struts*internal force of struts*x-direction component=0; Number of vertebrae*internal force of vertebrae*component in y direction+number of oblique cables*internal force of oblique cables*component in y direction+number of struts*internal force of struts*component in y direction=0; Number of ridge poles*ridge pole internal force*z-direction component + number of oblique cables*oblique cable internal force*z-direction component + number of struts*strut internal force*z-direction component + equivalent load at this node = 0.

6. The method according to claim 1, characterized in that In step 9, the method for calculating the internal force of the oblique cable in the central ring is: In the cable dome structure visualization inspection model, select each lower node in the central ring; search for several rods connected to the node; search for the other end number of the rod to obtain the corresponding rod number; search for the rectangular coordinates of the rod based on the rod number; calculate the direction vector of each rod; construct the force balance equations in the x-, y- and z-directions, and solve them simultaneously to obtain the internal force of the inclined cable. The specific formula is: Number of struts * internal force of struts * x-direction component + number of inclined cables * internal force of inclined cables * x-direction component = 0; Number of struts * internal force of struts * component in y direction + number of inclined cables * internal force of inclined cables * component in y direction = 0; Number of struts*internal force of struts*component in z direction+number of inclined cables*internal force of inclined cables*component in z direction+equivalent load of the node=0.

7. The method according to claim 1, characterized in that In step 9, the method for calculating the internal forces of the oblique cables and the hoop cables in the non-central ring is: In the cable dome structure visualization inspection model, select each lower node in the non-central ring; search for several rods connected to the node; search for the other end number of the rod to obtain the corresponding rod number; search for the rectangular coordinates of the rod based on the rod number; calculate the direction vector of each rod; construct the force balance equations in the x-, y- and z-directions, and solve them simultaneously to obtain the internal force of the inclined cable. The specific formula is: Number of struts * internal force of struts * x-direction component + number of oblique cables * internal force of oblique cables * x-direction component + number of hoop cables * internal force of hoop cables * x-direction component = 0; Number of struts * internal force of struts * component in y direction + number of oblique cables * internal force of oblique cables * component in y direction + number of hoop cables * internal force of hoop cables * component in y direction = 0; Number of struts*internal force of struts*component in z direction+number of inclined cables*internal force of inclined cables*component in z direction+number of hoop cables*internal force of hoop cables*component in z direction+equivalent load of this node=0.

8. The method according to claim 1, characterized in that: The key in the key-value pair represents the component type, and the value in the key-value pair represents the prestress value of the component at different rings.

9. A cable dome structure prestress calculation system, characterized in that: It includes a data receiving module, a data processing module and a result generating module; The data receiving module is used to receive basic structural data, known mechanical parameters and nodes that need to be calculated of the cable dome structure; The data processing module calculates the final design value of the cable prestress by the method according to any one of claims 1 to 8; The result generation module is used to send out the final design value of the cable prestressing.

10. A cable dome structure prestress calculation device, characterized in that: It includes a processor, a memory and a bus, wherein the memory stores instructions and data read by the processor, the processor is used to call the instructions and data in the memory to execute any method as claimed in claim 1-8, and the bus connects the functional components for transmitting information.

Images