Screw Motion via Matrix Algebra in Three-Dimensional Generalized Space

: This paper aims to investigate the screw motion in generalized space. For this purpose, ﬁrstly, the change in the screw coordinates is analyzed according to the motion of the reference frames. Moreover, the special cases of this change, such as pure rotation and translation, are discussed. Matrix multiplication and the properties of dual numbers are used to obtain dual orthogonal matrices, which are used to simplify the manipulation of screw motion in generalized space. In addition, the dual angular velocity matrix is calculated and shows that the exponential of this matrix can represent the screw displacement in the generalized space. Finally, to support our results, we give two examples of screw motion, the rotation part of which is elliptical and hyperbolic.


Introduction
A rigid body is a system of particles with a fixed distance from one another, so every displacement of a rigid body be an isometry. An isometry is called a rigid displacement when it has properties such as rotation and translation. The set of rigid displacements is called the Euclidean group [1,2].
The term "kinematics", which originates from Greek, means motion, and it is a branch of mechanics that deals with the analysis of the motion of particles and rigid bodies [3].
In order to represent rigid motion in Euclidean or Lorentzian space equipped with multiple coordinate frames, it is required to determine the "rotation matrix" and "translation axis" concepts. In particular, if the axis of rotation is parallel to the axis of translation, this transformation is called screw motion and the axis is called the screw axis. These concepts are used to create homogeneous transformation matrices, which represent the position and direction of a coordinate frame relative to the other. These transformations allow us to navigate from one coordinate frame to another [4,5].
To obtain frame M from frame F, it is necessary to apply a rotation determined by R, and translation (with respect to F) given as n, where F is a fixed frame and M is a moving frame. This transformation called coordinate transformation, which is denoted as D : F → M, and the coordinates of it are determined y = Rx + n, where y is the position of x after the displacement. In this notation, R is an n × n orthogonal matrix called a rotation matrix, and n is an n-dimensional vector called a translation. This transformation is denoted by D = (R, n) [3,5].
Screw theory was introduced by Chasles in 1830 in order to explain the concept of the twist motion of a rigid body. Then, in 1865, Plücker called this a screw. Later, in 1900, Ball explained the kinematics of a rigid body by using screw theory [6][7][8]. The screw theory is defined as the concept of rotation and translation about a common axis; it can be thought of as any helix of rotation. A common example is the thread of a screw (hence the name screw theory). In engineering, a screw thread is defined by a pitch, i.e., the ratio of translational motion to rotational motion. To describe the direction in which the screw points, an axis is needed along which it points. Finally, to describe a point on the screw, we need a magnitude, a measure of distance along the screw that allows us to find our position on the spiral. Thus, to describe rigid body motion using a screw, we need three components: a pitch, an axis, and an amount [9].
The derivative of motion represents the velocity of a point from frame F to frame M. Linear velocity is the instantaneous rate of change in the linear position of a point relative to some frame. The angular velocity describes the rotational motion of M with respect to F [3,5,10].
Dual numbers were initially introduced by Clifford in 1873. Their first applications to kinematics were in the study by Kotel'nikov in 1895 and in 1903. Bottema and Roth, in 1978, used dual numbers in theoretical kinematics. Then, dual quaternions were applied to spatial kinematics by Agrawal in 1987. Further information on dual quaternions can be found in the work of McCarthy from 1990 [3,5,[11][12][13][14][15].
Vector spaces with a quadratic form generate an associative algebra with unity called Clifford algebra. The quaternion algebra, a subalgebra of Clifford algebra, is isomorphic to the rotation group. The generalized quaternions H are a four-dimensional algebra that is associative but not commutative. The algebra of a generalized quaternion is a natural generalization of quaternion algebra H [16][17][18][19][20].
According to the definition, the generalization space is a Euclidean space if α = β = 1 and a semi-Euclidean space if α = 1 and β = −1. In this study, all situations of α and β are considered, except α = β = 0. Some concepts that are used in the definition of the scalar product and vector product in the generalization space are different from those in Euclidean and Lorentzian space.
Beggs gave a derivation for a screw matrix by using two different coordinate systems in 1965. In 2018, Sıddıka studied the screw theory in Lorentzian space L 3 . In 2021, Ata and Savcı showed that unit-generalized quaternions H(α, β) correspond to a rotation in generalized space E 3 (α, β). They obtained the screw axis of a displacement, the dual Rodrigues equation, and the Euler parameters of a spatial displacement, in generalized space. They stated that this motion gives the rotational motion in three-dimensional Euclidean space E 3 and in three-dimensional semi-Euclidean space for their specific choices of α, β ∈ R. In 2022, Savcı obtained the screw axis of a displacement, the dual Rodrigues equation, and the Euler parameters of a spatial displacement, in generalized space E 3 (α, β) [10,[20][21][22][23].
Since the rotational motion in Euclidean space is on the sphere, it is circular, whereas the rotational motion in Lorenz space is hyperbolic because it is on hyperboloids. With the appropriate selection of α and β in generalized space, rotational motions such as circular, hyperbolic, and elliptical can be obtained. Thus, this gives the opportunity to work in many areas, such as from the centrifugal force taken into account when calculating the slope of the road to the motion of an electric field [24] and the orbits of the planets [25,26]. Generalized space is also a scalar product space. For general information about scalar product spaces and their algebraic properties, see [16][17][18]20,[24][25][26][27][28][29].
Screw motion is an important part of rigid body motion, and it is vital in various fields, such as multi-body dynamics, mechanical design, computational geometry, and robot mechanics. From this point of view, it is important to reveal the wider application prospects of the screw motion. With our study, the details of which are given below, screw motion has been examined for all situations in generalized space.
In this paper, we explain the change in the screw coordinates according to the motion of the reference frames and consider the specific cases of changes in the coordinates. Then, we obtain Plücker coordinates for the screw axis in generalized space. Moreover, we obtain dual orthogonal matrices to simplify the screw coordinate transformations by using matrix multiplication and the properties of dual numbers. Finally, we obtain a dual angular velocity matrix and show that screw displacement can be represented as the exponential of this matrix in generalized space E 3 (α, β), which includes the Euclidean and semi-Euclidean space.

General Information
In this section, we will provide some basic information that is used in the manuscript. Definition 1. Let r = (r 1 , r 2 , r 3 ), s = (s 1 , s 2 , s 3 ) be two vectors in R 3 and α, β ∈ R. The generalized scalar product is defined by The vector space R 3 , with the generalized scalar product defined over it, is called a three-dimensional generalized space and is denoted by E 3 (α, β) = (R 3 , <, > G ). This space includes the Euclidean and semi-Euclidean spaces, i.e., The generalized vector product is defined as follows: Let {e 1 , e 2 , e 3 } be the orthonormal basis of this space; then, its Clifford algebra Cl(E 3 (α, β)) = Cl p,q has the basis {1, e 1 , e 2 , e 1 e 2 }, where p + q = 3, e 2 1 = −α, e 2 2 = −β and e 1 e 2 = −e 2 e 1 [18].

Definition 2. A dual number s has the form s
is the set of dual numbers, and it forms an Abelian ring with unity and divisors of zero [30].
On the other hand, the set . The Lorentzian scalar product of S and R is defined as follows: Dual space D 3 , together with the scalar product, is called the dual Lorentzian space and denoted by D 3 1 . For S = 0, the norm S of S is defined by S = S, S . A dual vector S = S + εS * is called a dual-space-like vector if S, S > 0 or S = 0, a dual time-like vector if S, S < 0, and a dual null vector if S, S = 0 for S = 0 [30,31]. The dual angle between two lines is the dual number θ = θ + εd, where θ is the angle between two lines about the common normal and d is the distance between them along their common normal. Some trigonometric and hyperbolic functions with dual variables can be written as follows: (2) These definitions arise from the Taylor series expansion of the function f (θ + εd) about θ, which yields [5,30,31] .
The set of all G-orthogonal matrices under matrix multiplication operation is called a rotation group in E 3 (α, β) [18,20].
A rotation about the origin can be given with the equation of Rx = y, which is a 3 × 3 G-orthogonal matrix and x ∈ E 3 (α, β). The generalized Cayley formula is defined as where C is a G-skew symmetric matrix. By using the G-Cayley formula, a G-orthogonal matrix can be obtained by a G-skew symmetric matrix C, which is defined in Definition 3 [20].

Definition 5.
The set of all generalized orthogonal matrices with |R| = 1 under the operation of matrix multiplication forms a group denoted as SO(α, β))(3), called a rotation group in E 3 (α, β) [18]. Definition 6. The set of displacements of the space E 3 (α, β) is an algebraic group. If T 1 : F → M 1 and T 2 : The matrix representation of the combination of these two displacements is where 0 is a row vector (0, 0, 0). This representation is called a homogeneous transformation. The set of homogeneous transformations in E 3 (α, β) is as follows: Although the displacement is not a linear transformation, the homogeneous transformation g = (R, n) is a linear transformation.
The tangent operators of G − H(4) must satisfy the relation where A = · RR −1 is the 3 × 3 angular velocity matrix and s = −An + · n is its threedimensional linear velocity vector. The tangent direction of the motion g(t) at t = t 0 is defined by · g(t 0 ).
Similarly, the condition for SO(α, β)(3) to be tangent operators is obtained from the relation R T R = in E 3 (α, β). Differentiating both sides, we obtain · R T R + R T · R = 0 which can be written as The last equation shows that R T · R = Φ is a skew-symmetric matrix that is called the angular velocity matrix of the rotation R(t) in E 3 (α, β). If we can calculate that The solution of this matrix-differential equation gives a one-parameter group of rotations as follows [23]: In the following section, the translation matrix between the coordinate transformations will be obtained and the changes in the screw motion depending on the coordinate transformations will be examined.

Coordinate Transformations of a Screw in E 3 (α, β)
The tangent operator of g is L = (A, s), where A is 3 × 3 G-skew symmetric matrix and s is a three-dimensional vector. Let the vector a correspond to matrix A. Then, the six-dimensional vector l = (a, s) is called a screw. In this section, we will analyze the change in the coordinates of the screw depending on the change in reference frame for the motion.
Given fixed F and moving frames M, let a motion be described by displacement g(t) : F → M, with parameter t. The displacements T : F → F ı and T : M → M ı define g ı (t) = Tg(t)T −1 . Then, the motion can be given as g ı (t) : F ı → M ı . The tangent operator L ı of g ı (t) is obtained from the tangent operator L of g(t) as follows: If we rewrite this equation more clearly for displacement T = (R, n), we obtain From the G-skew symmetric matrix A ı = RAR −1 , we can obtain the vector a ı . From the relationship between the generalized cross-product and G-skew symmetric matrices, This transformation can be written as Here, T is a G-skew symmetric matrix corresponding to vector n, and 0 is the 3 × 3 square zero matrix. Equation (4) describes how the screw changes with the change in reference frames.
In the following section, it will be examined how the screw motion, which changes depending on the coordinate transformations, changes the components that make up the screw, and their specific cases will be discussed.

The Forms of a Screw in E 3 (α, β)
Let us consider a specific case of the change in coordinates in Equation (4). It is the screw motion obtained as a pure translation g = (I, n) of the fixed and moving frames. Thus, l = (a, s) transforms to l ı = (a, n × G a + s).
This implies that the linear velocity s of the screw changes depending on the changes in the origin of the reference frames.
Linear velocity vector a of l = (a, s) can be split into two components: one that is parallel to a and another that is orthogonal to a. Let pa be a parallel component; then, p is obtained as follows: The orthogonal part is s − pa. Here, we can find a point k on the screw axis such that The solution of the above equation is Let us change the reference frames by the translation g = (I, −k); then, l = (a, s) is transformed to the screw l ı = (a, −k × G a + s). Since k × G a orthogonal to a, hence −k × G a + s = pa.
Therefore, in this reference frame, the screw has the form l ı = (a, pa). This implies that in the new frame, the angular and linear velocity vectors are in the same direction of a. It is clear that the body is moving in a screw motion at time t. The axis of this screw motion is line d = k + ta, called the instantaneous screw axis in E 3 (α, β). The general screw l = (a, s) can be written in the standard form using equation s = k × G a + pa l = (a, k × G a + pa) The line d = k + ta and constant p are the axis and the pitch of screw l, respectively. When p = 0 in Equation (7), we encounter a critical case of the screw where there is no velocity vector. Let the screw be rewritten: There is an angular velocity and no linear velocity relative to d = k + ta. Thus, the line l can state the screw, and the components of the l are known as the Plücker coordinates of the line.
When p = 0, and the screw l = (0, s), we obtain another special case where the screw degenerates into a line screw. This means that the motion with no angular velocity is pure translation.
In the generalized scalar product in Equation (1), we have four cases, depending on the selection of the signs of α and β, which are α,β > 0, α > 0 and β < 0, α < 0 and β > 0, α, β < 0. However, from the definitions of the scalar product, the last three cases are the same. Since all possible selections of α and β can be covered by two conditions α, β > 0 and α > 0, β < 0, in the following section, we will consider only these two cases.
In the following section, a dual orthogonal matrix of an orthogonal matrix will be obtained and this matrix will be used to describe screw motion.

Generalized Dual Orthogonal Matrices and Screw Motion in E 3 (α, β)
Now, we wish to simplify the screw coordinate transformations via matrix multiplication and the properties of dual numbers. Let l = (w, s) and l ı = (a ı , s ı ) be screws associated with the dual vectors a = a + εs and a ı = a ı + εs ı , respectively. Thus, the screw transformation in Equation (4) can be rewritten as and T is a G-skew symmetric matrix.

Theorem 1. Dual matrix R is a G-dual orthogonal matrix.
Proof. Let dual matrix R be given as in (9), with R as a G-orthogonal matrix and T as a G-skew symmetric matrix. Then, since Case I (α > 0 and β > 0): By using Equation (8), the G-orthogonal matrix related to the screw displacement has its angle θ about the axis y, and a distance n along axis y, and we can obtain In addition, the screw displacement has its angle φ about the axis x, and a distance n along axis x The dual orthogonal transformation Case II (α > 0 and β < 0): In this case, we have two subcases where the x-axis and y-axis are timelike or spacelike.
Case IIa Let x-axis and y-axis be timelike; then, we have Case 1.

Case IIb
Let the x-axis and y-axis be spacelike; then, where In the same way, the screw displacement has its angle φ about the axis y, and a distance n along axis y Here, Additionally, the dual orthogonal transformation D = Y X In Case I and Case II, D is used to refer to the coordinate frames in successive members of a manipulator. Example 1. Case 1. Let α = 1 and β = 2. Let the screw displacement have its elliptical angle θ = 0.6435011 about the axis y, and a distance n along axis y. Substituting Equation (9) into Equation (8), we obtain the G-dual orthogonal matrix as where G-skew symmetric matrix T corresponding to translation vector n = (0, 1, 0) is and for α = 1 and β = 2, the matrix is Note that R T R = and R T R = are the G-orthogonal matrix and G-dual orthogonal matrix, respectively. where G-skew symmetric matrix T corresponding to the translation vector n = (0, 1, 0) is and for α = 1 and β = −2, the matrix is Again, note that R T R = and R T R = are the G-orthogonal matrix and G-dual orthogonal matrix, respectively.
In the following section, the exponential equation of the obtained dual orthogonal matrix for rotation given in (2) will be obtained for the screw motion, and the screw components will be defined by this exponential equation.

Generalized Dual Exponential Matrix in E 3 (α, β)
Let R(t) be a G-dual orthogonal matrix parametrized with t; then, we have the equation R T R = . Differentiating both sides, we obtain · ( R) T R + R T · R = 0 which can be written as The last equation shows that R T · R = Φ is a skew-symmetric matrix in E 3 (α, β). Here, We obtain G-dual skew-symmetric matrix Ω, called the dual angular velocity matrix in E 3 (α, β). Let w be a dual vector corresponding to G-dual skew symmetric matrix Ω in E 3 (α, β). Then, it can be written as w = w + εw * . Let us split w * into components parallel and perpendicular to w, as in Section 3, written as w * = k × G w + pw. Here, Thus, dual vector w is written as where 1 + εp is a dual scalar. The dual vector w + εk × G w describes an instantaneous screw axis and p is the pitch of the screw motion defined by R(t).
Dual angular velocity Ω is a constant because R T R = . The solution of the matrix differential equation, R(t) is a screw motion; its axis w + εk × G w is obtained via Ω. We have two cases where k × G w is either equal to zero or not.
Now, assume that R(t) is the exponential Ω. Thus, Equation (16) defines the G-dual orthogonal matrix R(t) = R(t) + εpt ΩR(t). Consequently, this transformation describes a screw motion with rotation about the angle φ(t) = ωt, with axis w, (ω = w ) and the amount n(t) = ptw. Moreover, ω is the angular velocity, s = pω is its linear velocity, and p is the pitch of the motion.
The above equation defines a displacement affiliated with a screw motion with pitch p and axis w + εk × G w in E 3 (α, β).

Conclusions
Screw motion has been studied by many authors in Euclidean and Lorentzian spaces. In this paper, screw motion is studied in generalized space. Thus, the theory of screw motion in the three-dimensional Euclidean and Lorentzian spaces is generalized into generalized space. The results obtained in this study reveal that the theory of screw motion is independent of the choice of the metric tensor in the Euclidean and Lorentzian spaces.
As is known, screw motion is a combination of a rotation around an axis and translational along the same axis. The rotation part of the screw motion is circular rotation in Euclidean space and hyperbolic rotation in Lorenz space. However, the rotation part of the screw motion in generalized space can be circular rotation, hyperbolic rotation, elliptical rotation, etc., depending on the choice of alpha and beta. Thus, this gives the opportunity to study a space where rotation motion is richer. The first case of the example in this study is the screw motion whose rotation part is elliptical rotation, and the second case is the screw motion whose rotation part is hyperbolic rotation. As mentioned briefly in the Introduction, generalized space has many application areas, such as mechanics, electromagnetism, and astrophysics.
The study of screw motion in the more general space is more difficult and complex than in the Euclidean and Lorentzian. Nevertheless, this study will make an important contribution to future studies to better understand and interpret screw motion. Screw motion in other scalar product spaces is also an open problem for further study.