Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Tool Path Generation For CNC Milling Based On STL File: Nadia Sami Hasan Dr. Laith Abdullah Mohammed

Download as pdf or txt
Download as pdf or txt
You are on page 1of 12

Eng. &Tech.Journal, Vol.33,Part (A), No.

2, 2015

Tool Path Generation for CNC Milling Based on STL File

Nadia Sami Hasan


Production & Materials Engineering Department, University of Technology/Baghdad
Dr. Laith Abdullah Mohammed
Production & Materials Engineering Department, University of Technology/Baghdad
Email: dr.laith@uotechnology.edu.iq

Received on:6/4/2014 & Accepted on:8/1/2015

ABSTRACT
This paper proposes and develops algorithms to read (STL) files and extract
engineering entities required in CNC milling processes. The proposed algorithms are
dependent on some mathematical modeling and manipulations of the engineering models
by slicing an (STL) file to many slices and then building the required algorithms to adopt
these slices to generate machining paths for CNC milling machines as (G-Codes).
The proposed system is divided into three parts, in the first part, an algorithm proposed to
extract engineering object entities to some proposed models based on their (STL) files
using Matlab program. The proposed models include cube, cylinder, dome, cone and
cavity models. In the second part, a slicing algorithm is proposed to enable the slices
along the proposed models z-axis to find and navigate the required manufacturing data.
UGS program was used also to generate the tool paths and to simulate the machining
process and then generate NC part program of the proposed objects (G-Code). Finally,
the third part of this work includes comparing results produced based on both (STL) and
(UGS), to achieve the required aim of this paper by experimentally comparing objects
surface roughness, resulted values show that the maximum difference in average R a is
equal to (0.27µm).

Keywords: STL files, CNC Milling.

STL ‫ﺗﻮﻟﯿﺪ ﻣﺴﺎر اﻟﻌﺪة ﻟﻤﻜﺎﺋﻦ اﻟﺘﻔﺮﯾﺰ اﻟﻤﺒﺮﻣﺠﺔ ﺑﺄﻋﺘﻤﺎد ﻣﻠﻒ‬


:‫اﻟﺨﻼﺻﺔ‬
‫( واﺳﺘﻨﺒﺎط اﻟﺴﻤﺎت اﻟﮭﻨﺪﺳﯿﺔ اﻟﻼزﻣﺔ ﻟﻌﻤﻠﯿﺎت‬STL) ‫ﻓﻲ ھﺬا اﻟﺒﺤﺚ ﺗﻢ اﻓﺘﺮاض وﺗﻄﻮﯾﺮ ﺧﻮارزﻣﯿﺎت ﻟﻘﺮاءة ﻣﻠﻔﺎت‬
‫ اﻋﺘﻤﺪت اﻟﺨﻮارزﻣﯿﺎت اﻟﻤﻄﻮرة ﻋﻠﻰ اﻟﻨﻤﺬﺟﺔ اﻟﺮﯾﺎﺿﯿﺔ واﻟﻤﻌﺎﻟﺠﺔ ﻟﻨﻤﺎذج ھﻨﺪﺳﯿﺔ ﻣﻦ ﺧﻼل ﺗﻘﺴﯿﻢ‬.‫اﻟﺘﻔﺮﯾﺰ اﻟﻤﺒﺮﻣﺠﺔ‬
‫( اﻟﻰ ﻋﺪة ﺷﺮاﺋﺢ وﺑﻌﺪھﺎ ﺑﻨﺎء ﺧﻮارزﻣﯿﺎت ﻻﻋﺘﻤﺎد ھﺬه اﻟﺸﺮاﺋﺢ ﻓﻲ ﺗﻮﻟﯿﺪ ﻣﺴﺎرات اﻟﻘﻄﻊ ﻟﻤﻜﺎﺋﻦ اﻟﺘﻔﺮﯾﺰ‬STL) ‫ﻣﻠﻒ‬
.(G-codes) ‫اﻟﻤﺒﺮﻣﺠﺔ ﻛﺒﺮاﻣﺞ ﺗﺸﻐﯿﻞ‬
‫ ﻓﻲ اﻟﺠﺰء اﻷول ﺗﻢ اﻗﺘﺮاح ﺧﻮارزﻣﯿﺔ ﻻﺳﺘﻨﺒﺎط اﻟﺴﻤﺎت اﻟﮭﻨﺪﺳﯿﺔ‬.‫ﺗﻢ ﺗﻘﺴﯿﻢ اﻟﻨﻈﺎم اﻟﻤﻘﺘﺮح اﻟﻰ ﺛﻼﺛﺔ اﺟﺰاء‬
‫ ﺗﺸﻤﻞ اﻟﻨﻤﺎذج اﻟﻤﻘﺘﺮﺣﺔ‬.Matlab ‫( ﺑﺮﻧﺎﻣﺞ ﻣﻌﺪ ﻋﻠﻰ ﻧﻈﺎم‬STL) ‫ﻟﻤﺠﻤﻮﻋﺔ ﻣﻦ اﻟﻨﻤﺎذج اﻟﻤﻘﺘﺮﺣﺔ اﻋﺘﻤﺎدا ﻋﻠﻰ ﻣﻠﻔﺎﺗﮭﺎ‬
‫ ﻓﻲ اﻟﺠﺰء اﻟﺜﺎﻧﻲ ﺗﻢ اﻗﺘﺮاح ﺧﻮارزﻣﯿﺔ ﻟﻠﺘﻘﻄﯿﻊ ﻟﻠﺤﺼﻮل ﻋﻠﻰ‬.‫ اﻟﻤﺨﺮوط واﻟﺘﺠﻮﯾﻒ‬، ‫ اﻟﻘﺒﺔ‬، ‫ اﻷﺳﻄﻮاﻧﺔ‬، ‫اﻟﻤﻜﻌﺐ‬
(UGS) ‫ ﺗﻢ اﯾﻀﺎ اﺳﺘﺨﺪام ﺑﺮﻧﺎﻣﺞ‬.‫ ﻟﻠﻨﻤﺎذج اﻟﻤﻘﺘﺮﺣﺔ ﻟﻠﺤﺼﻮل ﻋﻠﻰ اﻟﺒﯿﺎﻧﺎت اﻟﺘﺼﻨﯿﻌﯿﺔ‬z ‫ﺷﺮاﺋﺢ ﻋﻠﻰ طﻮل ﻣﺤﻮر‬

350
Eng. &Tech.Journal, Vol.33,Part (A), No.2, 2015 Tool Path Generation for CNC Milling Based
on ST File

‫ اﻟﺠﺰء اﻟﺜﺎﻟﺚ ﺗﻀﻤﻦ ﻣﻘﺎرﻧﺔ‬.‫ﻟﺘﻮﻟﯿﺪ ﻣﺴﺎرات اﻟﻘﻄﻊ وﻣﺤﺎﻛﺎة ﻋﻤﻠﯿﺔ اﻟﺘﺸﻐﯿﻞ وﻣﻦ ﺛﻢ ﺗﻮﻟﯿﺪ ﺑﺮاﻣﺞ اﻟﺘﺸﻐﯿﻞ ﻟﮭﺬه اﻟﻨﻤﺎذج‬
‫( ﻣﻦ ﺧﻼل اﻟﻤﻘﺎرﻧﺔ اﻟﻌﻤﻠﯿﺔ ﻟﻘﯿﻢ‬UGS) ‫( واﻋﺘﻤﺎد ﻣﻠﻔﺎت ﺑﺮﻧﺎﻣﺞ‬STL) ‫اﻟﻨﺘﺎﺋﺞ اﻟﻤﺘﻮﻟﺪة ﻣﻦ ﻛﻞ ﻣﻦ اﻋﺘﻤﺎد ﻣﻠﻔﺎت‬
.(0.27µm) ‫ ﺣﯿﺚ ﺑﯿﻨﺖ اﻟﻨﺘﺎﺋﺞ أن أﻛﺒﺮ ﻓﺮق ﻓﻲ ﻣﻌﺪﻻت ﻗﯿﻢ اﻟﺨﺸﻮﻧﺔ اﻟﺴﻄﺤﯿﺔ ﯾﺴﺎوي‬،‫اﻟﺨﺸﻮﻧﺔ اﻟﺴﻄﺤﯿﺔ‬

INTRODUCTION

T oday manufacturing relies heavily on CNC machines. The operations of CNC


machines require modeling a part in the Computer Aided Design (CAD) system,
many CAD systems are being used and each of them uses a different kernel to
describe the geometry. It is not possible to read all of these different data formats using
one software package. A standard format that can be output by most commercial CAD
systems should be used as the input for Computer Aided Manufacturing (CAM) system
[1, 2].
An STL file is a type of standardized computer exchange file which contains a 3D
model. The representation of the surface(s) of the object(s) in the file is in the form of one
or more polygon meshes. The polygon meshes in an STL file are entirely composed of
triangular faces, edges and vertices. Further, the faces have assigned normal’s which
indicates their orientation (inside/outside).
The name STL is taken from its extension, .stl, originally because the files were
intended for the rapid prototyping process called Stereo lithography. The file format has
become a world standard for exchanging 3D polygon mesh type objects between
programs, and .stl’s are now used as input for virtually all rapid prototyping processes, as
well as some 3D machining [3].
There are many articles in literature that deal with different CAD formats used in CNC
milling operations. In this section the most relative articles related to the present work are
reviewed. The research of R. Lin and C. Ye (2012) [4] proposed a method to produce
accurate tool paths using STL file formats and machined using 5-axis machining. The
result shows that the trajectory contour error is improved about 5 times by using
conventional method and the proposed approach, respectively. The new method improves
not only the max. contour error, but part surface smoothness. A new tool-path generation
method is presented by M. Li, et al (2012) [5] by adopting the thought of generating
cutter-location (CL) data directly from corresponding cutter-contact (CC) data. A
computational model is proposed to calculate single-layer interference-free CL contour,
instead of computing a precise CL point by checking potential interferences from
candidate facet, vertices, and edges, respectively. So this method is more efficient and
easier for program implementation. Implementation tests prove that the new method is
effective and robust for STL models, and tool-path achieved is highly precise. G.
Kiswanto and M. Azka (2012) [6] research focuses on detection of primitive features,
based on faceted models, available in a part, this is done by applying part slicing and
grouping adjacent facets. The results showed that this method can identify the primitive
features automatically, accurately in all planes of tested part, this covered pocket,
cylindrical and profile features. To compute the vertex normal of triangular meshes more
accurately, the paper of Y. Chen, et al (2013) [7] presents an improved algorithm based
on angle and centroid weights. Firstly, algorithms are analyzed by comparing their
weighting characteristics such as angles, areas and centroids. Following that, an improved
algorithm is put forward based on angle and centroid weights. Finally, by taking the
deviation angle between the nominal normal vector and the estimated one as the error

351
Eng. &Tech.Journal, Vol.33,Part (A), No.2, 2015 Tool Path Generation for CNC Milling Based
on ST File

evaluation standard factor. Experimental results demonstrate that the algorithm is


effective.

The Proposed System


Figure (1) describes the main steps adopted in the algorithms which are achieved in
this work. The basic modeling primitives that are used in this study are parameterized as
follows:
1- A cube which is parameterized by its width, height, and length.
2- A cylinder which is parameterized by its radius and height.
3- A dome which is parameterized by its radius.
4- A cone which is parameterized by its radius and height.
5- Finally, A cavity which is parameterized by its radius.
AutoCAD program facilities are used to draw the object wireframe model as shown in
Figures (2 a) which represents the wireframe models for the cube starting from the origin
with 60 mm length, Figure (2 b) which represents the wireframe model for the cylinder
with a base circle center starting from the origin with radius equal to 25 mm and cylinder
height is equal to 60 mm, and Figure (2 c) which represents the wireframe model for the
dome starting from the origin with radius equal to 50 mm, and Figure (2 d) which
represents the wireframe model for the cone starting from the origin with radius equal to
25 mm and its height is equal to 60 mm, and finally Figure (2 e) which represents the
wireframe model for the cavity starting from the origin with radius equal to 50 mm.
These wireframe models were saved as an STL file extension format using “STLOUT”
AutoCAD command. Consequently, this file format will be modeled for manipulation
using MATLAB Package environment using an algorithm depending on the following
procedures:

Object Mathematical Representation


In this stage of modeling, the program is built using MATLAB R2010a as a first step
to generate the database of geometry as a set of points in three axes in Cartesian
coordinate system. The set of points that represent the shape of the objects shapes models
obtained from previous stage will be exported to the MATLAB environment to transfer
the objects as a set of points to any other CAD/CAM software without any distortion.
Figure (3) shows the models triangulation meshs for the proposed cube, cylinder, dome,
cone and cavity shape respectively.
Based on this figure, one can recognize 12 triangular meshes with 36 vertices which
represent the cube model, 144 triangular meshes with 432 vertices which represent the
cylinder model, 900 triangular meshes with 2700 vertices which represent the dome
model, 72 triangular meshes with 216 vertices which represent the cone model, and 900
triangular meshes with 2700 vertices which represent the cavity model.
The next step including the transforming process to the exported models to a set of
surfaces as shown in Figure (4) which represent the cube, the cylinder, the dome, the
cone, and the cavity triangular surface models respectively, according to a procedure
achieved using a developed algorithm in MATLAB environment.

352
Eng. &Tech.Journal, Vol.33,Part (A), No.2, 2015 Tool Path Generation for CNC Milling Based
on ST File

Triangular-Mesh Intersection Detection


The triangular-mesh intersection detection approach is commonly used in current 3D
graphics hardware and software, as it provides a precise intersection detection
mechanism. Essentially, the triangular-mesh technique involves the detection of an
intersection between a line and a triangle. However, this process needs to be applied for
every triangle that comprises the primitive.
The foundation of the triangular-mesh approach is the algorithm used to detect an
intersection between a line and a triangle in three dimensions. This process has been
studied in depth to develop fast and efficient techniques. Existing techniques can be
divided into two broad categories.
The first category of algorithms uses a two step process. These algorithms first locate
where the straight line intersects with the plane containing the triangle. Following this,
the algorithms detects if the intersection point lies within the triangle as shown in figure
(5). The equation for the plane containing this triangle is given by Equation (1).

Ax + By + Cz + D = 0
where: A = (y2 – y1) (z3 – z1) - (z2 – z1) (y3 – y1)
B = (z2 – z1) (x3 – x1) - (x2 – x1) (z3 – z1)
C = (x2 – x1) (y3 – y1) - (y2 – y1) (x3 – x1)
Hence;

D = -A × x1 - B ×y1 - C × z1 ….. (1)

Using Equation (1) for the plane in 3D space, it is possible to determine the
intersection point (x i, yi, zi) using Equation (2).
N = - (A ×xs + B×ys + C ×zs + D)
M = A (xe – xs) + B (ye – ys) + C (ze – zs)

N
S=
M

= S× − + ….. (2)

Using Equation (2) one can determine the point of intersection coordinates between the
line and the triangular plane, accordingly, Figure (6) represents the cube, the cylinder, the
dome, the cone, and the cavity models respectively including the normal vectors for each
triangular mesh.

The Slicing Process


The slicing planes, as shown in Figure (7), are planes parallel to the reference plane.
In this case, the XY plane is simply chosen as reference plane. Then the database is
developed for facets which are intersected by each slicing plane. Facets indexes are
stored in a vector data structure. The interval between slicing planes is determined based
on the specified accuracy and also the distance between Z min and Zmax value of the faceted
model. Once the information regarding the index of facet, index of vertex, and

353
Eng. &Tech.Journal, Vol.33,Part (A), No.2, 2015 Tool Path Generation for CNC Milling Based
on ST File

coordinates of vertex are identified and stored in the data array, then facets that are
intersected by the slicing plane are to be found. The search is done through the facet
index. The slicing algorithm can be explained as follows in the flowchart shown in Figure
(8).
In order to minimize searching and manipulating time and to avoid stair case, this
algorithm is performed by checking each triangular facet, which has a normal vector
parallel to the slicing plane, whether intersected by the slicing plane or not. If it is
intersected, the intersection point coordinates is stored in an array data. Intersection
checking with each facet was conducted for each slicing plane along the slicing axis.
Figure (9) represents the resulting slicing process achieved on faceted cube, cylinder,
dome, cone and cavity models respectively. Accordingly, these obtained figures represent
the array manufacturing data for the suggested models, which is then exported to the
UGS package as .dat file to achieve the machining process including tool path and g-code
generation. Figure (10) represent the machined objects using the proposed algorithm.

Measurements Of Surface Roughness


The average roughness (Ra) has been examined for the machined objects by using
surface roughness tester brand “Pocket Surf” USA made, this parameter (R a) is chosen
because it is simple and widely used as a parameter for the surface finishing. A full
experimental procedure is carried out for each machined object, using 1 mm long as a
measuring range that is based on the instrument probe movement, holding the object in
arrangement permits to achieve five lines measurement recorded along the machining
paths and the surface roughness has been measured from the object top towards its
bottom perpendicular to the direction of the tool movement. Table (1) lists the resulted
average surface roughness data. According to these data, it is obvious that the maximum
difference in average R a is equal to (0.27µm).

CONCLUSIONS
On the basis of this study, and observations recorded experimentally, the following
findings can be concluded:
1- STL file format consists of unordered list of triangular facets.
2- According to the conducted experiments and the plotted results, one can observe
that the objects produced using STL models are clearly have rough surfaces compared
with the surfaces produced by UGS with average percentage of 13µm.
3- The main cause of the high surface roughness values of the objects produced by
STL is the cutting tool length and this motive increases the vibration in the tool axis
leading to increase the magnitude of (R a).
4- The proposed system gives less number of machining blocks and minimum
machining time as compared with those resulted from the UGS system, this is because
the proposed algorithm is performed by checking each triangular facet, which has a
normal vector parallel to the slicing plane, whether intersected by the slicing plane or not.
If it is intersected, the intersection point coordinates are stored in a manufacturing array
data.
5- The points close to the object top are more prone to false feature extraction than
the points around the object surface area, this can be seen clearly in both dome and cone
shaped objects respectively.

354
Eng. &Tech.Journal, Vol.33,Part (A), No.2, 2015 Tool Path Generation for CNC Milling Based
on ST File

Table (1) the resulted average surface roughness data

Cube Cylinder Dome Cone Cavity


Propos Propose Propos UG Propose Propose
UGS UGS UGS UGS
ed d ed S d d
Average
1.7 1.73 1.6 1.74 1.41 1.56 1.67 1.94 1.46 1.52
Ra (µm)

Start

Stl File

Scan the Stl File


to Extract Facet Data

Scan Facet Data

Get Vertices and Normal Vectors Coordinates

Slicing Process

Slice 1
Slice Tracing
Slice 2
Finding Machining Points
Slice 3

Machining
Path file
Slice n

Figure (1) Flowchart describing the main steps of the algorithms achieved in this work

355
Eng. &Tech.Journal, Vol.33,Part (A), No.2, 2015 Tool Path Generation for CNC Milling Based
on ST File

(a) (b) (c)

(d) (e)

Figure (2) The Main Wireframe Models for the (a) Cube,
(b) Cylinder, (c) Dome, (d) Cone, and (e) Cavity

(a) (b) (c)

(d) (e)

Figure (3) The Models Triangular Mesh for the (a) Cube,
(b) Cylinder, (c) Dome, (d) Cone, and (e) Cavity

356
Eng. &Tech.Journal, Vol.33,Part (A), No.2, 2015 Tool Path Generation for CNC Milling Based
on ST File

(a) (b) (c)

(d) (e)
Figure (4) The Triangular Surface Models for the (a) Cube, (b) Cylinder,
(c) Dome, (d) Cone, and (e) Cavity

(x3, y3, z3)

End (E)
(xe, ye, ze) Intersection
(xi, yi, zi)

Start (S)
(xs, ys, zs)

(x1, y1, z1)


(x2, y2, z2)

Figure (5) Triangular Mesh Intersection in 3D Cartesian Space

357
Eng. &Tech.Journal, Vol.33,Part (A), No.2, 2015 Tool Path Generation for CNC Milling Based
on ST File

(a) (b) (c)

(d) (e)

Figure (6) The Normal Vectors for Each Triangular Mesh


for the (a) Cube, (b) Cylinder, (c) Dome, (d) Cone, and (e) Cavity

Faceted
Model

Slicing
Axis Slicing
Plane

Figure (7) The Slicing Process on the Faceted Models

358
Eng. &Tech.Journal, Vol.33,Part (A), No.2, 2015 Tool Path Generation for CNC Milling Based
on ST File

Start

Input STL Data Structure


(AutoCAD)

Define Slicing Interval


(Zslice)

For Every Facet

No
Zmin<Zslice<Zmax

Yes

Set of a Plane End


Points

Zslice = Zslice + Interval

Figure (8) The Developed Slicing Algorithm Flowchart

359
Eng. &Tech.Journal, Vol.33,Part (A), No.2, 2015 Tool Path Generation for CNC Milling Based
on ST File

Slicing Planes

(a) (b) (c)

(d) (e)

Figure (9) The Slicing Process Achieved on the 3D Models


for the (a) Cube, (b) Cylinder, (c) Dome, (d) Cone, and (e) Cavity

(a) (b) (c)

(d) (e)

Figure (10) The Machined Objects Using The Proposed Algorithm


for the (a) Cube, (b) Cylinder, (c) Dome, (d) Cone, and (e) Cavity

360
Eng. &Tech.Journal, Vol.33,Part (A), No.2, 2015 Tool Path Generation for CNC Milling Based
on ST File

REFERENCES:
[1] J. W. Gunnick, “Multi-Axis High Speed Milling: How to Speed up Prototyping &
Tooling Process by Using STL Technology”, Proceedings TCT Conference 1998,
Nottingham, pp.43-65,Oct. 1998.
[2] L. Lennings, “Selecting Either Layered Manufacturing or CNC Machining to Build
Your Prototype”, Rapid Prototyping and Manufacturing Conference,
Rosemont(Chicago), April 2000.
[3] M. Löwdin, “Time Reducing Actions With Freeform Fabrication and Rapid
Prototyping”, M.Sc. Thesis, Örebro University,Sweden, 2004.
[4] R. Lin and C. Ye, “Accurate Trajectory Control for Five-Axis Tool-Path Planning”,
Proceeding of the international Multi conference of engineers and computer scientists,
Vol. II, Hong Kong, March. 14-16, 2012. [www.ivsl.org].
[5] M. Li, L. Zhang, J. Mo, and Y. Lu, “Tool-path generation for sheet metal incremental
forming based on STL model with defects”, International Journal of Advanced
Manufacturing Technology, pp. 535–547, (2012). [www.ivsl.org].
[6] G. Kiswanto and M. Azka, “Automatic Part Primitive Feature Identification Based on
Faceted Models”, International Journal of Computer Science Issues, Vol. 9, Issue 5, No
2, September 2012. [www.ivsl.org].
[7] Y. Chen, J. Gao, H. Wen and X. Chen, “Estimation Normal Vector of Triangular
Mesh Vertex by Angle and Centroid Weights and its Application”, TELKOMNIKA, Vol.
11, No. 4, pp. 1841-1848, April 2013. [www.ivsl.org].

361

You might also like