A cable protective cover simulation analysis method and system

Abstract

The present invention provides a cable protective cover simulation analysis method and system, which relates to the field of data processing technology. The method includes: determining a cable protective cover model; discretizing the cable protective cover model into multiple units; calculating the stress corresponding to each unit; obtaining the response characteristics of the entire cable protective cover model when subjected to force according to the stress corresponding to each unit and the dynamic equilibrium equation; analyzing the performance of the protective cover during the force-bearing process according to the response characteristics to obtain the performance analysis results; optimizing the protective cover according to the performance analysis results to obtain the final protective cover. The present invention can more accurately simulate the performance of the cable protective cover in an actual working environment.

Description

Technical Field

[0001] The invention relates to the technical field of data processing, and in particular to a cable protective cover simulation analysis method and system. Background Art

[0002] With the development of science and technology, numerical simulation technology has gradually become an important tool for the design and optimization of cable sheaths. Through numerical simulation, the performance of cable sheaths can be accurately predicted in the design stage, thereby greatly improving design efficiency and product quality.

[0003] However, existing numerical simulation methods often only focus on a single physical field (such as stress field) and fail to fully consider the performance of cable sheaths under the coupling of multiple physical fields. Summary of the invention

[0004] The technical problem to be solved by the present invention is to provide a cable protective cover simulation analysis method and system, which can more accurately simulate the performance of the cable protective cover in an actual working environment.

[0005] In order to solve the above technical problems, the technical solution of the present invention is as follows:

[0006] In a first aspect, a cable protective cover simulation analysis method is provided, the method comprising:

[0007] Determine the cable sheath model;

[0008] Discretize the cable sheath model into multiple units;

[0009] Calculate the stress corresponding to each unit;

[0010] According to the stress and dynamic equilibrium equation corresponding to each unit, the response characteristics of the entire cable protection sheath model under stress are obtained;

[0011] According to the response characteristics, the performance of the protective cover during the stress process is analyzed to obtain the performance analysis results;

[0012] According to the performance analysis results, the protective cover is optimized to obtain the final protective cover.

[0013] Furthermore, the stress corresponding to each unit is calculated, including:

[0014] pass Calculate the stress corresponding to each unit;

[0015] Among them, σ represents stress; E represents elastic modulus; ε represents strain; T represents current temperature; T0 represents reference temperature; t represents time; α represents thermal expansion coefficient; β represents creep coefficient; τ represents creep time constant.

[0016] Furthermore, the calculation formula of the kinetic equilibrium equation is:

[0017]

[0018] Among them, K represents the stiffness matrix; u represents the displacement vector of the node; C represents the damping matrix; represents the velocity vector of the node displacement; M represents the mass matrix; The acceleration vector representing the node displacement; F ext represents the externally applied load vector; F int (u,t) represents the internally generated force vector.

[0019] Furthermore, according to the response characteristics, the performance of the protective cover during the stress process is analyzed to obtain performance analysis results, including:

[0020] Through the response characteristics, the stress distribution diagram of the protective cover during the stress process is obtained;

[0021] By comparing the stresses in different parts through the stress distribution diagram, the stress concentration area can be determined;

[0022] According to the stress concentration area, the displacement data of each node in the area is obtained, and the displacement of different nodes in the concentrated area is compared to determine the overall deformation of the protective cover;

[0023] According to the overall deformation, the energy change of the protective cover during the stress process can be obtained;

[0024] According to the energy change, the energy absorption and release are analyzed, and the buffering capacity and impact resistance of the protective cover are determined to obtain the performance analysis results.

[0025] Furthermore, according to the performance analysis results, the protective cover is optimized to obtain the final protective cover, including:

[0026] According to the performance analysis results, the performance shortcomings and risk areas of the protective cover are optimized;

[0027] In the optimization process, the design parameters of the protective cover are taken as particles in the particle swarm, and the position of each particle represents a design scheme of the protective cover;

[0028] According to the design scheme, a protective cover corresponding to a predetermined performance index is searched in the protective cover design space, and the speed and position of each particle are updated according to the performance evaluation result of each particle to obtain a global optimal solution;

[0029] According to the global optimal solution, the final protective cover is obtained.

[0030] Furthermore, the particle velocity update formula is:

[0031] v id (t+1)=w(t)·v id (t)+c1·r1·(p id -x id (t))+c2·r2·(gt d -x id (t))+c3·r3·(l id -x id (t))+α·(x id (t)-x id (t-1));

[0032] Among them, v id (t+1) represents the velocity of particle i at time t+1 in dimension d; w(t) represents the time-varying inertia weight; c1, c2, c3 represent the acceleration coefficients; p id represents the individual historical optimal position of particle i in dimension d; gt d represents the global optimal position of the entire particle swarm in dimension d; l id represents the local optimal position of particle i in dimension d; r1, r2, r3 represent random numbers in the range [0,1]; x id (t) represents the position of particle i at time t in dimension d; α represents the momentum coefficient; x id (t-1) represents the position of particle i at time t-1 in dimension d; i represents the index; t represents time.

[0033] Furthermore, the particle position update formula is:

[0034]

[0035] Among them, x id (t+1) represents the position of particle i at time t+1 in dimension d; x id (t) represents the position of particle i at time t in dimension d; v id (t+1) represents the velocity of particle i at time t+1 in dimension d; β and δ represent adjustment coefficients; x maxd represents the maximum possible position value in dimension d; γ represents a positive coefficient; e represents the base of the natural logarithm; sgn(v id (t))Speed ​​v id The sign function of (t); ∈ denotes an exponent.

[0036] In a second aspect, a cable protective cover simulation analysis system includes:

[0037] An acquisition module is used to determine the cable sheath model; discretize the cable sheath model into multiple units; calculate the stress corresponding to each unit, and combine the equations of each unit to form an overall stiffness matrix;

[0038] The processing module is used to obtain the dynamic equilibrium equation through the stiffness matrix to obtain the response characteristics of the entire cable protective cover model when subjected to force; according to the response characteristics, the performance of the protective cover during the force process is analyzed to obtain the performance analysis results; according to the performance analysis results, the protective cover is optimized to obtain the final protective cover.

[0039] According to a third aspect, a computing device includes:

[0040] one or more processors;

[0041] The storage device is used to store one or more programs. When the one or more programs are executed by the one or more processors, the one or more processors implement the method described.

[0042] In a fourth aspect, a computer-readable storage medium stores a program, and when the program is executed by a processor, the method described is implemented.

[0043] The above solution of the present invention includes at least the following beneficial effects:

[0044] By discretizing the cable sheath model into multiple units and calculating the stress of each unit, this method can provide more accurate stress distribution data; by analyzing the response characteristics of the entire cable sheath model when subjected to force, a comprehensive understanding of the performance of the sheath can be obtained. Based on the performance analysis results, the sheath can be optimized. This optimization can be structural or material, thereby improving the mechanical properties and durability of the sheath. The optimized sheath will be more adaptable to complex actual use environments and improve the safety and reliability of the cable. Through simulation analysis, potential problem areas, such as stress concentration points, can be discovered in the design stage, so that improvements can be made in the design stage to avoid problems such as breakage and damage in actual use. Using simulation analysis methods can predict performance in the product design stage, reduce the need for physical prototype testing, thereby saving development costs and time, and at the same time, the failure rate of products can be reduced through optimization. BRIEF DESCRIPTION OF THE DRAWINGS

[0045] Figure 1 It is a flow chart of a cable protective sheath simulation analysis method provided by an embodiment of the present invention.

[0046] Figure 2 It is a schematic diagram of a cable protective sheath simulation analysis system provided by an embodiment of the present invention.

[0047] Figure 3 It is a front view of a cable protective cover provided by an embodiment of the present invention.

[0048] Figure 4 The present invention provides a cable protective cover with one petal structure.

[0049] Figure 5 The present invention provides a cable protective cover with two petals.

[0050] Figure 6 The invention discloses a circular cross-sectional structure diagram of a cable protective cover provided by an embodiment of the invention.

[0051] Figure 7 The present invention provides an elliptical cross-sectional structure diagram of a cable protective sheath at a bend.

[0052] Figure 8 The present invention provides a cable protective sleeve bending inner diameter stress structure diagram.

[0053] Explanation of the reference numerals: 1. welding point; 2. sheath; 3. coaxial medium; 4. tightening ring; 5. conductor; 6. outer conductor; 7. coaxial outer sheath. DETAILED DESCRIPTION

[0054] The exemplary embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. Although the exemplary embodiments of the present disclosure are shown in the accompanying drawings, it should be understood that the present disclosure can be implemented in various forms and should not be limited by the embodiments set forth herein. On the contrary, these embodiments are provided to enable a more thorough understanding of the present disclosure and to fully convey the scope of the present disclosure to those skilled in the art.

[0055] like Figure 1 As shown, an embodiment of the present invention provides a cable protective cover simulation analysis method, the method comprising the following steps:

[0056] Step 1, determine the cable protective cover model;

[0057] Step 2, discretize the cable protective cover model into multiple units;

[0058] Step 3, calculate the stress corresponding to each unit;

[0059] Step 4, according to the stress corresponding to each unit and the dynamic equilibrium equation, the response characteristics of the entire cable protective cover model under stress are obtained;

[0060] Step 5, analyzing the performance of the protective cover during the stress process according to the response characteristics to obtain a performance analysis result;

[0061] Step 6: Optimize the protective cover according to the performance analysis result to obtain the final protective cover.

[0062] In an embodiment of the present invention, by obtaining a model of a cable protective sleeve, a finite element method is used to discretize a complex cable protective sleeve model into multiple small, simple units. For each small unit, a set of equations describing its mechanical behavior is established according to Hooke's law and the principle of elastic mechanics. These sets of equations are combined to form an overall stiffness matrix. This stiffness matrix describes the mechanical response of the entire cable protective sleeve model when subjected to external forces. Using the stiffness matrix, various response characteristics of the cable protective sleeve when subjected to force, such as displacement, stress distribution, etc., can be solved. These response characteristics are an important basis for evaluating the performance of the protective sleeve. Based on the obtained response characteristics, the performance of the cable protective sleeve during the force process is deeply analyzed. This includes evaluating key performance indicators such as its degree of deformation, stress distribution, stiffness and strength. Based on the performance analysis results, the cable protective sleeve is optimized. This includes adjusting material selection, changing geometric shapes, and strengthening specific areas to improve the overall performance of the protective sleeve.

[0063] By using the finite element method and the calculation of the stiffness matrix, the behavior of the cable protective cover under stress can be simulated more accurately, thereby improving the accuracy of the design. At the same time, this method can also reduce the number of experiments, shorten the design cycle, and improve design efficiency. Through in-depth analysis of the response characteristics, the potential weaknesses and deficiencies of the cable protective cover can be discovered. These deficiencies can be improved through optimized design, thereby improving the overall performance of the protective cover. The optimized design can reduce material waste and reduce production costs. At the same time, by discovering and solving potential problems in the design stage, the risk of product failure in actual use can be reduced.

[0064] In another preferred embodiment of the present invention, the above step 1 may include:

[0065] Collect the design drawings, specifications and related technical parameters of the cable protective cover to determine the geometry, size and material properties of the cable protective cover; use 3D modeling software (such as AutoCAD) to build the geometric model of the cable protective cover according to the design drawings and technical parameters. Define the corresponding physical properties of the cable protective cover, such as elastic modulus, Poisson's ratio, density, thermal expansion coefficient, etc., according to the material of the cable protective cover (such as rubber, plastic, metal, etc.). Compare the established geometric model with the original design drawings to ensure the accuracy of the model, and conduct a preliminary model check to ensure that there are no geometric errors or inconsistencies.

[0066] In another preferred embodiment of the present invention, the above step 2 may include:

[0067] According to the geometry and complexity of the cable sheath, select the meshing in the finite element method (FEM), use meshing software (such as ANSYS Meshing, HyperMesh, etc.) to mesh the cable sheath model, determine the density and distribution of the mesh according to the needs of the analysis to ensure the accuracy and efficiency of the analysis, and refine the mesh in key areas (such as stress concentration areas) to improve the accuracy of the analysis. According to the analysis requirements, select the appropriate unit type (such as linear unit, quadratic unit, etc.) for each mesh unit to ensure that the selected unit type can accurately simulate the behavior of the cable sheath during the actual stress process. Perform a quality check on the divided mesh to ensure that there is no deformed mesh or overly distorted mesh. Use tools such as mesh quality factor to evaluate the quality of the mesh, make necessary adjustments, and export the divided mesh data to a format that can be recognized by the numerical simulation software (such as .msh, .cdb, etc.).

[0068] In a preferred embodiment of the present invention, the above step 3, calculating the stress corresponding to each unit, includes:

[0069] pass Calculate the stress corresponding to each unit;

[0070] Among them, σ represents stress; E represents elastic modulus; ε represents strain; T represents current temperature; T0 represents reference temperature; t represents time; α represents thermal expansion coefficient; β represents creep coefficient; τ represents creep time constant.

[0071] In an embodiment of the present invention, E(ε, T, t)·ε reflects the direct connection between stress and strain. The elastic modulus E, itself changes with the change of strain ε, temperature T and time t, and determines the degree of response of the material to the strain. When the material is deformed by external force, corresponding stress will be generated inside to resist this deformation. α(T-T0)·E reveals the mechanism of temperature influence on stress. The thermal expansion coefficient α indicates the degree of expansion or contraction of the material when the temperature changes. When the temperature of the material changes from the reference temperature T0 to the current temperature T, additional stress will be generated inside the material due to the thermal expansion and contraction effect. This stress is proportional to the elastic modulus E and the temperature change. It reflects the creep behavior of materials under long-term stress. The creep coefficient β describes the characteristic of the material gradually deforming over time under constant stress. The creep time constant τ reflects the rate at which this deformation occurs. Over time, even if the stress remains unchanged, the material will produce additional deformation and stress due to the creep effect.

[0072] The effects of temperature and time can more accurately predict the performance of materials in actual working environments, thereby improving the accuracy and reliability of designs. Parameters can compare the performance of different materials and select the most suitable materials for specific application environments. By accurately predicting the creep behavior of materials under long-term stress, measures can be taken in the design stage to reduce the impact of creep on the structural integrity of the product, thereby extending the service life of the product. Accurately predicting the stress response of materials helps to reduce damage caused by stress concentration during product use, thereby reducing the cost of repair and replacement of parts. Considering the impact of temperature on material stress, early warning and prevention of potential structural failure risks under extreme temperature conditions can be used to enhance product safety.

[0073] In a preferred embodiment of the present invention, the calculation formula of the kinetic equilibrium equation is:

[0074]

[0075] Among them, K represents the stiffness matrix; u represents the displacement vector of the node; C represents the damping matrix; represents the velocity vector of the node displacement; M represents the mass matrix; The acceleration vector representing the node displacement; F ext represents the externally applied load vector; F int (u, t) represents the internally generated force vector.

[0076] In the embodiment of the present invention, the stiffness matrix K reflects the connection stiffness and mutual influence between the nodes of the structure. When the structure is displaced u, the stiffness matrix K determines the restoring force (i.e., elastic force) generated by the deformation of the structure, which restores the structure to its original undeformed state. Therefore, [K]u represents the internal elastic force generated by the deformation of the structure. The damping matrix C describes the energy dissipation characteristics of the structure during dynamic motion. When the structure moves at a speed When vibrating, damping results in a gradual loss of energy, helping to reduce the amplitude of the structure's vibrations. It represents the resistance due to damping, which affects the decay rate of structural vibration. The mass matrix M represents the mass distribution of each node of the structure. During motion, according to Newton's second law, an inertial force will be generated in the opposite direction of acceleration. Represents the inertial force due to acceleration or deceleration of the structure. ext Represents the forces acting on the structure from the outside, which are static (such as gravity) or dynamic (such as earthquake waves, wind, etc.). External loads are the main driving force that causes structural deformation and vibration. int(u, t) represents the forces generated by internal factors of the structure, which depend on the current displacement u of the structure and the time t. Internal forces include forces due to material nonlinearity, thermal expansion, prestress, etc.

[0077] This equation can be used to accurately predict the behavior of a structure under dynamic loads, including vibration frequency, amplitude, and phase. Understanding the dynamic response of a structure helps engineers optimize the design to ensure the safety and comfort of the structure under the expected loads. This equation can be used to evaluate the performance of a structure in extreme events (such as earthquakes and storms) so that corresponding risk mitigation strategies can be developed. By comparing the effects of different materials on the dynamic response of a structure, the most appropriate material can be selected to build the structure. Dynamic analysis of a structure can help identify potential failure points so that preventive measures can be taken to extend the service life of the structure. By optimizing the design and material selection, costs can be reduced while meeting performance requirements.

[0078] In another preferred embodiment of the present invention, the above step 4, obtaining the response characteristics of the entire cable protective sheath model under stress according to the stress and dynamic equilibrium equation corresponding to each unit, may include:

[0079] According to the geometry, material properties and connection mode of each unit, the stiffness matrix [K] of the entire cable sheath model is assembled. This matrix describes the relationship between the forces generated by elastic deformation between the nodes in the model. Similarly, according to the physical properties of the model (such as the damping coefficient), the damping matrix [C] is assembled. This matrix describes the relationship between the damping forces generated by the velocity difference between the nodes in the model. According to the mass distribution of the model, the mass matrix [M] is assembled. This matrix describes the relationship between the inertial forces generated by acceleration at the nodes in the model.

[0080] Determine the externally applied load vector F ext and the internally generated force vector F int (u, t), externally applied load vector F ext Including tension, pressure, gravity, etc., defined according to actual conditions; the internally generated force vector F int (u, t) includes thermal stress due to temperature change, force due to material creep, etc. These forces are related to the node displacement u and time t.

[0081] The stiffness matrix [K], damping matrix [C], mass matrix [M], external load vector F ext and the internal force vector F int Substituting (u, t) into the dynamic equilibrium equation, we can obtain a set of equations describing the stress behavior of the cable protective sheath.

[0082] The Newmark method is used to solve the dynamic equilibrium equations to discretize the continuous time domain into a series of time points, and to solve the approximate solution of the equations at each time point. Through iterative calculation, the node displacement vector u and the velocity vector of the node displacement are gradually solved. The acceleration vector ü of the node displacement changes with time. According to the solution results, the response characteristics of the cable protective cover under stress are analyzed, including maximum displacement, stress distribution, vibration frequency, etc. The performance of the cable protective cover is evaluated, such as whether the strength is sufficient, whether the stiffness is reasonable, whether the stability meets the requirements, etc. If the performance is insufficient or there are potential problems, the model can be optimized.

[0083] Among them, the Newmark method is used to solve the dynamic equilibrium equation to discretize the continuous time domain into a series of time points and solve the approximate solution of the equation system at each time point, including:

[0084] Set the time step Δt, total time T, and parameters γ and β in the Newmark method; initialize the displacement vector u0, velocity vector v0, and acceleration vector a0 (the displacement and velocity at the initial moment are known, and the acceleration can be solved by the dynamic equilibrium equation under the initial conditions).

[0085] For each time point t i =i×Δt (i=0, 1, 2, ..., N), where N=T / Δt; Calculate t according to the prediction formula of the Newmark method i+1 The predicted displacement u′ at time i+1 and predicted velocity v′ i+1 ,in:

[0086] u′ i+1 =u i +Δtv i +0.5Δt 2 ×(1-2β)×a i ;

[0087] v′ i+1 =v i +Δt×(1-γ)×a i ;

[0088] According to the dynamic equilibrium equation and the formula of Newmark method, the equivalent stiffness matrix [K eq ] and the equivalent load vector F eq :

[0089] [K e q]=[K]+(1 / βΔt 2 )[M]+(γ / βΔt)[C];

[0090] Fe q=F ext (t i+1 )+F int (u′ i+1 , t i+1 )+[M](a i ×(1 / β-1)+Δt / βv i +Δt 2 / 2βu i )+[C](Δt / βv i +Δt 2 / 2βa i );

[0091] Among them, F ext (t i+1 ) and F int (u′ i+1 , t i+1 ) indicates that at t i+1 The external load vector and internal force vector at the moment. Solve the linear equations [K eq ]Δu=F eq , and obtain the displacement increment Δu.

[0092] According to the displacement increment Δu, update t i+1 The displacement at time u i+1 , speed v i+1 and acceleration a i+1 ,in:

[0093] u i+1 =u′ i+1 +Δu;

[0094] v i+1 =v′ i+1 +yΔta i+1 ;

[0095] a i+1 =(F ext (t i+1 )+F int (u i+1 , t i+1 )-[K]u i+1 -[C]v i+1 ) / [M];

[0096] At each time point or a specific output time point, output the displacement, velocity, acceleration and other results, and end the loop when the time t reaches the total time T.

[0097] In a preferred embodiment of the present invention, the above step 5, analyzing the performance of the protective cover during the force-bearing process according to the response characteristics to obtain the performance analysis result, includes:

[0098] Step 51, obtaining a stress distribution diagram of the protective cover during the stress process through the response characteristics;

[0099] Step 52, comparing stresses at different locations through the stress distribution diagram to determine areas of stress concentration;

[0100] Step 53, according to the stress concentration area, obtain the displacement data of each node in the area, and compare the displacement of different nodes in the concentrated area to determine the overall deformation of the protective cover;

[0101] Step 54, obtaining the energy change of the protective cover during the force-bearing process according to the overall deformation;

[0102] Step 55, analyzing the energy absorption and release according to the energy change, determining the buffering capacity and impact resistance of the protective cover, and obtaining the performance analysis result.

[0103] In an embodiment of the present invention, the response characteristics include stress, strain, and displacement, etc., and the stress data is imported into the visualization software. According to the size of the stress data, different colors or contour maps are used on the protective cover model to represent the stress distribution. The generated stress distribution map can clearly show which parts of the protective cover are subjected to higher stress during the force process. By observing the stress distribution map, it can be clearly seen that the stress values ​​in some areas are much higher than those in the surrounding areas. These high-stress areas are stress concentration areas, which are potential damage points in practical applications. Record the location and range of these areas, and extract the displacement data of the nodes in these stress concentration areas in the finite element analysis. Compare the displacements of different nodes, and pay attention to local deformation, such as abnormal movement of a node relative to other nodes. The protective cover will undergo energy absorption and release processes when external forces are applied and unloaded. By continuously monitoring the energy changes in this process, a curve of energy changes over time is drawn. This curve helps to understand how the protective cover effectively absorbs and releases energy at different force stages. Combined with deformation analysis and energy change trends, the buffering capacity and impact resistance of the protective cover can be comprehensively evaluated. By analyzing stress distribution and deformation, structural weaknesses can be discovered and improved in a timely manner, thereby improving the safety of the protective cover in actual use. Based on the performance analysis results, designers can optimize the structure of the protective cover to improve its mechanical properties and durability. Through accurate performance analysis, over-design can be avoided, thereby reducing material and manufacturing costs while meeting performance requirements.

[0104] When applied in practice, step 51, in this step, a test device or software is used to simulate or actually measure the stress distribution of the protective cover when it is subjected to force. This involves applying a certain force to the protective cover and capturing and recording stress data at each point through a sensor. These data are then organized into a graphical form, i.e., a stress distribution diagram, which can intuitively show the stress magnitude of the protective cover at different locations.

[0105] Step 52, after obtaining the stress distribution diagram, carefully observe the stress values ​​of each part in the diagram. By comparing the stress magnitudes of different parts, it can be clearly seen which areas have higher stress values, and these areas are where the stress is concentrated. These areas are usually the places where the protective cover is most likely to be damaged or fail during the stress process, so special attention should be paid.

[0106] Step 53, after determining the stress concentration area, it is necessary to further obtain the displacement data of each node in these areas. This can be obtained by experimental measurement or numerical simulation. By comparing the displacement of different nodes, the overall deformation of the protective cover during the stress process can be understood. For example, the displacement of some nodes may be large, indicating that these parts have undergone large deformation.

[0107] Step 54, by analyzing the overall deformation of the protective cover, the energy change inside the protective cover during the force application process can be inferred. For example, the area with greater deformation may absorb more energy. This analysis of energy change helps to understand the impact resistance and buffering effect of the protective cover.

[0108] Step 55, the last step is to evaluate the cushioning capacity and impact resistance of the protective cover based on the previous analysis results. If the protective cover can effectively absorb and disperse energy during the force process, then its cushioning capacity and impact resistance are good. This can be quantitatively evaluated by comparing the energy absorption and release of different protective covers under the same force conditions. In this way, a comprehensive analysis result of the performance of the protective cover can be obtained.

[0109] In a preferred embodiment of the present invention, the above step 6, optimizing the protective cover according to the performance analysis result to obtain the final protective cover, includes:

[0110] Step 61, optimizing the performance shortcomings and risk areas of the protective cover according to the performance analysis results;

[0111] Step 62, during the optimization process, the design parameters of the protective cover are represented by particles in the particle group, and the position of each particle represents a design scheme of the protective cover;

[0112] Step 63, according to the design scheme, searching for a protective cover corresponding to a predetermined performance index in the protective cover design space, and updating the speed and position of each particle according to the performance evaluation result of each particle to obtain a global optimal solution;

[0113] Step 64, obtaining the final protective cover according to the global optimal solution.

[0114] In an embodiment of the present invention, by analyzing the performance of the protective cover in simulation or actual testing, its performance shortcomings and risk areas are identified. These shortcomings include problems such as stress concentration, excessive deformation or insufficient energy absorption. In the optimization process, the design parameters of the protective cover (such as material thickness, structural shape, rib position, etc.) are regarded as particles in the particle swarm. The position of each particle represents a protective cover design scheme. A certain number of particles are randomly generated as an initial design scheme set. According to predetermined performance indicators (such as maximum stress, deformation range, energy absorption capacity, etc.), the design scheme of each particle is evaluated. The speed and position of the particles are updated, and particles with better performance will obtain a greater speed, thereby performing a wider search in the design space. At the same time, information will be exchanged between particles to accelerate the search for the global optimal solution. By continuously iterating the above process, the particle swarm will gradually converge to the vicinity of the global optimal solution. This optimal solution represents the best protective cover design scheme that meets the predetermined performance indicators. According to the design parameters corresponding to the global optimal solution, the final protective cover design scheme is obtained.

[0115] By optimizing the performance shortcomings and risk areas, the cushioning and impact resistance of the protective cover can be improved, thereby more effectively protecting the wrapped objects or people. The particle swarm optimization algorithm can automatically search for the optimal solution in the design space, avoiding the tedious and time-consuming traditional trial and error method and improving the design efficiency. This method can adapt to different design requirements and performance indicators, and provide flexible and efficient solutions for different types of protective cover designs. By optimizing the design, it is possible to reduce material usage or simplify the structure while meeting performance requirements, thereby reducing production costs.

[0116] In a preferred embodiment of the present invention, the particle velocity update formula is:

[0117] v id (t+1)=w(t)·v id (t)+c1·r1·(p id -x id (t))+c2·r2·(gt d -x id (t))+c3·r3·(l id -x id (t))+α·(x id (t)-x id (t-1));

[0118] Among them, v id (t+1) represents the velocity of particle i at time t+1 in dimension d; w(t) represents the time-varying inertia weight; c1, c2, c3 represent the acceleration coefficients; p id represents the individual historical optimal position of particle i in dimension d; gt d represents the global optimal position of the entire particle swarm in dimension d; l id represents the local optimal position of particle i in dimension d; r1, r2, r3 represent random numbers in the range [0, 1]; x id (t) represents the position of particle i at time t in dimension d; α represents the momentum coefficient; x id (t-1) represents the position of particle i at time t-1 in dimension d; i represents the index; t represents time.

[0119] In the embodiment of the present invention, w(t)·v id (t) represents the tendency of the particle to maintain its current velocity, where w(t) is the time-varying inertia weight, which determines how much the particle maintains its current velocity. Larger values ​​of w(t) mean that the particle is more inclined to maintain its current velocity, while smaller values ​​make it easier for the particle to change the direction of its velocity. id -x id (t)) reflects the particle's memory of its best historical position. id is the individual best historical position of particle i in dimension d. The particle will move towards its best historical position in the hope of finding a better solution. c1 is the acceleration coefficient, which controls the influence, and r1 is a random number used to increase the diversity of the search. c2·r2·(gt d -x id (t)) reflects the information sharing between particles. gt d represents the global best position of the entire particle swarm in dimension d, that is, the best solution found among all particles. Particles will be attracted by the global best position to accelerate convergence to the global optimal solution. c2 is the acceleration coefficient that controls this part of the influence, and r2 is a random number. c3·r3·(l id -x id (t)) allows the particle to move to the optimal position l among its neighbors or a local subgroup id Learning. It helps to increase the diversity of the search and may help avoid premature convergence to the local area of ​​the global optimum. c3 is the corresponding acceleration coefficient and r3 is a random number. id represents the best position of particle i in dimension d within its neighborhood. Allowing particles to be influenced by the best position in their neighborhood helps to perform a more detailed search in the search space. id (t)-x id(t-1)) takes into account the rate of change of the particle's position, i.e. momentum. It helps the particle maintain a dynamic tendency to continue moving in the direction it was previously moving. The momentum coefficient α regulates the strength of this tendency.

[0120] By adjusting the inertia weight w(t), acceleration coefficients c1, c2, c3, and momentum coefficient α, a balance can be found between global search and local search. This helps to strike a balance between exploring new areas and using known information, thereby increasing the probability of finding the global optimal solution. The introduction of random numbers r1, r2, and r3 increases the diversity of the search process and helps avoid falling into the local optimum too early. At the same time, by using the individual optimal p id , global optimal gt d and the local optimum l id The time-varying inertia weight w(t) can be adaptively adjusted according to the algorithm's iterative process, so that the algorithm focuses more on global search in the early stage and more on local fine search in the later stage. At the same time, it can flexibly respond to different optimization problems and scenarios.

[0121] In a preferred embodiment of the present invention, the particle position update formula is:

[0122]

[0123] Among them, x id (t+1) represents the position of particle i at time t+1 in dimension d; x id (t) represents the position of particle i at time t in dimension d; v id (t+1) represents the velocity of particle i at time t+1 in dimension d; β and δ represent adjustment coefficients; represents the maximum possible position value in dimension d; γ represents a positive coefficient; e represents the base of the natural logarithm; sgn(v id (t))Speed ​​v id The sign function of (t); ∈ denotes an exponent.

[0124] In the embodiment of the present invention, the particle is updated according to its velocity v id (t+1) from the current position x id (t) Move to a new location. This part represents the distance between the particle and the boundary of the search space, and is adjusted accordingly. , the value of this item increases, which helps push particles to the interior of the search space and prevents particles from getting stuck on the boundary too early. is a function that increases with time, which means that the influence of boundary effects may gradually increase with time. id (t))·|vid (t)| ∈ This part represents a decay term related to the current velocity of the particle. id (t))·|v id (t)| ∈ is the speed v id (t) is a sign function used to maintain the direction of velocity. id (t)| ∈ The magnitude of the velocity is taken into account and adjusted by the exponent ∈. When ∈ is greater than 1, particles moving at high speeds experience greater attenuation, while when ∈ is less than 1, particles moving at low speeds experience greater attenuation. The parameter δ controls the overall strength of the velocity attenuation.

[0125] By v id (t+1) term, particles can move according to their updated speed, which helps particles to explore extensively in the search space, thereby enhancing the global search capability of the algorithm. Taking into account the distance between the particle and the upper bound of the search space helps to avoid particles from gathering near the boundary too early, thereby maintaining the diversity of the solution. At the same time, the function that changes over time The influence of the boundary effect can be adjusted dynamically, so that different emphases can be placed at different stages of the search. id (t))·|v id (t)| ∈ The speed decay mechanism helps control the movement speed of particles and prevents particles from passing through possible optimal areas too quickly. This mechanism is particularly important in the later stages of the search because it helps the algorithm perform more refined local searches. Multiple parameters in the formula (such as β, δ, ∈, and γ) can be adjusted according to specific problems and are highly adaptive. It helps prevent the algorithm from converging to a local optimal solution too early. By comprehensively considering the particle's speed, position, and relationship to the search boundary, a balance can be achieved between global exploration and local development, thereby reducing the risk of premature convergence.

[0126] like Figure 3 and Figure 4 As shown, when applied specifically, the cable protective sheath specifically includes a conductor 5, a coaxial medium 3 is arranged on the outside of the conductor 5, an outer conductor 6 is arranged on the outside of the coaxial medium 3, a coaxial outer sheath 7 is arranged on the outside of the outer conductor 6, a sheath 2 is arranged on the outside of the coaxial outer sheath 7, a plurality of grooves are arranged on the sheath 2, and a tightening ring 4 is arranged in the groove.

[0127] The outer diameter of the sheath 2 has multiple grooves for tools (binding ropes, metal hoops, etc., in any way) to wrap and fix. The groove spacing is inconsistent. From the welding position to the farther away, the groove spacing decreases. The inner diameter of the sheath has a protruding structure corresponding to the outer diameter groove position. After being assembled and fixed with the cable, the protruding structure of the protective sheath contacts the outer surface of the cable to form a bite connection. The area between the two adjacent protruding structures will leave a certain gap on the cable contact surface. The inner diameter gap area of ​​the protective sheath does not contact the outer surface of the cable. When the cable is bent under force, the inner curved surface (compression deformation) and the outer curved surface (tensile deformation) provide a release space for the deformation of the inner curved surface and the outer curved surface of the cable, reducing axial stretching.

[0128] This structure makes the rigidity of the protective sleeve the greatest at the welding point 1, and the rigidity decreases as it moves farther away (controlled by the unequal spacing of the inner diameter protrusion structure, with a small spacing at the welding point and a large spacing as it moves farther away). When the cable assembly is subjected to external bending stress, the bending radius becomes smaller as the welding point moves farther away, and the bending radius becomes larger as it moves closer to the welding point, effectively preventing the welding point from being damaged by direct force, thereby improving the reliability of the cable.

[0129] The sheath 2 can also be a closed integral structure, or an integral structure with a slit, or a two-petal or multi-petal split structure. When in use, the protective sleeve is assembled onto the cable, and then a tool (binding rope, metal hoop, etc., the method is not limited) is used to wrap and fix it in the groove of the outer diameter of the protective sleeve, and the non-metallic material is used for protection, and the rigid joint part is gradually transitioned to the soft cable. The sheath is made of plastic, which has a certain rigidity and strong weather resistance. Multiple gaps and micro-shrinkage fastening are used alternately. The local contraction of the sheath is used to compress the cable skin to shrink, and the sheath and the cable are "tightened" to achieve axial fastening protection. In addition, the gap can be used to release the compression deformation of the inner fillet material of the outer skin when the cable is bent. The spacing can be gradual, combined with the wall thickness of the sheath to achieve a gradual change from rigid to soft to protect the root of the interconnection.

[0130] like Figures 6 to 8 As shown, for RF coaxial cable, the characteristic impedance Where Z0 represents characteristic impedance; ε r = represents the dielectric constant; D represents the inner diameter of the outer conductor; d represents the outer diameter of the inner conductor. When the cable is bent, the inner circle of the bend is squeezed and the outer circle is stretched. The bend will become elliptical (the coaxial circle interface at the unbent part gradually becomes elliptical at the bend and then returns to a circle). The characteristic impedance will also change accordingly. The smaller the bending radius, the more obvious the change. At the same time, during the cable bending process, the stress on the coaxial cable medium will cause the dielectric constant ε rThe greater the stress, the greater the change, the characteristic impedance Z0 will also change accordingly, and the smaller the bending radius, the more obvious the change. The change of the cable characteristic impedance means that the impedance of the cable is discontinuous, and signal reflection will occur at the point where the characteristic impedance changes, increasing the transmission loss of the signal and affecting the port voltage standing wave ratio index and insertion loss index of the cable assembly.

[0131] When the coaxial cable itself receives axial tensile (compressive) force, the outer sheath is required to be deformed first and then transmitted to the cable. Only when the sheath is deformed too much, the force is transmitted to the metal wire braid layer (and / or metal tape winding layer), causing the metal wire braid layer (and / or metal tape winding layer) to shrink (expand), the inner diameter becomes smaller (or larger), and then causes the coaxial cable dielectric layer to be thinner (larger) (the inner conductor is generally difficult to change), thus reducing the change in the ratio of the inner and outer diameters of the dielectric, reducing the change in the characteristic impedance of the coaxial cable, and reducing the discontinuity of the characteristic impedance. Therefore, it is necessary to fix the limit of the outer sheath of the cable.

[0132] Therefore, there is a minimum bending radius limit for coaxial cables when used, otherwise the electrical performance will deteriorate sharply and cannot be used normally. The higher the frequency of the transmitted signal, the more obvious the impact, especially in the millimeter wave band.

[0133] Therefore, after the coaxial cable connector is interconnected with the flexible cable, the sheath 2 is partially in contact with the micro-surface of the coaxial cable, and a tightening ring 4 is added to the outer side of the sheath corresponding to the contact point, so that the cable sheath at the contact position is slightly contracted and interlocked with each other (avoiding the influence on the internal cross-sectional shape and size of the coaxial cable and affecting the electrical performance), and the cable sheath can be fixed. The existence of the elastic sheath can not only limit the force on the welding part during axial stretching, but also limit the bending of the cable. At the same time, the sheath is gradually thinner, the limiting force is gradually reduced, and the allowed bending is gradually increased to the minimum bending radius allowed by the cable itself. The sheath 2 is fixed in contact with the local micro-surface of the coaxial cable, and a micro-gap is left in the non-contact part, allowing the cable sheath and the internal medium to produce a slight deformation when the cable is bent, releasing the bending deformation stress, and the micro-gap limits the degree of bending deformation of the cable, playing a role in preserving the bending shape of the coaxial cable (reducing the ellipticity of the cross section at the bending part), while reducing the axial stretching and compressive stress during bending, and protecting the cable itself and the welding part.

[0134] By adjusting the micro-contact surface area and force, as well as the micro-contact surface length and gap size, the comprehensive optimization of welding point (interconnection point, which can be crimped and other fixing methods) protection, bending protection and microwave performance can be achieved. Close to the welding point, the sheath thickness is thick, the limiting force is large, and the allowed deformation is minimal, which can fully reduce the deformation stress (bending and stretching) of the welding point.

[0135] like Figure 2 As shown, an embodiment of the present invention further provides a cable protective sheath simulation analysis system 20, comprising:

[0136] The acquisition module 21 is used to determine the cable protective sheath model; discretize the cable protective sheath model into multiple units; calculate the stress corresponding to each unit, and combine the equations of each unit to form an overall stiffness matrix;

[0137] The processing module 22 is used to obtain the dynamic equilibrium equation through the stiffness matrix to obtain the response characteristics of the entire cable protective cover model when subjected to force; according to the response characteristics, the performance of the protective cover during the force process is analyzed to obtain the performance analysis results; according to the performance analysis results, the protective cover is optimized to obtain the final protective cover.

[0138] It should be noted that the system is a system corresponding to the above method, and all implementation methods in the above method embodiment are applicable to this embodiment and can achieve the same technical effect.

[0139] The embodiment of the present invention further provides a computing device, comprising: a processor, a memory storing a computer program, wherein when the computer program is executed by the processor, the method described above is executed. All implementations in the above method embodiment are applicable to this embodiment and can achieve the same technical effect.

[0140] The embodiment of the present invention also provides a computer-readable storage medium storing instructions, which, when executed on a computer, enable the computer to execute the method described above. All implementations in the above method embodiment are applicable to this embodiment and can achieve the same technical effect.

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

[0142] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working processes of the systems, devices and units described above can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.

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

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

[0145] In addition, each functional unit in each embodiment of the present invention may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit.

[0146] If the functions are implemented in the form of software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, or the part that contributes to the prior art, or the part of the technical solution, can be embodied in the form of a software product. The computer software product is stored in a storage medium, including several instructions for a computer device (which can be a personal computer, server, or network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes: various media that can store program codes, such as USB flash drives, mobile hard disks, ROM, RAM, magnetic disks, or optical disks.

[0147] In addition, it should be noted that in the apparatus and method of the present invention, it is obvious that each component or each step can be decomposed and / or recombined. These decompositions and / or recombinations should be regarded as equivalent schemes of the present invention. Moreover, the steps of performing the above-mentioned series of processing can naturally be performed in chronological order according to the order of description, but it is not necessary to perform them in chronological order, and some steps can be performed in parallel or independently of each other. For those of ordinary skill in the art, it is understood that all or any steps or components of the method and apparatus of the present invention can be implemented in any computing device (including processors, storage media, etc.) or a network of computing devices in hardware, firmware, software or a combination thereof, which can be achieved by those of ordinary skill in the art using basic programming skills after reading the description of the present invention.

[0148] Therefore, the purpose of the present invention can also be achieved by running a program or a group of programs on any computing device. The computing device can be a well-known general device. Therefore, the purpose of the present invention can also be achieved by simply providing a program product containing a program code that implements the method or device. That is to say, such a program product also constitutes the present invention, and the storage medium storing such a program product also constitutes the present invention. Obviously, the storage medium can be any well-known storage medium or any storage medium developed in the future. It should also be pointed out that in the device and method of the present invention, it is obvious that each component or each step can be decomposed and / or recombined. These decompositions and / or recombinations should be regarded as equivalent schemes of the present invention. In addition, the steps of performing the above-mentioned series of processing can naturally be performed in chronological order according to the order of description, but it is not necessary to perform them in chronological order. Some steps can be performed in parallel or independently of each other.

[0149] The above is a preferred embodiment of the present invention. It should be pointed out that for ordinary technicians in this technical field, several improvements and modifications can be made without departing from the principles of the present invention. These improvements and modifications should also be regarded as the scope of protection of the present invention.

Claims

1. A cable protective sheath simulation analysis method, characterized in that: The method comprises: Determine the cable sheath model; Discretize the cable sheath model into multiple units; Calculate the stress corresponding to each unit; According to the stress and dynamic equilibrium equation corresponding to each unit, the response characteristics of the entire cable protection sheath model under stress are obtained; According to the response characteristics, the performance of the protective cover during the stress process is analyzed to obtain the performance analysis results; According to the performance analysis results, the protective cover is optimized to obtain the final protective cover; Calculate the stress corresponding to each element separately, including: pass Calculate the stress corresponding to each unit; in, Indicates stress; represents the elastic modulus; Indicates strain; Indicates the current temperature; Indicates the reference temperature; Indicates time; represents the coefficient of thermal expansion; represents the creep coefficient; represents the creep time constant; The calculation formula of the kinetic equilibrium equation is: ； in, represents the stiffness matrix; represents the displacement vector of the node; represents the damping matrix; The velocity vector representing the node displacement; represents the mass matrix; The acceleration vector representing the node displacement; represents the externally applied load vector; represents the internally generated force vector; According to the response characteristics, the performance of the protective cover during the stress process is analyzed to obtain the performance analysis results, including: Through the response characteristics, the stress distribution diagram of the protective cover during the stress process is obtained; By comparing the stresses in different parts through the stress distribution diagram, the stress concentration area can be determined; According to the stress concentration area, the displacement data of each node in the area is obtained, and the displacement of different nodes in the concentrated area is compared to determine the overall deformation of the protective cover; According to the overall deformation, the energy change of the protective cover during the stress process can be obtained; According to the energy changes, analyze the energy absorption and release, determine the cushioning capacity and impact resistance of the protective cover, and obtain the performance analysis results; According to the performance analysis results, the protective cover is optimized to obtain the final protective cover, including: According to the performance analysis results, the performance shortcomings and risk areas of the protective cover are optimized; In the optimization process, the design parameters of the protective cover are taken as particles in the particle swarm, and the position of each particle represents a design scheme of the protective cover; According to the design scheme, a protective cover corresponding to a predetermined performance index is searched in the protective cover design space, and the speed and position of each particle are updated according to the performance evaluation result of each particle to obtain a global optimal solution; According to the global optimal solution, the final protective cover is obtained; The particle velocity update formula is: ； in, Represents particles Time in dimension d speed at 1000 s; represents the time-varying inertia weight; , , represents the acceleration factor; Represents particles The individual's historical optimal position on dimension d; represents the global optimal position of the entire particle swarm in dimension d; Represents particles The local optimal position in dimension d; , , Represents a random number in the range [0,1]; Represents particles The position at time t in dimension d; represents the momentum coefficient; Represents particles In Dimension Time location; represents index; t represents time; The particle position update formula is: ； in, Represents particles The position at time t+1 in dimension d; Represents particles Position at time t in dimension d; Represents particles The velocity at time t+1 in dimension d; Adjustment factor; Represents the maximum possible position value in dimension d; represents a positive coefficient; represents the base of natural logarithms; speed The symbolic function of An index.

2. A cable protective sheath simulation analysis system, characterized in that: The system is used to execute the method according to claim 1, comprising: An acquisition module is used to determine the cable sheath model; discretize the cable sheath model into multiple units; calculate the stress corresponding to each unit, and combine the equations of each unit to form an overall stiffness matrix; The processing module is used to obtain the dynamic equilibrium equation through the stiffness matrix to obtain the response characteristics of the entire cable protective cover model when subjected to force; according to the response characteristics, the performance of the protective cover during the force process is analyzed to obtain the performance analysis results; according to the performance analysis results, the protective cover is optimized to obtain the final protective cover.

3. A computing device, characterized in that: include: one or more processors; A storage device, used for storing one or more programs, when the one or more programs are executed by the one or more processors, the one or more processors implement the method as claimed in claim 1.

4. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a program, which implements the method according to claim 1 when executed by a processor.

Images