Abstract

This paper presents a novel two degrees-of-freedom planar parallel manipulator (PPM) designed for infinite-axis 3D printing, alongside tools for facilitating future design iterations. Unlike traditional gantry-supported designs used in infinite-axis 3D printing, which impose significant mass movement requirements, the examined new design prioritizes reducing overall weight to enhance speed potential at the cost of a reduced work area. In this innovative approach, the PPM effectively reduces weight by decoupling the motion of the hot end from that of the motor. Motors are attached to the frame, controlling a system of pulleys, and connecting arms to drive the hot-end's motion. Due to the length of the arms, the hot end will be unable to fully explore the entire printing plane. Verification of the angled PPM for 3D printing involved developing kinematic and dynamic equations, conducting finite element analysis on critical components, and testing a completed prototype. A metaheuristic optimization method was employed to derive optimal design parameters, focusing on optimizing the arm length of the connectors while maximizing dynamic performance. Considerations included the usable workspace and the angle between the connecting arm and end-effector. The final prototype validated the stability and rigidity of the PPM during movement, indicating its viability for 3D printing. The results presented in this paper demonstrate the capabilities of using an angled PPM in infinite 3D printing, providing fundamental knowledge crucial for future designs involving this innovative mechanism.

Graphical Abstract Figure
Graphical Abstract Figure
Close modal

1 Introduction

The development of new 3D printing technologies has been rapidly growing in the past few decades, largely due to the expiration of the patents surrounding the technology. Currently, ASTM recognizes seven distinct categories of 3D printing technologies that include binder jetting, directed energy deposition, material extrusion, material jetting, powered bed fusion, sheet lamination, and vat photopolymerization [1]. These categories were created to differentiate 3D printers based on the processes they use to create a 3D object. Each process uses different techniques that allow for the utilization of a wide range of materials and precision. The aerospace, automotive, food, healthcare, construction, textile/fashion, and electronic industries have all been disrupted by the advent of 3D printing [2]. Some marvels of 3D printing include custom medical implants [3], fabrication of soft robotics [4], automotive parts [5], 3D printed houses [6], 3D printed jewelry [7], and even 3D printed rockets [8].

In 2017, the company BlackBelt released the first commercially available 3D printer featuring an angled gantry design mounted on a rotating bed, as depicted in Figs. 1 and 2 [9]. This new design, referred to as a conveyor belt 3D printer, introduced the ability to continuously print objects of infinite length. The process of printing at an angle is functionally identical to traditional extrusion 3D printers. Instead of building layers up vertically, layers are stacked on top of each other at an angle as depicted in Fig. 2. Fast forward to 2021, and there are only two commercially available conveyor belt 3D printers—the BlackBelt and the Creality CR-30 [10]. Both printers rely on the same hot end mounted gantry design. The hot end being mounted to a gantry necessitates moving a significant mass during printing, posing limitations on printing speed. The lack of diversity with the design of commercially available belt 3D printers, points to a substantial opportunity for innovation.

Fig. 1
Blackbelt printer (Reprinted with permission from [9])
Fig. 1
Blackbelt printer (Reprinted with permission from [9])
Close modal
Fig. 2
Blackbelt printing a crankshaft (Reprinted with permission from [9])
Fig. 2
Blackbelt printing a crankshaft (Reprinted with permission from [9])
Close modal

Motivated by the aim to create a faster alternative, this paper proposes a new implementation of a planar parallel manipulator (PPM) as a weight reduction design to replace the conventional gantry design. This innovative approach aims to decrease the moving weight by eliminating the need for a cross beam and fixing the axis motors stationary to the frame. A system of pulleys and arms will facilitate the movement of the hot end. The primary objective of this paper is to validate the feasibility of utilizing an angled PPM and introduce a methodology for generating an optimized design. However, comparisons between types of conveyor belt 3D printers will not be investigated, as this is beyond the scope of the current research.

1.1 Literature Review.

Given that belt 3D printing, also known as infinite-axis 3D printing, is still in its infancy with very few commercial offerings, there are no peer-reviewed articles covering the design and abilities of these printers. Most of the development of these printers has been achieved by 3D printing enthusiasts who publish their final designs in an open-source format [1113]. The most relevant academic research available on this topic only details the capabilities of traditional selective compliance assembly robot arm (SCARA) and delta-style 3D printers, as those types of printers use parallel manipulators as a means of achieving motion. These articles provide insight into the governing principles of parallel manipulators used for 3D printing. Unfortunately, research surrounding the use of parallel manipulators for 3D printing focuses on SCARA/delta printers. There is no examination of a hybrid approach where bed movement is used [1418].

However, there are a series of papers published that detail the use of a PPM for the translational motion of an industrial-sized milling machine [1922]. The PPM detailed in these papers consists of a symmetric design whereby the end-effector is attached to two arms, which are attached to a carriage as shown in Fig. 3. The arms are attached using revolute joints and the carriage is attached to the frame via a prismatic joint. In the papers, this design was mounted in an upright orientation. This design can be angled and fitted with a hot end, enabling its adaptation for creating an infinite 3D printer.

Fig. 3
PPM design schematic
Fig. 3
PPM design schematic
Close modal

2 Design

To create a functioning infinite 3D printer, the PPM shown in Fig. 3 will need to be attached to the frame of the printer at an angle. In this design, the PPM will be attached at 45 deg with respect to the build plane. Figure 4 shows a computer-aided design (CAD) representation of the final design. The kinematic and dynamic derivations of this PPM are detailed in the following sections.

Fig. 4
CAD render of PPM infinite 3D printer
Fig. 4
CAD render of PPM infinite 3D printer
Close modal

2.1 Kinematics.

It is important to note that the kinematics are fully described regardless of the orientation of the frame. Therefore, previous work analyzing the kinematics of this PPM still applies. The kinematics of the PPM have been derived in full by Liu et al. and restated with respect to Fig. 3 [19]. The center of the end-effector, C, can be written as
(1)
The coordinate points of Bi (i = 1, 2) can be written as
(2)
The coordinate points of Pi (i = 1, 2) can be written as
(3)
The inverse kinematics for the manipulator can be written as [19]
(4)
(5)
where 2r represents the length of the end-effector and 2R is the width of the frame.

2.2 Inverse Dynamics.

Dynamic equations for an upright PPM have been previously established [19,20]. In these works, the authors independently derive these equations using two distinct methods—the Lagrangian and the virtual work principle. However, these derivations cannot be directly applied when the PPM is at an angle. The problem is fundamentally different now that gravity is influencing the mechanism outside the plane of motion.

To address this angled scenario, Kane's method will be applied to create the equations of motion. Kane's method offers the advantage of systematically analyzing a system without introducing excessive complexity, while also allowing for modifications even after the derivation is complete. This unique approach enables the PPM to be initially analyzed as an unconstrained system. Subsequently, additional constraints can be applied post-derivation to formulate the entire scenario, resulting in the finalized dynamic equations. For this analysis, the simplified unconstrained PPM in Fig. 5 will be used as a reference.

Fig. 5
Unconstrained schematic of the PPM
Fig. 5
Unconstrained schematic of the PPM
Close modal

2.2.1 Partial Velocities.

It is important to note that seen in Fig. 5, there exists a revolute joint at all points P that allows the arms to rotate about the n^3-axis and P1 is free to move along the n^2-axis. The generalized speeds can be written as
(6)

The transformation matrix between the N and A frame of reference is given in Table 1, while the transform between N and B is given in Table 2. It is assumed that the points on the system are massless and the bodies, A and B, have the same mass (m) and length (L). The center of mass for body A is located at Lcm. The center of mass for body B is located at LLcm. The partial velocities Vr and the partial angular velocities ωr for the PPM are recorded in Table 3.

Table 1

Transformation matrix between N and A

n^1n^2n^3
a^1cos(q2)sin(q2)0
a^2sin(q2)cos(q2)0
a^3001
n^1n^2n^3
a^1cos(q2)sin(q2)0
a^2sin(q2)cos(q2)0
a^3001
Table 2

Transformation matrix between N and B

n^1n^2n^3
b^1cos(q3)sin(q3)0
b^2sin(q3)cos(q3)0
b^3001
n^1n^2n^3
b^1cos(q3)sin(q3)0
b^2sin(q3)cos(q3)0
b^3001
Table 3

Partial angular velocities and partial velocities for Fig. 5 

r = 1r = 2r = 3
ωrA0n^30
ωrB00b^3
VrP1n^200
VrA¯n^2Lcma^20
VrP2n^2La^20
VrB¯n^2La^2(LLcm)b^2
VrP3n^2La^2Lb^2
r = 1r = 2r = 3
ωrA0n^30
ωrB00b^3
VrP1n^200
VrA¯n^2Lcma^20
VrP2n^2La^20
VrB¯n^2La^2(LLcm)b^2
VrP3n^2La^2Lb^2

2.2.2 Translational and Angular Accelerations.

The translational acceleration a and angular acceleration α for the elements on the PPM can be constructed by relating the accelerations of two points fixed on a rigid body. The accelerations for the PPM are given as follows:
(7)
(8)
(9)
(10)
(11)
(12)
(13)

2.2.3 Generalized Inertial Forces.

The generalized inertial forces can be derived using Kane's method and the previously derived accelerations and velocities [23]. The generalized inertial forces Fr* (r = 1, 2, 3) for body β in the system is
(14)
(15)
(16)
where Fr* is in the reference frame N for a rigid body β, a is the acceleration of the center of mass of β in N, M is the total mass of β, I is the central dyadic of β, R* is the inertial force for β in N, and T* is the inertia torque for β in N. The generalized inertial force for the system is
(17)
where N is the total number of bodies in the system.

2.2.4 Generalized Active Forces.

The final component of Kane's method is to include the generalized active forces on the system. For this scenario, the only active force that needs to be accounted for is the force due to gravity, keeping in mind the PPM is angled, and the active force used to move P1. The equation for the generalized active forces is
(18)

where k is the total number of active forces acting on the system.

2.2.5 Equation of Motion for Unconstrained Planar Parallel Manipulator.

With the generalized forces and active forces accounted for, the equations of motion for the unconstrained mechanism can be defined as
(19)

2.2.6 Equation of Motion for Constrained Planar Parallel Manipulator.

The equations of motion for the unconstrained system can be built upon to account for the addition of a constraint acting on the system. Figure 6 shows the new constraint placed upon P3. This constrains the system so that P3 can only move in the n2-axis like P1. To derive the new equations of motion, a new generalized velocity u4 is selected to represent the velocity of q4.
(20)
Fig. 6
Constrained schematic of the PPM
Fig. 6
Constrained schematic of the PPM
Close modal
The new constraint dictates that
(21)
Given the constraint applies to the velocity of point 3, the following independent constraint equations can be formed:
(22)
(23)
Given below are the new generalized inertia and active forces caused by the new constraint
(24)
(25)
where m4 is the mass of P4. Note that the new constraint has given rise to two dependent generalized speeds. Equations (22) and (23) can be rewritten in terms of the independent generalized speeds. This will give form to
(26)
where n = 4 is the number of independent generalized speeds (prior to constraints), and m = 2 is the number of independent constraints applied. From there, the new equation of motion for the constrained system can be written as
(27)
where r = 1,…, nm.

2.3 Performance Indicators.

The performance of the planar parallel manipulator is assessed using two key indicators in this study: the percentage of not usable vertical workspace (PNUVW) and the global conditioning index (GCI). PNUVW measures the efficiency of the mechanism in utilizing the available space, with a lower value indicating better performance. On the other hand, GCI serves as an indicator of dynamic performance, where a higher value is indicative of superior results.

2.3.1 Percentage of Not Usable Vertical Workspace.

The PNUVW serves as a performance indicator, gauging the efficiency of the manipulator in utilizing the available space. A value closer to one is undesirable, as it suggests that the overall machine size is disproportionately larger than the usable space. Ideally, a more efficient use of space is sought, emphasizing the importance of minimizing PNUVW for optimal manipulator performance. To calculate PNUVW, the total usable work area in the y-axis (yw) must be determined with the following equation:
(28)
where h is the total height of the PPM, and d is the distance from the carriage tower to the workspace where the end-effector cannot reach. The value of d can be written as
(29)
where xw is the workspace of the PPM in the x-direction. The PNUVW can be represented as the following ratio:
(30)

2.3.2 Global Conditioning Index.

The global conditioning index is based on the analysis of the CI over an entire workspace. The CI is a well-established all-around indexing method used to evaluate parallel manipulators [2426]. A higher CI is more desirable as that indicates the manipulator has a high degree of accuracy, stiffness, and dexterity. The CI can be determined based on the reciprocal of the conditioning number.
(31)
where λ1 and λ2 are the minimum and maximum singular values of the Jacobian matrix.
The Jacobian matrix was derived by Xu et al. [16].
(32)
The GCI will be
(33)
where W is the workspace of the end-effector.

3 Design Analysis/Optimization

This section will focus on the analysis and optimization of the PPM. For the analysis portion, a finite element method is utilized to understand the deflection and stress that the arms will be experiencing under static loading conditions assuming larger than expected load values to reflect the maximum potential loads during motion. For the optimization portion, a metaheuristic approach will take place to determine a range of desirable design variables that result in optimized performance indicators.

3.1 Finite Element Analysis Complexity Reduction.

The PPM arms are attached to the prismatic joints via carriage and revolute joints. When the PPM is at an angle, the revolute joints will experience an axial load, and some amount of deformation is expected to occur in the arms. Finite element analysis (FEA) was performed under static loading conditions of the PPM to determine the bearing loading conditions and the deflection in the arms. The static load values envelop possible max dynamical loading conditions by overestimating the static loading variables.

To perform a quick examination of the mechanism, the system can be reduced to the representative version as shown in Fig. 7. Figure 7 shows that the system can be modeled as a two-element truss problem experiencing an out-of-plane force applied at the node attaching both elements together.

Fig. 7
Simplified PPM for FEA
Fig. 7
Simplified PPM for FEA
Close modal
Each element has six degrees-of-freedom at each end for a total of 12 for the given conditions. The displacement of the element can be modeled using the following equation:
(34)
where K is the stiffness matrix, U is the displacement vector, and F is the force vector.
Given the total degrees-of-freedom for the elements is 12, it can be expected that the element K matrix be a 12 × 12 matrix. The global K matrix is expected to be 18 × 18 matrix because of the arrangement of all the elements. This means that both U and F will be a vector with a size of 18. For this problem, there are two fixed ends where the elements are restricted to experience no deformation at those nodes. Because of these boundary conditions, the global matrix can be reduced to a 6 × 6 matrix, and the force vector can also be reduced to a size of 6. Now the deformation vector can be calculated by the following equation:
(35)

The displacement vector solved for in Eq. (35) can be back substituted into Eq. (34) after reassembling the full-sized vector and be used to solve for the reaction forces.

3.2 Finite Element Analysis Via ansys.

The above method presents a simplified scenario that can be analyzed without any supporting software. However, given the reduction in complexity, the accuracy of the result is also slightly diminished. ansys is used to generate more accurate results because it has the capability to create a mesh of an object and break the object into thousands of different elements/nodes. ansys is then able to numerically solve the system of equations under given boundary conditions. For an accurate representation of the problem, the CAD model of the PPM was imported into ansys and the material properties for the arms and end-effector were specified as carbon fiber and polyethylene terephthalate respectively. All six arms were selected to be a fixed support on the face they would be attached to the carriage. Finally, a gravitational force was implemented at a 45-deg rotation, and a point mass representing that hot end was defined. The point load was estimated at double the actual mass to account for any additional loads caused by the acceleration of the machine. Figures 8 and 9 show the results.

Fig. 8
FEA reaction results
Fig. 8
FEA reaction results
Close modal
Fig. 9
FEA deflection results
Fig. 9
FEA deflection results
Close modal

From the simulation, it can be shown that the reaction at the fixed support is about 4 N from Fig. 8, and the maximum deflection is about 0.6 mm from Fig. 9. The 0.6 mm deflection can be considered negligible considering the change in length would result in a 0.2% difference. This indicates that the stiffness of the mechanism would impart a negligible amount of slack to the system resulting in no noticeable impact on 3D printing accuracy.

From the simulation, it shows that the bearings are expected to experience about 4 N of combined loading. The concern is that most of the loading experienced is axial loading and not normal loading. The bearings that are used have an inner diameter of 3 mm, an outer diameter of 10 mm, and a thickness of 4 mm. The bearing's spec sheet specifies a maximum allowable dynamic load of 628 N and a maximum static load of 216 N [27]. Importantly, the maximum dynamic and static load conditions exceed the expected static loading conditions by 157 and 54 times, respectively. Hence, it is reasonable to assume that the selected ball bearing will be suitable for this application, even when subjected to dynamic conditions typical in 3D printing and experiencing combined axial loading.

3.3 Metaheuristics Optimization.

Using the performance indicators of GCI and PNUVW, an optimum configuration of selected design variables can be determined using metaheuristics. To determine the optimum values for the length of the arms (L) and the width of the frame (R), the multi-objective genetic algorithm is employed and implemented with matlab. For this implementation, the cost functions used were the performance indicators. The bounds for the design variables are
(36)
The nonlinear constraints used are
(37)
(38)
where α is the angle in degrees made between the end-effector and the arm when the end-effector is at the edge of the workspace. The value of α can be written as
(39)

Equation (37) is used to ensure that the length of the arms can span the entire width of the PPM. Equation (38) is used to ensure that the PPM does not approach a singularity configuration when at the edge of the workspace. The specified parameters passed into the algorithm were r = 25 mm, xw = 220, and h = 700.

3.3.1 Optimization Results.

The Pareto front in Fig. 10, and the plot compared α to the PNUVW in Fig. 11 were created using the scenario described above and the multi-objective genetic algorithm from matlab's optimization toolbox.

Fig. 10
Pareto front of GCI and PNUVW
Fig. 10
Pareto front of GCI and PNUVW
Close modal
Fig. 11
Comparison between α and PNUVW
Fig. 11
Comparison between α and PNUVW
Close modal

The Pareto front in Fig. 10 represents the most efficient solutions given the competing optimum objectives. Figure 10 should be used as a reference for determining what parameters to use as there is a tradeoff between performance and spatial efficiency. Figure 11 was included to visualize the optimum results in terms of α. This figure is important for understanding how close the PPM is to achieving singularity when at the boundary of the workspace.

With the results of the optimization presented, there is no single solution for what values result in the best possible PPM. Instead, design desires need to be taken into consideration before an optimum design can be selected. For creating an infinite-angled 3D printer, the entire footprint of the machine is an important factor to consider. Machines that are very big are cumbersome and can get in the way. Therefore, a recommendation is made to select a solution that prioritizes the efficient use of space. This would lead to the selection of the first point on the Pareto front. Comparing this selection to Fig. 11, the parameters would result in α=6.2deg. This value is comparatively small, so the next point on the Pareto front is recommended. This selection would result in α=11.54deg, a more desirable result. The parameters resulting from the optimization would be L = 320 mm and R = 221.17 mm. This configuration would result in the best possible performance of the machine.

4 Prototype and Testing

This section details the construction of the first prototype along with the methodology used during the testing of the machine. The prototype discussed in this section had been built prior to the exploration of optimized parameters. The design variables used for the construction of the prototype can be found in Table 4.

Table 4

Prototype design variables

R =190 mm
L =380 mm
r =25 mm
h =700 mm
xw =220 mm
yw =235 mm
R =190 mm
L =380 mm
r =25 mm
h =700 mm
xw =220 mm
yw =235 mm

4.1 Prototype Construction.

The preliminary design for the angled PPM was modeled using Inventor. The model includes the ability to change the configuration of the manipulator to see how the PPM can move in the given space. Figures 12 and 13 show an example of the CAD model displaying different end-effector configurations. This was helpful in determining component size and potential interference that may occur prior to building the prototype.

Fig. 12
CAD model example configuration right side
Fig. 12
CAD model example configuration right side
Close modal
Fig. 13
CAD model example configuration left side
Fig. 13
CAD model example configuration left side
Close modal

The frame was constructed out of 20 × 20 mm aluminum extrusion and was attached together using a combination of 3D printed parts and steel brackets. The frame for the PPM mechanism was attached at a 45 deg angle to the base frame of the 3D printer using an angled bracket and additional side plates to increase frame rigidity.

The PPM consisted of linear rails, bearings, carbon fiber rods, stepper motors, a pulley system, and 3D printed parts. The arms of the mechanism featured carbon fiber rods with a 3D printed attachment where a bearing was press fitted into. Additionally, the prototype featured a tensioning system between two of the arms on either side of the end-effector. This system was included to reduce any sag in the end-effector. The controller for the system was a repurposed Anet A8 control board. The assembled prototype is shown in Fig. 14.

Fig. 14
Assembled prototype
Fig. 14
Assembled prototype
Close modal

4.2 Prototype Firmware/Software.

The firmware used to operate the machine is based on a fork of Marlin 1.1.8, an open-source firmware widely adaptable and utilized by various 3D printers. The firmware employed on the angled PPM printer is a derivative of Stuart's firmware [28], with modifications made to accommodate the angled attachment of the PPM and the extensive capabilities of the rotating bed.

To operate the machine, G-code is necessary to communicate to the firmware and instruct the end-effector's movements. The BlackBelt Cura (BBC) program was used to generate G-code used for the test. BBC is a free program provided by the BlackBelt company that is used to convert 3D models into the code used by angled 3D printers.

4.3 Testing Methodology.

With the completion of the prototype, testing was done to examine the dynamic performance of the machine. The following test was conducted to closely examine the acceleration of the carriages and the end-effector as the machine is programmed to run a series of ellipses of decreasing size. With the machine operating at high-speeds, the acceleration values collected will provide insight into the stability and rigidity of the system. If the high-speed testing results in irregular and large peaks in acceleration in unexpected directions, this will indicate poor performance for a 3D printer. Large spikes in acceleration would induce vibrations that could cause visible defects on the printed part. Ideal performance would exhibit acceleration only in the axis of motion, with clearly defined acceleration peaks and little noise. This would represent consistent circular motion where the end-effector is constantly experiencing a change of acceleration.

To test the dynamic performance of the angled PPM, an accelerometer was attached to three different test points: (1) the end-effector, (2) the right carriage, and (3) the left carriage. The accelerometer used was an ADXL335, and the data were collected using an Arduino Uno. According to the datasheet, the sensor is capable of measuring acceleration with a minimum full-scale range of ±3 g [29]. The accelerometer was calibrated using the table the printer was sitting on as the reference point. To calibrate a specific direction, the accelerometer was laid flat on the table, aligning the positive direction of the calibrating axis downwards. The readout of the accelerometer was then tuned to output 1 g. This process was repeated for each axis.

The accelerometer was attached flush to the surface of the test point so that the positive x-direction was aligned along the angled plane to the ground, the positive y-direction oriented to the right of the frame, and the positive z-direction was perpendicular to the angled frame downward. To initiate movement for testing, G-code was generated to move the end-effector in a series of decreasingly sized ellipses. Additionally, the motor controlling the rotating bed was unplugged to eliminate the possibility of induced vibrations from an external source. Due to resource limitations, acceleration values could only be collected from one test point at a time. A total of 15 tests were conducted, with each test point undergoing five tests. The results of the testing were then averaged for each position, and the final plot of the accelerations can be found in Sec. 5.

5 Results and Discussion

Prototype testing results are shown in Figs. 1517. These findings align with the expected outcomes. Figure 15 illustrates consistent levels of acceleration for the z-axis throughout the testing duration, due to its orthogonality to the plane of motion. Moreover, the findings depicted in Fig. 15 indicate that the end-effector exhibits stiffness, displaying minimal oscillation during rapid movements. This characteristic is highly desirable for 3D printing processes, as it contributes to enhanced precision and accuracy in the final printed results.

Fig. 15
End-effector acceleration over time
Fig. 15
End-effector acceleration over time
Close modal
Fig. 16
Left carriage acceleration over time
Fig. 16
Left carriage acceleration over time
Close modal
Fig. 17
Right carriage acceleration over time
Fig. 17
Right carriage acceleration over time
Close modal

The outcomes depicted in Figs. 16 and 17 align with expectations but unveil an intriguing phenomenon. In both the left and right carriages, the acceleration in the Y and Z directions should ideally manifest as constant flat lines, given that the carriage can only travel along the frame. However, Fig. 16 reveals occasional peaks in Y-axis accelerations, introducing irregularities in comparison to the smoother profile observed in Fig. 17. This suggests that the carriage deviates from its linear path along the frame toward the outside or center of the frame.

This deviation is likely attributed to the linear rails incorporating roller bearings. As the end-effector moves laterally and undergoes directional changes, inertia exerts a force on the PPM, compelling it to move away from the anchor point. The minimal play present between the bearings allows the carriage to experience momentary changes in accelerations. To mitigate this phenomenon, an alternative attachment for the prismatic joint is proposed. Utilizing V-slot roller wheels could offer a more effective solution by functioning as grippers, enhancing adherence to the frame, and minimizing undesired lateral movements.

Throughout the testing phase, the stepper motors emitted a discernible whine, accompanied by an audible “stepping” sound that persisted the entire time. The absence of an optimized firmware package to efficiently drive the machine might have influenced the performance of the PPM. The evident lack of smooth movement is present at all testing points, manifested as noisy and jagged accelerations. This erratic movement implies a stop-and-go pattern as the printer traveled. Addressing these issues could be achieved through the implementation of updated software.

6 Concluding Remarks and Future Work

This paper introduces a novel belt 3D printer design featuring an angled PPM instead of the traditional gantry design. The new design significantly reduces the moving weight, reducing inertia, and facilitating higher accelerations and dynamic performance. However, it is essential to acknowledge that this design alteration introduces additional unused space, resulting in a larger footprint of the 3D printer. To address the tradeoff between spatial efficiency and dynamic performance, an optimization using the multi-objective genetic algorithm was employed. The outcomes of this optimization effort yielded a Pareto front, delineating optimal design parameters based on key performance indicators—specifically, the GCI for dynamic performance and the PNUVW for spatial utilization. A static FEA evaluation conducted under conservative loading assumptions reveals the impact of the design on bearing load and arm deflection within the PPM. Importantly, findings suggest that these factors are unlikely to compromise the performance of the angled PPM. Finally, prototype testing substantiates the stability and rigidity of the angled PPM during motion, providing empirical support for the proposition that this innovative approach holds promise for advancing belt 3D printing technology. This comprehensive investigation highlights the potential viability of the angled PPM as a progressive and suitable solution in the realm of belt 3D printing.

The work done in this paper highlights the novelty in kinematic and dynamic formulation and optimization of planar parallel manipulator needed for the design of a 3D printer. While these findings mark a significant contribution, the journey toward a competitive alternative to the prevalent Cartesian gantry-designed infinite 3D printers remain ongoing. Future endeavors will focus on fine-tuning the control system, enhancing overall construction, and implementing optimal design parameters. Through the refinement of the machine, direct comparisons with various belt printers can be conducted, allowing for a comprehensive evaluation of the effectiveness of the angled PPM design. In essence, this work lays the groundwork for continued exploration into advanced 3D printing technologies.

Acknowledgment

Saint Martin's University and the Hal and Inge Marcus School of Engineering provided access to the Robotics and Manufacturing labs. This access was instrumental in enabling the fabrication and testing of a live prototype.

Conflict of Interest

There are no conflicts of interest.

Data Availability Statement

The datasets generated and supporting the findings of this article are obtainable from the corresponding author upon reasonable request.

References

1.
ASTM International
,
2012
,
Standard Terminology for Additive Manufacturing Technologies: Designation F2792-12a
,
ASTM International
,
West Conshohocken, PA
.
2.
Shahrubudin
,
N.
,
Lee
,
T. C.
, and
Ramlan
,
R.
,
2019
, “
An Overview on 3D Printing Technology: Technological, Materials, and Applications
,”
Procedia Manuf.
,
35
, pp.
1286
1296
.
3.
Yan
,
Q.
,
Dong
,
H.
,
Su
,
J.
,
Han
,
J.
,
Song
,
B.
,
Wei
,
Q.
, and
Shi
,
Y.
,
2018
, “
A Review of 3D Printing Technology for Medical Applications
,”
Engineering (Beijing, China)
,
4
(
5
), pp.
729
742
.
4.
Yap
,
Y. L.
,
Sing
,
S. L.
, and
Yeong
,
W. Y.
,
2020
, “
A Review of 3D Printing Processes and Materials for Soft Robotics
,”
Rapid Prototyp. J.
,
26
(
8
), pp.
1345
1361
.
5.
Nichols
,
M. R.
,
2019
, “
How Does the Automotive Industry Benefit From 3D Metal Printing?
,”
Met. Powder Rep.
,
74
(
5
), pp.
257
258
.
6.
Contour Crafting Corporation
,
2017
, “
Home
https://www.contourcrafting.com/.
7.
Formlabs
,
2019
, “
How 3D Printing is Disrupting the Jewelry,
” Industry. Formlabs, https://formlabs.com/blog/3d-printed-jewelry/.
8.
Relativity
,
2020
, “
Building Humanity's Multiplanetary Future
https://www.relativityspace.com/.
9.
Blackbelt 3D
,
2021
, “
The Printer,
https://blackbelt-3d.com/.
10.
Creality
,
2022
, “
Creality 3DPrintMill(CR-30) Belt 3D Printing for Everyone,
https://www.kickstarter.com/projects/3dprintmill/creality-infinite-z-axis-build-volume-printer-3dprintmill.
11.
NAK3DDesigns
,
2021
, “
Nak3ddesigns/White-Knight: White Knight Belt Printer,
https://github.com/NAK3DDesigns/White-Knight.
12.
Swaleh Owais
,
2018
, “
Automatic Infinite 3D Printer,
https://hackaday.io/project/114738-automatic-infinite-3d-printer.
13.
AdamfromPowerbelt3D
,
2022
, “
Ender EZ Belt 3D Printer Conversion,
https://www.printables.com/model/255684-ender-ez-belt-3d-printer-conversion.
14.
Carabin
,
G.
,
Scalera
,
L.
,
Wongratanaphisan
,
T.
, and
Vidoni
,
R.
,
2021
, “
An Energy-Efficient Approach for 3D Printing With a Linear Delta Robot Equipped With Optimal Springs
,”
Rob. Comput. Integr. Manuf.
,
67
, p.
102045
.
15.
Rodriguez
,
E.
,
Alvares
,
A. J.
, and
Jaimes
,
C. I. R.
,
2019
, “
Conceptual Design and Dimensional Optimization of the Linear Delta Robot With Single Legs for Additive Manufacturing
,”
Proc. Inst. Mech. Eng. Part I J. Syst. Control Eng.
,
233
(
7
), pp.
855
869
.
16.
Xu
,
J.
,
Feng
,
X.
,
Cen
,
J.
, and
Zhang
,
S.
,
2020
, “
Precision Forward Design for 3D Printing Using Kinematic Sensitivity Via Jacobian Matrix Considering Uncertainty
,”
Int. J. Adv. Manuf. Technol.
,
110
(
11–12
), pp.
3257
3271
.
17.
Rehman
,
S. U.
,
Raza
,
M.
, and
Khan
,
A.
,
2019
, “
Delta 3D Printer: Metal Printing
,”
J. Electr. Eng. Electron. Control Comput. Sci.
,
5
(
17
), pp.
19
24
. https://jeeeccs.net/index.php/journal/article/viewFile/130/102
18.
Borrelli
,
V.
,
2018
, “
Kinematic and Dynamic Analysis of a Machine for Additive Manufacturing
,”
Master’s thesis
,
Politecnico Di Torino
.
19.
Liu
,
X.-J.
,
Wang
,
Q.-M.
, and
Wang
,
J.
,
2005
, “
Kinematics, Dynamics, and Dimensional Synthesis of a Novel 2-DOF Translational Manipulator
,”
J. Intell. Rob. Syst.
,
41
(
4
), pp.
205
224
.
20.
Hao
,
Q.
,
Guan
,
L.
, and
Wang
,
J.
,
2012
, “
Dynamic Performance Evaluation of a 2-DOF Planar Parallel Mechanism
,”
Int. J. Adv. Rob. Syst.
,
9
(
6
), p.
250
.
21.
Wu
,
J.
,
Wang
,
D.
, and
Wang
,
L.
,
2015
, “
A Control Strategy of a Two Degrees-of-Freedom Heavy Duty Parallel Manipulator
,”
ASME J. Dyn. Sys. Meas. Control
,
137
(
6
), p.
061007
.
22.
Wu
,
J.
,
Yu
,
G.
,
Gao
,
Y.
, and
Wang
,
L.
,
2018
, “
Mechatronics Modeling and Vibration Analysis of a 2-DOF Parallel Manipulator in a 5-DOF Hybrid Machine Tool
,”
Mech. Mach. Theory
,
121
, pp.
430
445
.
23.
Kurfess
,
T. R.
,
2015
,
Robotics and Automation Handbook
,
CRC Press
,
Boca Raton, FL
.
24.
Angeles
,
J.
, and
López-Cajún
,
C.
,
1992
, “
Kinematic Isotropy and the Conditioning Index of Serial Robotic Manipulators
,”
Int. J. Rob. Res.
,
11
(
6
), pp.
560
571
.
25.
Liu
,
X.-J.
,
Jin
,
Z.-L.
, and
Gao
,
F.
,
2000
, “
Optimum Design of 3-DOF Spherical Parallel Manipulators With Respect to the Conditioning and Stiffness Indices
,”
Mech. Mach. Theory
,
35
(
9
), pp.
1257
1267
.
26.
Salisbury
,
J. K.
, and
Graig
,
J. J.
,
1982
, “
Articulated Hands: Force Control and Kinematic Issues
,”
Int. J. Rob. Res.
,
1
(
1
), pp.
4
17
.
27.
Bearings Direct
, “
Miniature Radial Ball Bearing,
https://bearingsdirect.com/600mini.spec.htm
28.
Stuart, Rob, Allison
,
2019
, “
Deltesian Printer,
” GitHub Repository, https://github.com/bornity/deltesian-part-files.
29.