Table of Contents Author Guidelines Submit a Manuscript
Advances in Materials Science and Engineering
Volume 2011 (2011), Article ID 137407, 5 pages
Research Article

Precise Hole Drilling in PMMA Using 1064 nm Diode Laser CNC Machine

1Institute of Laser for Postgraduste Studies, University of Baghdad, P.O. Box 47314, Jadriha, Baghdad, Iraq
2Laser and Optoelctronics Engineering Department, Al-Nahrain university, Jadriha, Baghdad, Iraq

Received 30 December 2010; Accepted 21 February 2011

Academic Editor: J. Dutta Majumdar

Copyright © 2011 Jinan A. Abdulnabi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


This paper represents the outcome of efforts that intended to achieve laser hole drilling execution in polymethylmethacrylate (PMMA) of 2.5 mm thickness using 1064 nm diode laser of 5 W output power. Different laser beam powers, exposure time, and positions of the laser spot were taken into consideration with respect to the workpiece. The workpieces were tested in the existence of low-pressure assist gas (20–60 mmHg of air). The experimental results were supported by the predicted results of the analytical model.

1. Introduction

The conduction of heat in a three-dimensional solid is given by the solution of the following equation:𝜌𝐶𝑃𝜕𝑇=𝜕𝜕𝑡𝐾𝜕𝑥𝜕𝑇+𝜕𝜕𝑥𝐾𝜕𝑦𝜕𝑇+𝜕𝜕𝑦𝐾𝜕𝑧𝜕𝑇𝜕𝑧+𝐴(𝑥,𝑦,𝑧,𝑡),(1) where, the thermal conductivity 𝐾 (Wcm−1 K−1), the density 𝜌 (g cm−3), and the specific heat 𝐶𝑃 (Jkg−1 K−1) are dependent on the temperature and position. The rate of the applied heat to the solid is 𝐴(𝑥,𝑦,𝑧,𝑡) per unit time per unit volume, and 𝑡 is the time [1].

Using the cylindrical coordinates 𝑟 and 𝑧 (Figure 1), the temperature distribution is [1, 2] =𝑇(𝑟,𝑧,𝑡)𝑃𝜀2𝜋a𝐾0𝐽0(𝑚𝑟)𝐽1(×𝑧𝑚𝑎)exp(𝑚𝑧)erfc2(𝑘𝑡)1/2𝑚(𝑘𝑡)1/2𝑧𝑒(𝑚𝑧)erfc2(𝑘𝑡)1/2+𝑚(𝑘𝑡)1/2𝑑𝑚𝑚+𝑇0,(2) where, 𝑟 is the radial coordinate (hole radius), 𝑧 is the axial coordinate (thermal penetration depth), 𝐽𝑜 and 𝐽1 are Bessel functions of the first kind, 𝑃 is the constant power during a laser pulse, 𝑎 is the radius of the laser spot at the surface, 𝐾 is the thermal conductivity of the material, 𝑘 is the thermal diffusivity, 𝜀 is the fraction of incident radiation absorbed, 𝑚 is an integer that represents the limit of integration, 𝑡 is the exposure time, and 𝑇𝑜 is the initial temperature.

Figure 1: Coordinate system [1].

The numerical solution of (2) for determining the temperature distribution as a function of time at any point inside the material along the 𝑍-axis for 𝑟=0 is given by [1, 2] =𝑇(0,𝑧,𝑡)2𝑃𝜀(𝑘𝑡)1/2𝜋𝑎2𝐾𝑧𝑖erfc2(𝑘𝑡)1/2𝑧𝑖erfc2+𝑎21/22(𝑘𝑡)1/2+𝑇0.(3) The dimensionless variables for temperature and time, respectively are defined as [1, 2] 𝜃=𝑇𝑎𝐾𝜋2𝑃𝜀,(4)Γ=(𝑘𝑡)1/2𝑎.(5)

Hence, (3) becomes𝑧𝜃(0,𝑧,𝜏)=Γ𝑖erfc𝑧𝑎Γ𝑖erfc2+𝑎21/2𝑎Γ+𝑇0.(6) For steady-state condition at any depth below or at the center of the focal spot, 𝑡 in (6) [1, 2] 𝑇(0,0,𝑡)=2𝑃𝜀(𝑘𝑡)1/2𝜋𝑎2𝐾1𝜋1/2𝑎𝑖erfc2(𝐾𝑡)1/2+𝑇0(7)1𝜃(0,𝑧,)=𝑎𝑧2+𝑎21/2𝑧.(8) This implies that the maximum surface temperature attainable is given by [1] 𝑇(0,0,)=𝑃𝜀𝜋𝑎𝐾.(9) Figure 2 illustrates a plot diagram for (6) as the variation of the dimensionless temperature 𝜃 with 𝜏=Γ2=4𝑘𝑡/𝑎2 and depth below the focal spot. The temperature rises rapidly at first and approaches 75% of its steady-state value within 𝜏=4 or 𝑡=𝑎2/𝑘, then the change in temperature with time proceeds at a progressively decreasing rate. The time 𝑡=𝑎2/𝐾 can, therefore, be considered as the thermal time constant [1].

Figure 2: Dimensionless temperature 𝜃 versus dimensionless time for various depths when heating with continuous disk (—) and Gaussian (- -) sources [1].

2. Experimental Work

A special nozzle for the laser head was designed and constructed for achieving the optimum performance as shown in Figure 3.

Figure 3: The designed nozzle.

In the absence of the assist gas nozzle, the output power of the laser head was 2.75 W as examined by a power meter type (Gentec TPM 300 CE) as an average value of the output power. In this work, the nozzle was designed in a way that acts as a normal assist gas nozzle and as an assist gas nozzle with variable orifice diameters (with a changeable nozzle tips of 0.4, 0.6, 0.8, 1, 1.2, 1.5, and 3.0 mm in diameters), as shown in Figure 4, that were used as an a apertures (placed at the waist of the collimating lens [3, 4]) to allow choosing various values of the laser beam output power and suppress the higher-order modes.

Figure 4: The assist gas nozzle and the seven different tips.

All the workpieces were illuminated by the laser beam for many different time periods in order to reach the vaporization temperature and executing drilling process. The drilling process was examined for different laser spot positions (at the surface, at the mid surface, and at the lower surface), without the use of assist gas, and with the assist gas (the used pressures were 20, 30, 40, 50, and 60 mmHg).

Figure 5 shows the values of the aspect ratio (depth to diameter) and the taper ratio (outlet to inlet diameter) using laser power of 2.45 W, the focal position at the surface of the workpiece, and the assist gas pressure of 20 mmHg for 11 exposures time.

Figure 5: Different holes drilled by 𝑃=2.45W, focal position at the surface, and with 20 mmHg of assist gas.

Figures 8, 9, and 10 illustrate the temperature distributions for different working conditions.

3. Simulation Results

The maximum temperature in the center of the focus at the surface was measured from (7). This temperature was substitute in (2) for each used output laser power with the rest values mentioned in this equation (𝑎,𝐾,𝑘,𝑡,, etc.) and plotting the temperature distribution. The value of 𝑡 that substituted in (2) was measured as the thermal time constant (𝑡=𝑎2/𝑘) for 𝑧 below the focal spot. This t represents the maximum exposure time needed for the material to reach the vaporization temperature over which the heat diffused inside the material and not along the depth. This value of 𝑡 was used for the rest stages of the drilling process. The hole depth was measured using the plot diagram and substituted in (8) for determining 𝜃. This 𝜃 was substituted in (4) for determining the maximum temperature at this depth, which will be considered as the new surface. The newly measured temperature was substituted in (9) for measuring the new laser power at this depth. The above-mentioned steps were followed again many times until the hole covered the whole thickness of the workpiece. Knowing that each time the new thickness was added to the depth of the previous stage and then substituted in (2).

4. Results and Discussion

The analytical steps were followed for laser output powers 2.45 W, 1.82 W, and 0.96 W using the nozzle orifices of 3 mm, 1.5 mm, and 1.2 mm, the spot radius measured as 0.5 mm, The thermal time constant 𝑡 was 2.25 s, 𝐾 of the workpiece material (PMMA) 0.2 × 10−3 W mm−1 °K−1, 𝑘 of the workpiece material is 0.11 mm2 s−1, and the reflectivity 𝑅 of the workpiece material was assumed as 0.01. Therefore, the emissivity of the surface of the workpiece 𝜀=1𝑅=0.99, and the exposure time versus the surface temperature was plotted for each of the used powers, 2.45 W, 1.82 W, and 0.96 W as shown in Figure 6.

Figure 6: The laser beam exposure time versus the surface temperature.

A MATLAB package was used for presenting the case under study. For one stage drilling process, it was found that whatever the increase in the exposure time the hole did not exceed certain depth of the whole thickness of the workpiece and the heat was diffused (dissipated) inside the material and not along the depth (Figure 7).

Figure 7: The hole depth versus the hole radius.
Figure 8: The first-stage temperature distribution; exposure time (𝑠)=2.25, total power (𝑊)=2.45, and surfacetemperature=1.5364e+004. The second-stage temperature distribution; exposure time (s)=2.25, total power (𝑊)=0.3528, and surfacetemperature=2.2124e+003.
Figure 9: The first-stage temperature distribution; exposure time (s)=2.25, total power (𝑊)=1.82, and surfacetemperature=1.1413e+004. The second-stage temperature distribution; exposure time (𝑠)=2.25, total power (𝑊)=0.2953, and surfacetemperature=1.8518e+003. The third-stage temperature distribution; exposure time (𝑠)=2.25, Total power (𝑊)=0.19, and surfacetemperature=1.2559e+003.
Figure 10: The first-stage temperature distribution; exposure time (𝑠)=2.25, total power (𝑊)=0.96, and surfacetemperature=6.0201e+003. The second-stage temperature distribution; exposure time (𝑠)=2.25, total power (𝑊)=0.1782, and surfacetemperature=1.1175e+003. The third-stage temperature distribution; exposure time (𝑠)=2.25, total power (𝑊)=0.19, and surfacetemperature=1.2559e+003. The forth-stage temperature distribution; exposure time (𝑠)=2.25, total power (𝑊)=0.0989, and surfacetemperature=620.196.

The aspect ratio (ratio of depth to diameter) and the taper ratio (ratio of outlet to inlet diameter) of the drilled holes by using laser beam powers; 2.45 W, 1.82 W, and 0.96 W as measured from Figures 8, 9, and 10 are listed in Table 1.

Table 1: The theoretical values of the aspect and taper ratios at different powers.

These values of 𝐴𝑟 and 𝑇𝑟 clarify that irradiating PMMA workpiece of 2.5 mm thickness by 1064 nm CW diode laser of 1 W output power and exposure time of 2.25 s leads to achieving an acceptable quality of the hole drilling process {that matches the most improved values represented by the highest aspect ratio (20–30) and lowest taper ratio (1) [5, 6].

5. Conclusion

The experimental results for different working conditions that were illustrated in Table 2 show that without the use of the assist gas the best holes can be achieved by focusing almost of 1 W laser output power on the surface of the chosen workpiece material for about 2.5 s which matches with the best result of the analytical model as shown in the last row of the table.

Table 2: The different experimental values of 𝐴𝑟 and 𝑇𝑟 at different working conditions.

Moreover, the use of the low-pressure assist gas of about 60 mmHg with laser output power of 2.45 W and focusing the spot on the surface of the workpiece enhanced the drilling process and reduced the required exposure time. Therefore, this study concluded that for executing laser hole drilling process in black acrylic (PMMA) material of 2.5 mm thickness using CW diode laser of low laser output power around 1 W with the absence of the assist gas, the focal position represents the most affecting parameter for getting best results (highest aspect ratio and lowest taper ratio) while when using higher laser output power of 2.45 W with the existence of the assist gas, the most affecting parameters are the assist gas pressure and focal position.


  1. W. W. Duley, CO2 Laser Effects and Applications, Academic Press, New York, NY, USA, 1976.
  2. H. S. Carslaw and J. C. Jaegar, Conduction of Heat in Solids, Oxford Press, Oxford, UK, 1969.
  3. R. Menzel, Photonics, Linear and Nonlinear Interaction of Light and Matter, Springer, Berlin, Germany, 2001.
  4. J. F. Ready, LIA Handbook of Laser Materials Processing, Laser Institute of America, Magnolia Publishing, Orlando, Fla, USA, 2001.
  5. H. El-Hofy, Advanced Machining Processes, Nontraditional and Hybrid Machining Processes, Production Engineering Department, Alexandria University, Alexandria, Egypt, 2005.
  6. W. M. Steen, Laser Material Processing, Springer, London, UK, 2nd edition, 1998.