Dynamical Analysis and Simulation Validation of Incompletely Restrained Cable-Suspended Swinging System Driven by Two Cables

The flexibility of the suspension multicables and driven length difference between two cables cause the translation and rotation of the platform in the incompletely restrained cable-suspended system driven by two cables (IRCSWs2), which are theoretically investigated in this paper. The suspension cables are spatially discretized using the assumed modes method (AMM) and the equations of motion are derived from Lagrange equations of the first kind. Considering all the geometric matching conditions are approximately linear with external actuator, the differential algebraic equations (DAEs) are transformed to a system of ordinary differential equations (ODEs). Using linear boundary conditions of the suspension cable, the current method can obtain not only the accurate longitudinal displacements of cable and posture of the platform, but also the tension between the platform and cables, and the current method is verified by ADAMS simulation.


Introduction
The dynamics of cable-driven parallel mechanisms have been a topic of interest since the introduction of the first designs.Cable-driven parallel mechanisms, due to cable light weight and the ability to resist relatively large axial load, have a good choice for a variety of applications, such as telescope [1], haptic interface [2], elevators [3,4], rescue robotics [5], heavy load transportation [6], and mine hoisting [7].Cable-driven parallel mechanisms are suitable for lifting, positioning large heavy load, and digging a vertical shaft.
The cable-driven mechanisms are generally divided into two categories: the full constrained cable-driven mechanism [8,9] and cable-suspended mechanism [10][11][12].The moving platform of the full constrained cable-driven mechanism can fully control all degrees of freedom; however, the suspended cable mechanism uses external force to maintain all cables in tension.Many cable-driven parallel mechanisms have been mostly considered quasi-static devices; for instance, classic determination method of the static workspace is addressed in the related literature [13].Dynamic models of cablesuspended mechanism are built and studied on the optimal tension distribution among the cable [14] and the deflection of the suspended platform under the varying cable length.Dynamic behaviors of cable-suspended mechanism are also studied in [15,16] and others.
All the works as mentioned above did not consider the influence of the longitudinal vibration of cable on cablesuspended mechanisms; however, they are subjected to the longitudinal vibration caused by the external excitation for high flexibility of cable.The longitudinal vibrations of cables have been studied extensively for decades.Particularly, in the typical sinking winch mechanism [17][18][19], simulation platforms of a ship and a portable crane, the mechanism used to sink the vertical shaft is composed of winches, head sheaves, a construction platform, and suspension cables.In this mechanism, the platform is hung by suspension cables.Deb [20] analyzed the lateral vibration of a string fixed at both ends excited by the lateral impact of an elastic load at any point using Heaviside's expansion theorem.Subsequently, many scholars have focused on the vibrations of cables characterized by time-varying length.Nguyen and Hong [21] analyzed the longitudinal and transverse vibrations of an axially moving string system at right side exerted by a hydraulic actuator equipped with a damper using Hamilton's principle.Sandilo and van Horssen [22] studied the free lateral responses of vertically translating media with variable length, velocity, and tension, subject to general initial conditions.The cable-suspended mechanism consists of not only distributed-parameter components, such as the suspension cable, but also lumped-parameter components, such as the head sheave and the suspended platform [23].A reduction method is developed by Zhu et al. [24,25] to deal with such structural system into distributed-and lumped-parameter components, and the distributed-parameter components can be discretized using the AMM and the complicated natural matching conditions can serve as the constraints of the generalized coordinates [26].
Indeed, cable-driven parallel mechanisms have the potential to produce very fast motions and the control of such motions requires a proper understanding of the vibration of the mechanical system [22,27].However, swing conditions often occur on these platforms.In order to improve reliability of equipment on these swing conditions, the swing experiments are designed for testing the behaviors of equipment and elevators.An important feature of the system is that the platform is shaken by the suspension cables.In addition, the tensions in the suspension cables should be properly allocated to meet the safety and performance requirements.Hence, it is necessary to evaluate the design by calculating the dynamic responses of the platform and the cables.

Description of Incompletely
Restrained Cable-Suspended Swinging System (IRCSWs2) According to the schematic, the 3D model of the novel IRCSWs2 is designed in Figure 1.The IRCSWs2 is composed of two cables, pulleys, actuators, base frame, and suspended platform.Four connecting joints are symmetrically distributed around the suspended platform.Each cable is driven by an actuator in the middle of each cable, and it can be driven back and forth with a periodic motion.The system is equivalent to the suspended platform hung with four suspension cables.
For simplicity and without loss of generality, IRCSWs2 (Figure 1) can be described as the th suspension system model shown in Figures 2(a) and 2(b).Table 1 explains the meaning of the symbols adopted in Figure 2.This model is composed of four suspension cables of displacement   (t) at time instant t and the total length   ; each suspension cable is passed over a set of head sheaves, and the suspended platform is attached to the lower ends of the four suspension cables shown in Figure 2(a).Based on invariability of suspension cable's length, it can be straightened in vertical line shown in Figure 2(b) and the corresponding sheave is located with  The suspended platform is driven by four suspension cables with the upper ends fixed at actuators.It takes the actuator displacement as input and then takes the pose of the suspended platform and tensions and longitudinal vibration in cables as output.The input-output mathematics model of such a system is derived.When one actuator is working

Suspension cable
Sheave Figure 2: The th suspension system model.independently, the suspended platform is swinging motion about a single coordinate axis.
In this paper, the AMM and Lagrange equations of the first kind are combined and applied to IRCSWs2.The suspension cables and the platform components are connected by mechanical joints that impose restrictions on their relative motion.Finite rotations introduce geometric nonlinearities; hence, the systems are inherently nonlinear.Mechanical joints result in algebraic constraints leading to a set of governing equations that combines differential and algebraic equations.Several challenging problems, such as the tensions between four suspension cables and the suspended platform, the deformations of the suspension cables, and the pose of the suspended platform, are treated.Four suspension cables are spatially discretized using the AMM and the equations of motion are derived from Lagrange equations of the first kind, while the geometric matching conditions at the interfaces of the cables are accounted for by the Lagrangian multiplier.The performance of dynamically periodic motions with small amplitude will be investigated.Thus, all the geometric matching conditions can be approximately linear, and the resulting spatially discretized equations, which are differential algebraic equations (DAEs), are transformed to ordinary differential equations (ODEs) and solved by an ODE solver.

Theoretical Model of IRCSWs2
3.1.Spatial Discretization.Ignoring the lateral vibration of the suspension cables and the influences of the frictional force, the IRCSWs2 can be modeled as four cables with fixed length and the pulleys can be set as lumped-parameter on cables.And the kinetic energy   associated with the longitudinal vibrations is where   () is the mass distribution function of the suspension cable and it is defined as   () =  + ∑   ( −   ),  = 1 ∼   (  is the number of the sheave of th suspension cable);  and (⋅), respectively, represent the unit mass and the Dirac delta function;   and I  are the mass and inertia tensor matrix of the suspended platform; r  and  are the position and the angular velocity vectors of the suspended platform.
The relationship of the body-fixed angular velocity vector  can be determined by resolving the Euler rates into the bodyfixed coordinate fame.The operator / is given by Furthermore, the angular and linear velocity vectors of the suspended platform can be written as where The inertia tensor of the platform with respect to the origin of   - is denoted as The total potential energy   of all the cables is where E is Young's modulus of the cable; A is the cross section area of cable;  is gravity acceleration;   (, ) is the static tensions in the cables at position  due to the gravitational acceleration  and given as where  0 is the force of th suspension cable caused by the suspended platform and ℎ() is unit step function; if  ≥ 0, then ℎ() = 1; else if  < 0, then ℎ() = 0.The strain   can be approximated as   =   /.The first term in (8) represents the potential energy caused by the tension and strain in the suspension cable, while the second term represents the potential energy associated with gravity of the suspended platform and suspension cables.The geometric boundary conditions for the suspension cable and the suspended platform are obtained as where in which the rotation matrix R can be written in terms of the roll-pitch-yaw Euler angles as follows: where c(⋅)and s(⋅) denote shorthand writings for sin(⋅) and cos(⋅) functions, respectively.And , , and  are the rotation angles around the fixed axes , , and , respectively.The paper is focused on the influence of small cables' axial motion due to the longitudinal variation of cable length on their dynamic behaviors.With the linear vibration approximations, that is, || ≪ 1, sin  ≈ , cos  ≈ 1, ( = , , ), and second-order item omitted, the nonlinear geometric matching conditions at the interface between the suspension cable and the suspended platform are rewritten as The transformation  = ( −   ())/  is used to transform the time-varying spatial domain [  (),   () +   ], which is converted to a fixed domain [0, 1] for .By using the separation of variables method, the longitudinal displacements can be approximated by expansions of a complete set of trial functions and expressed as where  represent the number of assumption mode methods;  , (), in which  = 1 ∼ 4, are the generalized coordinates;  , are the trial functions and should satisfy the homogeneous boundary conditions in (13) and can be expressed as Substituting ( 1), (8), and (13) into Lagrange equations of the first kind yields the equations of motion with holonomic constraints where M is mass matrices; Q is force vector; q = (q  1 , . . ., q  4 , r   ,   )  (q  = [ ,1 , . . .,  , ]  ,  = 1∼ 4) is the vector of generalized coordinates; and it should be noted that these (4 + 6) generalized coordinates are not independent; however, the geometric matching conditions (13) yield the holonomic constraints g of the generalized coordinates, where g = ( 1 ,  2 ,  3 ,  4 )  is a vector including all the constraint conditions in (13); G = g/q is the Jacobian of the constraint equations, which is a 4 × (4 + 6) matrix; and  is the array of 4 Lagrange's multipliers; denote the tensions between the suspension cable and the suspended platform.
For selection of the independent generalized coordinates, a transformation matrix T  is introduced, where p = ( where Suspended platform Pulley Cable 1 Base frame Cable 2

DAEs to ODEs for Solution.
Considering the constraint conditions in ( 13) are linear algebraic equations ( 17), a system of DAEs could be transformed to ODEs.Equation ( 17) could be expressed in the form By reduction method, the independent generalized coordinate vector p 1 is selected and the DAEs of constrained mechanical system can be transformed to ODEs problem.Thus the cable tensions are obtained by evaluating ( 23) The parameters G, g  (), p 1 , ṗ 1 , Φ, and ü () in ( 21) and ( 22) are shown in Appendix.And the initial displacement and velocity of the suspension cable are necessary for the numerical solution using Runge-Kutta method.

Simulation Model of IRCSWs2
An ADAMS simulation is built in order to validate the performance, which is shown in Figure 3.The method and procedure are presented by referring to literature [28].
The key parameters involved in the dynamic simulation of IRCSWs2 are listed in Notations.
The distance between the cable attachment points to the origin of   - is 2.07 m.Two translational actuators are built to control the suspension cable with periodic motion.The actuator movement profile is shown in Figure 4.The actuator displacement is sine function with amplitude 0.2 m and frequency 0.1 Hz, while the jerk is constant during startup.The equivalent mass of sheaves is assumed to be constant.

Shock and Vibration
1st and 2nd suspension cable 3rd and 4th suspension cable

The Convergence Analysis of IRCSWs2.
It is necessary to know the number of terms about assumption mode method which decides the reliability of the result and computation time.Above all, the convergence of the solution for governing equation is examined by varying the number of included modes in Figure 5.The results indicate that the model ( = 10) converges very fast under lower order modes.In order to decrease computing time, the number of the trial function used in the AMM is  = 10.

Examples and Discussion
. The results of theoretical calculation and ADAMS simulation are denoted by black lines and red lines in Figures 6, 7, and 8, respectively.Two translation joints are used to control the status of the cables.Swing motions were successfully demonstrated with the results.All results present approximately sinusoidal variation with a sinusoidal excitation in Figure 6.The axial tension has finite impulse response that is caused by sudden switch at 1 second just as shown in Figure 7.The displacements along the  axis,  axis and  axis are nearly zero, which validated the assumptions under small excitation.
The two results of theoretical calculation and ADAMS simulation are in excellent agreement (Figure 7).The cable tension ranges between 28 kN and 33 kN which indicates that the cable tension is small compared to the tensile capacity of the cable (safety factor larger than 4) and that auxiliary protections are usually designed in the practice to prevent failures.The results show that the longitudinal vibration is only slightly smaller than ADAMS results due to neglecting the rigid body displacement of the suspended platform in  and  direction (Figures 6 and 8); however, during a cycle time of 10 s, their variable trends have no great difference.The higher mode of high-frequency variation is caused by the constant jerk at the beginning of the motion.Due to the low driving frequency and damping effect, this highfrequency variation is quick to disappear.Therefore, the Shock and Vibration theoretical model and the ADAMS simulation model are both reasonable, and these two models could be adopted independently to study the behavior of IRCSWs2.
For the normal case of system resonance and low levels of excitation, the phase trajectory of the generalized coordinate  11 when the initial stage is zeros is shown in Figure 9. Circle family of phase trajectory is changing when the excitation is different.The status of the system accelerations over an extended period and the system behavior are on steady state of maintaining swing motion of suspended platform with system dissipation.In the upper half plane (d 11 > 0), the motion direction of system status is along positive direction of horizontal axis; in the down half plane (d 11 < 0), the phase trajectory moves to left similarity.The rate of phase trajectory at horizontal axis is always parallel to the vertical axis of the graph.

Conclusions
The dynamical characteristic of IRCSWs2 is investigated in aspects of both theoretical model and the ADAMS simulation model.The following conclusions could be drawn: (1) The focus of this paper is to introduce and analyze the motion of this novel IRCSWs2 by the AMM and Lagrange equations of the first kind.The performance of dynamically periodic motions with small amplitude is investigated.Thus, all the geometric matching conditions can be approximately linear, and the resulting spatially discretized equations are transformed to ordinary differential equations (ODEs) by reduction method.
(2) In order to validate the theoretical model, the corresponding ADAMS simulation is built.The results of theoretical model are consistent with the ADAMS simulation by comparison.
(3) The translation displacement of suspended platform is almost near zero, and rotations of angles about, respectively, and -axes present sinusoidal variation.The rotation about  axes is tiny enough, which can be neglected.Tensions and longitudinal displacement also present sinusoidal variation, which is useful to design and optimize the cable-suspended system.

Theoretical model ADAMS model
Theoretical model ADAMS model    where p can be obtained by substituting the solution p 1 into the geometric constraint equation (A.5) and q can be assembled as q = T  (Φ () p 1 + u ()) . (A.8) ) and 2(b) describe the longitudinal vibrations along the axial transport motion.The displacement   represents the longitudinal displacement of the particle of th suspension cable at position  at time t.

Figure 5 :
Figure 5: Longitudinal response with different numbers of modes.

Figure 6 :
Figure 6: The pose of the suspended platform.

Figure 9 :
Figure 9: The phase trajectory of  11 : (a) no longitudinal excitation; (b) the excitation of the four suspension cables is given by () = 0.01 sin().
The matrices M and Q are expressed as