A force density-based adaptive morphing method for long-span suspension antenna

Abstract

The application provides a large-span suspension antenna adaptive shape adjustment method based on force density, which comprises the following steps: given material parameters, cross-sectional dimensions and antenna structure facade parameters of antenna suspension cables and down lead; grid division is performed on the main cable structure according to a given grid size or division number, a node coordinate matrix and a node connection relationship matrix are constructed; a self-weight load vector corresponding to the node number is constructed, and a concentrated force item corresponding to an insulator chain is additionally added; the part of the main cable structure not including the down lead is constructed into a final load vector; grid division is performed on the down lead to form unit and node information and a load vector; initial values of horizontal tension of all cable sections of the antenna are given to construct a force density vector; all node force balance equations are constructed to obtain updated positions of the nodes; and it is judged whether the antenna satisfies the design requirement of the antenna shape. The whole process of the embodiment of the application does not need manual intervention, the design efficiency is improved under the premise of ensuring accuracy.

Description

Technical Field

[0001] This invention relates to the field of suspension antenna morphology design technology, and in particular to an adaptive morphology adjustment method for a long-span suspension antenna based on force density. Background Technology

[0002] Long-wave antennas play a crucial role in underwater communications, and long-span suspension cable (wire) structures, acting as radiators, are a typical type of long-wave antenna. They operate in the frequency range of 3–300 kHz, with wavelengths in vacuum ranging from 1000 to 100,000 meters. These antennas are typically deployed in plains or valleys, and their dimensions are generally several kilometers in size. The feed end is usually located at the lower end of the suspension cable, connected to the middle of the cable by a lower lead. Considering the cable sag under load, the lower lead is usually adjusted by a regulating rope to prevent it from touching the ground. The design of a large suspension antenna requires finding a suitable antenna shape that meets the tension requirements of the suspension end and the lower lead adjustment end, while also satisfying structural strength requirements and providing support for system construction. Shape design is an essential and crucial aspect of suspension antenna design and construction.

[0003] Traditional suspension antenna design methods typically employ finite element analysis (FEM), which involves manually adjusting the virtual thermal expansion coefficient of the structure step-by-step to ensure the boundary conditions meet design requirements. The internal forces and morphology of the current structure are then extracted as the final design result. The drawback of this method is that it requires manual intervention in the parameter adjustment process, resulting in low efficiency and sometimes leading to design failures. Summary of the Invention

[0004] The technical problem to be solved by this invention is how to effectively improve the efficiency of the design of large suspension antennas and avoid the low design efficiency caused by traditional methods that rely on manual experience and trial calculations. In view of this, this invention provides an adaptive shape adjustment method for large-span suspension antennas based on force density.

[0005] The technical solution adopted in this invention is an adaptive morphology adjustment method for a long-span suspension antenna based on force density, comprising:

[0006] Step 1: Given the material parameters, cross-sectional dimensions, and antenna structure elevation parameters of the antenna suspension cable and lower lead wire;

[0007] Step 2: According to the given mesh size or number of segments, divide the main cable structure into a mesh, and construct the node coordinate matrix and the node connection matrix;

[0008] Step 3: Based on the main cable structure, construct the self-weight load vector corresponding to the node number according to the principle of equal distribution of self-weight, and add the concentrated force term corresponding to the insulator chain;

[0009] Step 4: For the portion of the main cable structure excluding the lower lead wire, construct the final load vector based on the horizontal tension using a single variable method.

[0010] Step 5: Construct a full antenna structure for the lower lead portion of the main cable structure, and perform meshing on the lower lead to form cell and node information and load vector;

[0011] Step 6: Based on the design tension requirements of the suspension cable end and the lower lead adjustment rope end, give the initial value of the horizontal tension of all cable segments of the antenna and construct the force density vector;

[0012] Step 7: Construct the force equilibrium equations for all nodes and solve them to obtain the updated positions of the nodes;

[0013] Step 8: Determine whether the error between the positions of all antenna nodes and the positions calculated in the previous calculation is less than the set value. If it is, output the antenna shape that meets the design requirements and end the calculation; otherwise, update the force density vector and load vector according to the target tension requirement and return to step 7.

[0014] In one implementation, step 1 includes:

[0015] Step 101: Given the spatial coordinates of the two suspension ends, the feed end of the lower lead, and the balance end of the lower lead of the main cable of the suspension antenna; initially, give the coordinates of the bifurcation point of the lower lead;

[0016] Step 102: Given the material parameters of all cables involved in the main cable and the lower lead, including the density vector, elastic modulus vector, and breaking tensile force vector; given the cross-sectional parameters of all cables; wherein, the number of elements in the above four types of vectors depends on the number of cable segment types that constitute the actual antenna;

[0017] Step 103: Given the weight and position of each insulator chain in the main cable and the lower lead wire, it is reflected in the analysis model as a concentrated load;

[0018] Step 104: Based on the actual counterweight, set the target tension for the traction force at the suspension cable traction end and the lower lead wire balance end respectively.

[0019] In one implementation, step 3 includes;

[0020] Step 301: Divide the self-weight load of each unit of the main cable equally into the front and rear nodes as nodal loads;

[0021] Step 302: Superimpose the self-weight loads assigned to each node;

[0022] Step 303: Add the weight of the insulator chain to the corresponding nodes to construct the final load vector.

[0023] In one implementation, step 4 includes:

[0024] Step 401: Determine the range of horizontal tension based on the required tension at the suspension end;

[0025] Step 402: Calculate the antenna shape corresponding to the lower limit of the horizontal tension to obtain the corresponding suspension end tension;

[0026] Step 403: Calculate the antenna shape corresponding to the lower limit of the horizontal tension to obtain the corresponding suspension end tension;

[0027] Step 404: Using the golden section method commonly used in univariate optimization, update the lower bound of the horizontal tension until the error between the tension at the suspension point and the target tension is less than the given tolerance. The suspension cable shape at this point is the final shape of the current step.

[0028] In one embodiment, step 4 further includes:

[0029] Step 405(a): For the current horizontal tension, write the nodal equilibrium equation for a node based on the force density method;

[0030] Step 405(b): Construct the balance equations for all nodes;

[0031] Step 405(c): Solve the equilibrium equations to obtain the new node y-coordinates. At this point, the position of each node is uniquely determined, and the node coordinate matrix is ​​updated.

[0032] Step 405(d): Determine whether the error between the position of all antenna nodes and the position calculated in the previous step is less than the set value. If the error is satisfied, the calculation of the current horizontal tension ends and the tension at the current suspension end is obtained. Otherwise, update the force density vector and load vector according to the target tension requirement and return to step 405(a).

[0033] In one implementation, step 5 includes:

[0034] Step 501: Determine the initial position of the entire antenna structure based on the current suspension cable configuration and the location of the lower lead bifurcation point;

[0035] Step 502: Based on the given unit size or number of segments, divide the lower lead wire and the lower lead wire adjustment rope into units to form unit and node information.

[0036] Step 503: Based on the unit and node information, determine the node coordinate matrix and node connection matrix of the overall structure.

[0037] In one implementation, step 6 includes:

[0038] Step 601: Based on the requirements and current shape of the suspension end and the lower lead wire balance end, determine the horizontal tension of the suspension end and the lower lead wire adjustment end, and then calculate the remaining horizontal tension values ​​based on the force balance of the lower lead wire bifurcation point and the suspension point.

[0039] Step 602: Calculate the force density vector of all cable segments based on all horizontal tensions.

[0040] In one embodiment, the displacement sensor and the shape memory alloy spring are configured in a non-contact manner.

[0041] Another aspect of the present invention provides an electronic device comprising: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the computer program, when executed by the processor, implements the steps of the force density-based adaptive morphology adjustment method for a long-span suspension antenna as described in any of the preceding claims.

[0042] Another aspect of the present invention provides a computer storage medium storing a computer program that, when executed by a processor, implements the steps of the adaptive morphology adjustment method for a long-span suspension antenna based on force density as described in any of the preceding claims.

[0043] This invention proposes an adaptive morphology adjustment method for long-span suspension antennas based on force density, requiring no manual intervention throughout the entire process. Compared to existing methods relying on experience and simulation software calculations, this method significantly improves design efficiency while maintaining accuracy. Attached Figure Description

[0044] Figure 1 This is a schematic flowchart of an adaptive morphology adjustment method for a long-span suspension antenna based on force density according to an embodiment of the present invention.

[0045] Figure 2 This is a schematic diagram of a suspension antenna structure according to an embodiment of the present invention;

[0046] Figure 3 This is a flowchart of antenna morphology analysis under a given horizontal tension according to an embodiment of the present invention;

[0047] Figure 4 A morphological diagram obtained from a specific implementation example according to an embodiment of the present invention;

[0048] Figure 5 The deformation cloud diagram is a finite element analysis result of the morphological design according to an embodiment of the present invention.

[0049] Figure 6 This is a schematic diagram of an electronic device structure according to an embodiment of the present invention. Detailed Implementation

[0050] To further illustrate the technical means and effects of the present invention in achieving its intended purpose, the present invention will be described in detail below with reference to the accompanying drawings and preferred embodiments.

[0051] In the accompanying drawings, the thickness, size, and shape of the objects have been slightly exaggerated for ease of illustration. The drawings are for illustrative purposes only and are not drawn to scale.

[0052] It should also be understood that the terms "comprising," "including," "having," "containing," and / or "comprising," when used in this specification, indicate the presence of the stated features, integrals, steps, operations, elements, and / or components, but do not exclude the presence or addition of one or more other features, integrals, steps, operations, elements, components, and / or combinations thereof. Furthermore, when expressions such as "at least one of..." appear after a list of listed features, they modify the entire listed feature, not individual elements in the list. Additionally, when describing embodiments of this application, the word "may" is used to mean "one or more embodiments of this application." And the term "exemplary" is intended to refer to an example or illustration.

[0053] As used herein, the terms “basically,” “approximately,” and similar terms are used as terms of approximation rather than terms of degree, and are intended to describe inherent biases in measured or calculated values ​​that will be recognized by those skilled in the art.

[0054] Unless otherwise specified, all terms used herein (including technical and scientific terms) shall have the same meaning as commonly understood by one of ordinary skill in the art to which this application pertains. It should also be understood that terms (e.g., those defined in common dictionaries) shall be interpreted as having the meaning consistent with their meaning in the context of the relevant art and shall not be interpreted in an idealized or overly formal sense unless expressly so specified herein.

[0055] It should be noted that, unless otherwise specified, the embodiments and features described in this application can be combined with each other. This application will now be described in detail with reference to the accompanying drawings and embodiments.

[0056] The first embodiment of the present invention provides an adaptive morphology adjustment method for a long-span suspension antenna based on force density, such as... Figure 1 As shown, it includes the following steps:

[0057] Step 1: Given the material parameters, cross-sectional dimensions, and antenna structure elevation parameters of the antenna suspension cable and lower lead wire;

[0058] Step 2: According to the given mesh size or number of segments, divide the main cable structure into a mesh, and construct the node coordinate matrix and the node connection matrix;

[0059] Step 3: Based on the main cable structure, construct the self-weight load vector corresponding to the node number according to the principle of equal distribution of self-weight, and add the concentrated force term corresponding to the insulator chain;

[0060] Step 4: For the portion of the main cable structure excluding the lower lead wire, construct the final load vector based on the horizontal tension using a single variable method.

[0061] Step 5: Construct a full antenna structure for the lower lead portion of the main cable structure, and perform meshing on the lower lead to form cell and node information and load vector;

[0062] Step 6: Based on the design tension requirements of the suspension cable end and the lower lead adjustment rope end, give the initial value of the horizontal tension of all cable segments of the antenna and construct the force density vector;

[0063] Step 7: Construct the force equilibrium equations for all nodes and solve them to obtain the updated positions of the nodes;

[0064] Step 8: Determine whether the error between the positions of all antenna nodes and the positions calculated in the previous calculation is less than the set value. If it is, output the antenna shape that meets the design requirements and end the calculation; otherwise, update the force density vector and load vector according to the target tension requirement and return to step 7.

[0065] refer to Figures 1 to 5 The method provided by the present invention will be described in detail step by step below.

[0066] In this embodiment, step 1 may include the following:

[0067] Step 101: Given the spatial coordinates of the two suspension ends, the feed end of the lower lead, and the balance end of the lower lead of the main cable of the suspension antenna; to facilitate the shape-finding design, initially give the coordinates of the bifurcation point of the lower lead;

[0068] Step 102: Given the material parameters of all cables involved in the main cable and the lower lead wire, including the density vector ρ, the elastic modulus vector E, and the breaking tensile force vector F. f Given the cross-sectional parameters A of all cables, the number of elements in the above four types of vectors depends on the number of cable segment types that make up the actual antenna.

[0069] Step 103: Given the weight and position of each insulator chain in the main cable and the lower lead wire, it is reflected in the analysis model as a concentrated load;

[0070] Step 104: Based on the actual counterweight, specify the target tension for the traction force at the suspension cable traction end and the lower lead wire balance end, respectively. and

[0071] In this embodiment, step 3 may include the following:

[0072] Step 301: Divide the self-weight load of each unit of the main cable equally into the front and rear nodes as nodal loads;

[0073] Step 302: Superimpose the self-weight loads assigned to each node;

[0074] Step 303: Add the weight of the insulator chain to the corresponding nodes to construct the final load vector.

[0075] In this embodiment, step 4 may include the following:

[0076] Step 401: Determine the range of horizontal tension based on the required tension at the suspension end.

[0077] Step 402: Calculate the lower limit of horizontal tension The corresponding antenna configuration yields the corresponding suspension end tension.

[0078] Step 403: Calculate the lower limit of horizontal tension The corresponding antenna configuration yields the corresponding suspension end tension.

[0079] Step 404: Using the golden section method commonly used in univariate optimization, update the lower bound of the horizontal tension until the error between the tension at the suspension point and the target tension is less than the given tolerance. The suspension cable shape at this point is the final shape of the current step.

[0080] The calculation of the antenna configuration involved in steps 402-404 above is performed as follows:

[0081] Step 405(a): For the current horizontal tension, write the nodal equilibrium equations based on the force density method. For node q, its preceding and succeeding nodes are q-1 and q+1, respectively. The force equilibrium equations for this node are as follows:

[0082]

[0083] In the formula, subscript 1 represents the units connected to nodes q-1 and q, and subscript 2 represents the units connected to nodes q+1 and q. The specific node connection relationships are stored in the node connection relationship matrix R. zs Middle. F x1 and F x2 These represent the horizontal tensions of two units, where x and y are the coordinates of the nodes. and This represents the force density of the corresponding element. The third term on the left side of the equation is the self-weight load, and the fourth term is the concentrated load at that node.

[0084] Step 405(b): Write the above equilibrium equations for all nodes, and rearrange them with the y-coordinates of all nodes as unknowns to obtain the following system of linear equations.

[0085] K m*m Y m*1 =P m*1

[0086] Where K is the coefficient matrix formed by all force density vector sets, Y and P are the unknown vector and load vector, respectively, and m is the number of nodes of a single suspension cable.

[0087] Step 405(c): Solve the above system of linear equations to obtain the new y-coordinates of the nodes, while the x-coordinates of the nodes remain unchanged. At this point, the position of each node is uniquely determined. Update the node coordinate matrix L. zs .

[0088] Step 405(d): Determine whether the error between the position of all antenna nodes and the position calculated in the previous step is less than the set value. If the error is satisfied, the calculation of the current horizontal tension ends and the tension at the current suspension end is obtained. Otherwise, update the force density vector and load vector according to the target tension requirement and return to step 405(a).

[0089] In this embodiment, step 5 may include the following:

[0090] Step 501: Based on the suspension cable shape and the location of the lower lead bifurcation point from the previous step, give the initial position of the entire antenna structure;

[0091] Step 502: Based on the given unit size or number of segments, divide the lower lead wire and the lower lead wire adjustment rope into units to form unit and node information.

[0092] Step 503: Form the node coordinate matrix L of the overall structure total And node connection matrix R total .

[0093] In this embodiment, step 6 may include the following:

[0094] Step 601: Based on the requirements and current shape of the suspension end and the lower lead wire balance end, calculate the horizontal tension of the suspension end and the lower lead wire adjustment end, and then calculate the other three horizontal tension values ​​based on the force balance of the lower lead wire bifurcation point and the suspension point.

[0095] Step 602: Calculate the force density vector of all cable segments based on all horizontal tensions.

[0096] The specific operation of step 7 is described in step 405. The difference between this step and step 405 is that step 7 involves shape finding for all structures of the antenna (main cable and bottom lead), while step 405 only involves the main cable.

[0097] Compared with the prior art, this embodiment has at least the following advantages:

[0098] The force density-based adaptive morphological design method provided in this invention requires no manual intervention throughout the entire process. Compared to existing methods that rely on experience and simulation software calculations, it significantly improves design efficiency while ensuring accuracy.

[0099] A second embodiment of the present invention provides an electronic device, such as... Figure 6 As shown, it can be understood as a physical device, including a processor and a memory storing processor-executable instructions. When the instructions are executed by the processor, the following operations are performed:

[0100] Step 1: Given the material parameters, cross-sectional dimensions, and antenna structure elevation parameters of the antenna suspension cable and lower lead wire;

[0101] Step 2: According to the given mesh size or number of segments, divide the main cable structure into a mesh, and construct the node coordinate matrix and the node connection matrix;

[0102] Step 3: Based on the main cable structure, construct the self-weight load vector corresponding to the node number according to the principle of equal distribution of self-weight, and add the concentrated force term corresponding to the insulator chain;

[0103] Step 4: For the portion of the main cable structure excluding the lower lead wire, construct the final load vector based on the horizontal tension using a single variable method.

[0104] Step 5: Construct a full antenna structure for the lower lead portion of the main cable structure, and perform meshing on the lower lead to form cell and node information and load vector;

[0105] Step 6: Based on the design tension requirements of the suspension cable end and the lower lead adjustment rope end, give the initial value of the horizontal tension of all cable segments of the antenna and construct the force density vector;

[0106] Step 7: Construct the force equilibrium equations for all nodes and solve them to obtain the updated positions of the nodes;

[0107] Step 8: Determine whether the error between the positions of all antenna nodes and the positions calculated in the previous calculation is less than the set value. If it is, output the antenna shape that meets the design requirements and end the calculation; otherwise, update the force density vector and load vector according to the target tension requirement and return to step 7.

[0108] In the third embodiment of the present invention, the process of the adaptive morphology design method for a long-span suspension antenna based on force density is the same as that of the first and second embodiments. The difference lies in the engineering implementation: this embodiment can be implemented using software plus necessary general-purpose hardware platforms. While hardware implementation is also possible, the former is often a better approach. Based on this understanding, the method of the present invention can be embodied in the form of a computer software product stored in a storage medium (such as ROM / RAM, magnetic disk, optical disk), including several instructions to cause a device to execute the method of the embodiments of the present invention.

[0109] Through the description of specific embodiments, a more in-depth and specific understanding should be gained of the technical means and effects adopted by the present invention to achieve the intended purpose. However, the accompanying drawings are only provided for reference and illustration and are not intended to limit the present invention.

Claims

1. A method for adaptive morphological adjustment of a long-span suspension antenna based on force density, characterized in that, include: Step 1: Given the material parameters, cross-sectional dimensions, and antenna structure elevation parameters of the antenna suspension cable and lower lead wire; Step 2: According to the given mesh size or number of segments, divide the main cable structure into a mesh, and construct the node coordinate matrix and the node connection matrix; Step 3: Based on the main cable structure, construct the self-weight load vector corresponding to the node number according to the principle of equal distribution of self-weight, and add the concentrated force term corresponding to the insulator chain; Step 4: For the portion of the main cable structure excluding the lower lead wire, construct the final load vector based on the horizontal tension using a single variable method. Step 5: Construct a full antenna structure for the lower lead portion of the main cable structure, and perform meshing on the lower lead to form cell and node information and load vector; Step 6: Based on the design tension requirements of the suspension cable end and the lower lead adjustment rope end, give the initial value of the horizontal tension of all cable segments of the antenna and construct the force density vector; Step 7: Construct the force equilibrium equations for all nodes and solve them to obtain the updated positions of the nodes; Step 8: Determine whether the error between the positions of all antenna nodes and the positions calculated in the previous calculation is less than the set value. If it is satisfied, output the antenna shape that meets the design requirements, and the calculation ends. Otherwise, update the force density vector and load vector according to the target tensile force requirement, and return to step 7.

2. The adaptive morphology adjustment method for a long-span suspension antenna based on force density according to claim 1, characterized in that, Step 1 includes: Step 101: Given the spatial coordinates of the two suspension ends, the feed end of the lower lead, and the balance end of the lower lead of the main cable of the suspension antenna; initially, give the coordinates of the bifurcation point of the lower lead; Step 102: Given the material parameters of all cables involved in the main cable and the lower lead, including the density vector, elastic modulus vector, and breaking tensile force vector; given the cross-sectional parameter vectors of all cables; wherein the number of elements contained in each of the density vector, elastic modulus vector, breaking tensile force vector, and cross-sectional parameter vector depends on the number of cable segment types constituting the actual antenna. Step 103: Given the weight and position of each insulator chain in the main cable and the lower lead wire, it is reflected in the analysis model as a concentrated load; Step 104: Based on the actual counterweight, set the target tension for the traction force at the suspension cable traction end and the lower lead wire balance end respectively.

3. The adaptive morphology adjustment method for a long-span suspension antenna based on force density according to claim 1, characterized in that, Step 3 includes: Step 301: Divide the self-weight load of each unit of the main cable equally into the front and rear nodes as nodal loads; Step 302: Superimpose the self-weight loads assigned to each node; Step 303: Add the weight of the insulator chain to the corresponding nodes to construct the final load vector.

4. The adaptive morphology adjustment method for a long-span suspension antenna based on force density according to claim 3, characterized in that, Step 4 includes: Step 401: Determine the range of horizontal tension based on the required tension at the suspension end; Step 402: Calculate the antenna shape corresponding to the lower limit of the horizontal tension to obtain the corresponding suspension end tension; Step 403: Calculate the antenna shape corresponding to the lower limit of the horizontal tension to obtain the corresponding suspension end tension; Step 404: Using the golden section method commonly used in univariate optimization, update the lower bound of the horizontal tension until the error between the tension at the suspension point and the target tension is less than the given tolerance. The suspension cable shape at this point is the final shape of the current step.

5. The adaptive morphology adjustment method for a long-span suspension antenna based on force density according to claim 4, characterized in that, Step 4 further includes: Step 405(a): For the current horizontal tension, write the nodal equilibrium equation for a node based on the force density method; Step 405(b): Construct the balance equations for all nodes; Step 405(c): Solve the equilibrium equations to obtain the new node y-coordinates. At this point, the position of each node is uniquely determined, and the node coordinate matrix is ​​updated. Step 405(d): Determine whether the error between the position of all antenna nodes and the position calculated in the previous step is less than the set value. If the error is satisfied, the calculation of the current horizontal tension ends and the tension at the current suspension end is obtained. Otherwise, update the force density vector and load vector according to the target tension requirement and return to step 405(a).

6. The adaptive morphology adjustment method for a long-span suspension antenna based on force density according to claim 1, characterized in that, Step 5 includes: Step 501: Determine the initial position of the entire antenna structure based on the current suspension cable configuration and the location of the lower lead bifurcation point; Step 502: Based on the given unit size or number of segments, divide the lower lead wire and lower lead wire adjustment rope into units to form unit and node information; Step 503: Based on the unit and node information, determine the node coordinate matrix and node connection matrix of the overall structure.

7. The adaptive morphology adjustment method for a long-span suspension antenna based on force density according to claim 1, characterized in that, Step 6 includes: Step 601: Based on the requirements and current shape of the suspension end and the lower lead wire balance end, determine the horizontal tension of the suspension end and the lower lead wire adjustment end, and then calculate the remaining horizontal tension values ​​based on the force balance of the lower lead wire bifurcation point and the suspension point. Step 602: Calculate the force density vector of all cable segments based on all horizontal tensions.

8. An electronic device, characterized in that, The electronic device includes a memory and a processor. The processor is used to execute the adaptive morphology adjustment method for a long-span suspension antenna based on force density as described in any one of claims 1 to 7.

9. A storage medium, characterized in that, The storage medium stores a computer program, which, when executed by a processor, implements the adaptive morphological adjustment method for a long-span suspension antenna based on force density as described in any one of claims 1 to 7.

Images