Nonlinear mechanics of flexible cables in space robotic arms subject to complex physical environment

A nonlinear model of a special cable in space robotic arms is developed in space environment. The mechanic effects of control cables in powerful robots can often be neglected. However, in complex space multi-physics environments, involving ultra-low temperature, radiation, and other extreme conditions of outer space, the externally mounted cables (protected by shielding layers) can induce strong nonlinear interference to robot arms; and this can induce further small-range slow rotations or oscillations of the flexible joint of robots at a specific posture, which consequently affect the precision and operation performance of end effectors. Effective mathematical models on nonlinear mechanics of strong cables under multi-physics environments and their effects on weak robots have not been well developed yet. Complex key factors, such as low gravity, nonlinear friction, and unexpected curved surface constraints, have not been extensively investigated either. In this study, considering all these key factors, a Kirchhoff nonlinear mechanical model of cables in complex space environments is developed, and a relatively improved algorithm based on a trust-region strategy is proposed for solving this nonlinear model, based on which the geometry and terminal force of the modeled robot cable can be obtained. The validity and accuracy of the proposed algorithm and theoretical calculation results are verified via experiments. The theoretical findings revealed in this study are significant to future research on the slow rotations and oscillations of weak robot joints in space exploration with robotic arms.


Introduction
The complex nonlinear mechanical characteristics of control cables could be critical to the stability and accuracy of the joint operation of robot arms. Flexible cables, which are employed for various functions including power supply and information transmission, are indispensable parts of robotic systems. In general, the internal and external flexible cables employed in robot arms are relatively very weak compared to expected robotic payload and output torques. Their nonlinear mechanical properties can usually be negligible consequently. However, considering space environments, including ultra-low temperature, radiation, and other extreme environment of outer space, coin-cided with resource constraints, the precision manipulator joint is often relatively weak with limited output torque. Moreover, in order to preserve communication security and provide a reliable supply of power, flexible cables are often wrapped with thick nonmetallic shielding layers which can become exceptionally rigid under the conditions of long-term radiation or a cold environment. Consequently, externally mounted cables protected by shielding layers can induce strong nonlinear interference effects on robot arms. Owing to this interference, the flexible joints of robot arms could slowly rotate or oscillate within a small range when the robot arm is in a certain posture; this significantly affects the precision and operation performance of the end effectors of robots. It is thus very relevant to study the nonlinear mechanics of flexible cables in space robots under the conditions of multi-physics coupling dynamics. Robots employed for planetary explorations are often exposed to special gravity fields and placed in extremely limited space. Gravity, friction, and surface geometric constraints are often the key factors that need to be studied to determine the mechanical properties of cables.
In recent years, flexible bodies including cables have been extensively investigated in various fields. Payandeh et al. [1] revealed the geometric shape of a surgical suture by using particle-spring model, which accurately represents the tensile deformation of a flexible body that is prone to elastic deformations. Servin et al. [2] proposed a connecting rod chain model, i.e., a series of connecting rods, to characterize the cable shape and study the pulling deformation of the cable without considering twisting effects. Lazarus et al. [3] employed quaternions to develop a flexible cable model. In addition, an asymptotic numerical method has been employed to obtain the cable shape; however, this method only considered cable twisting. The above-mentioned methods focused on the shape of the cable. In contrast, a nonlinear mechanical model of an elastic rod can accurately reveal the cable shape and mechanical properties at discrete points along the cable. The Cosserat [4,5] and Kirchhoff [6][7][8][9][10][11][12][13][14] theories are mainly employed in this type of model. Owing to its simplicity, the Kirchhoff cable mechanics model is currently widely used [15][16][17][18][19][20][21]. By implementing the Kirchhoff theory and elastic thin rod static model [22], Liu et al. [23][24][25] developed a cable mechanical model with smooth plane contact to improve the virtual reality assembly technology for cables. Although their model is useful, they ignored the effects of gravity on the cable; moreover, the contact between the cable and curved surface was not considered. Consequently, their model is still insufficient for practical application. A model to describe the contact of DNA molecules with curved surfaces has been proposed in [26]. However, this method does not consider various effects such as gravity, friction, and length, and the flexible body was required to remain in full contact with the curved surface; thus, the method cannot be employed when the flexible body and curved surface are partially restrained. In addition, it can only be employed to solve the initial value problem and cannot be employed to estimate the actual working conditions. Based on previously reported theoretical results [22][23][24][25][26], the effects of gravity and friction on the cable are considered to propose a nonlinear mechanical model of a cable under localized constraints of a curved surface in this paper. When implementing multi-physics coupling methods, the gravity force distribution on the heavy cable, which is determined by the safety shield, is an important factor that should not be neglected. Owing to the limited space resources, it is necessary to consider the friction and geometric constraints between the cable and curved surface. Taking into account the above issues, it is clear that the mechanical constraints on the cable are very complex, the dimensions of the equations increase in size, and the solution of the cable nonlinear mechanics model is even more challenging. In order to ensure global convergence of the iterative calculations, a more reliable solution algorithm is yet to be developed. Therefore, by considering the gravity, friction, and curved surface constraints, and using a planetary exploration robot as the main object of the model, we employ the Kirchhoff theory to study the joint stability and develop a nonlinear mechanical model of a flexible cable that has one end fixed and other end free. Using the finite difference method, the differential equations are discretized into a set of nonlinear equations, and then transformed into a global optimization problem. In order to ensure global convergence of the iterative calculations, an improved algorithm is proposed based on the trust-region strategy. The large-scale nonlinear optimization problem is solved, and the mechanical parameters of each discrete point of the cable are obtained. For the calculation results, Euler parameters are employed to obtain the coordinates of each discrete point. In addition, a cubic spline curve is used to fit the discrete points in order to obtain the geometric shape of the cable in the world coordinate system. The contributions of this study lie in the following points: 1. First, we introduce and analyze a planetary exploration robot subjected to multi-physics coupling dynamics and consider the space robot joint stability in the presence of strong nonlinear cable interference; 2. Then, we develop a mathematical model of the flexible cable in space robotic arms; 3. Using the developed methods, the mechanical properties and shape of a cable in space robotic arms constrained to a curved surface are revealed. Furthermore, the nonmetallic shielding layers age and harden with the increase in the space radiation time.

2-DOF observation robot joint: slow rotation phenomenon
As is shown in Fig. 2a, for a single joint, the cable is fixed to the upper and lower links of the joint in the 2-DOF space observation robotic arm, with one end of the joint being installed on the fixed link, whereas the other end is connected to the free link. When the robot moves within a certain range, the cable follows the free end of the joint.
In addition, there are many research results on the error elimination control of flexible structures. A neural network controller is developed to suppress the vibration of a flexible robot system with input deadzone in [34], and the proposed controller has good performance in suppressing vibration. Also, through a boundary controller, a control problem about input backlash in a flexible robotic manipulator is well solved in [35]. A cooperative control with good control effect is proposed for a nonuniform gantry crane (consisting of a flexible cable and other components) with constrained tension in [29]. And an optimal feedback linearization approach is used to find the optimal path between the initial and the target points in the case of a cable robot under the maximum dynamic load. These are very good research results in controlling cable robots in [28,30,31]. Different from the above error elimination algorithms, high-precision position sensors cannot be used at present (due to the strict size constraints and limitations of high and low temperature space environments), the joint rotation can only be open-loop controlled by a stepper motor. Therefore, the problem proposed can only be solved by developing a theoretical cable model instead of using a control algorithm to eliminate the error.
As shown in Figures 2b, c, the joint is driven by a single stepper motor through a planetary reducer. The inherent gap θ j1 (less than 1 degree) of this joint is mainly caused by the planetary reducer. The axis of rotation of the robotic joint is parallel to the direction of the microgravity field.
Defining T m as the load torque of the joint: where l c is the arm of F c with respect to the joint axis. The internal damping torque T c of the joint includes the bearing friction torque, internal cable resistance torque of the joint, and damping torque induced by the planetary reducer. There is an inherent gap θ j1 at the joint, and when the joint rotates, the cable force F c changes nonlinearly. Slow rotation of a robotic joint may occur when the load torque T m is close to 0, as is illustrated in Fig. 2b. Because the presence of the gap is evident at a particular position when the load torque T m is close to 0. In Fig. 2b, Position 1 is the target stop position of the robotic joint for end-effector operation. Owing to the cable force F c , the joint slowly rotates until it reaches θ j2 , which corresponds to Position 2. Owing to the internal damping of the joint, the joint gap θ j2 is smaller than the inherent gap θ j1 . Position 2 is the final position after the slow rotation of the joint. Moreover, at Position 2, the torque of the cable force F c , offsets the internal damping torque T c of the joint with respect to the joint axis; this places the robotic system in equilibrium (T m = 0). When the load torque T m is much greater than 0, under the effect of the load torque T m , the joint will reach Position 2 in a flash, and the slow rotation of a robotic joint will not occur. When F c tends 0, the joint will remain motionless under the effect of the internal damping torque T c , the joint will stay at Position 1, and this phenomenon will also not occur.
In [27], although constraints were imposed at both ends, only the influence of the disturbance torque of a single cable on a stabilized platform was studied. However, there are more complex constraints on cables used in the manufacture of robots purposed for outer space missions; furthermore, the robotic joint is influenced by the mechanical characteristics of multiple cables. The strong interference effects induced by the cables are most significant when the joint is in a particular position. The free end of the robot joint is equipped with camera sensors, and the slow rotation of the joint can significantly reduce the quality of the images recorded by the cameras. It is worth noting that there is currently no extensive research being carried out in this area. The mechanical properties of cables implemented in robots are the main reasons for these phenomena. In order to study the above phenomena, we develop and solve a nonlinear mechanical model of a robot-purposed cable under the conditions of a complex outer space environment.

Basic assumptions
We made the following assumptions for the development of a simple cable mechanical model in the 2-DOF space observation robotic arm [22][23][24][25]: 1. The cross section of the cable is a rigid circular section, and the cable centerline is always perpendicular to the cable cross section. 2. The centerline of the cable is a smooth curve exceeding 2nd order; the length of the cable is significantly larger than the diameter of its cross section. 3. The cable is isotropic; there is a linear relationship between stress and strain. 4. The deformations of the cable attributed to stretching and shearing are ignored.

Coordinate system transformation
Three different coordinate systems are employed for the developed cable mechanical model ( Fig. 3):

Coordinate system O-XY Z
The coordinate system O-XY Z is employed as the world coordinate system. The starting point of the cable centerline corresponds to the origin O of the coordinate system. In this coordinate system, an arc coordinate s is introduced along the cable centerline; any point P at the cable centerline can be determined via the arc coordinate s. The vector of the point P, which is set with respect to the origin O, is denoted as r; the vector that corresponds to the point P 0 is r 0 .

Coordinate system P-N BT
The frame {P |N , B, T }, which can be regarded as a rigid body, follows the movement of the point P. The parameters affecting the point P are as follows: tangent vector T(s) = dr/ds, curvature κ = |dT/ds|, main normal vector N(s) = dT/(κds) (κ = 0), vice normal vector B(s) = T(s) × N(s), and torsion τ = |dB/ds| . The cable centerline can be obtained via Eq. (2); ε is the integration variable. 3 Coordinate system transformation for the cable model

Coordinate system P-x yz
The coordinate system P-x yz is the spindle coordinate system for the cable section (cable section and coordinate system P-x yz are fixed). The z-axis coincides with the T -axis, while the plane N P B coincides with the plane x Py. The angle between the N -and x-axes is equal to that between the B-and y-axes, and is denoted as α, which represents the twist angle. The twist angle α reveals the torsional characteristics of the cable. The use of the transformation matrix D can overcome the pose problem between the P-N BT and P-x yz coordinate systems: Using the curvature κ, and torsion τ and twist angle α occurring at the point P, we define the curvaturetwisting vector ω of the point P as where (ω x , ω y , ω z ) is the curvature-twisting vector ω projected onto the coordinate system P-x yz.

Static equilibrium equations for the cable
Taking into account that the motion speed of the robot is relatively small and stable, and that θ j1 < 1 • , we can statically analyze the cable; this means that we can ignore the dynamical processes. As is shown in Fig. 3, in the world coordinate system O-XY Z, the microarc P 0 P of the cable is considered in the static analysis. A gravitational force per unit length f acts on the microarc P 0 P, which is oriented opposite to the Z -axis. Assuming that the arc coordinate of the point P 0 is s 0 , the arc coordinate of the point P is s, with s = s 0 + s and s > 0. When the microarc P 0 P is in equilibrium, the torque generated by the gravitational force per unit length f is the second-order small value of s, which is negligible because the shape of the cross section of the cable is round. Thus, based on Yanzhu Liu's study [22], we obtain the following: Then, we divide Eqs. (6) and (7) by s ( s → 0). Therefore, dF ds By changing the reference coordinate system for Eqs. (8) and (9) to the coordinate system P-x yz, we obtaiñ dF ds Equations (10) and (11) are projected onto the coordinate system P-x yz. Thus, we obtain The wavy lines in Eqs. (10) and (11) indicate that the reference coordinate system is P-x yz, which is different from that employed for Eqs. (8) and (9). Assuming that the cable has no initial curvature or twist, and that the constitutive relationship of the cable is known, we obtain (14) where A and B represent the flexural stiffnesses about the x-and y-axes of the cable cross section, respectively, and C is the torsional stiffness of the cable cross section about the z-axis. (15) Assuming that the cable is uniform and isotropic, E and G represent the Young's modulus and shear modulus of the cable, respectively. I x and I y are the moments of inertia for the cable cross section with respect to the x-and y-axes, and I z is the polar moment of inertia for the cable cross section with respect to the z-axis.
For a cable with a diameter d c , As is shown in Fig. 4, q k (k = 1, 2, 3, 4) are the Euler parameters for the cable cross section (the cable cross section is perpendicular to the z-axis.). The Euler parameters can be determined by Eq. (17). The coordinate system P-x yz is rotated about the axis e by an angle φ with respect to the coordinate system O-XY Z.
where e k (k = 1, 2, 3) correspond to the direction cosines of the limited axis of rotation e with respect to the world coordinate system O-XY Z; φ is the limited rotation angle. The derivatives of q k (k = 1, 2, 3, 4) with respect to the arc coordinate s are employed as new variables, as follows: Next, we define a vector c = c 1 c 2 c 3 c 4 c 5 c 6 c 7 T .

Discretization of the equations
Equation (22) can be solved by using the finite difference method. We assume that the length of the cable before and after deformation is the same (L 0 ) (Fig. 5). The cable is uniformly divided into n parts with lengths ofl i =l 0 = L 0 /n. Therefore, the cable has (n + 1) discrete points and n calculation units. The calculation unit between the discrete points P i and P i+1 is denoted When n → ∞, For Part i (ith calculation unit), Using Eq. (22), we determine that, for each calculation unit (Part i), w i = u i1l0 u i2l0 u i3l0 u i4l0 u i5l0 u i6 l 0 u i7 T . Therefore, where Next, we define a vector When the cable is in equilibrium, we obtain The rotation matrix of the coordinate system P i+1 -x yz at the discrete point P i+1 , with respect to the world coordinate system O-XY Z, is The rotation matrix P i+1 O R is an orthogonal matrix, i.e., its inverse matrix is equal to its transposed matrix.
The distribution forces on the calculation unit Part i are concentrated at the discrete point P i+1 . The gravitational force per unit length f at the discrete point P i+1 is projected as (0, 0, f 0 ) in the world coordinate system O-XY Z. The projection of f onto the coordinate system P i+1 -x yz is: If we substitute Eq. (30) into Eqs. (32), (33) is obtained.
The point P 1 coincides with the origin O of the world coordinate system. In the O-XY Z coordinate system, P 1 P i+1 = r i+1 = r X,i+1 , r Y,i+1 , r Z ,i+1 , and the value of P 1 P n+1 = (r X , r Y , r Z ) should be known. Using Eqs. (2) and (30), we obtain the following [22,23,25]: By implementing the cable length condition [23][24][25], we obtain The boundary conditions of the cable (Euler parameters of the discrete points P 1 and P n+1 ) are u n1 corresponds to the calculation unit Part n, and the squared sum of the Euler parameters at the discrete point P n+1 is repeatedly defined; hence, u n1 is replaced by c 0 in Eq. (39). If the cable is in equilibrium, Eq. (39) is valid.
Equation (39) determines the parameters for all discrete points in Eq. (40). In an ideal case, the optimal solution of the problem is determined by minF(x) = 0, where x is a 7 (n + 1)dimensional vector, x 0 represents the iteration initial value [Eq. (41)], and x k represents the k th iteration point. The output value of F(x) is compared to the numerical iteration termination criterion for each iteration.
In [24], the Levenberg-Marquardt algorithm is proposed to solve the cable model under the condition of constraints imposed at both ends; however, this algorithm may cannot guarantee numerical convergence of the cable model under complex conditions. Thus, in this study, the algorithm is improved by using the trustregion method to ensure that the theoretical model of the cable in contact with a curved surface has good numerical stability. A flowchart for the numerical iterative algorithm is presented in Fig. 6.
The convergence of the calculation results is influenced by the iterative initial value, and the iterative initial value x 0 is selected as a column vector whose most components are 1. Define g(x k ) as the gradient of the function F(x) at x k , J(x k ) is the Jacobian matrix for v at x k , I is the unit matrix, d k is the iterative step of x k , and μ k is the adjustment parameter for the iteration step d k . Therefore, According to Eqs. (44) and (45), the parameter μ k is adjusted by using the trust-region method; r k is the radius of the trust region of the parameter μ k , q is the scaling factor of the parameter μ k , where q = 2 ∼ 10. The critical values for r k and scaling factor q can be, The step size factor for the iterative calculations is obtained by using the Armijo criterion, where β ∈ (0, 1) and σ ∈ (0, 0.5).

Calculation results
The Denavit-Hartenberg coordinate systems for the 2-DOF observation robot are shown in Fig. 7 (unit of length: mm). The robot has two DOFs; it has an orientation joint and a pitch joint, which can realize the movement of the cameras carrying cables within a certain range. The origin of the world coordinate system O-XY Z is set to correspond to the starting point O of the cable, and O 1 and O 2 are located at the geometric centers of the orientation and pitch joints, respectively. The experiment shows that the phenomenon of slow rotation may occur at the orientation joint, whereas it does not occur at the pitch joint. Therefore, in this study, a single joint, i.e., the orientation joint, is considered, while the pitch angle of the pitch joint is set to zero. For simplicity, only the right cable of the robot was employed for the analysis, and the experiments are currently being conducted using ordinary cables instead of real cables covered by the nonmetallic shielding layers. The z-axis direction of the coordinate system P-x yz of the discrete point P 1 is opposite to the Y -axis direc- tion, and the z-axis direction at the discrete point P n+1 is the same as the Z -axis. The X -and Z -coordinates of the end point P n+1 of the right cable are, respectively, 20 and 25 mm smaller than those of the point O 2 when γ = 0. The coordinates of the point P n+1 in the world coordinate system O-XY Z are ⎧ ⎨ ⎩ r X,n+1 = 160 − 145 sin γ + 80 cos γ r Y,n+1 = 80 sin γ + 145 cos γ r Z ,n+1 = 75 where γ is the rotation angle of the orientation joint in the 2-DOF space observation robotic arm. The material parameters of the cable are shown in Table 1. According to Eq. (39), the cable mechanical model considers the cable shape and mechanical properties. With the change in the rotation angle γ of the orientation joint, the morphological curves of the robot cable are projected onto the X O Z, Y O Z, and X OY planes. Figure 8 shows the projection of the calculated trajectory of the robot cable onto each coordinate plane, and Fig. 9 shows the comparison between the theoretical shape of the robot cable and the actual trajectory. The trends of the cable theoretical shape and the actual trajectory are consistent. In addition, when the actual trajectory of the robot cable in the XOY plane is collected (by a camera sensor), the cable is blocked by the robotic arm. To obtain a more complete actual trajectory of the robot cable, the optical axis of a camera sensor used is not perpendicular to the XOY plane (there is a small angle of inclination), resulting in a slightly large error in the comparison of the cable shape in the XOY plane. Figure 10 shows the projection of the internal force onto the end of the cable in the P-x yz coordinate system.

Curved surface contact
Taking into account the complexity of the contact between the cable and curved surface, in order to simplify the analysis, the diameter of the cable is ignored, and the surface is considered to be spherical. Therefore, a nonlinear model of the cable with spherical contact is developed. As is shown in Fig. 11, the coordinates of the center of the sphere O 0 are (X 0 , Y 0 , Z 0 ), the radius is R 0 , and the equation for the spherical surface G is For the actual solution, we first assume that the curved surface does not exist; by using Eq. (39), we can identify the coordinates of all discrete points. Then, we Next, let us assume that the number of elements in the point set U is n c (n c ≥ 1). The subscript size (i + 1) of P i+1 ∈ U is used as the standard for sorting. Then, we select the intermediate point P d in the point set U . If the number of cases is even, we choose a larger discrete subscript to ensure that there is a complete calculation unit in U . When the serial numbers of the cable discrete points in the curved surface are sorted from small to large, t 1 is the subscript size of the first discrete point. d is the subscript size of the middle point of the cable discrete points in the curved surface. When n c is an even number, d = t 1 + n c /2, and when n c is an odd number, d = t 1 + (n c − 1)/2. Only the elements that belong to the point set U are likely to be in contact with the curved surface. When  Fig. 11 Constraint of the cable which is in contact with a curved surface the cable is in contact with the curved surface, the tangent unit vector p 1 at the cable contact point is perpendicular to the unit vector of the sphere radius p 2 , i.e., p 1 · p 2 = 0 [26]. In order to ensure convergence of the solution of the nonlinear mechanical model, the discrete point P d is regarded as an element at the curved surface G. For the discrete point P d , Eq. (55) can be employed; otherwise, Eq. (56) can be employed for P i+1 ∈ U . When P i+1 ∈ U , the curved surface may exert a contact distribution force f sp,i+1 on the cable along the same direction p 2 as that of the sphere radius at the point P i+1 . f sp,i+1 = 0 when the discrete point P i+1 is on the curved surface G, whereas f sp,i+1 = 0 when it is not on the curved surface G [25]. By combining Eqs. (55) and (56), we obtain n c equations, where f sp,i+1 (t 1 ≤ i +1 ≤ t 1 +n c −1) are n c Eq.(39) requires the addition of n c unknown varibles and n c equations (But in this paper, we have added n unknown variables and n equations in order to easily modify the calculation code under different calculation conditions. When P i+1 / ∈ U , the additional constraint added at the discrete point P i+1 is: f sp,i+1 − 0 = 0.). If P i+1 ∈ U , the force distribution in Eq. (33) needs to be modified according to Eq. (57).
In addition, the distribution forces at the calculation unit Part i are concentrated at the discrete point P i+1 . As is shown in Fig. 12, v t 0 is the direction of movement of the end point P n+1 of the cable (v t 0 is also perpendicular to the Z -axis in the coordinate system O-XY Z), and the cable segment at discrete point P 1 is fixed. When the rotation angle of the orientation joint is relatively small, the movement of the cable at the P 1 H segment is also relatively small, so the movement at the discrete point P i+1 can be simplified as the rotation of the point P i+1 around the point H . Therefore, the relative motion trend v t at the discrete point P i+1 is per- Fig. 12 Frictional constraint of the cable on a curved surface pendicular to p 1 and p 2 , which means that the direction of v t can be determined;f i+1 is the friction force distribution at the discrete point P i+1 (Part i) that is oriented in the direction opposite that of v t . By neglecting the cable axial friction, we obtain μ c = 0.35, which is the friction coefficient between the cable and surface.
The projection off i+1 onto the world coordinate system O-XY Z is given as It is worth noting that v provides (7n + 11 + n c )-dimensional nonlinear overdetermined equations, which can be employed to obtain 7 (n + 1)+n c parameters (But in the calculation code of this paper, ν is (7n + 11 + n)-dimensional and we get 7(n + 1) + n parameters.). According to the algorithm illustrated in Fig. 6, the nonlinear mechanical model of a cable with curved surface contact can be solved. When a spline curve is used to fit the discrete points, the simulated cable may slightly overlap the curved surface. Thus, the interpolation points at the curved surface have to be modified by using Eq. (58); Z is the corrected coordinate value.
With the change in the rotation angle γ of the orientation joint, the shape of the cable in contact with the curved surface can be obtained. Figure 13 shows the projections of the calculated trajectory of the cable, which is in contact with a curved surface, onto each coordinate plane. Figure 14 shows the comparison between the theoretical shape of the robot cable and the actual trajectory in contact with a curved surface. When γ = π/4, π/3 and 5π/12, the robot cable comes into contact with the curved surface (spherical obstacle). The trends of the cable theoretical shape and the actual trajectory are consistent. Figure 15 shows the projection of the internal force on the end of the cable, which is in contact with the curved surface, in the coordinate system P-x yz. When γ = π/3, due to the contact force of the curved surface, the end-point force (in Fig. 15) under the constraint of a curved surface is smaller than the end-point force under the restraint of both ends (in Fig. 10). This is consistent with the actual situation, and to a certain extent, it also verifies the correctness of the theoretical results.

Conclusion
Under the complex conditions of multi-physics outer space environment, cables are the main sources of unexpected nonlinear disturbance effects on space robot joints. This induces slow rotation or oscillation of the joint at a certain position. In this study, by considering the gravitational force distribution, friction force distribution, and contact with the curved surface of a robot cable, we employ the mechanical parameters for each discrete point of the cable as variables and develop a discrete nonlinear mechanical model of a robot cable.
The following conclusions can be drawn: 1. In this paper, the nonlinear mechanical model of a cable in the space robotic arm is transformed into a large-scale nonlinear optimization problem, and the mechanical characteristics of each discrete point of the cable are systematically revealed including the end-point force and geometry of the cable. This model could be employed to facilitate stable operation of robotic arms designed for future outer space missions; 2. A solution algorithm, which ensures global convergence of the algorithm and is based on a trust-region adjustment strategy, is proposed to solve the nonlinear mechanical model of a cable. The rationality and accuracy of the proposed nonlinear mechanical cable model are verified by performing a morphological comparison; 3. At present, we use ordinary cables for our experiments. The next step is to use real space cables for experimental research. At the same time, we will further improve the experimental environment and make it as consistent as possible with the outer space environment. In addition, it is also important for future studies to theoretically determine the specific positions at which a robotic joint may exhibit slow rotation or oscillation, and to perform parallel light pipe experiments to estimate these positions on a micro-scale.

Compliance with ethical standards
Conflict of interest The authors declare that they have no conflict of interest.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http:// creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.