Geometric Nonlinear Analysis of the Catenary Cable Element Based on UPFs of ANSYS

: The catenary cable element has more advantages than other nonlinear truss elements but is less used in commercial programs. In this paper, the initial geometric conﬁguration of the element is solved iteratively by the dichotomous method. Then the Updated Lagrangian (UL) Formulation for the two-node catenary cable element is combined with the element secondary development tool provided by ANSYS platform-User Programmable Features (UPFs) to develop a three-dimensional cable element—user101. The algorithm and procedure of this paper are veriﬁed through examples and a real bridge. The study shows that the developed cable element user101 is more accurate and faster than the ANSYS self-contained element. The method can effectively use the computational theory of nonlinear cable elements with catenary geometry and combine it with ANSYS commercial program, which saves computation time without reducing accuracy and has good practicality.


Introduction
As a prestressed tensioning system with tension cables as the main load-bearing members, the cable structure is widely used in suspension bridges [1,2], cable-stayed bridges [3][4][5], housing structures [6,7] or arch bridges with cable-stayed buckling construction [8,9] because of its simple forces and beautiful shape.However, because it is a highly flexible material, it will produce large rotation and displacement under external load.So, the geometric nonlinear effects must be considered in the accurate calculation of the cable structure.The calculation methods for cable structure can be divided into the analytical method for solving linear equations and the finite element method.
The analytical method is only applicable to the simple cable structure, while the finite element method is generally used for complex and variable cable structures.In 1965, Ernst modeled the cable as a two-node truss element to account for the non-linear effects of the cable element.He equated the modulus of elasticity of the truss element with taking into account the cable's sag effects, employing the Ernst formula [10].Still, this method was only for short cables with high stresses, and the longer the cable length, the less accurate it was.Later in 1971, Knudson improved this method by equating the cables to multiple truss elements [11], but the method required a sufficient number of truss elements to achieve accuracy.With the advancement of computer technology, some researchers consider the influence of the sag effects of the cable by modeling the cable as a multi-node curvilinear element based on the Lagrangian shape function [12][13][14], which started from the geometry of the cables and improved the computational accuracy of the geometric nonlinearity of the cable structure, but with many nodes, the degrees of freedom increased, and it was not easy to converge.In addition, to avoid the detrimental effect of numerous nodes, researchers have also proposed two-node curved cable elements [15,16], but still in the form of an infinite approximation to the shape of the cable structure.
On the other hand, in order to be closer to the geometry of the cable suspension, researchers have equated the cable to a segment of parabolic form and derived an explicit expression for its tangential stiffness matrix [17][18][19][20][21], which is more accurate and straightforward to calculate, but ignores the sagging effect of the cable and the error increases with the increase of the span.For this case, many researchers have developed two-node suspension curve elements by solving the equilibrium equations satisfying the elastic suspension curve of the cable element [22][23][24][25][26][27][28][29][30][31][32][33][34], which is widely used in the analysis of nonlinearities in cable structures due to the small number of nodes, high accuracy, and fast efficiency.In addition, Chen et al., developed a suspension link element with rigid arm for accurately simulating the rigid arm connection of cable-stayed bridges [35]; Thongchai et al., proposed a scalable cable element based on the variational approach, which considered one end to be free to slide horizontally [36]; Crusells-Girona et al., proposed a cable element based on a hybrid variational formulation of finite deformation curve coordinates, in which the material nonlinearity and geometric nonlinearity are considered to obtain high accuracy with few elements [37]; Andreu et al., proposed a cable network element based on the formula of deformable suspension curve element, which improves the calculation accuracy and efficiency [38]; Vanja and Ivica proposed a three-dimensional suspension curve element, which can perform static and dynamic analysis and can be directly applied to the dynamic calculation of suspension bridges under live load [39]; Chung et al., proposed a threedimensional elastic catenary cable element considering slip, and verified the correctness of the element by comparing experimental and theoretical data [40]; Riabi and Shooshtari proposed a stiffness matrix and nodal force equation for a three-dimensional cable structure under self-weight based on the suspension curve equation, thus developing a cable element considering material and geometric nonlinearity [41]; Rodrigo et al., investigated the effects of material and geometric nonlinearities on the cable element by introducing temperature effects [42]; Rezaiee-Pajand proposed a three-dimensional thermoelastic catenary cable element without rotating coordinates, which eliminated the solution of the inverse matrix of the catenary equation and improved the computational efficiency [43].
In addition to the above finite element methods, several researchers have adopted new methods to consider the nonlinearity of the cable structure.Croce proposed an exact expression for the equivalent stiffness of the cable element based on the principle of virtual work, which can be used as an approximate expression for the calculation of geometric nonlinearity [44]; Lee and Park developed a three-dimensional catenary cable element based on the Absolute Nodal Coordinates Formula (ANCF) and verified the accuracy of the element by comparing the static and dynamic calculation results [45]; Huttner et al., analyzed the computational efficiency and degree of convergence of nine different dynamic relaxation methods in the geometric nonlinearity of the cable element [46]; Bouaanani et al., investigated the thermoelastic response of the catenary cable element by the finite difference approach [47,48].
In summary, the advanced nonlinear theory of catenary cable element has been quite abundant, but most researchers still use the method of preparing finite element program to apply it to engineering examples, and the preparation of the program will lead to a series of shortcomings such as unfriendly pre-and post-processing interface, unable to store a large amount of computational data, low computational efficiency and not having a visualization function.ANSYS, an extensive general purpose finite element program, provides users with a secondary element development platform-UPFs, which facilitates researchers to develop units by combining proposed more advanced finite element algorithms with ANSYS secondary development tools.In order to improve the accuracy and efficiency of ANSYS to calculate the nonlinearity of complex structures, Deng et al., then developed planar beam elements using the more advanced CR algorithm [49], and the authors also developed planar truss elements by using the variational approach [50].When the pre-stress of the cable structure is high and the network topology is involved, the accurate initialization solution is essential for the convergence of the nonlinear structure.For this reason, in this paper, the initial configuration of the catenary cable element is first derived by dichotomous iteration, and then the algorithm of the element is derived based on the explicit expression of the tangent stiffness matrix of the two-node catenary cable element in the UL column and combined with the secondary development tools (UPFs) in the ANSYS platform.The two-node three-dimensional elastic catenary cable element (user101) is developed.The element only needs to obtain its single stiffness matrix, which is automatically formed in the ANSYS program for the total stiffness and operates the same way as the other elements.The method makes use of the powerful pre-and post-processing and visualization capabilities of the ANSYS program to combine a more accurate algorithm with a commercial program and has been validated by comparison with examples and natural bridges to show that the method improves the efficiency and accuracy of the calculation and has solid practical value.

Tangential Stiffness Matrix of Two-Node Catenary Cable Element Based on UL
The stress-strain relationship of the cable conforms to Hooke's law.

2.
The cable cannot be subjected to compression and bending but only to tension, being ideally flexible.

3.
The cross-section of the cable does not change under external loading.

4.
The self-weight constant load set of the cable is consistent along the length of the cable.

Tangential Stiffness Matrix of the Catenary Cable Element
Research manuscripts reporting large datasets that are deposited in a publicly available database should specify where the data have been deposited and provide the relevant accession numbers.If the accession numbers have not yet been obtained at the time of submission, please state that they will be provided during review.They must be provided prior to publication.In order to highlight the advancement and accuracy of developing the cable element, this paper adopts the tangential stiffness matrix of the cable element derived based on the displacement pattern consistent with the cable properties and basic assumptions [26], as shown in Figure 1, in the right-angle coordinate system OXYZ, with the suspension chain lead element as ij and the coordinates of the branch point i as (x i , y i , z i ) and the coordinates of pivot point j are (x j , y j , z j ).u i , v i , w i are node i displacements, and u j , v j , w j are node j displacements.2. The cable cannot be subjected to compression and bending but only to tension, being ideally flexible.3. The cross-section of the cable does not change under external loading.4. The self-weight constant load set of the cable is consistent along the length of the cable.

Tangential Stiffness Matrix of the Catenary Cable Element
Research manuscripts reporting large datasets that are deposited in a publicly available database should specify where the data have been deposited and provide the relevant accession numbers.If the accession numbers have not yet been obtained at the time of submission, please state that they will be provided during review.They must be provided prior to publication.In order to highlight the advancement and accuracy of developing the cable element, this paper adopts the tangential stiffness matrix of the cable element derived based on the displacement pattern consistent with the cable properties and basic assumptions [26], as shown in Figure 1, in the right-angle coordinate system OXYZ, with the suspension chain lead element as ij and the coordinates of the branch point i as (xi, yi, zi) and the coordinates of pivot point j are (xj, yj, zj).ui, vi, wi are node i displacements, and uj, vj, wj are node j displacements.Displacement pattern of the catenary cable element: where {u} = {u v w} T , {u ij } = {u i v i w i u j v j w j } We have L is the chord length of the cable element; β = ql/2H; q is the mean load set along the cable curve; l is the length of the horizontal projection of the cable; H is the horizontal force of the cable.
The tangential stiffness matrix of the cable element is then obtained from the UL columnar imaginary work increment equation.
where [k 0 ] is the elastic stiffness matrix and [k σ ] is the geometric stiffness matrix, expressed as follows [26]: where , 5. The tangential stiffness matrix [k] in the local coordinate system is converted to the tangential stiffness matrix [K] in the general coordinate system to obtain the equilibrium equation for the cable element in the general coordinate system [51]: where {∆U} is the displacement increment in the general coordinate system of the demand solution, {P} is the external load vector in integral coordinates, {F} is the force at the end of the cable in necessary coordinates, and the expression for the force at the end of the rope in its local coordinate system is as follows [26]: The iterative solution of the tangential stiffness matrix of the cable element can be divided into two types: (1) solving for the stress-free length s 0 of the cable according to the initial shape of the cable element; (2) solving for the constraint equation of the cable element by the known s 0 to obtain β, which can be substituted into Equation ( 3) to obtain the tangential stiffness matrix of the cable element.

Iterative Solution of a Catenary Cable Element with Known Pretension T j or T i
As can be seen from the tangential stiffness matrix expression of the catenary cable element in the previous section, the most critical thing is to solve for the tangential stiffness matrix β, and the need to determine the value of β requires the solution of the horizontal force H.In the initial tensioning of the cables of a cable-stayed bridge or the buckled anchor cables of a cable-stayed arch bridge with buckled construction, the cable end force T j or T i is often known, so the initial state horizontal force H 0 and the unstressed cable length s 0 can be obtained by seeking the cable end.For the catenary cable element shown in Figure 2, its differential cable equation under uniformly distributed vertical loads along the ties is: where D 1 = sin h −1 β 0 c/l sin hβ 0 + β 0 , β 0 = ql/2H 0 , H 0 is the horizontal force in the initial state of the cable, c and l are the vertical and horizontal lengths in the initial tensioned state of the cable.The formula for calculating the cable length is: Appl.Sci.2022, 12, x FOR PEER REVIEW 5 of 23 in Figure 2, its differential cable equation under uniformly distributed vertical loads along the ties is: where  = sinh / +  , β 0 = ql/2H 0 , H 0 is the horizontal force in the initial state of the cable, c and l are the vertical and horizontal lengths in the initial tensioned state of the cable.The formula for calculating the cable length is: The conditions for the equilibrium of forces along the x-axis lead to the following: From the above equation, the tensile forces on the ties at points i and j, respectively, are: The derivation is now made using the known Ti as an example and is obtained from Equation (11): From Equation ( 12), we have: Substituting Equation ( 12) into Equation ( 13) gives: The conditions for the equilibrium of forces along the x-axis lead to the following: From the above equation, the tensile forces on the ties at points i and j, respectively, are: The derivation is now made using the known T i as an example and is obtained from Equation (11): From Equation ( 12), we have: Substituting Equation ( 12) into Equation ( 13) gives: Let: a = ql/H 0 , the left-hand side of Equation ( 14) be f (a), then: When a ≥ ql/T i 0 (that is T i 0 ≥ H), f(a) is a single increasing function (f'(a) ≥ 0).Therefore, the equation f(a) = c + T i 0 /q has a single root and it can be obtained by the dichotomy method.For this purpose, it is first necessary to determine two lower and upper limits, a 1 and a 2 , that will encompass the value of a.In this paper, the initial lower and upper limits a 1 and a 2 are determined in the following way.
Because f(a) ≥ T i 0 cosha/q; so, c + T i 0 /q ≥ coshaT i 0 /q; which is: Make: Then, because, f (a) < l a T 0 i a ql 2 sin ha + T 0 i q exp(a), so, f (a) < T i 0 exp(a)/q, which is: c + T i 0 /q < T i 0 exp(a)/q, which is: a > ln(ql/T i 0 + 1), in addition, a ≥ ql/T i 0 , therefore: From Equations ( 17) and ( 18) are the initial upper and lower limits found.And by using the dichotomy method to solve the equation and thus obtain a value of a with sufficient accuracy.Following from a = ql/H 0 , H 0 is obtained, and D 1 is received from Equation (12), which in turn gives the cable equation and the unstressed cable length s 0 in the initial state of the tension cable.

Iterative Solution of the Tangential Stiffness Matrix of a Catenary Cable Element with
Known S 0 β can be obtained by solving the constraint equation for the unit conforming to the hanging chain trail if s 0 is known.This means that: Newton's method can be used to substitute the solved β i into Equation (3) obtain the cable's tangential stiffness matrix [k i ].

User Element Secondary Development Process
The user element is most notably the development of the nested programs UECxxx.F and UELxxx.F; the introduction of these two subroutines can be found in the literature [52] and will not be repeated here.The compilation process of the subroutine UECxxx.F, which provides the element property parameters, can be found in the literature [50], and this paper focuses on the compilation process of the core subroutine UELxxx.F: (1) Enter the variables to be used in the cable element; (2) Obtain basic information about the cable element: modulus of elasticity E, area A, Poisson's ratio µ, density den, initial strain ε, temperature ∆t; (3) Call the function to calculate the initial chord length L 0 of the cable element and the transformation matrix [T 0 ], horizontal and projection lengths l 0 and c 0 based on the initial position of the cable element; (4) Solve for β 0 with sufficient accuracy according to Equations ( 17) and ( 18) where the end force T i or T j is known and substituted into Equation ( 9) to derive the unstressed length s 0 of the cable and the initial state horizontal force H 0 ; (5) Calculation of the tangential stiffness matrix that is the elastic stiffness matrix [k 0 ] and the geometric stiffness matrix [k σ ], for the cable element in the structural coordinate system: when the first iteration is performed, the elastic stiffness matrix is calculated as in the linear analysis and the geometric stiffness matrix is 0; if it is not the first iteration, it is calculated according to Equation ( 3

Verification and Numerical Examples
In order to verify the correctness and superiority of the development element, the following calculation example is used.

Single Cable Structures Subjected to Concentrated Loads
Figure 4 shows a single cable structure (f0/L = 0.1) with a load set of q = 46.12N/m, an

Verification and Numerical Examples
In order to verify the correctness and superiority of the development element, the following calculation example is used.

Single Cable Structures Subjected to Concentrated Loads
Figure 4 shows a single cable structure (f 0 /L = 0.1) with a load set of q = 46.12N/m, an area of 548.4 mm 2 , a self-weight drape of f 0 = 30.48,a modulus of elasticity of the cable E = 131 × 105 kPa and a span L = 304.8m, calculate the displacement at point C at the position shown when subjected to a concentrated load of P = 35.586kN.In the analysis, the cable structure was divided into two cable elements, AC and CB, and the load was applied in one step.The calculation results are shown in Table 1, and the load-displacement curve at point C under 20 load steps is shown in Figure 5.

Verification and Numerical Examples
In order to verify the correctness and superiority of the development element, the following calculation example is used.

Single Cable Structures Subjected to Concentrated Loads
Figure 4 shows a single cable structure (f0/L = 0.1) with a load set of q = 46.12N/m, an area of 548.4 mm 2 , a self-weight drape of f0 = 30.48,a modulus of elasticity of the cable E = 131 × 105 kPa and a span L = 304.8m, calculate the displacement at point C at the position shown when subjected to a concentrated load of P = 35.586kN.In the analysis, the cable structure was divided into two cable elements, AC and CB, and the load was applied in one step.The calculation results are shown in Table 1, and the load-displacement curve at point C under 20 load steps is shown in Figure 5.As can be seen from the calculation results in Table 1 and Figure 5, the results of this paper using the developed two-node catenary cable element user101 are very close to those found in other literature, and a total of two elements were divided as proof of the correctness and feasibility of the developed element.

Verification and Analysis of Tensioned Cable Structures
A calculation model for the construction phase of an asymmetric cable-stayed bridge is shown in Figure 6, with a cable-stayed cable with rigid arms on each side of the tower, elements 5 and 6, for a total of eight nodes and nine elements. 1 -4 in Figure 6 indicate main beam elements 1-4, 5 -6 in Figure 6 indicate cable-stayed cable elements 5-6, 7  in Figure 6 indicates lower tower column element 7. 8 -9 in Figure 6 indicate upper tower column elements 8-9.The modulus of elasticity E = 1.097×10 7 MPa and other basic calculation parameters are shown in Table 2.   rigid arm elements are shown in Table 3.The load set of the diagonal cables q = 4.37 kN/m and the relative positions of the rigid arm elements are shown in Table 3.The calculation is divided into two working conditions: (1) The two cables are synchronized for initial tensioning at the end of the tower, with the initial tension of 27,154 kN for the left-hand cable, that is, element No. 5, and 33,037 kN for the right-hand cable, that is, element No. 6, anchored after tensioning is completed; (2) A concentrated force couple of 8 × 10 6 kN•m is applied simultaneously at each of the two cantilevered ends of the main beam.
Figure 7 and Table 4 show that the displacement results of the catenary cable element user101 developed in this paper based on ANSYS secondary development and the cable element developed by the program prepared in the literature [34] differ by 0.06%.The bending moment results only differ by 0.02% under the action of working condition one.The displacement results only differ by 0.08%, and the bending moment results only differ by 0.24% under working condition two, thus verifying the accuracy of the cable element developed in this paper.The accuracy of the cable element developed in this paper is thus confirmed.

Engineering Example Validation and Analysis
The Shatuo Bridge in Guizhou Province, China, is a reinforced concrete cantilevered arch bridge with a net span of 240 m, a net vector height of f = 40 m, a net vector-to-span ratio of ƒ/L = 1/6, and an arch axis factor of m = 1.85.The arch ring section form is a single-box double-chamber structure, box width 10 m, height 4.5 m, where the thickness of the top plate from the foot of the arch is 80 cm to L/4 section of 35 cm, and then to the top of the arch L/2 section of 50 cm gradual; bottom plate thickness from the foot of the arch is 80 cm to L/4 section of 35 cm gradual, L/4 section to the top of the arch bottom plate thickness; web thickness from the foot of the arch 65 cm to L/4 section of 40 cm at the foot of the arch and then tapered to 50 cm at the L/2 section of the top plate, as shown in Figure 8.The main arch circle is divided longitudinally into 37 sections and is constructed by hanging basket cantilever casting.The site layout of the buckling anchor cables for the Shatuo Special Bridge is shown in Figure 9. 4        In order to verify the correctness and practicality of the developed elements, this paper uses the measured data from the Shatuo Special Bridge for comparison with the ANSYS theoretical model data.APDL parametric modeling is used for the Shatuo Bridge; as shown in Figure 10b, the whole bridge is divided into 1808 elements and 1011 nodes, of which the 3D beam element Beam188 simulates the central arch circle, abutment pier, and buckling tower, the buckling anchor cables are simulated by the development element user101 and truss element link10 respectively.The buckling anchor cables are simulated by the virtual rigid beam Beam44 between the main arch circle and the abutment pier.Boundary conditions: Solidified form at the foot of both arches, anchor cable anchorage position is solidified.Loads: 90 t for the hanging baskets plus formwork and 40 t for the hangers plus formwork.Both were acting as concentrated loads on the arch ring.In order to verify the correctness of the development element user101, the finite element model was corrected by reference [55], and the updated values of each material parameter are shown in Table 5.Table 6 shows the number, area and initial tensioning force value of the buckled anchor cable of Shatuo Bridge.In order to verify the correctness and practicality of the developed elements, this paper uses the measured data from the Shatuo Special Bridge for comparison with the AN-SYS theoretical model data.APDL parametric modeling is used for the Shatuo Bridge; as shown in Figure 10b, the whole bridge is divided into 1808 elements and 1011 nodes, of which the 3D beam element Beam188 simulates the central arch circle, abutment pier, and buckling tower, the buckling anchor cables are simulated by the development element  For verifying the advancement and accuracy of the development element user101, the measured values of the maximum vertical downward deflection, the measured values of the maximum cable force, and the measured values of the maximum tensile stress in the top and bottom slab of the arch were compared with the calculated values simulated by element user101 and element link10 respectively.The following data and graphs were obtained.The bridge was divided into 63 construction phases.The newly developed element user101 was used to simulate the buckled anchor cables in just 37 min, while the ANSYS element link10 was used to simulate a total of 43 min, an improvement of 6 min in calculation efficiency.
Table 7 and Figure 11 show that the cable structure simulated by the user101 element can be matched well with the experimental results.And the relative error of the newly developed element user101 in calculating the main arch deflection of each casting section is 8.4% at most, while the maximum relative error of element link10 is 13.6%.This is because the link10 element, as a two-node truss element, does not take into account the sag effect of the cable, while the sag effect of the buckled anchor cable will cause the main arch to deflect more.Table 8 and Figure 12 show that the relative error between the maximum buckling force and the measured value of the buckling anchor cable calculated using the newly developed element user101 is 5.6%.In comparison, the maximum relative error of element link10 is 5.7%.The maximum buckling force simulated by element link10 is smaller than the measured value and is smaller than the maximum buckling force simulated by element link10, which is due to the increase in the sag effect caused by the larger buckling Table 8 and Figure 12 show that the relative error between the maximum buckling force and the measured value of the buckling anchor cable calculated using the newly developed element user101 is 5.6%.In comparison, the maximum relative error of element link10 is 5.7%.The maximum buckling force simulated by element link10 is smaller than the measured value and is smaller than the maximum buckling force simulated by element link10, which is due to the increase in the sag effect caused by the larger buckling length, resulting in a smaller buckling force and an increase in the lower deflection of the arch.length, resulting in a smaller buckling force and an increase in the lower deflection of the arch.Figure 13 shows that the maximum tensile stress in the top slab is more significant, and the maximum tensile stress in the bottom slab is less than that in the link10 element and is more in line with the measured converted tensile stress values due to the buckling Figure 13 shows that the maximum tensile stress in the top slab is more significant, and the maximum tensile stress in the bottom slab is less than that in the link10 element and is more in line with the measured converted tensile stress values due to the buckling effect of the buckling anchor cable causing a greater deflection in the cast section of the arch ring and a smaller buckling anchor cable force.effect of the buckling anchor cable causing a greater deflection in the cast section of the arch ring and a smaller buckling anchor cable force.In summary, by comparing the deflection values of the above arch casting sections, the maximum cable force of the buckled anchor cable, and the maximum tensile stress of the top and bottom slabs, it can be concluded that the use of the user101 element to simulate the buckled anchor cable is more matching with the experimental results, with higher accuracy and faster calculation efficiency than the link10 element.

Conclusions
With the emergence and development of more advanced nonlinear mechanics theories, researchers have realized combining theory and experiment through computer programming.Still, due to the shortcomings of programming, it cannot be applied to large and complex structures, while ANSYS, as a large general-purpose software, can solve this problem well.This research mainly combines the nonlinear mechanic's theories with the UPFs technology of ANSYS, integrates the advantages of both, and applies the more advanced nonlinear mechanics theories in ANSYS for computational analysis.
1.For the nonlinear analysis of cable structures, higher computational accuracy and efficiency are always the goals to be pursued, and when the pretension of the cable structure is high and involves network topology, the accurate initialization solution is essential for the convergence of the nonlinear structure, and this paper adopts the dichotomous method to iterate the initial configuration of the cable unit with pretension to accelerate the convergence and thus improve the computational efficiency.2. Through the comparison and verification between the example and the engineering model, it can be seen that, according to the characteristics of large span cable-stayed bridges or arch bridges with cable-stayed suspension construction, the initial geometric configuration of the cable is first solved iteratively by the dichotomous method under the known initial tension of the cable, and then the equilibrium equation of the cable element is solved iteratively according to the nonlinear geometric algorithm of the two-node catenary cable element, so the catenary cable element user101 developed in this paper can be directly applied to the ANSYS model for analysis and calculation; In summary, by comparing the deflection values of the above arch casting sections, the maximum cable force of the buckled anchor cable, and the maximum tensile stress of the top and bottom slabs, it can be concluded that the use of the user101 element to simulate the buckled anchor cable is more matching with the experimental results, with higher accuracy and faster calculation efficiency than the link10 element.

Conclusions
With the emergence and development of more advanced nonlinear mechanics theories, researchers have realized combining theory and experiment through computer programming.Still, due to the shortcomings of programming, it cannot be applied to large and complex structures, while ANSYS, as a large general-purpose software, can solve this problem well.This research mainly combines the nonlinear mechanic's theories with the UPFs technology of ANSYS, integrates the advantages of both, and applies the more advanced nonlinear mechanics theories in ANSYS for computational analysis.

1.
For the nonlinear analysis of cable structures, higher computational accuracy and efficiency are always the goals to be pursued, and when the pretension of the cable structure is high and involves network topology, the accurate initialization solution is essential for the convergence of the nonlinear structure, and this paper adopts the dichotomous method to iterate the initial configuration of the cable unit with pretension to accelerate the convergence and thus improve the computational efficiency.

2.
Through the comparison and verification between the example and the engineering model, it can be seen that, according to the characteristics of large span cable-stayed bridges or arch bridges with cable-stayed suspension construction, the initial geometric configuration of the cable is first solved iteratively by the dichotomous method under the known initial tension of the cable, and then the equilibrium equation of the cable element is solved iteratively according to the nonlinear geometric algorithm of the two-node catenary cable element, so the catenary cable element user101 developed in this paper can be directly applied to the ANSYS model for analysis and calculation; 3.
A new type of element, user101, developed by using the nonlinear finite element algorithm of the catenary cable element combined with the ANSYS secondary development platform (UPFs), has a certain improvement in nonlinear calculation accuracy and efficiency compared to the existing commercial program ANSYS developed truss

Figure 2 .
Figure 2. Local coordinate system of the cable.

Figure 2 .
Figure 2. Local coordinate system of the cable.

Figure 4 .
Figure 4. Single cable structures subjected to concentrated loads.

Figure 4 .
Figure 4. Single cable structures subjected to concentrated loads.

Table 1 .
Comparison of displacement at point C (unit: m).

Figure 5 .
Figure 5. Calculated results for a single cable structure: (a) Vertical displacement, units in m; (b) Lateral displacement, units in m; (c) Axial force diagram, units in kN; (d) Load-displacement curve.

Figure 6 .
Figure 6.Calculation model diagram.Units in m.

23 4. 3 .
-5 and 8 -9 in the figure indicate No. 4-5 piers and No. 8-9 piers. 6-7 in the figure indicate No. 6-9 archs.The Scene picture of Shatuo Bridge is shown in Figure 10a.Appl.Sci.2022, 12, x FOR PEER REVIEW 15 of Engineering Example Validation and Analysis The Shatuo Bridge in Guizhou Province, China, is a reinforced concrete cantilevered arch bridge with a net span of 240 m, a net vector height of f = 40 m, a net vector-to-span ratio of ƒ/L = 1/6, and an arch axis factor of m = 1.85.The arch ring section form is a singlebox double-chamber structure, box width 10 m, height 4.5 m, where the thickness of the top plate from the foot of the arch is 80 cm to L/4 section of 35 cm, and then to the top of the arch L/2 section of 50 cm gradual; bottom plate thickness from the foot of the arch is 80 cm to L/4 section of 35 cm gradual, L/4 section to the top of the arch bottom plate thickness; web thickness from the foot of the arch 65 cm to L/4 section of 40 cm at the foot of the arch and then tapered to 50 cm at the L/2 section of the top plate, as shown in Figure 8.The main arch circle is divided longitudinally into 37 sections and is constructed by hanging basket cantilever casting.The site layout of the buckling anchor cables for the Shatuo Special Bridge is shown in Figure 9. ④-⑤ and ⑧-⑨ in the figure indicate No. 4-5 piers and No. 8-9 piers.⑥-⑦ in the figure indicate No. 6-9 archs.The Scene picture of Shatuo Bridge is shown in Figure 10a.

Figure 8 .
Figure 8. Diagram of the cross-sectional form of the arch ring: (a) section of arch foot; (b) L/4 cross section; (c) vault cross section.Units in cm.

Figure 9 .
Figure 9. Site layout of the Shatuo Bridge buckling anchor cable.Units in cm.

Figure 8 . 23 4. 3 .
Figure 8. Diagram of the cross-sectional form of the arch ring: (a) section of arch foot; (b) L/4 cross section; (c) vault cross section.Units in cm.

Figure 8 .
Figure 8. Diagram of the cross-sectional form of the arch ring: (a) section of arch foot; (b) L/4 cross section; (c) vault cross section.Units in cm.

Figure 9 .
Figure 9. Site layout of the Shatuo Bridge buckling anchor cable.Units in cm.

Figure 9 .
Figure 9. Site layout of the Shatuo Bridge buckling anchor cable.Units in cm.

*Figure 11 .
Figure 11.Comparison of maximum deflection values for each section of the arch ring: (a) Comparison of different elements with experimental values; (b) Comparison of the absolute error of different elements.

Figure 11 .
Figure 11.Comparison of maximum deflection values for each section of the arch ring: (a) Comparison of different elements with experimental values; (b) Comparison of the absolute error of different elements.

Figure 12 .
Figure 12.Comparison of maximum buckle force values: (a) Comparison of different elements with experimental values; (b) Comparison of the absolute error of different elements

Figure 12 .
Figure 12.Comparison of maximum buckle force values: (a) Comparison of different elements with experimental values; (b) Comparison of the absolute error of different elements.

Figure 13 .
Figure 13.Maximum tensile stresses of the main arch ring: (a) Comparison of maximum tensile stresses in the top slab for different elements and experimental values; (b) Comparison of maximum tensile stresses in the bottom slab.

Figure 13 .
Figure 13.Maximum tensile stresses of the main arch ring: (a) Comparison of maximum tensile stresses in the top slab for different elements and experimental values; (b) Comparison of maximum tensile stresses in the bottom slab.

Table 2 .
Basic calculation parameters.Figure 6. Calculation model diagram.Units in m.
The calculation is divided into two working conditions:

Table 4 .
Deformation value of square frame under concentrated tension.
* (1) ∆v 1 denotes the vertical displacement of node 1, M 2 and M 3 are the end section bending moments of element 2 and element 3, respectively; (2) BDCMS considers only the drape effect of the cable using the Ernst formula.Appl.Sci.2022, 12, x FOR PEER REVIEW 13 of 23

Table 6 .
The initial tension and area of the anchor cable.

Table 7 .
Maximum vertical deflection of each segment of the arch ring.
* Numbers in parentheses indicate the absolute error percentage with respect to experimental results.

Table 7 .
Maximum vertical deflection of each segment of the arch ring.

Table 8 .
Comparison of calculated and measured values of the maximum cable force of the arch ring buckles.
* Numbers in parentheses indicate the absolute error percentage with respect to experimental results.

Table 8 .
Comparison of calculated and measured values of the maximum cable force of the arch ring buckles.
* Numbers in parentheses indicate the absolute error percentage with respect to experimental results.