Next Article in Journal
Enhanced Three Layer Hybrid Clustering Mechanism for Energy Efficient Routing in IoT
Previous Article in Journal
Quasi-All-Passive Thermal Control System Design and On-Orbit Validation of Luojia 1-01 Satellite
Previous Article in Special Issue
HEALPix-IA: A Global Registration Algorithm for Initial Alignment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Design and DOF Analysis of a Novel Compliant Parallel Mechanism for Large Load

1
Key Laboratory of Electronic Equipment Structure Design, Ministry of Education of China, Xidian University, Xi’an 710071, China
2
College of Information Science and Technology, Shihezi University, North 4th Road, Shihezi 832003, China
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(4), 828; https://doi.org/10.3390/s19040828
Submission received: 20 December 2018 / Revised: 9 February 2019 / Accepted: 13 February 2019 / Published: 17 February 2019
(This article belongs to the Special Issue Sensors for MEMS and Microsystems)

Abstract

:
The degree of freedom (DOF) and motion characteristics of a kind of compliant spherical joint were analyzed based on the screw theory, and a new design scheme for force-inversion of the compliant spherical joint was proposed in this paper. A novel type of six DOF compliant parallel mechanism (CPM) was designed based on this scheme to provide a large load capacity and achieve micrometer-level positioning accuracy. The compliance matrix of the new type of CPM was obtained through matrix transformation and was then decomposed into its generalized eigenvalues. Then, the DOF of the mechanism was numerically analyzed based on the symbolic formulation. The finite element analysis model of the compliant parallel mechanism was established. The static load analysis was used to verify the large load capacity of the mobile platform. By comparing the deformation obtained by the compliance matrix numerical method with the deformation obtained by the finite element method, the correctness of the compliance matrix and the number of the DOF of the CPM was verified.

1. Introduction

The motion of a complaint mechanism is caused by the elastic deformation of flexure elements when bearing loads [1,2,3,4]. Compared with the traditional rigid-body mechanism, it has the advantages of no friction, no lubrication required, compact structure, integrated molding, no assembly and high precision. Complaint mechanisms have been widely used in various precision instruments [5,6] including nanomanipulator [7], micropositioner [8], optical scanning mirrors [9], scanning transmission X-ray microscopy [10] and sensors [11], etc. However, the compliant mechanism usually fails to bear large loads. In the current research, there are few published results taking the load into account. In the process of judging the DOF of the compliant mechanism, the boundary between the DOF and the constraint is different from that of the rigid member, which can be directly calculated by the formula of the DOF; therefore, the DOF determination of the compliant mechanism is challenging in the analysis and synthesis for compliant mechanisms.
The DOF of a compliant mechanism is decided by geometric parameters, material properties and load. Currently, there is less research on the load capacity of the compliant mechanism, and the compliant mechanism is weak in bearing large loads due to the structural characteristics of the compliant hinge in the mechanism. Shi et al. presented workspace optimization of MEMS flexure-based hexapod nanopositioner previously built by the National Institute of Standards and Technology (NIST) and presented an analytical formulation and a search algorithm to determine the workspace of the flexure-based parallel mechanisms. A novel adaptive genetic algorithm was developed to conduct the single and bi-objective optimization for maximum translational and rotational workspace [12]; Shi et al. covered the kinematic modeling of a flexure-based, hexapod nanopositioner and calculated the actuation data for a set of commands for decoupled and coupled motions and obtained the Jacobian matrix of the mobile platform for the controller to calibrate the precision of the nanopositioner; however, it can hardly bear a considerable load [13]. Brouwer presented a precision MEMS-based six degrees-of-freedom (DOFs) manipulator. The purpose of the manipulator is to position a small sample (10 μm × 20 μm × 0.2 μm) in a transmission electron microscope. A parallel kinematic mechanism with slanted leaf-springs was used to convert the motion of six in-plane electrostatic comb-drives into six DOFs at the end-effector, which can only carry a small sample of 10~50 micrograms [14,15].
At the beginning of the design, it is necessary to analyze the motion characteristics, in which the DOF analysis is of the most importance. The DOF of a traditional rigid-body mechanism can be directly calculated by Grubler–Kutzbach’s formula, but this fails to judge the DOF of compliant mechanisms due to the flexure joints. Howell et al. first proposed the pseudo-rigid-body model method to analyze the DOF of a planar compliant mechanism [16]. Deshmukh et al. used a pseudo-rigid-body model (PRBM) method to design a flexure-based compliant parallel (4-bar) mechanism for a linear translational motion actuated via a precision slide, compared the results with those from finite element analysis, and verified the PRBM theory through experiments [17]. Su analyzed the motion characteristics of the general complaint mechanism based on screw theory and provided an important method for guiding the qualitative design of flexure mechanisms [18,19,20]. However, these methods proposed by Su are qualitative in process of judging the DOF and do not consider the influence of the size factor on the number of DOF. Li et al. defined a new generalized mechanism and derived the three main formulas for calculating the maximum and minimum DOF with screw theory and verified the factors affecting the DOF [21].
Based on the screw theory, the degree of freedom of the compliant spherical joint was firstly analyzed in this paper. A novel six-DOF CPM able to bear significant loads based on the compliant spherical joint was designed [22,23,24,25]. This type of mechanism is a mechatronic device with six DOFs and high-precision positioning in the space, which enables the mechanism to accurately detect the position of six DOFs in space after each branch installation is picked up. Thus the mechanism is widely applied in sensors [26], assembly of MEMS devices and micro-systems, and precision measurement of micro-systems. Instead of the present pseudo-rigid-body model method [27] and the DOF criteria based on screw theory [28], the method of the eigenwrench and the eigentwist decomposition in the mechanism compliance matrix is used to judge the DOF of the compliant mechanism. Firstly, by using the matrix transformation, the compliance matrix of the whole mechanism is established [29,30,31]. Then, the DOF of the mechanism is numerically analyzed by solving the generalized eigenvalues and eigenvectors of the compliance matrix. Finally, the compliance matrix of the parallel mechanism and the number of DOF were validated by comparing the deformation obtained by the compliance matrix numerical method with that from finite element simulation.

2. Analysis of the Compliant Spherical Joint Based on Screw Theory

2.1. Compliance Matrix of an RCCS (Right Circular Cross-Section) Flexural Joint

Under the condition that the ratio of the length to the radius of the elongated cylindrical rod is less than five, the slender rod joint has a total of three DOFs, which are rotated by the three axes of the coordinate system ( θ x , θ y and θ z ).
An RCCS flexural joint [32,33,34,35,36] is shown in Figure 1a. The ratio of the minimum cross-sectional diameter t-min of the compliant joint to the maximum cross-sectional diameter t-max is less than or equal to 1/10, and the distance of the joint between ends is 2r. The joint can achieve three DOFs rotation in the three-dimensional space shown in Figure 1a, namely bending along the y-axis and z-axis and axial torsion along the x-axis. This type of flexural joint has higher precision than the slender compliant spherical joint [37].
The two rods in Figure 1 are the two ends of the two different parts connected to it, because the deformation of the two parts is negligible relative to the compliant joint, so they are not taken into account in the numerical analysis. In the finite element analysis, in view of the actual force of the RCCS compliant spherical joint in the single branch of the mechanism, the boundary condition is gravity, the fixed constraint is applied at one end of the compliant spherical joint, and the other end of the round end face is applied with a direction-constant force or a constant amount of torque.
When the finite element analysis software is used for the RCCS compliant spherical joint, the linear analysis module is adopted. The reason for not using a non-linear analysis module is the verification of compliance matrix of the mechanism is compared and analyzed by the numerical method based on the compliance matrix and deformation obtained by the finite element method. The force matrix added during the compliance matrix numerical method is added in one step, and the compliance matrix is not changed in the process (even if the force matrix is applied in several steps, the compliance matrix does not change for the numerical method). For the finite element method, the compliance/stiffness matrix is dynamically changing after each loading of the generalized force. Therefore, in the process of finite element analysis, the force matrix is not used in several steps in the nonlinear analysis module. The material of the mechanism is 60 steel of structural steel, and the elastic modulus of the material is E = 200 GPa. Poissons’ ratio is ʋ = 0.3, density is 7.85 g/cm3; element type is solid186; automatic meshing is adopted; and mesh density Relevance = 20. After mesh, the number of nodes is 7281 and the number of elements is 3588. As for mesh quality, the element quality average is 0.61, the aspect ratio average is 3.19 and the Jacobian ratio average is 1.02; the boundary conditions include a fixed constraint applied on the bottom surface of the base, gravity, a constant magnitude force along the x-axis applied at the center point of the mobile platform, and the simulation analysis is completed in one step; after finite element analysis, the compliant spherical joint has three DOFs, which are rotations around the three coordinate axes (θx, θy, and θz). The angle of rotation of the RCCS compliant spherical joint around the x, z, and y-axes is −0.88°~+ 0.88°, −0.88°~+0.88° and −0.57°~+0.57°, respectively.
According to the mechanics of elasticity, the compliance matrix of an RCCS flexural joint can be written as [38]:
C hinge = ( 0 0 0 c θ x M x 0 0 0 0 c θ y F x 0 c θ y M y 0 0 c θ z F y 0 0 0 c θ z M z c x F y 0 0 0 0 0 0 c y F y 0 0 0 c y M z 0 0 c z F y 0 c z M y 0 )
where c θ x M x = 32 π G ( 0 2 r 1 t ( x ) 4 d x ) , c θ y M y = c θ z M z = 64 π E ( 0 2 r 1 t ( x ) 4 d x ) , c x F y = 4 π E ( 0 2 r 1 t ( x ) 2 d x ) , c y F y = c z F z = 64 π E ( 0 2 r x 2 t ( x ) 4 d x ) + 7 E 6 G c x F y , c y M z = c z M y = 64 π E ( 0 2 r x t ( x ) 4 d x ) , c y M z = c θ z F y , and c z M y = c θ y F z . where E is the Youngs modulus and G is the shear modulus. t ( x ) can be expressed as:
{ t ( x ) = t + 2 [ r r 2 ( r x ) 2 ] , x ( 0 , r ] t ( x ) = t + 2 [ r r 2 ( x r ) 2 ] , x ( r , 2 r ]
Then, the corresponding material parameters are substituted into Equation (1) to obtain a three DOFs RCCS flexural joint with a 6 × 6 compliance matrix.

2.2. The Compliance Matrix of RCCS Serial Chain

As shown in Figure 2, one end of an RCCS flexural joint was fixed as the reference object 1, and the other end of the RCCS flexural joint can be regarded as the mobile body 3. A serial compliant mechanism consisting of two RCCS flexural joints connected through the rigid body 2 can be obtained [39,40]. The parts 1, 2 and 3 are rigid bodies in this case.
[ C i ] denotes the compliance matrix of the i-th flexure member. Within the elastic limit range, the deformation of a serial chain is a superposition of each flexure element’s deformation in the same coordinate system. For a serial mechanism, the relationship between the overall compliance matrix and other flexure member compliance matrices can be expressed as follows:
[ C ] = i = 1 n [ A d i ] [ C i ] [ A d i ] 1
where [ A d i ] is the transformation matrix transforming the compliance matrix of the i-th flexure element to the functional body coordinate system.
As shown in Figure 2, o-xyz is considered as the global coordinate system; the coordinate system o-x2y2z2 is fixed to the moving body. The compliance matrix of the flexure element under the coordinate systems o-x1y1z1 and o-x2y2z2 can be converted to the coordinate system o-xyz through the transformation matrices [ A d 1 ] and [ A d 2 ] . Where R 1 = R 2 = I , I is a 3 × 3 unit matrix, d 1 = ( l , 0 , 0 ) T , d 2 = ( 0 , 0 , 0 ) T . [ A d 1 ] and [ A d 2 ] can be written as:
[ A d 1 ] = [ R 1 0 D 1 R 1 R 1 ] , [ A d 2 ] = [ R 2 0 D 2 R 2 R 2 ]
According to Equation (3), the compliance matrix of a serial chain based on compliant joints can be expressed as:
C leg = [ A d 1 ] [ C 1 ] [ A d 1 ] 1 + [ A d 2 ] [ C 2 ] [ A d 2 ] 1

2.3. A combined Design of Serial Compliant Joints with a Function of Interchanging Press and Pull

When a slender rod is subjected to pressure, it exhibits a property completely different from the strength failure. As shown in Figure 3, the lower end of the slender rod is fixed, and the upper end is free. When the pressure is beyond the limit value of F c r , the slender bar can return to its original state after the force is removed. When the pressure F gradually exceeds this limit value, the slender bar will maintain the balance of the curve shape and cannot restore the original shape. The slender bar loses its balance of a linear shape and transitions to a curve balance called buckling phenomenon [41,42,43].
As shown in Figure 3b, the RCCS flexural joint is similar to the compression bar and tends to exhibit the above-mentioned buckling phenomenon. After the compression bar is destabilized, a small increase in pressure will cause a significant increase in bending deformation, at which point the compression bar has lost its ability to carry the load. Instability causes the failure of the compression bar, which can cause damage to the entire machine or the corresponding structure, resulting in irreparable damage.
To avoid the damage of the mechanism caused by the instability of the RCCS flexural joint due to excessive pressure, this paper adopted a new structure to change the pressure on the RCCS flexural joint to the tension shown in Figure 4a,b.
A novel serial RCCS flexural joint chain structure based on an RCCS flexural joint is illustrated in Figure 4a. Part 1 is the fixed end; Part 2 is the rigid rod connecting two RCCS flexural joints. Part 3 is the free end bearing various forces or torque; and Part 4 is the RCCS flexural joint. As shown in Figure 4b, when the two U-shaped free ends are subjected to the axial pressure F/2 along the coordinate system, the two ends of the RCCS flexural joint will be subjected to two forces of the same value F and opposite directions. In this way, the pressure applied to the RCCS flexural joint is converted into tension, so that the structure can bear a much larger axial pressure without an instability problem.
The new serial RCCS flexural joint chain compliance matrix that can bear a large pressure load is solved in the same way as the previous compliance matrix solution. According to Equations (1)–(4), substituting l = 230 mm into the Equation (5) obtains a new serial RCCS flexural joints branch compliance matrix C newleg :
C newleg = [ A d 1 ] [ C 1 ] [ A d 1 ] 1 + [ A d 2 ] [ C 2 ] [ A d 2 ] 1

2.4. Compliance Matrix of CPM Based on the New Serial RCCS Flexural Joint Branch Chain

As shown in Figure 5a, the CPM in this paper is mainly composed of six new serial RCCS flexural joint branch chains as shown in Figure 4a. The angle between the center of the adjacent compliant spherical joint in each of branches near the mobile platform is 60° and the angle between the single and vertical branch is α = 12° [44].
From Equation (6), the compliance matrix of a single new serial RCCS flexural joint branch can be calculated. According to the relationship between the compliance matrix and stiffness matrix K newleg = C newleg 1 , the stiffness matrix of a new serial RCCS flexural joint can be written as:
K newleg = C newleg 1
To facilitate the study of the motion of the platform on the entire CPM and the deformation and displacement of the mobile platform during the stress process, it is necessary to transform the stiffness matrix of each branch into the fixed coordinate system of the mobile platform. As shown in Figure 5b, the single-branched coordinate system o-x1y1z1 needs to be matrix transformed into the upper platform fixed coordinate system o-xyz.
The coordinate transformation matrix of each single-branched local coordinate system oi-xiyizi (i = 1, 2, …, 6) to the coordinate system o-xyz can be expressed as:
R i = R x ( α i ) R y ( β i ) R z ( γ i ) , i = 1 , 2 , , 6
d i = ( r 1 c ρ i , h , r 1 c ζ i ) T , i = 1 , 2 , , 6
where c ρ = cos ( ρ i ) ,   s ζ = sin ( ζ i ) , ρ and ζ respectively indicate the angle between the line which is composed of the origin of each local coordinate system projected on the xz plane of the global coordinate system; the origin of the global coordinate system; and the x-axis and z-axis of the global coordinate system; R x ( α i ) , R y ( β i ) and R z ( γ i ) represent the transformation matrix around the x-axis, y-axis and z-axis, respectively. They can be written as:
R x ( α i ) = [ 1 0 0 0 c α i s α i 0 s α i c α i ] , R y ( β i ) = [ c β i 0 s β i 0 1 0 s β i 0 c β i ] , R z ( γ i ) = [ c γ i s γ i 0 s γ i c γ i 0 0 0 1 ]
where c α i = cos ( α i ) ,   s α i = sin ( α i ) , c β i = cos ( β i ) ,   s β i = sin ( β i ) , c γ i = cos ( γ i ) ,   s γ i = sin ( γ i ) .
The geometric parameters of the CPM and the corresponding material property parameters are shown in Table 1. The Z-Y-X Euler angle coordinate transformation is applied. γi represents the angle at which the i-th branch rotates about the z-axis in the local coordinate system; βi represents the angle at which the i-th branch rotates about the y-axis in the local coordinate system; αi represents the angle at which the i-th branch rotates about the x-axis in the local coordinate system, i = 1, 2, …, 6. r 1 represents the radius of the dotted circle consisting of the each single-branch local coordinate system origin as shown in Figure 5c. The material and geometric parameters are designed according to the pre-planning according to the medium-scale prototype, and the material cost and fabricating cost are taken into account. The compliant spherical joint symmetry center section’s size of the circle of 2 mm is much smaller than the size of other components in the mechanism. The reason why it does not choose a smaller size is also based on the consideration of whether it can withstand large loads and the results of finite element simulation.
According to Equation (7), the stiffness matrix of the large load CPM can be calculated as:
K i = K newleg K = i = 1 6 [ A d i ] [ K i ] [ A d i ] 1 = ( 0 0 2.7298 × 10 6 9.2440 × 10 3 0 0 0 0 0 0 4.0877 × 10 5 9.2440 × 10 3 2.7298 × 10 6 0 0 0 0 2.8298 × 10 6 8.0849 × 10 8 0 0 0 0 0 0 2.1867 × 10 5 0 0 0 0 0 0 8.0849 × 10 8 2.7298 × 10 6 0 0 )
The compliance matrix of the large load CPM can be written as:
C = K 1 = ( 0 0 1.2621 × 10 4 4.2738 × 10 7 0 0 0 0 0 0 4.5732 × 10 6 0 1.2621 × 10 4 0 0 0 0 4.2738 × 10 7 0.0374 0 0 0 0 1.2621 × 10 4 0 2.4464 × 10 6 0 0 0 0 0 0 0.0374 1.2621 × 10 4 0 0 )

3. The DOF Analysis of CPM Based on Compliant Branch

The compliant mechanism usually consists of a series of compliant components. In the process of the DOF analysis, the difference between DOFs and constraints are indiscernible. It is different from rigid members which can be directly calculated using the analytical formula of DOFs. The DOF determination of the compliant mechanism is of difficulty in the actual process. In this paper, the eigenwrench and eigentwist decomposition method in the mechanism compliance matrix were mainly used to judge the DOF of the compliance mechanism.
Because the deformation of the compliant mechanism satisfies the small deformation hypothesis according to elastic theory, the relationship between the wrench W ^ and twist T ^ can be expressed as:
T ^ = C W ^ , W ^ = K T ^
T ^ = ( θ T δ T ) T = ( θ x θ y θ z δ x δ y δ z ) T
W ^ = ( f T m T ) T = ( f x f y f z m x m y m z ) T
where θ ,   δ ,   m ,   f are all composed of 3 × 1 vector elements, which represent rotational deformation, moving deformation, force, and torque respectively.
The compliance matrix C and stiffness matrix K are both 6 × 6 square arrays. If the stiffness matrix K is nonsingular, the compliance matrix C can be written as C = K 1 .
According to the calculation of the reciprocal product, the wrench W ^ work was done along the twist T ^ , which can be written as:
T ^ W ^ = T ^ Δ W ^ = f T · δ + m T · θ
where is the reciprocal product of the wrench and twist, and ∆ is the swap operator defined as:
Δ = ( 0 3 × 3 I 3 × 3 I 3 × 3 0 3 × 3 )
In a coordinate system Q , the wrench and twist are represented by W ^ and T ^ , respectively. In the new coordinate system Q , the wrench and twist are represented by W ^ and T ^ , respectively. The relationship between them can be expressed as:
T ^ = Ad g T ^ , W ^ = Ad g W ^
where Ad g is the adjoint transformation matrix:
Ad g = ( R 0 D R R )
As shown in Equation (18), the matrix R represents the 3 × 3 rotation matrix of the new coordinate system Q with respect to the coordinate system Q , D is the skew-symmetric matrix defined by the translational vector d, and 0 represents the 3 × 3 zero matrix.
It can be seen from the above equation that the compliance matrix and stiffness matrix in the new coordinate system can be obtained by the following equation
C = Ad g C Ad g 1 , W = Ad g W Ad g 1

3.1. Generalized Eigenvalue Decomposition of Compliance Matrix

For a compliant mechanism, its DOF space and constrained space can be determined by the eigenwrench and eigentwist decomposition of its compliance matrix [C] [45], which can be decomposed as follows:
C W ^ = c δ T 1 W ^
C T 1 T ^ = c θ T ^
where T 1 = ( 0 0 I 0 ) . I represents the 3 × 3 identity matrix, 0 is a zero matrix of 3 × 3 , and c δ and c θ are called the translational and rotational eigencompliances, respectively. Equations (20) and (21) can be derived from two constraint minimization problems.
The two generalized eigenvalue problems 20 and 21 can be solved with linear algebra. It is well known that each problem yields up to three nonzero eigenvalues c δ i and c δ i . Their corresponding eigenvectors are called the eigenwrench W ^ δ i and eigentwist T ^ δ i , which are written as:
W ^ δ i = ( f i m i ) , T ^ δ i = T 1 W ^ δ i = ( 0 f i ) , ( i = 1 , 2 , 3 )
T ^ θ i = ( θ i δ i ) , W ^ θ i = T 1 T ^ θ i = ( 0 θ i ) , ( i = 1 , 2 , 3 )
The above two generalized eigenvalue decompositions can be summarized as the following equation:
C ( W ^ δ 1 W ^ δ 2 W ^ δ 3 W ^ θ 1 W ^ θ 2 W ^ θ 3 ) = ( T ^ δ 1 T ^ δ 2 T ^ δ 3 T ^ θ 1 T ^ θ 2 T ^ θ 3 ) ( c δ 1 c δ 2 c δ 3 c θ 1 c θ 2 c θ 3 )
where W ^ δ i and T ^ δ i are the wrench–compliance axes and twist–compliance axes respectively. Physically, the six eigencompliances c δ i and c δ i are the magnitudes of translation and rotation along these compliance axes assuming that the directional vector of these axes are normalized.

3.2. Judging the DOF of the CPM by the Symbolic Formulation

The compliance matrix for the novel CPM can be obtained by C = K 1 , which can be written as:
C = K 1 = ( 0 0 1.2621 × 10 4 4.2738 × 10 7 0 0 0 0 0 0 4.5732 × 10 6 0 1.2621 × 10 4 0 0 0 0 4.2738 × 10 7 0.0374 0 0 0 0 1.2621 × 10 4 0 2.4464 × 10 6 0 0 0 0 0 0 0.0374 1.2621 × 10 4 0 0 )
where the units of the upper left, upper right, and lower left 3 × 3 blocks, are rad/N, rad/Nmm and mm/N, respectively.
According to Equations (20) and (21), the translational eigencompliances and rotational eigencompliances can be expressed as:
c δ = ( 1.0818 × 10 4 , 2.4464 × 10 6 , 1.0818 × 10 4 )   mm / N c θ = ( 4.2738 × 10 7 , 4.5732 × 10 2 , 4.2738 × 10 7 )   rad / Nmm
According to Equation (26), since the units are inconsistent, the DOF of the mechanism cannot be directly determined. A characteristic length should be selected in order to determine the DOF. Since the mechanism is a closed-chain type, selection criterion 1 is used to select the characteristic length value l c = r 1 = 72.55   mm . According to Reference [38], the characteristic length lc only affects the rotational eigencompliances, instead of the translational eigencompliances, c ˜ δ = c δ ,   c ˜ θ = l c 2 c δ . The rotation eigencompliances in Equation (26) are multiplied by l c 2 . Then the translational eigencompliances and rotational eigencompliances are respectively expressed as follows:
c ˜ δ = ( 1.0818 × 10 4 ,   2.4464 × 10 6 ,   1.0818 × 10 4 )   mm / N c ˜ θ = ( 2.1944 × 10 3 ,   2.3481 × 10 2 ,   2.1944 × 10 3 )   mm / N
According to Equation (27), the eigencompliance values are listed in ascending order as follows:
( c 1 = 2.3481 × 10 2 ) > ( c 2 = 2.1944 × 10 3 ) ( c 3 = 2.1944 × 10 3 ) > ( c 4 = 1.0818 × 10 4 ) ( c 5 = 1.0818 × 10 4 ) > ( c 6 = 2.4464 × 10 6 )
According to Equation (28), the calculation of the compliance ratio can be expressed as:
C R 1 = 1 ,   C R 2 =   9.3455 ×   10 2 ,   C R 3 =   9.3455 ×   10 2 C R 4 =   4.6071 ×   10 3 ,   C R 5 =   4.6071 ×   10 3 ,   C R 6 =   1.0419 ×   10 4
According to Equation (29), if the threshold is chosen to be ε = 0.001, obviously there are m relatively small eigencompliances, that is m = 1. According to DOF criterion 1, the DOF of the compliance mechanism is N = 6 − m = 5.
According to Equations (22) and (23), the eigenwrench and eigentwist of the mechanism can be written as:
( W ^ δ W ^ θ ) = ( 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 295.3095 1 0 0.0044 0 0 0 0 1 0 295.3095 1 0.0082 228.4822 0 1 )
( T ^ δ T ^ θ ) = ( 0 0 0 1 0 0.0044 0 0 0 0 1 0 0 0 0 227.8769 0 1 1 0 0 6.7474 × 10 4 0 296.1 0 1 0 0 0 0 0 0 1 296.2400 0 1.300 )
where the eigentwists represent the motion of the compliant mechanism when exerted to the eigenwrenches. If the eigencompliance is larger, the corresponding eigentwist represents the mobility of the flexure mechanism along the direction. The direction depends on the vector θ of the eigentwist when the normal value of θ is not equal to zero and otherwise depends on the vector δ. As for the CPM shown in Figure 5a, it can be seen from Equation (31) that the eigentwist corresponding to the five larger eigencompliances are T ^ δ i and T ^ θ i , so the mechanism has five DOFs: three rotational DOFs about the x, y and z-axes, and two translational DOFs about x and z-axes. According to the screw theory, T = ( S | S 0 ) = ( S | c × S + p S ) , where vector S denotes the direction of the twist line and is independent of origin; the vector S0 is origin-dependent. c denotes the radius vector of the twist line to the origin of the global coordinate system, and p is called the pitch of a screw given by p = S · S 0 S · S and c = S × S 0 S · S . The vector c = ( 0 0 0 ) T is located at the origin of the global coordinate, so its translation along and rotation around the global coordinate system x, y and z-axes can be expressed as: T δ x = ( 0 0 0 1 0 0 ) T , T δ y = ( 0 0 0 0 1 0 ) T , T δ z = ( 0 0 0 0 0 1 ) T , T θ x = ( 1 0 0 0 0 0 ) T , T θ y = ( 0 1 0 0 0 0 ) T , T θ z = ( 0 0 1 0 0 0 ) T .
For example, T δ x = ( 0 0 0 1 0 0 ) T denotes the translation along the x-axis of the global coordinate system. By calculation, p = 0 and c = ( 0 0 0 ) T can be obtained which denotes the pitch of a screw and is located at the origin of the global coordinate, respectively. As for T δ y = ( 0 0 0 0 1 0 ) T , since C R 6 = 1.0419 × 10 4 < ε , there is no DOF along the y-axis of the global coordinate system.
Moreover, T θ z = ( 0.0044 0 1 296.1 0 1.3 ) T denotes the rotation around the z-axis of the global coordinate system. By calculation, p = 1 and c = ( 296.1 0 2.3 ) T can be obtained, which denotes the pitch of a screw and the radius vector of the twist line to the origin of the global coordinate system, respectively.
Since each branch of the parallel mechanism does not yet include an actuating unit, there is one less DOF for the mobile platform along the y-axis. After adding the actuating unit, according to the parallel robot actuating mode and the symmetry of the parallel mechanism, the mobile platform can perform the displacement in the y-axis.
The characteristic length lc = r1 = 72.55 mm selected according to the selection criterion 1 in Reference [40] is not an exact value and is selected according to the simplified formula. The value of the characteristic length does not need high accuracy, and the characteristic length value within a certain range does not affect the judgment of the number DOF of the mechanism based on the symbolic formulation method, as shown in Figure 6.

4. Numerical and Finite Element Method Normal Variable Analysis of CPM

Two methods are used to obtain the deformation data of the CPM, one is the numerical method by applying a certain load, and the other is the finite element simulation with commercial software. The deformation data obtained by the two methods were analyzed and compared in this section. Since the CPM can achieve a high-precision motion, the load is static or quasi-static with a low speed, and the mobile platform is very slow, so there is no need to analyze other performances such as velocity and acceleration, and this paper mainly considers its static stiffness. One can transform the coordinate matrix transformation based on Equation (32):
K = i = 1 6 [ A d i ] [ K i ] [ A d i ] 1
As shown in Equation (32), K is the stiffness matrix of the mechanism obtained by the coordinate transformation method. According to Equation (12), the stiffness matrix can be inversely transformed to obtain the compliance matrix of the mechanism.
Applying a load matrix CW to the origin of the fixed coordinate system at the center of the upper surface of the mobile platform, the deformation of the compliant mechanism can be obtained numerically, as follows:
W = ( F x F y F z M x M y M z ) T
T = C W
For the pure force along the x -axis, the deformation is mainly the translation along the x -axis. The translation and rotation on the other axes are so small that they can be ignored. Therefore, for pure force, the amount of translation along the x-axis direction is mainly acquired.
In this paper, the material of the compliant spherical joint is high-quality carbon structural steel 60, with high strength, hardness and elasticity, and it has a tensile strength of σ b 675   MPa . The basic configuration of the finite element static analysis is consistent with the compliant spherical joint; According to the FEA result, when the force applied to the mobile platform along the x-axis reaches 32.59 N, the tension of the compliant spherical joint will exceed 675 MPa. Therefore, within the range of yield strength of the selected material, by applying a force in the range of 0 N ~ 32.59 N along the x-axis to the mobile platform, a total of 33 sets of numerical and finite element methods are used to obtain the deformation of the mobile platform, as shown in Figure 7.
In Figure 7, Δδx denotes the displacement deformation along the x-axis and μ represents the displacement deformation difference obtained by the two methods along the x-axis.
As shown in Figure 7a, in the range of yield strength, the magnitude of the force acting on the mobile platform in the x-direction is proportional to the displacement deformation of the mobile platform along this direction. It can be seen from Figure 7b, that the larger the force is, the larger the error between the numerical method and the finite element method is; moreover, the error increases linearly, and the maximum error of the analytical method relative to the numerical method is 6.89%.
The above data is an example of displacement deformation caused by applying a force along the x-axis to the mobile platform. The force applied to the z-axis and the torque around the x, y, and z-axes will correspondingly produce a displacement and angular deformation. After the branches are installed with the driving, according to the parallel robot driving mode and the symmetry of the parallel mechanism, the mobile platform can perform the y-direction displacement movement along the y-axis.
Under the load condition, the displacement obtained by the numerical method and analytical method has some deviations on some data points, but the error is less than 6.89% relative to the numerical method, which is within the acceptable range. Therefore, the results of the numerical and finite element methods on the displacement and angular deformation of the six DOFs parallel mechanism under load conditions agree with each other and satisfy the linear relationship.
Additionally, the load capacity of the CPM was analyzed by the FEA method in this paper. As shown in Figure 8, under the condition of existing gravity, a load was applied vertically downward along the y coordinate axis to the origin of the fixed coordinate system at the center of the mobile platform in the CPM. When the load was fixed, the coordinate system reached 3826.82 N; the equivalent stress of each component of the mechanism reached the tensile strength limit of 675 MPa for 60 steel materials at each compliant spherical joint. However, the load capacity of the CPM was weakened under the non-initial position, because after the motion of the mobile platform, the compliant spherical joint was not only subjected to the axial force but also to the bending moment. According to the finite element analysis, when the mobile platform moves 1.20 mm along the x or z-axis, the CPM mobile platform can withstand 733.82 N without exceeding the material elastic limit. According to the finite element simulation, the closer to the boundary of the workspace of the CPM, the smaller the bearing capacity of the CPM, and the bearing capacity can be maximized at its initial position. It can be seen that the CPM has a certain ability to withstand large loads in its workspace. This static simulation thus shows that the CPM is able to bear large loads.
Table 2 shows the corresponding deformation values obtained by finite element simulation analysis of some force/torque values taken randomly within the linear elastic range of the mechanism material. The deformation values measured in the table are the value outputs by the finite element software.

5. Conclusions

A novel CPM able to effectively improve the load capacity was designed by employing a force-inversed concept of the compliant spherical joint in this paper. This type of mechanism is recommended for high-precision, multi-degree-of-freedom and large-load applications, and is widely applied in sensors andassembly of micro-systems. The DOF was analyzed with the screw theory and symbolic formulation method. The DOF and structural stiffness were verified by finite element simulation. Some meaningful conclusions can be drawn as follows.
(1)
The existing six DOFs CPM mobile platforms fail to bear a considerable load because of the bucking of the flexural joints. However, the proposed novel mechanism will be able to bear large loads, as much as a 32.59 N force along the x and z-axes.
(2)
The DOF of the mechanism was numerically analyzed by the generalized eigenvalue and eigenvector decomposition of the compliance matrix for the compliance matrix symbolic formulation of the whole mechanism, and the number of the DOF is verified by the numerical and finite element method. Under the same load, the deformation error obtained by the two methods is less than 6.89% compared with the finite element method; in addition, the load and deformation increases linearly within the yield strength range of the material.

Author Contributions

The authors contributions are described as follows: conceptualization, X.D. and D.Z.; methodology, X.W.; the FEA model, Y.L.; formal analysis, X.W.; writing—original draft preparation, X.W.; writing—review and editing, W.D. and X.D.; supervision, D.Z.; project administration, Y.L.; funding acquisition, X.D.

Funding

This research was funded by Advanced Armament Research Project of 13th Five-year Plan under Grant. 30508040102 and the National Natural Science Foundation of China under Grant. 51679398.

Acknowledgments

Acknowledge to Guimin Chen and Ying Zhang for their instructive discussions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Smith, S.T. Flexures: Elements of Elastic Mechanisms; Gordon and Breach Science Publishers: New York, NY, USA, 2000. [Google Scholar]
  2. Howell, L.L. Compliant Mechanisms; John Wiley & Sons: New York, NY, USA, 2013. [Google Scholar]
  3. Lobontiu, N. Compliant Mechanisms: Design of Flexure Hinges; CRC Press: Boca Raton, FL, USA, 2008. [Google Scholar]
  4. Patil, C.B. Robust Design of Selectively Compliant Flexure-Based Precision Mechanisms. Ph.D. Thesis, The University of Texas at Austin, Austin, TX, USA, 2008. [Google Scholar]
  5. Li, S.Z.; Yu, J.J.; Pei, X.; Su, H.-J.; Hopkins, J.B.; Culpepper, M.L. Type Synthesis Principle and Practice of Flexure Systems in the Framework of Screw Theory: Part III–Numerations and Synthesis of Flexure Mechanisms. In Proceedings of the ASME 2010 International Design Engineering Technical Conferences and Computers and Information in Engineering Conference, Montreal, QC, Canada, 15–18 August 2010; pp. 619–628. [Google Scholar]
  6. Smith, S.T. Foundations of Ultra-Precision Mechanism Design; Gordon and Breach Science Publishers: Philadelphia, PA, USA, 2003. [Google Scholar]
  7. Culpepper, M.L.; Anderson, G. Design of a Low-Cost Nano-Manipulator Which Utilizes a Monolithic, Spatial Compliant Mechanism. Precis. Eng. 2004, 28, 469–482. [Google Scholar] [CrossRef]
  8. Dagalakis, N.G.; Amatucci, E. Kinematic Modeling of a 6 Degree of Freedom Tri-Stage Micro-Positioner. In Proceedings of the American Society for Precision Engineering 16th Annual Meeting, Crystal City, VA, USA, 10–15 November 2001; pp. 200–203. [Google Scholar]
  9. Aguirre, A.D.; Hertz, P.R.; Chen, Y.; Fujimoto, J.G.; Piyawattanametha, W.; Fan, L.; Wu, M.C. Two-Axis MEMS Scanning Catheter for Ultrahigh Resolution Three-Dimensional and En Face Imaging. Opt. Express. 2007, 15, 2445–2453. [Google Scholar] [CrossRef] [PubMed]
  10. Henein, S.; Frommherz, U.; Betemps, R.; Kalt, H.; Ellenberger, U.; Flechsig, U.; Raabe, J. Mechanical Design of a Spherical Grating Monochromator for the Microspectroscopy Beamline Pollux at the Swiss Light Source. In Proceedings of the AIP Conference, Daegu, Korea, 28 May–2 June 2006; pp. 643–646. [Google Scholar]
  11. Wang, Z.; Chen, L.; Sun, L. An Integrated Parallel Micromanipulator with Flexure Hinges for Optical Fiber Alignment. In Proceedings of the IEEE International Conference on Mechatronics and Automation, Harbin, China, 5–8 August 2007; pp. 2530–2534. [Google Scholar]
  12. Shi, H.L.; Duan, X.; Su, H.-J. Optimization of the workspace of a MEMS hexapod nanopositioner using an adaptive genetic algorithm. In Proceedings of the IEEE International Conference on Robotics and Automation (ICRA), Hong Kong, China, 31 May–7 June 2014; pp. 4043–4048. [Google Scholar]
  13. Shi, H.; Su, H.-J.; Dagalakis, N.; Kramar, J.A. Kinematic modeling and calibration of a flexure based hexapod nanopositioner. Precis. Eng. 2013, 37, 117–128. [Google Scholar] [CrossRef]
  14. Brouwer, D.M.; Soemers, H. 6-DoFs MEMS-based precision manipulators. Mikroniek 2008, 48, 5–13. [Google Scholar]
  15. Brouwer, D.M.; Jong, B.R.D.; Soemers, H.M.J.R. Design and modeling of a six DOFs MEMS-based precision manipulator. Precis. Eng. 2010, 34, 307–319. [Google Scholar] [CrossRef]
  16. Howell, L.L.; Midha, A. Determination of the degrees of freedom of compliant mechanisms using the pseudo-rigid-body model concept. In Proceedings of the Ninth World Congress on the Theory of Machines and Mechanisms, Milan, Italy, 1–3 September 1995; pp. 1537–1541. [Google Scholar]
  17. Deshmukh, B.; Pardeshi, S.; Mistry, R. Development of a Four Bar Compliant Mechanism using Pseudo Rigid Body Model (PRBM). Procedia Mater. Sci. 2014, 6, 1034–1039. [Google Scholar] [CrossRef]
  18. Su, H.J. Mobility Analysis of Flexure Mechanisms via Screw Algebra. J. Mech. Rob. 2011, 3, 041010. [Google Scholar] [CrossRef]
  19. Su, H.J.; Yue, C. Type synthesis of freedom and constraint elements for design of flexure mechanisms. Mech. Sci. 2013, 4, 263–277. [Google Scholar] [CrossRef]
  20. Su, H.J.; Tari, H. On Line Screw Systems and Their Application to Flexure Synthesis. J. Mech. Rob. 2011, 3, 011009. [Google Scholar] [CrossRef]
  21. Li, T.; Deng, H.; Zhang, L. Mobility Analysis of Generalized Mechanisms via Screw Algebra. Mech. Mach. Sci. 2017, 408, 581–596. [Google Scholar]
  22. Yu, J.; Li, S.; Su, H.J.; Culpepper, M.L. Screw theory based methodology for the deterministic type synthesis of flexure mechanisms. J. Mech. Rob. 2011, 3, 03100814. [Google Scholar] [CrossRef]
  23. Su, H.J.; Dorozhkin, D.V.; Vance, J.M. A Screw Theory Approach for the Conceptual Design of Flexible Joints for Compliant Mechanisms. J. Mech. Rob. 2009, 1, 0410098. [Google Scholar] [CrossRef]
  24. Wei, D.; Sun, L.N.; Du, Z.J. Design of a precision compliant parallel positioner driven by dual piezoelectric actuators. Sens. Actuators A Phys. 2007, 135, 250–256. [Google Scholar]
  25. Su, H.J.; Tari, H. Realizing Orthogonal Motions with Wire Flexures Connected in Parallel. J. Mech. Des. 2010, 132, 121002. [Google Scholar] [CrossRef]
  26. Liang, Q.; Zhang, D.; Wang, Y.; Ge, Y. Design and Analysis of a Novel Six-Component F/T Sensor based on CPM for Passive Compliant Assembly. Meas. Sci. Rev. 2013, 13, 253–264. [Google Scholar] [CrossRef]
  27. Panda, B.; Dutta, A. Design of a partially compliant crank rocker mechanism using Ionic Polymer Metal Composite for path generation. Mater. Des. 2010, 31, 2471–2477. [Google Scholar] [CrossRef]
  28. Huang, Z.; Liu, J.F.; Zeng, D.X. A general method for analyzing the degree of freedom based on screw theory. Sci. China 2009, 39, 84–93. (In Chinese) [Google Scholar]
  29. Selig, J.M.; Ding, X. A screw theory of static beams. In Proceedings of the IEEE/RSJ International Conference on Intelligent Robots and Systems, Maui, HI, USA, 29 October–3 November 2001; pp. 312–317. [Google Scholar]
  30. Su, H.-J.; Shi, H.; Yu, J.J. Analytical Compliance Analysis and Synthesis of Flexure Mechanisms. In Proceedings of the ASME 2011 International Design Engineering Technical Conferences and Computers and Information in Engineering Conference, ASME, Washington, DC, USA, 29–31 August 2011; pp. 169–180. [Google Scholar]
  31. Shi, H.L.; She, Y.; Duan, X.C. A Structural Model for Loading Analysis of a Hexapod Platform. In Proceedings of the 2015 world congress on Aeronautics, Nano, Bio, Robotics, and Energy (ANBRE15), Incheon, Korea, 25–28 August 2015. [Google Scholar]
  32. Jia, M.; Jia, R.P.; Yu, J.J. A Compliance-Based Parameterization Approach for Type Synthesis of Flexure Mechanisms. J. Mech. Rob. 2015, 7, 031014. [Google Scholar] [CrossRef]
  33. Lobontiu, N.; Paine, J.S.N.; Garcia, E.; Goldfarb, M. Corner-Filleted Flexure Hinges. J. Mech. Des. 2000, 123, 346–352. [Google Scholar] [CrossRef]
  34. Furqan, M.; Alam, M.N. Finite Element Analysis of a Stewart Platform using Flexible Joints. In Proceedings of the 1st International and 16th National Conference on Machines and Mechanisms (iNaCoMM2013), IIT Roorkee, India, 18–20 December 2013. [Google Scholar]
  35. Chen, G.M.; Jia, J.Y.; Li, Z.W. Right-circular corner-filleted flexure hinges. In Proceedings of the 2005 IEEE International Conference on Automation Science and Engineering, Edmonton, AL, Canada, 1–2 August 2005; pp. 249–253. [Google Scholar]
  36. Lobontiu, N.; Garcia, E.; Hardau, M.; Bal, N. Stiffness characterization of corner-filleted flexure hinges. Rev. Sci. Instrum. 2004, 75, 4896–4905. [Google Scholar] [CrossRef]
  37. Dong, W.; Du, Z.; Sun, L. Conceptional design and kinematics modeling of a wide-range flexure hinge -based parallel manipulator. In Proceedings of the IEEE International Conference on Robotics and Automation, Barcelona, Spain, 18–22 April 2005. [Google Scholar]
  38. Lobontiu, N.; Paine, J.S.N. Design of circular cross-section corner-filleted flexure hinges for three-dimensional compliant mechanisms. J. Mech. Des. 2002, 124, 479–484. [Google Scholar] [CrossRef]
  39. Shi, H. Modeling and Analysis of Compliant Mechanisms for Designing Nano-positioners. Ph.D. Thesis, Ohio State University, Columbus, OH, USA, 2013. [Google Scholar]
  40. Zhang, Y. Kinematics analysis and synthesis of space agencies and compliant mechanisms. Ph.D. Thesis, Beijing University of Posts and Telecommunications, Beijing, China, 2015. [Google Scholar]
  41. Miller, J.T.; Su, T.; Pabon, J.; Wicks, N.; Bertoldi, K.; Reis, P.M. Buckling of a thin elastic rod inside a horizontal cylindrical constraint. Extreme Mech. Lett. 2015, 3, 36–44. [Google Scholar] [CrossRef]
  42. Su, T.X.; Wicks, N.; Pabon, J.; Bertoldi, K. Mechanism by which a frictionally confined rod loses stability under initial velocity and position perturbations. Int. J. Solids Struct. 2013, 50, 2468–2476. [Google Scholar] [CrossRef]
  43. Shams, A.; Aureli, M.; Porfiri, M. Nonlinear buckling of a spherical shell embedded in an elastic medium with imperfect interface. Int. J. Solids Struct. 2013, 50, 2310–2327. [Google Scholar] [CrossRef]
  44. Su, H.-J.; Shi, H.; Yu, J.J. A Symbolic Formulation for Analytical Compliance Analysis and Synthesis of Flexure Mechanisms. J. Mech. Des. 2012, 134, 51009. [Google Scholar] [CrossRef] [PubMed]
  45. Zhang, Y.; Su, H.-J.; Liao, Q. Mobility criteria of compliant mechanisms based on decomposition of compliance matrices. Mech. Mach. Theory 2014, 79, 80–93. [Google Scholar] [CrossRef]
Figure 1. The model and geometry of an RCCS compliant spherical joint: (a) model of an RCCS; (b) geometry of an RCCS compliant spherical joint.
Figure 1. The model and geometry of an RCCS compliant spherical joint: (a) model of an RCCS; (b) geometry of an RCCS compliant spherical joint.
Sensors 19 00828 g001
Figure 2. A serial chain of two RCCS compliant spherical joints.
Figure 2. A serial chain of two RCCS compliant spherical joints.
Sensors 19 00828 g002
Figure 3. Buckling phenomenon of the slender rod: (a) the slender rod buckling phenomenon; (b) pressure on an RCCS flexural joints and (c) pull on RCCS flexural joints.
Figure 3. Buckling phenomenon of the slender rod: (a) the slender rod buckling phenomenon; (b) pressure on an RCCS flexural joints and (c) pull on RCCS flexural joints.
Sensors 19 00828 g003
Figure 4. A novel serial RCCS flexural joint chain structure based on an RCCS flexural joint: (a) a serial chain of two compliant spherical joints with the force-inverted structure; (b) compliant spherical joint force-inverted structure.
Figure 4. A novel serial RCCS flexural joint chain structure based on an RCCS flexural joint: (a) a serial chain of two compliant spherical joints with the force-inverted structure; (b) compliant spherical joint force-inverted structure.
Sensors 19 00828 g004
Figure 5. Views of the CPM: (a) three-dimensional model of the CPM; (b) left view of the CPM; (c) top view of the CPM.
Figure 5. Views of the CPM: (a) three-dimensional model of the CPM; (b) left view of the CPM; (c) top view of the CPM.
Sensors 19 00828 g005
Figure 6. The relationship between the number of DOF and the value of the characteristic length.
Figure 6. The relationship between the number of DOF and the value of the characteristic length.
Sensors 19 00828 g006
Figure 7. Deformation of the mobile platform that is obtained by analytical and numerical methods: (a) displacement along the x-axis under axis under load; (b) displacement difference along the x load.
Figure 7. Deformation of the mobile platform that is obtained by analytical and numerical methods: (a) displacement along the x-axis under axis under load; (b) displacement difference along the x load.
Sensors 19 00828 g007
Figure 8. The load nephograms of the CPM with bearing gravity.
Figure 8. The load nephograms of the CPM with bearing gravity.
Sensors 19 00828 g008
Table 1. Parameters of the 6-DOF CPM.
Table 1. Parameters of the 6-DOF CPM.
VariableValueUnit
E200GPa
ν0.3mm
t2mm
r5mm
l240mm
r172.55mm
α10°
α20°
α30°
α40°
α50°
α60°
β130°
β290°
β3150°
β4210°
β5270°
β6330°
γ1102°
γ2102°
γ3102°
γ4102°
γ5102°
γ6102°
ρ130°
ρ290°
ρ3150°
ρ4210°
ρ5270°
ρ6330°
ζ1150°
ζ2180°
ζ3210°
ζ4240°
ζ5270°
ζ6300°
h35.78mm
Table 2. Partial force/torque values corresponding to deformation values.
Table 2. Partial force/torque values corresponding to deformation values.
Fz (N)Δδz (mm)Mx (Nmm)Δθx (°)My (Nmm)Δθy (°)Mz (Nmm)Δθz (°)
5.540.225536.680.171474.420.444281.420.13
9.910.406619.620.201539.200.464958.330.15
14.090.576876.590.211675.830.505547.470.17
14.450.577302.790.221915.120.575839.590.18
17.790.717394.970.231945.700.586254.180.19
22.210.898110.670.252396.660.727216.390.22
24.090.978962.770.272435.860.737534.120.23
25.711.039247.490.282667.900.807888.290.24
26.441.069405.110.292966.740.908155.980.25
29.621.039751.750.303047.800.918976.380.27
31.031.2510,003.810.313070.200.929150.620.28
32.341.3010,125.490.313593.101.089909.630.30

Share and Cite

MDPI and ACS Style

Wu, X.; Lu, Y.; Duan, X.; Zhang, D.; Deng, W. Design and DOF Analysis of a Novel Compliant Parallel Mechanism for Large Load. Sensors 2019, 19, 828. https://doi.org/10.3390/s19040828

AMA Style

Wu X, Lu Y, Duan X, Zhang D, Deng W. Design and DOF Analysis of a Novel Compliant Parallel Mechanism for Large Load. Sensors. 2019; 19(4):828. https://doi.org/10.3390/s19040828

Chicago/Turabian Style

Wu, Xiaochuan, Yi Lu, Xuechao Duan, Dan Zhang, and Wenyao Deng. 2019. "Design and DOF Analysis of a Novel Compliant Parallel Mechanism for Large Load" Sensors 19, no. 4: 828. https://doi.org/10.3390/s19040828

APA Style

Wu, X., Lu, Y., Duan, X., Zhang, D., & Deng, W. (2019). Design and DOF Analysis of a Novel Compliant Parallel Mechanism for Large Load. Sensors, 19(4), 828. https://doi.org/10.3390/s19040828

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop