Tension evaluation system and method based on static theory

Abstract

The invention relates to the technical field of tension evaluation, in particular to a tension evaluation system and method based on a static theory, and the system comprises a node deformation extraction module, a component direction reconstruction module, a force direction feature collection module, a node grade division module and a tension path evaluation module. According to the method, the stress deformation of the component at the node is comprehensively evaluated, the mutual relation between the node position and the component stress is accurately captured, the tension distribution in the structure can be effectively extracted and quantified, and the data such as the axis deviation angle, the displacement change and the attitude angle of the component at the node are integrated, so that the accuracy of the structure is improved. Errors caused by neglecting tiny deformation or complex component interaction in traditional static analysis are avoided, the reliability of the evaluation process and the system response speed are improved through judgment of the acting force direction of the incident component and elimination of a force direction redundant path, the accuracy of data integration is improved in details, and the evaluation accuracy is improved. And the precision and practicability of the evaluation result are ensured.

Description

Technical Field

[0001] This invention relates to the field of tension assessment technology, and in particular to a tension assessment system and method based on static theory. Background Technology

[0002] The field of tension assessment technology involves the quantitative detection and evaluation of tensile forces acting on objects. Its core aspects include accurately measuring the tension borne by components such as wire ropes, cables, and suspension cables in static or dynamic states through mechanical sensor structures, mechanical response devices, and theoretical modeling methods. The entire technical field primarily employs mechanical, photoelectric, and strain gauge tension sensing devices combined with static measurement methods or dynamic response analysis methods to acquire and analyze tension values. Within this field, different measurement schemes select appropriate response mechanisms and establish mechanical models based on the component's usage scenario and stress characteristics to calculate, estimate, and analyze tension values, forming a relatively systematic tension assessment system.

[0003] Among them, the tension assessment system based on static theory refers to the system that calculates and assesses the tensile force on structural members by using Newton's static equilibrium principle. It is used to obtain the internal tension value of the member under static conditions without external dynamic interference. It usually uses known member boundary conditions and structural geometric parameters as the basis for calculation, combines force equilibrium theory to derive the magnitude of the tensile force on the target member, and completes the indirect estimation of tension by measuring the deformation angle or position change of the member under stress.

[0004] Traditional assessment systems rely on static equilibrium derivations based on known boundary conditions and structural geometric parameters. In practical applications, they cannot accurately reflect the complexity of minute deformations and nodal interactions caused by stress on components in real time, leading to significant deviations in tension assessment results. This is especially true in complex structures with large stress variations, where they cannot accurately capture structural details and force transmission paths. Traditional systems fail to fully consider axial offsets and directional changes caused by component deformation and cannot effectively eliminate redundant tension paths, affecting the reliability and accuracy of assessment results. In complex structural analyses requiring high precision, they struggle to cope with variable stress conditions, resulting in errors in assessing structural safety and impacting the accuracy of decision-making. Summary of the Invention

[0005] The purpose of this invention is to address the shortcomings of existing technologies by proposing a tension assessment system and method based on static theory.

[0006] To achieve the above objectives, the present invention adopts the following technical solution: a tension assessment system based on static theory, the system comprising: The node deformation extraction module obtains the connection relationship between components and nodes in the target structure, integrates the component axis offset angle, triaxial displacement value and attitude angle change value of each component at the node, and generates a node component deformation information group. The component orientation reconstruction module calls the node component deformation information group, performs a magnitude judgment on the component axis offset angle and component axis deformation limit for each component, assigns the offset vector and node coordinates to the corresponding labels, and generates the node component update force set; The force feature collection module calls the node components to update the force set, compares the direction angles of the action vectors of the incident components under each node, takes the distribution ratio of the number of component action vectors in the direction of the main axis of the structure as the spatial density value of the direction vector, calculates the mean and difference values ​​of the magnitude of the incident force, and generates a cluster of node direction feature parameters. The node level classification module calls the node orientation feature parameter cluster to classify nodes that meet the three conditions of small included angle, high mean, and high density as main participating nodes, and the rest as edge nodes, and generates a structural node participation identification table. The tension path evaluation module calls the structural node participation identification table and the node component update force vector set, constructs a force transmission path map according to the node, marks the tension redundant path items, removes them from the path, and performs spatial vector summation on the force vectors of non-redundant paths to synthesize the total vector of the node force transmission path and generate the node tension result set.

[0007] The present invention improves upon this invention by including the following: the node component deformation information group includes the relative displacement difference between nodes, the node connection number, and the component deformation constraint state; the node component updated force vector set includes the node force transmission vector, the node force reconstruction number set, and the node connection component projection vector group; the node direction feature parameter cluster includes the node direction distribution clustering results, the force direction spatial coverage, and the angle difference statistical index; the structural node participation identification table includes the main node identification set, the edge node identification set, and the node level classification standard number; and the node tension result set includes the force transmission path spatial vector set, the vector synthesis modulus value, the redundant path elimination record table, and the node force result scalar result.

[0008] The present invention is improved in that the node deformation extraction module includes: The component coordinate acquisition submodule obtains the connection relationship between components and nodes in the target structure, including the connection number of the component and node, and acquires the three-axis coordinate values ​​of the component endpoints before and after being subjected to force. The three-axis coordinate values ​​before and after being subjected to force are calculated separately in the global coordinate system. The difference is calculated to correspond to the spatial displacement change of each component endpoint under the force state, and the difference is classified and organized according to the component number and node number to generate endpoint spatial displacement information. The attitude change calculation submodule extracts the displacement change direction of each component at the node based on the endpoint spatial displacement information, calculates the three-axis angle difference with the component's attitude angle before being subjected to force, calls the node normal direction vector, performs differential calculation on the angle between the attitude angle change value of each component after being subjected to force and the node normal direction, records the angle change value into the angle change set according to the component number, and generates the component attitude change parameter group. The offset angle generation submodule selects the angle between the axial direction vector and the node normal direction vector of each component before it is subjected to force according to the component attitude change parameter group. The angle change value is calculated as the component axis offset angle. The offset angle is indexed and matched with the three-axis displacement value and attitude angle change value, and integrated into a node-level component behavior set to generate a node component deformation information group.

[0009] The present invention is improved in that the component orientation reconstruction module includes: The offset angle judgment submodule calls the node component deformation information group to obtain the component axis offset angle and the corresponding component number attribute of each component. It calculates the difference between the component axis offset angle and the component axis deformation limit. Based on the relationship between the difference and zero, it judges whether the offset condition is met. It records the component number of the offset angle value that exceeds the component axis deformation limit and generates a set of offset component numbers. The force direction substitution calculation submodule extracts the triaxial displacement values ​​of the corresponding components based on the offset component number set, constructs an offset vector with the triaxial displacement values ​​as components, replaces the original axial direction vector of the component with the offset vector, adds component number identifier and node affiliation label to each set of offset vectors, establishes the current force transmission direction set of the component, and generates a sequence of substitution force direction vectors. The component vector aggregation submodule constructs a multi-level index field by adding component number and node number to each vector according to the alternative force vector sequence, and classifies and integrates the offset component vectors according to the node dimension to construct the direction input set in the component set corresponding to the node, and generates the node component update force set.

[0010] The present invention is improved in that the force orientation feature collection module includes: The angle parameter calculation submodule calls the node component to update the force set, extracts all component action vectors under each node and classifies them according to the node number, performs unit vectorization on any two component action vectors under the same node, performs cosine angle calculation, records the minimum angle value and assigns the angle value to the node angle set, and generates a node angle feature set. The orientation density extraction submodule extracts all component action vectors in the node based on the node angle feature set, calculates the angle between each vector direction and the main axis direction of the structure, determines whether the angle is lower than the main axis angle tolerance value, and counts the ratio between the number of component vectors that meet the angle tolerance condition and the total number of incident vectors of the node to generate a node orientation density ratio group. The force value statistics and merging submodule obtains the magnitude of the force vector of all components under each node based on the node direction density ratio group. After grouping based on the node number, it performs mean and range calculations on the force magnitude of each group. It integrates the minimum included angle value, direction density ratio and force value statistics of each node into a multi-field array structure to generate a cluster of node direction feature parameters.

[0011] The present invention is improved in that the specific steps for determining whether the included angle is lower than the principal axis included angle tolerance value are as follows: the included angle value between the action vector of each component and the direction of the principal axis of the structure is calculated and compared with the principal axis included angle tolerance value, which is the sum of the average included angle value between the action vector of each component and the direction of the principal axis and the standard deviation at the principal axis included angle tolerance value node.

[0012] The present invention is improved in that the node level classification module includes: The included angle threshold judgment submodule calls the node direction feature parameter cluster, extracts the minimum included angle value of each node, compares the included angle value with the component convergence angle threshold item by item, and includes the node number that meets the condition of being less than the threshold into the candidate set, generating a set of included angle compliant nodes. The orientation density filtering submodule extracts the orientation vector spatial density value based on the set of compliant nodes with the included angle, compares it with the set main orientation component ratio threshold, filters the node numbers with a ratio not lower than the threshold, and generates an effective orientation density node set. The tension strength comparison submodule obtains the average incident force modulus value of the corresponding node based on the effective node set of the directional density, performs difference calculation with the tension benchmark value of the component, and calculates the node participation strength index value of each node by combining the included angle value and the directional density value. It then compares the index value with the participation level threshold, classifies nodes with node participation strength index values ​​less than the participation level threshold as edge nodes, and classifies nodes with index values ​​greater than or equal to the threshold as main participating nodes, thus establishing a structural node participation identification table.

[0013] The present invention is improved in that the formula for calculating the node participation strength index value of each node is specifically as follows: ; in, Indicates the first The participation strength index value of each node, Indicates the first The normalized value of the mean tension modulus of the incident component at each node. This represents the normalized value of the component's tension reference value. Indicates the first The spatial density value of the direction vector of each node. Indicates the first The normalized value of the minimum included angle of each node. Indicates the first Normalized standard deviation of the tension values ​​of the incident component at each node The directional density attenuation coefficient, is the base of the natural logarithm.

[0014] The present invention is improved in that the tension path evaluation module includes: The path graph generation submodule calls the structural node participation identifier table and the node component update force set, filters the main participating nodes, extracts the corresponding incident component action vector, groups the incident vectors according to the node number and records the direction attribute, constructs the force path set graph of each group of node corresponding to the incident component vector set, and generates the node path construction graph group. The redundant path removal submodule constructs a graph group based on the path of the upper node, obtains the angle between any two vectors in the action vector group of each node, filters vector pairs with an angle smaller than the path redundancy judgment angle threshold, calculates the difference between the vector magnitude ratio and the force direction overlap judgment ratio threshold, calculates the redundant force direction compression coefficient of the node, marks the node vector pairs with the redundant force direction compression coefficient greater than the redundancy judgment compression threshold as tension redundant path items, and removes them from the path set to obtain the redundant path pruning result set; The node resultant force extraction submodule extracts the remaining non-redundant action vectors in each node based on the redundant path trimming result set, and sums the components in three-dimensional space according to the vector coordinate values ​​to synthesize the total force vector of the node on the action plane. The magnitude of the total force vector is extracted as the tension output index to establish the node tension result set.

[0015] A tension assessment method based on static theory, used to implement the aforementioned tension assessment system based on static theory, includes the following steps: S1: Obtain the connection relationship between components and nodes in the target structure, integrate the component axis offset angle, triaxial displacement value and attitude angle change value of each component at the node, and generate a node component deformation information group; S2: Call the node component deformation information group, perform a magnitude judgment on the component axis offset angle and component axis deformation limit for each component, assign the offset vector and node coordinates to the corresponding markers, and generate a node component update force set; S3: Call the updated force set of the node component, compare the direction angles of the action vectors of the incident components under each node, take the distribution ratio of the number of component action vectors in the direction of the main axis of the structure as the spatial density value of the direction vector, calculate the mean and difference of the magnitude of the incident force, and generate the node direction feature parameter cluster. S4: Call the node orientation feature parameter cluster, classify nodes that meet the three conditions of small included angle, high mean, and high density as main participating nodes, and otherwise as edge nodes, and generate a structural node participation identification table; S5: Call the structural node participation identification table and update the force vector set of the node components, construct the force transmission path map according to the node, mark the tension redundant path items, remove them from the path, and perform spatial vector summation on the force vectors of the non-redundant paths to synthesize the total vector of the node force transmission path and generate the node tension result set.

[0016] Compared with the prior art, the advantages and positive effects of the present invention are as follows: In this invention, by comprehensively evaluating the stress and deformation of components at nodes, the relationship between node position and component stress is accurately captured, effectively extracting and quantifying the tension distribution inside the structure, and achieving optimized analysis of force transmission paths. By integrating data such as the axial offset angle, displacement change, and attitude angle of components at nodes, more accurate force vector data is formed, avoiding errors caused by neglecting subtle deformations or complex component interactions in traditional static analysis. The judgment of the force direction of incident components and the elimination of redundant force paths ensure accurate evaluation of tension paths. Furthermore, the reliability of the evaluation process and the system response speed are improved through precise node participation division. By improving the accuracy of data integration in detail, the tension evaluation method is optimized, ensuring the accuracy and practicality of the evaluation results. Attached Figure Description

[0017] Figure 1 This is a system flowchart of the present invention; Figure 2 This is a flowchart of the node deformation extraction module of the present invention; Figure 3 This is a flowchart of the component orientation reconstruction module of the present invention; Figure 4 This is a flowchart of the force orientation feature aggregation module of the present invention; Figure 5 This is a flowchart of the node level division module of the present invention; Figure 6 This is a flowchart of the tension path evaluation module of the present invention. Detailed Implementation

[0018] To make the objectives, technical solutions, and advantages of this invention clearer, the invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the invention.

[0019] In the description of this invention, it should be understood that the terms "length," "width," "upper," "lower," "front," "rear," "left," "right," "vertical," "horizontal," "top," "bottom," "inner," and "outer," etc., indicating orientation or positional relationships, are based on the orientation or positional relationships shown in the accompanying drawings and are only for the convenience of describing the invention and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, or be constructed and operated in a specific orientation, and therefore should not be construed as a limitation of the invention. Furthermore, in the description of this invention, "a plurality of" means two or more, unless otherwise explicitly specified.

[0020] Please see Figure 1 The present invention provides a technical solution: a tension assessment system based on static theory, the system including a node deformation extraction module, a component orientation reconstruction module, a force orientation feature collection module, a node level classification module, and a tension path assessment module; The node deformation extraction module obtains the connection relationship between components and nodes in the target structure, collects the triaxial coordinate values ​​and attitude angle changes of each component endpoint before and after being subjected to force, and records the change in the angle between the axis direction vector and the node normal before being subjected to force as the component axis offset angle. It integrates the component axis offset angle, triaxial displacement value and attitude angle change value of each component at the node to generate a node component deformation information group. The component orientation reconstruction module calls the node component deformation information group and performs a size judgment on the component axis offset angle and component axis deformation limit for each component. If it is greater than the component axis deformation limit, an offset vector is constructed using the triaxial displacement value and replaced with the original axis direction vector as the current force input of the component. At the same time, the offset vector and node coordinates are assigned and marked, and integrated into a component action vector set to generate the node component update force set. The allowable range of axial deformation of a component refers to the range of axial deflection angles that a component can undergo under normal working conditions. It is usually defined according to the structural material mechanics code and generally does not exceed 1° to 5°. It can be set according to the maximum allowable value of the rotation angle before yielding of different materials. The force feature collection module calls the node components to update the force set, compares the direction angles of the action vectors of the incident components under each node pairwise, obtains the minimum angle value and the angle set, takes the distribution ratio of the number of component action vectors in the direction of the main axis of the structure as the spatial density value of the direction vector, calculates the mean and difference of the magnitude of the incident force, and classifies the three data into a node feature vector set to generate a cluster of node direction feature parameters. The spatial density value of the direction vector refers to the proportion of the incident vector of a component at a certain node in a three-dimensional coordinate system that is distributed in the direction of the principal axis of the structure (usually the direction of the main load transmission). It can be obtained by normalizing the count of the direction vector in the solid angle region of ±5° of the principal axis in spherical coordinates. The node level classification module calls the node direction feature parameter cluster, judges the minimum included angle value of each node with the component convergence angle threshold, and combines the direction vector spatial density value with the main direction component proportion threshold to determine whether it is lower than the participation judgment lower limit. At the same time, it compares the mean value with the component tension benchmark value. Nodes that meet the three conditions of small included angle, high mean value and high density are classified as main participating nodes, and the rest are edge nodes, generating a structural node participation identification table. The reference value of component tension is the lower limit value for judging the participation of tension in the static analysis of the structure. It is usually set according to the initial tensile force calculation value of the structural component or the standard design load, and is often expressed in N or kN.

[0021] The tension path assessment module calls the structural node participation identification table and the node component update force vector set, filters the main participating nodes and corresponding incident component action vectors, constructs the force transmission path map according to the nodes, compares the pairwise ratios of the action vectors with the same node under the path map whose directional angle is less than the path redundancy judgment angle threshold, if the ratio difference is less than the force direction overlap judgment ratio threshold, the smaller vector is marked as a tension redundancy path item and removed from the path, obtains the tension flow path set after direction trimming and node screening, and performs spatial vector summation on the force vectors of non-redundant paths to synthesize the total vector of the node force transmission path, extracts the modulus as the tension value, and generates the node tension result set; The force overlap determination threshold refers to the threshold set for the difference in the force magnitude ratio between two approximately directional forces. If the difference between the two is within 10% and the directional angles are close, they are considered as duplicate paths.

[0022] The node component deformation information group includes the relative displacement difference of nodes, the node number of the component connection, and the deformation constraint state of the component. The node component updated force vector set includes the node force transmission vector, the node force reconstruction number set, and the node connection component projection vector group. The node direction feature parameter cluster includes the node direction distribution clustering results, the spatial coverage of the force direction, and the statistical index of the angle difference. The structural node participation identification table includes the main node identification set, the edge node identification set, and the node level classification standard number. The node tension result set includes the force transmission path spatial vector set, the vector synthesis modulus value, the redundant path elimination record table, and the node force result scalar result.

[0023] Please see Figure 2 The node deformation extraction module includes: The component coordinate acquisition submodule obtains the connection relationship between components and nodes in the target structure, including the connection number of the component and node, and acquires the three-axis coordinate values ​​of the component endpoints before and after being subjected to force. The three-axis coordinate values ​​before and after being subjected to force are calculated separately in the global coordinate system. The difference is calculated to correspond to the spatial displacement change of each component endpoint under the force state, and the difference is classified and organized according to the component number and node number to generate endpoint spatial displacement information. The pre-set structural data information table is retrieved. This table records in detail the unique number of each component and the node numbers connected to its two ends. For example, component C1 connects to node N1 and node N2, and component C2 connects to node N1 and node N3. Subsequently, before the structure is subjected to force, the spatial coordinates of each node are measured using a 3D laser scanner or displacement sensors deployed at the endpoints of the components to obtain their initial triaxial coordinate values ​​in the global coordinate system. For example, the initial coordinates of node N1 are... (unit: The initial coordinates of node N2 are: (unit: When the structure is subjected to external loads, the spatial coordinates of each node are measured again using the same method to obtain the triaxial coordinate values ​​after the load is applied. For example, the coordinates of node N1 after the load is applied become... (unit: The coordinate differences between the three axes before and after the force is applied are calculated in the global coordinate system. Specifically, the coordinate components after the force is applied are subtracted from the corresponding coordinate components before the force is applied. Taking the calculation node N1 as an example, its... The displacement change in the axial direction is , The displacement change in the axial direction is , The displacement change in the axial direction is The difference vector corresponds to the spatial displacement change of each component endpoint under stress. This refers to the spatial displacement of node N1, which is then categorized and organized according to component number and node number. The spatial displacement of endpoint N1 of component C1 is recorded as follows: This integrates all such information from the endpoints of the components to generate endpoint spatial displacement information.

[0024] The attitude change calculation submodule extracts the displacement change direction of each component at the node based on the spatial displacement information of the endpoints, calculates the three-axis angle difference with the attitude angle of the component before being subjected to force, calls the node normal direction vector, performs differential calculation on the angle between the attitude angle change value of each component after being subjected to force and the node normal direction, records the angle change value into the angle change set according to the component number, and generates the component attitude change parameter group. Extract the direction of displacement change of each component at the node. The direction is determined by the difference of the three-axis coordinates in the spatial displacement information of the endpoints. The displacement vector formed is characterized by the displacement vector of component C1 at node N1. For example, the three-axis angle difference is calculated between the component's attitude angle before being subjected to force and the component's attitude angle before being subjected to force. The attitude angle of component C1 before being subjected to force passes through its two endpoints N1. and N2 The initial coordinates are calculated, and the initial axis direction vector is: The attitude angle corresponding to the global coordinate system is After being subjected to force, the displacement vector of its endpoint N1 This results in a slight attitude shift, triggering the use of the node normal vector. The node normal vector is a unit vector representing the characteristics of the node's connection plane, calculated based on the initial geometric configuration of the multiple components connected to the node. For example, the normal vector of node N1 is set to... For each component, the angle between the change in attitude angle after being subjected to force and the angle between the angle and the normal direction of the node is calculated using a differential method. Specifically, the displacement vector is calculated first. With node normal vector The included angle Then calculate the initial axis direction vector of component C1. With node normal vector The included angle In this example Then calculate and The difference is the angular change of the component's attitude relative to the node normal. The angular change value is recorded in the angular change set according to the component number, for example, recorded as (C1, N1, ...). ),in This is the calculated angle difference. Perform this operation on all components to generate a set of component attitude change parameters.

[0025] The offset angle generation submodule selects the axial direction vector of each component before it is subjected to force and calculates the angle between the node normal direction vector and the component attitude change parameter group. The calculated angle change value is used as the component axial offset angle. The offset angle is indexed and matched with the three-axis displacement value and attitude angle change value, and integrated into a node-level component behavior set to generate a node component deformation information group. The angle between the axial direction vector and the node normal direction vector of each component before loading is selected and calculated. Taking component C1 as an example, its axial direction vector before loading is... The normal direction vector of the node N1 it connects to is The angle between the two can be calculated using the dot product formula. Subsequently, based on the component's attitude after being subjected to force (determined by the original attitude and attitude change parameters), the angle between the component's axial direction vector and the normal direction vector at the same node after being subjected to force is calculated. Calculate the change in the included angle. As the component axis offset angle, this The value is the axis offset angle of component C1. The offset angle is matched with the triaxial displacement value and attitude angle change value by index matching. Specifically, the axis offset angle of component C1 at node N1 is... Triaxial displacement values and attitude angle change value By associating component number C1 and node number N1, they are integrated into a single record. This record is the set of behaviors of component C1 at node N1 level. Perform this operation on all components under all nodes to generate a node component deformation information group.

[0026] Please see Figure 3 The component orientation reconstruction module includes: The offset angle judgment submodule calls the node component deformation information group to obtain the component axis offset angle and the corresponding component number attribute of each component. It calculates the difference between the component axis offset angle and the component axis deformation limit, and judges whether the offset condition is met based on the relationship between the difference and zero. It records the component number of the component whose offset angle value exceeds the component axis deformation limit and generates a set of offset component numbers. Extract the axial offset angle of component C1 at node N1 from the information group. The offset angle of component C2 on the axis at node N1 is . The difference between the component axis offset angle and the component axis deformation limit is calculated. This component axis deformation limit is determined through a series of structural mechanics analysis experiments based on the component's material properties, cross-sectional dimensions, and connection method. For example, for Q345B steel components in the current structure, tensile-torsional coupled loading experiments were conducted on 100 similar components, recording the critical rotation angle data as they transitioned from the elastic to the plastic stage. After removing the highest and lowest 5% of values, the average of the remaining data was taken, and the final deformation limit was set as follows: The offset angle of component C1 With tolerance Perform the difference calculation to obtain The offset condition is determined based on the relationship between the difference and zero. If the difference is greater than zero, component C1 is determined to satisfy the offset condition. For component C2, the difference is... If the value is less than zero, it is determined that the offset condition is not met. The component number of the component whose offset angle value exceeds the deformation limit of the component axis is recorded. The component number C1 is recorded. The judgment results of all components are summarized to generate a set of offset component numbers.

[0027] The force direction substitution calculation submodule extracts the triaxial displacement values ​​of the corresponding components based on the set of offset component numbers, constructs an offset vector with the triaxial displacement values ​​as components, replaces the original axial direction vector of the component with the offset vector, adds component number identifier and node affiliation label to each set of offset vectors, establishes the current force transmission direction set of the component, and generates a sequence of substitution force direction vectors. Obtain component number C1 from the offset component number set, and then extract the displacement value of C1 at node N1 from the endpoint spatial displacement information. And construct an offset vector using the displacement values ​​in the three axes as components, that is, construct a vector The offset vector is replaced with the original axial direction vector of the component. The original axial direction vector of component C1 is... Now, the direction of force transmission at node N1 is changed from Replace with For components not appearing in the offset component number set, such as C2, their force transmission direction remains the original axial direction vector. A component number identifier and node affiliation label are added to each set of offset vectors. Labeled as (component: C1, node: N1), establish the current force transmission direction set of the component. This set contains the latest force transmission direction vectors of all components at each node, and generate a sequence of alternative force direction vectors.

[0028] The component vector collection submodule constructs a multi-level index field by adding component number and node number to each vector according to the sequence of substitute force vectors, and classifies and integrates the offset component vectors according to the node dimension to construct the direction input set in the component set corresponding to the node, and generates the node component update force set. Vectors from the sequence The components and their identifiers (component: C1, node: N1) are organized into a single data object, with node number N1 serving as the primary index and component number C1 as the secondary index. The offset component vectors are then categorized and integrated according to the node dimension. Node N1 is also set to connect components C2 and C3, where C2 has not been offset and its force transmission vector is the original axis vector. C3 is shifted, and its substitution force vector is Then these three vectors , , The input directions are aggregated under the name of node N1, and the direction input set in the corresponding component set of the node is constructed. Specifically, it is a mapping with node N1 as the key and a data list containing three vectors as the value. This operation is repeated for all nodes in the structure to generate the updated force vector set of the node component.

[0029] Please see Figure 4 The force-direction feature aggregation module includes: The angle parameter calculation submodule calls the node component to update the force set, extracts all component action vectors under each node and classifies them according to the node number, performs unit vectorization on any two component action vectors under the same node, performs cosine angle calculation, records the minimum angle value and assigns the angle value to the node angle set, and generates the node angle feature set. Extract all component action vectors under each node and categorize them according to node number. Taking node N1 as an example, extract its three action vectors. , , For any two components acting at the same node, perform unit vectorization processing. For example, for and Unitize to obtain and And perform cosine angle calculation, calculate and The included angle , and The included angle ,as well as and The included angle Then compare the values ​​of these three included angles. Record the minimum included angle value and add it to the node angle set. Then, calculate the... It is stored in association with node number N1, and the same calculation process is performed on all nodes to generate a set of node angle features.

[0030] The orientation density extraction submodule extracts all component action vectors in a node based on the node angle feature set, calculates the angle between each vector direction and the main axis direction of the structure, determines whether the angle is lower than the main axis angle tolerance value, and calculates the ratio between the number of component vectors that meet the angle tolerance condition and the total number of incident vectors of the node to generate a node orientation density ratio group. The specific steps to determine whether the included angle is lower than the principal axis included angle tolerance value are as follows: calculate the included angle value between the action vector of each component and the direction of the principal axis of the structure, and compare it with the principal axis included angle tolerance value. The sum of the average included angle value and the standard deviation of the action vector of each component and the direction of the principal axis under the principal axis included angle tolerance value node; Call the three action vectors of node N1 , , The angle between each vector direction and the principal axis direction of the structure is calculated. The principal axis direction is predefined based on the overall geometry of the structure and its main load-bearing directions. For example, for a truss bridge, the principal axis direction can be defined as a vector along the bridge length. Here, let the principal axis direction vector be... Calculate separately , , and The included angle is obtained To determine whether the included angle is lower than the principal axis included angle tolerance value, the specific steps are as follows: Calculate the included angle between the action vector of each component and the principal axis direction of the structure, and compare it with the principal axis included angle tolerance value. The sum of the average included angle value and standard deviation of the action vector of each component and the principal axis direction under the principal axis included angle tolerance node is used. First, calculate the average value of the three included angles under node N1. Then calculate its standard deviation. Then the upper limit of the tolerance value of the included spindle angle is The lower limit is The ratio between the number of component vectors satisfying the included angle tolerance condition and the total number of incident vectors at the nodes is calculated. and Located within the tolerance range, and If the vector is located outside the interval, then the number of vectors satisfying the condition is 2, and the total number of incident vectors at node N1 is 3. The calculated ratio is... This is the directional density ratio of node N1. This ratio is associated with the node number N1, and the calculation is performed on all nodes to generate a group of node directional density ratios.

[0031] The force value statistics and merging submodule obtains the magnitude of the force vector of all components under each node based on the node direction density ratio group. After grouping based on the node number, it performs mean and range calculations on the force magnitude of each group. It integrates the minimum included angle value, direction density ratio and force value statistics of each node into a multi-field array structure to generate a cluster of node direction feature parameters. The magnitudes of the force vectors acting on all components at each node are obtained. These magnitudes are calculated by measuring strain using strain gauges or fiber optic grating sensors mounted on the components, and then combining this with the component's material elastic modulus and cross-sectional area. For example, the magnitudes of the force vectors acting on the three components C1, C2, and C3 at node N1 are as follows: , , After grouping by node number, the mean and range of force magnitudes in each group are calculated. For node N1, the mean force magnitude is... The range is The minimum included angle value, direction density ratio, and force direction value statistics for each node are uniformly integrated into a multi-field array structure. For example, for node N1, the minimum included angle value calculated in the previous paragraphs is integrated. directional density ratio Mean value of force modulus and range , forming a record Perform this operation on all nodes to generate a cluster of node orientation feature parameters.

[0032] Please see Figure 5 The node level classification module includes: The included angle threshold judgment submodule calls the node direction feature parameter cluster, extracts the minimum included angle value of each node, compares the included angle value with the component convergence angle threshold item by item, and includes the node number that meets the condition of being less than the threshold into the candidate set, generating a set of included angle compliant nodes. Extract the minimum included angle value of node N1. The minimum included angle value of node N2 The minimum included angle value of node N3 The included angle value is compared with the component convergence angle threshold item by item. The component convergence angle threshold is set with reference to the requirements for node construction in relevant structural design codes. Combined with the finite element simulation results, for a certain type of node, when the included angle between components is less than a certain value, the stress concentration in the node domain will increase sharply. Here, stress analysis of node models with 10 different convergence angles is carried out. It is found that when the included angle is less than 20°, the maximum principal stress increases by more than 30%. Therefore, the component convergence angle threshold is set to 20°. The 25° of node N1 is compared with 20°. 25°>20°, which does not meet the condition of being less than the threshold. The 18° of node N2 is compared with 20°. 18°<20°, which meets the condition. The 30° of node N3 is compared with 20°. 30°>20°, which does not meet the condition. The node numbers that meet the condition of being less than the threshold are included in the candidate set. Therefore, the number "N2" of node N2 is included in this set. Finally, after judging all nodes, a set of nodes with compliant included angles is generated.

[0033] The orientation density filtering submodule extracts the spatial density value of the orientation vector based on the set of nodes with compliant included angles, compares it with the set ratio threshold of the main orientation component, filters the node numbers with a ratio not lower than the threshold, and generates a set of valid orientation density nodes. The spatial density values ​​of the directional vectors are extracted. The set of nodes with compliant included angles includes nodes N2 and N5. From the node directional density ratio group, the directional density value of N2 is extracted to be 0.85, and the directional density value of N5 is 0.60. These are compared with a set threshold for the proportion of main directional components. This threshold is based on the analysis of the force transmission path of the entire structure under the main load conditions. By calculating the distribution of all node directional density values, it was found that the directional density values ​​of nodes on the main force transmission path are generally concentrated above 0.75. To screen these key nodes, the threshold for the proportion of main directional components is set to the upper quartile of the distribution of all node directional density values, calculated to be 0.75. The directional density value of node N2 (0.85) is compared with the threshold of 0.75; 0.85 is not lower than 0.75, thus meeting the condition. The directional density value of node N5 (0.60) is compared with the threshold of 0.75; 0.60 is lower than 0.75, thus not meeting the condition. Nodes with ratios not lower than the threshold are selected, thus node N2 is selected, finally obtaining the effective node set for directional density.

[0034] The tension strength comparison submodule obtains the average incident force modulus value of the corresponding node based on the effective node set of directional density, performs difference calculation with the component tension reference value, and combines the included angle value and directional density value using the following formula: ; The node participation strength index value of each node is obtained by calculation and compared with the participation level threshold. Nodes with a node participation strength index value less than the participation level threshold are classified as edge nodes, and nodes with an index value greater than or equal to the threshold are classified as main participating nodes. A structural node participation identification table is established. in, Indicates the first The participation strength index value of each node, Indicates the first The normalized value of the average tension modulus of the incident component at each node is obtained by dividing the original average modulus by the maximum average modulus of all nodes in the set. The normalized value representing the reference value of the component tension is obtained by dividing the reference value by the upper limit of the nodal tension modulus. Indicates the first The spatial density value of the direction vectors of each node is calculated by dividing the number of incident vectors that satisfy the directional angle tolerance by the total number of incident vectors. Indicates the first The normalized value of the minimum included angle of each node is obtained by dividing the minimum included angle of that node by the maximum included angle limit defined by the system, thereby eliminating the influence of angle units on the calculation. Indicates the first The normalized standard deviation of the tension values ​​of the incident components at each node is obtained by calculating the standard deviation of the tension modulus values ​​of all incident components at that node and then normalizing it. This is the directional density attenuation coefficient, which adjusts the degree of suppression on density weights based on the structural directional aggregation characteristics. The base of the natural logarithm is used, and the participation level threshold is set based on taking all... After the values ​​form a distribution set, a threshold is set based on the quantile of the distribution median, and the upper quartile is used as the dividing point for participation or non-participation, so as to balance the judgment criteria of node participation coverage and effective participation.

[0035] Obtain the average incident force modulus value of the corresponding node, calculate the difference with the component tension reference value, and combine the included angle value and directional density value using the following formula: ; The node participation strength index value for each node is calculated, and the parameters in the formula are defined as follows. Indicates the first The participation strength index value of each node, It is the first The normalized value of the mean tension modulus of the incident component at each node. It is the normalized value of the component tension reference value. It is the first The spatial density value of the direction vector of each node. It is the first The normalized value of the minimum included angle of each node. It is the first Normalized standard deviation of the tension values ​​of the incident component at each node The directional density attenuation coefficient, The formula uses the natural logarithm as its base. The calculation logic is that the numerator evaluates the deviation between the node's average tension and the benchmark, and corrects it using the directional density. The denominator comprehensively considers the compactness of the node's geometric convergence (minimum included angle) and the dispersion of the force distribution (standard deviation). Overall, the formula quantifies the comprehensive importance of a node in terms of mechanical performance and geometric structure. Taking node N2 as an example, the detailed calculation process is as follows: First, the parameters are obtained and assigned. The effective node set of the directional density includes three nodes: N2, N4, and N8. Their original data are shown in Table 1. Table 1 Original Node Parameters

[0036] Table 1 lists the original parameters of each node required to calculate the node participation strength index value. The parameters are then calculated and set. The maximum value of the average of all moduli in the node set is 400. ,therefore (Normalized mean of node N2) = The reference value for component tension is the safe bearing capacity calculated based on the allowable stress of the material and the average cross-sectional area of ​​the component, and is set at 300. The upper limit of the nodal tension modulus is taken as 500 according to the design specifications. ,but (Normalized value of the benchmark) = , (Directional density value of node N2) = The system defines the maximum included angle as 90°. (Normalized minimum included angle value of node N2) = The standard deviation of the tension modulus values ​​of all incident components at node N2 is The maximum standard deviation of all nodes is ,but (Normalized standard deviation of node N2) directional density attenuation coefficient The setting is based on the fact that when the structural directional aggregation characteristics are obvious, that is, most components tend to converge toward the principal axis, the density value is generally high, and in this case, it should be reduced. To reduce its inhibitory effect, and vice versa. This section analyzes historical data with similar structures. The value of is positively correlated with the slenderness ratio of the structure. Given a large current slenderness ratio, the setting... Substitute the above parameter values ​​into the formula to perform the calculation: ; The advantage of the formula lies in the introduction of directional density. exponential decay term This suppresses the index values ​​of nodes with highly concentrated directions, avoiding false strong points caused by parallel local components. Simultaneously, the logarithmic term in the denominator... It can smoothly handle the influence of tension standard deviation, enhancing the model's robustness to force fluctuations. Subsequently, its value is calculated for all effective nodes. Value, set to get set The participation level threshold is set based on taking all After constructing the distribution set of values, a threshold is set based on the quantiles of the distribution's median. The upper quartile is used as the dividing point for participation or non-participation. The upper quartile (75th percentile) of this set is calculated to be 0.0889. The values ​​of each node are then... The value is compared with the threshold of 0.0889, and the value of node N2 is... Classified as an edge node, node N8 The nodes were classified as primary participating nodes, and the results indicate that node N8 has a higher overall importance under the current operating conditions. A structural node participation identification table was then established.

[0037] Please see Figure 6 The tension path evaluation module includes: The path graph generation submodule calls the structural node participation identifier table and the node component update force set, filters the main participating nodes, extracts the corresponding incident component action vector, groups the incident vectors by node number and records the direction attribute, constructs the force path set graph of each group of node corresponding to the incident component vector set, and generates the node path construction graph group. Filter the node numbers marked as "major participating nodes" from the identifier table, for example, filter nodes N8 and N15. Then, based on these node numbers, extract all incident member action vectors corresponding to these nodes from the node member update force vector set. For example, extract the action vector set under node N8. and the set of action vectors under node N15 The incident vectors are grouped by node number and their directional attributes are recorded. That is, the above vector sets are associated with node numbers N8 and N15 respectively. The directional attribute of each vector is its three-dimensional coordinate component. The set of incident component vectors corresponding to each node is constructed into a force path set graph. In terms of data structure, this graph is a hash table or dictionary with the main participating node number as the index and its value being a list of all incident force vectors of that node, thus generating a node path construction graph group.

[0038] The redundant path removal submodule constructs a graph group based on the path of the upper node, obtains the angle between any two vectors in the action vector group of each node, filters vector pairs whose angle is less than the path redundancy judgment angle threshold, and calculates the difference between the vector magnitude ratio and the force direction overlap judgment ratio threshold, using the formula: ; The redundant force compression coefficient of the node is obtained by calculation. The node vector pairs with the redundant force compression coefficient greater than the redundancy judgment compression threshold are marked as tension redundant path items and removed from the path set to obtain the redundant path pruning result set. in, Indicates the first The redundancy force of each node is directed towards the compression coefficient. Indicates the first The number of vector pairs in a node that satisfy the redundancy criterion is derived from the statistics of vector combinations within that node whose included angle is below the included angle threshold. Indicates the first The node in the node The spatial angle between vectors is obtained by calculating the vector dot product. , They represent the first The node in the node The two force magnitudes of the vector are obtained from the vector magnitude in three-dimensional coordinate space. This represents a very small positive constant to prevent the denominator from being zero. Indicates the first The node in the node The direction aggregation weights for vectors are derived from the normalized distribution density of the incident vector at that node along the principal axis of the structure. The redundancy compression threshold is set by considering all nodes... The distribution analysis of the set of values ​​is performed, and the sum of the mean and standard deviation is taken as the compression threshold. All vector pairs that exceed the threshold are marked as redundant terms. Vector pairs with an included angle less than the path redundancy threshold are selected, and the difference between the vector magnitude ratio and the force overlap ratio threshold is calculated using the formula: ; The redundancy force compression coefficient of the node is calculated, and the parameters in the formula are defined as follows. It is the first The redundancy force of each node is directed towards the compression coefficient. It is the first The number of vector pairs whose included angle within a node is less than the path redundancy determination angle threshold. It is the first The node in the node For the spatial angle between vectors, , It is the magnitude of the vector pair. It is a very small positive constant to prevent the denominator from being zero. This is the directional aggregation weight of the vector pair. The calculation logic of the formula lies in summing up the redundant contributions of all approximately parallel vector pairs within the node through the summation term, where the first term... At the angle It increases sharply as it approaches 0, thus amplifying the influence of collinear vector pairs, the second term This is a normalized force difference term. When the magnitudes of two forces are similar, this term approaches 0, suppressing redundant calculations for vector pairs with significantly different magnitudes. (The third term...) This introduces the consideration of global directionality. Taking the main participating node N8 as an example, it has three vectors. Their moduli are respectively , , The path redundancy determination angle threshold is set based on experimental data. Through analysis of a large number of nodes, it was found that when the component angle is less than 10°, their contribution directions to the resultant force of the node highly overlap. Therefore, the threshold is set to 10°. Calculations show that... and The included angle Less than 10° and The included angle is 35°. and The included angle is 40°, and all are greater than 10°. Therefore, the number of vector pairs in node N8 that satisfy the redundancy criterion is [missing information]. That is, vector pairs ,constant Set as directional aggregation weight This is derived from the normalized value of the distribution density of the incident vector at this node along the principal axis of the structure, i.e., the node direction density value calculated earlier. Substitute into the formula to calculate the redundancy force compressibility coefficient of N8: ; The advantage of this formula lies in the fact that it doesn't simply eliminate paths based on the angle size in a "one-size-fits-all" manner. Instead, it integrates three dimensions: directional similarity, force balance, and global importance (directional aggregation weight). This allows for a refined and quantitative evaluation of redundant paths. The redundancy judgment compression threshold is set based on all nodes. The set of values ​​is subjected to distribution analysis, and the sum of its mean and standard deviation is taken as the compression threshold. After calculation for all main participating nodes, the following is obtained: The average value is 0.0030, and the standard deviation is 0.0020. Therefore, the threshold is... The redundancy of the nodes is directed towards the compression coefficient. Compared to the threshold of 0.0050, since The vector pairs in node N8 Marked as a tension redundant path term, and remove the vector with the smaller magnitude from the path set, i.e., remove. This yields the redundant path pruning result set.

[0039] The node resultant force extraction submodule extracts the remaining non-redundant action vectors in each node based on the redundant path pruning result set, and adds the components in three-dimensional space according to the vector coordinate values ​​to synthesize the total force vector of the node on the action plane. The magnitude of the total force vector is extracted as the tension output index to establish the node tension result set. Extract the remaining non-redundant action vectors from each node. Taking node N8 as an example, after redundant path removal, its remaining non-redundant action vectors are: and The set coordinate components are respectively and All units are Then, the components of the vector coordinates are summed in three-dimensional space. Specifically, the calculation involves adding the corresponding coordinate components of all remaining vectors. Component resultant force , Component resultant force , Component resultant force The total force vector of the composite node on the plane of action, i.e. The magnitude of the total force vector is extracted as the tension output index, and the magnitude is calculated. The modulus value is 293.9. As the final tension output of node N8, and associated with node number N8, the same operation is performed on all nodes in the trimming result set to establish the result set of node tension result forces.

[0040] A tension assessment method based on static theory includes the following steps: S1: Obtain the connection relationship between components and nodes in the target structure, integrate the component axis offset angle, triaxial displacement value and attitude angle change value of each component at the node, and generate a node component deformation information group; S2: Call the node component deformation information group, perform a magnitude judgment on the component axis offset angle and component axis deformation limit for each component, assign the offset vector and node coordinates to the corresponding labels, and generate the node component update force set; S3: Call the node component to update the force set, compare the direction angles of the action vectors of the incident components under each node, take the distribution ratio of the number of component action vectors in the direction of the main axis of the structure as the spatial density value of the direction vector, calculate the mean and difference of the magnitude of the incident force, and generate the node direction feature parameter cluster. S4: Call the node orientation feature parameter cluster, classify nodes that meet the three conditions of small included angle, high mean, and high density as main participating nodes, and otherwise as edge nodes, and generate a structural node participation identification table. S5: Call the structural node participation identification table and update the force vector set of node components, construct the force transmission path map according to the node, mark the smaller vector as tension redundant path item, remove it from the path, and perform spatial vector summation on the force vectors of non-redundant paths to synthesize the total vector of the node force transmission path and generate the node tension result set.

[0041] The above are merely preferred embodiments of the present invention and are not intended to limit the present invention in any other way. Any person skilled in the art may make changes or modifications to the above-disclosed technical content to create equivalent embodiments that can be applied to other fields. However, any simple modifications, equivalent changes, and modifications made to the above embodiments based on the technical essence of the present invention without departing from the scope of the present invention shall still fall within the protection scope of the present invention.

Claims

1. A tension assessment system based on statics theory, characterized in that, The system comprises: The node deformation extraction module obtains the connection relationship between the members and nodes in the target structure, integrates the member axis offset angle, three-axis displacement value and attitude angle change value of each member at the node, and generates a node member deformation information group; The member direction reconstruction module calls the node member deformation information group, performs size judgment on the member axis offset angle and member axis deformation tolerance of each member, attributes the offset vector and node coordinates, generates a node member update force vector set, and calls the node member update force vector set; The force vector feature collection module calls the node member update force vector set, compares the direction included angle of the action vector of each incident member under each node, takes the quantity distribution ratio of the member action vector in the structure main shaft direction as the direction vector space density value, calculates the mean value and difference value of the incident force vector size, generates a node direction feature parameter cluster, and calls the node direction feature parameter cluster. The node level division module classifies the nodes that meet the three conditions of small included angle, high mean value and high density as main participating nodes, and otherwise as edge nodes, generates a structure node participation identification table, and calls the structure node participation identification table and the node member update force vector set. The tension path evaluation module constructs a force vector transmission path atlas according to the node division, marks the tension redundant path item, eliminates it from the path, and adds the force vectors of the non-redundant path to synthesize the total vector of the node force transmission path, generates a node tension resultant force result set.

2. The tension evaluation system based on statics theory according to claim 1, characterized in that, The node member deformation information group comprises a node relative displacement difference value, a member connection node number and a member deformation constraint state, the node member update force vector set comprises a node action force transmission vector, a node force vector reconstruction number set and a node connection member projection vector group, the node direction feature parameter cluster comprises a node direction distribution clustering result, a force vector direction space coverage and an included angle difference statistical index, the structure node participation identification table comprises a main force node identification set, an edge node identification set and a node level division standard number, and the node tension resultant force result set comprises a force transmission path space vector set, a vector synthesis module length value, a redundant path elimination record table and a node force value scalar result.

3. The tension evaluation system based on statics theory of claim 2, wherein, The node deformation extraction module comprises: The member coordinate acquisition submodule obtains the connection relationship between the members and nodes in the target structure, including the connection number of the members and nodes, acquires the three-axis coordinate values before and after the force of the member end point, performs coordinate difference calculation on the three-axis coordinate values before and after the force in the global coordinate system, calculates the spatial displacement change amplitude of each member end point in the stressed state, and classifies and arranges the coordinate values according to the member number and node number, and generates an end point spatial displacement information; The attitude change calculation submodule extracts the displacement change direction of each member at the node based on the end point spatial displacement information, performs three-axis angle difference calculation on the attitude angle before the force of the member, calls the node normal direction vector, performs difference calculation on the attitude angle change value of each member after the force and the node normal direction included angle, records the angle change value into an angle change set according to the member number, and generates a member attitude change parameter group. The offset angle generation submodule selects the axial direction vector of each component before being stressed and the node normal direction vector, calculates the included angle, calculates the included angle change value as the component axial offset angle, indexes and matches the offset angle, the three-axis displacement value and the attitude angle change value, integrates the component behavior set of the node level, and generates the node component deformation information group.

4. The tension evaluation system based on statics theory of claim 3, wherein, The component direction reconstruction module comprises: The offset angle judgment submodule calls the node component deformation information group, obtains the component axial offset angle of each component and the number attribute of the corresponding component, calculates the difference between the component axial offset angle and the component axial deformation tolerance, judges whether the offset condition is met based on the size relationship between the difference and zero, records the component number of the component axial offset angle value exceeding the component axial deformation tolerance, generates the offset component number set, and extracts the three-axis displacement value of the corresponding component based on the offset component number set. The force direction replacement calculation submodule constructs an offset vector by taking the three-axis direction displacement value as a component, replaces the offset vector with the original axial direction vector of the component, adds a component number identifier and a node attribution label to each offset vector, establishes the current force transmission direction set of the component, and generates the replacement force direction vector sequence. The component vector collection submodule constructs a multi-level index field according to the replacement force direction vector sequence, and classifies and integrates the offset component vectors according to the node dimension, constructs the direction input set in the node corresponding component set, and generates the node component updated force direction set.

5. The tension evaluation system based on statics theory according to claim 4, characterized in that, The force direction feature collection module comprises: The included angle parameter calculation submodule calls the node component updated force direction set, extracts all component action vectors under each node, classifies the action vectors according to the node number, unit vectorizes any two component action vectors under the same node, performs cosine included angle calculation, records the minimum included angle value, and integrates the included angle value into the node angle set, generates the node included angle feature set, and extracts all component action vectors in the node based on the node included angle feature set. The direction density extraction submodule calculates the included angle between the direction of each vector and the main shaft direction of the structure, judges whether the included angle is lower than the main shaft included angle tolerance value, and calculates the ratio value between the number of component vectors satisfying the included angle tolerance condition and the total incident vector number of the node, generates the node direction density ratio group, and extracts the modulus value of all component force vectors under each node according to the node direction density ratio group. The force value statistical merging submodule groups the modulus values of all component force vectors under each node based on the node number, calculates the mean value and the range of each group of force direction modulus values, integrates the minimum included angle value, the direction density ratio and the force value statistical result of each node into a multi-field array structure, generates the node direction feature parameter cluster, and extracts the modulus value of all component force vectors under each node according to the node direction density ratio group.

6. The tension evaluation system based on statics theory of claim 5, wherein, The specific steps of judging whether the included angle is lower than the main shaft included angle tolerance value are calculating the included angle value between the action vector of each component and the main shaft direction of the structure, and comparing the included angle value with the main shaft included angle tolerance value, wherein the main shaft included angle tolerance value is the sum of the average included angle value and the standard deviation between the action vector of each component under the node and the main shaft direction.

7. The tension evaluation system based on statics theory of claim 6, wherein, The node level division module comprises: The included angle threshold judgment submodule calls the node direction characteristic parameter cluster, extracts the minimum included angle value of each node, compares and judges the included angle value with the component convergence angle threshold item by item, and classifies the node number meeting the condition of being less than the threshold into a candidate set to generate an included angle compliant node set; The direction density screening submodule extracts the direction vector space density value according to the included angle compliant node set, compares and judges the ratio with the set main direction component proportion threshold value, screens the node number with a ratio not lower than the threshold value, and generates a direction density effective node set; The tension strength comparison submodule obtains the incident component force direction module value mean value of the corresponding node according to the direction density effective node set, performs a difference operation with the component tension reference value, and combines the included angle value and the direction density value to obtain the node participation strength index value of each node, and compares the node participation strength index value with the participation level threshold value, classifies the nodes with a node participation strength index value less than the participation level threshold value as edge nodes, and classifies the nodes with an index value greater than or equal to the threshold value as main participation nodes, and establishes a structure node participation identification table.

8. The tension evaluation system based on statics theory of claim 7, wherein, The formula for obtaining the node participation strength index value of each node is specifically: ； wherein, represents a participation intensity index value of the th node, represents a normalized value of the mean value of the incident member tension modulus of the th node, represents a normalized value of the member tension reference value, represents a directional vector space density value of the th node, represents a normalized value of the minimum angle value of the th node, represents a normalized value of the standard deviation of the incident member tension value of the th node, is a directional density decay coefficient, is a base number of a natural logarithm.

9. The tension evaluation system based on statics theory of claim 8, wherein, The tension path evaluation module comprises: The path atlas generation submodule calls the structure node participation identification table and the node component updated force direction set, screens the main participation nodes, extracts the corresponding incident component action vector, groups the incident vectors according to the node number and records the direction attribute, constructs the incident component vector set of each group of nodes into a force direction path set atlas, generates a node path construction graph group; The redundant path elimination submodule obtains the included angle value between any two vectors in the action vector group of each node according to the node path construction graph group, screens the vector pairs with an included angle less than the path redundancy judgment included angle threshold value, calculates the difference between the vector module value ratio and the force direction overlap judgment ratio threshold value, obtains the redundant force direction compression coefficient of the node by operation, marks the node vector pairs with a redundant force direction compression coefficient greater than the redundant judgment compression threshold value as tension redundant path items, and eliminates them from the path set to obtain a redundant path pruning result set; The node resultant force extraction submodule extracts the remaining non-redundant action vector of each node based on the redundant path pruning result set, and performs component addition in a three-dimensional space according to the vector coordinate value to synthesize the total force vector of the node on the action plane, extracts the module value of the total force vector as the tension output index, and establishes a node tension resultant force result set.

10. A tension assessment method based on statics theory, characterized in that, The method is used for realizing the tension evaluation system based on the static force theory in any one of claims 1-9, comprising the following steps: S1: obtaining the connection relationship between the components and the nodes in the target structure, integrating the component axis line offset angle, three-axis displacement value and attitude angle change value of each component at the node to generate a node component deformation information group; S2: calling the node component deformation information group, performing size judgment on the component axis line offset angle and the component axis line deformation tolerance of each component, performing attribution marking on the offset vector and the node coordinates to generate a node component updated force direction set; S3: Call the node component update force vector set, compare the direction angle of the action vector of each node incident component, take the quantity distribution ratio of the component action vector in the direction of the structure main shaft as the direction vector space density value, calculate the mean value and difference value of the incident force vector size, and generate the node direction characteristic parameter cluster; S4: Call the node direction characteristic parameter cluster, classify the nodes that meet the three conditions of small angle, high mean value and high density as main participating nodes, and vice versa as edge nodes, and generate the structure node participation identification table; S5: Call the structure node participation identification table and node component update force vector set, construct a force vector transmission path atlas according to the node division, mark the tension redundant path item, eliminate it from the path, and perform spatial vector addition on the force vector of the non-redundant path, synthesize the total vector of the node force transmission path, and generate the node tension resultant force result set.

Images