0
Research Papers

Numerical Analysis of Ultrashort Pulse Laser-Material Interaction Using ABAQUS

[+] Author and Article Information
Dongkyun Lee

Department of Mechanical Engineering, University of Michigan, Ann Arbor, MI 48105-2125dongkyun@umich.edu

Elijah Kannatey-Asibu

Department of Mechanical Engineering, University of Michigan, Ann Arbor, MI 48105-2125asibu@umich.edu

J. Manuf. Sci. Eng 131(2), 021005 (Mar 06, 2009) (15 pages) doi:10.1115/1.3075869 History: Received April 12, 2007; Revised December 17, 2008; Published March 06, 2009

Ultrafast lasers of subpicosecond pulse duration have thus far been investigated for ablation, drilling, and cutting processes. Ultrafast lasers also have the potential for laser welding of small components of the order of microns and for laser shock peening to enhance the peening depth. In this paper, the potential for welding is investigated using the two-temperature model to analyze laser-matter interaction for ultrafast lasers. The model is implemented in a general-purpose commercial finite element method package, ABAQUS , to enable broad based application of the two-temperature model in practical engineering problems. The implementation is validated by comparison with linear solutions obtained using separation of variables. It is then used to investigate the potential for microwelding using ultrafast laser pulse for both low and high fluence laser input.

FIGURES IN THIS ARTICLE
<>
Copyright © 2009 by American Society of Mechanical Engineers
Your Session has timed out. Please sign back in to continue.

References

Figures

Grahic Jump Location
Figure 1

The concept of dual domain configuration for the TTM in ABAQUS

Grahic Jump Location
Figure 2

(a) Temperature information sharing via common memory block in the user subroutine, USDFLD , shown as USD in the figure: The ◇ marks indicate that the subroutine is called at the integration points of 4-node planar elements, Ωen and Ωlm; (b) effective heat capacity with the latent heat (Area B) for phase change in USDFLD

Grahic Jump Location
Figure 3

Four-node linear planar finite element: (a) actual coordinates and (b) normalized coordinates

Grahic Jump Location
Figure 4

A dual domain setup of the TTM implementation for ABAQUS (Lx=200 nm, Δx=1 nm)

Grahic Jump Location
Figure 5

(a) Lattice temperature histories at the top surface for heat sources of different temporal profiles and (b) corresponding temperature distributions at select times; (c) lattice temperature histories at the top surface for heat sources of different spatial distributions and (d) corresponding temperature distributions at select times; the legend label “FEM” stands for results from ABAQUS and “Analytic” for results from analytical series solutions

Grahic Jump Location
Figure 6

Normalized homogeneous electron temperature histories for a laser pulse of tp=96 fs at the top and back surface for different material thicknesses of (a) 100 nm and (b) 200 nm. Data set “QT” taken from Ref. 30. In the legend, dt stands for Δt, and tp for tp.

Grahic Jump Location
Figure 7

Temperature distributions inside a material (Lx=100 nm) for a laser pulse of tp=100 fs at select times: (a) electron temperature and (b) lattice temperature. Data set “QT” taken from Ref. 30. In the legend, dt stands for Δt and tp for tp.

Grahic Jump Location
Figure 8

(a) Homogeneous and (b) hybrid element configurations (Lx=200 nm)

Grahic Jump Location
Figure 9

Comparison of FEM results for homogeneous and hybrid elements: (a) electron and lattice temperature histories at select locations and (b) electron and lattice temperature distributions at select times. Marks in (b) represent the locations of nodes.

Grahic Jump Location
Figure 10

(a) Lattice temperatures at the top surface and ablation depth histories and (b) lattice temperature distributions for t=2.0 and 22.0 ps, J=800 mJ/cm2. In the legends, dT stands for ΔTm, dt for Δt.

Grahic Jump Location
Figure 11

(a) Ablation depth with respect to the input fluence. Data set “CLB” taken from Ref. 34 and “Exp” from Ref. 14. (b) Corresponding ablation starting and ending times and average ablation rate with respect to the input fluence.

Grahic Jump Location
Figure 12

(a) Molten pool depth change history traced with melting point for select fluences and (b) estimation on molten pool thickness traced with melting point (Tm in the legend) and the liquidus (Tm+dT in the legend), with respect to the input fluence

Tables

Errata

Discussions

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In