A finite element procedure using a semi-implicit time-integration scheme has been developed for a cyclic thermoviscoplastic constitutive model for Pb-Sn solder and OFHC copper, two common metallic constituents in electronic packaging applications. The scheme has been implemented in the commercial finite element (FE) code ABAQUS (1995) via the user-defined material subroutine, UMAT. Several single-element simulations are conducted to compare with previous test results, which include monotonic tensile tests, creep tests, and a two-step ratchetting test for 62Sn36Pb2Ag solder; a nonproportional axial-torsional test and a thermomechanical fatigue (TMF) test for OFHC copper. At the constitutive level, we also provide an adaptive time stepping algorithm, which can be used to improve the overall computation efficiency and accuracy especially in large-scale FE analyses. We also compare the computational efforts of fully backward Euler and the proposed methods. The implementation of the FE procedure provides a guideline to apply user-defined material constitutive relations in FE analyses and to perform more sophisticated thermomechanical simulations. Such work can facilitate enhanced understanding thermomechanical reliability issue of solder and copper interconnects in electronic packaging applications.

1.
ABAQUS, 1995, Hibbitt, Karlsson & Sorenson, Inc., Pawtucket, RI, Verson 5.4.
2.
Argyris
 
J. H.
,
Vaz
 
L. E.
, and
William
 
K. J.
,
1978
, “
Improved Solution Methods for Inelastic Rate Problems
,”
Computer Methods in Applied Mechanics and Engineering
, Vol.
16
, pp.
231
277
.
3.
Busso
 
E. P.
,
Kitano
 
M.
, and
Kumazawa
 
T.
,
1994
, “
A Forward Gradient Time Integration Procedure for an Internal Variable Constitutive Model of Sn-Pb Solder
,”
International Journal for Numerical Methods in Engineering
, Vol.
37
, pp.
539
558
.
4.
Cormeau
 
I.
,
1975
, “
Numerical Stability in Quasi-Static Elasto/Visco-Plasticity
,”
International Journal for Numerical Methods in Engineering
, Vol.
9
, pp.
109
127
.
5.
Dombrovsky
 
L. A.
,
1992
, “
Incremental Constitutive Equations for Miller and Bodner-Partom Viscoplastic Models
,”
Computer & Structures
, Vol.
44
, No.
5
, pp.
1065
1072
.
6.
Freed, A. D., 1988, “Thermoviscoplasticity Model With Application to Copper,” NASA Tech. Paper 2845.
7.
Fu, C., Ume, I.C., and McDowell, D. L., 1996, “Thermomechanical Stress Analyses of Plated-Through Holes in PWB Using Internal State Variable Constitutive Models,” Sensing, Modeling and Simulation in Emerging Electronic Packaging, Proc. 1996 ASME International Mechanical Engineering Conference and Exposition, pp. 65–72, Atlanta, GA.
8.
Hughes, T. J. R., 1984, “Numerical Implementation of Constitutive Models: Rate-Independent Deviatoric Plasticity,” Theoretical Foundation for Large-Scale Computations of Nonlinear Material Behavior, Martinus Nijhoff Publisher, Boston, MA, pp. 29–57.
9.
Kumar
 
V.
,
Morjaria
 
M.
, and
Mukherjee
 
S.
,
1980
, “
Numerical Integration of Some Stiff Constitutive Models of Inelastic Deformation
,”
ASME Journal of Engineering Materials and Technology
, Vol.
102
, pp.
92
96
.
10.
Lamba
 
H. S.
, and
Sidebottom
 
O. M.
,
1978
a, “
Cyclic Plasticity for Nonproportional paths: Part 1—Cyclic Hardening, Erasure of Memory, and subsequent Strain Hardening Experiments
,”
ASME Journal of Engineering Materials and Technology
, Vol.
100
, pp.
96
103
.
11.
Lamba
 
H. S.
, and
Sidebottom
 
O. M.
,
1978
b, “
Cyclic Plasticity for Nonproportional Paths: Part 2—Comparison With Predictions of Three Incremental Plasticity Models
,”
ASME J. Engng. Matls. Tech.
, Vol.
100
, pp.
104
111
.
12.
Lush
 
A. M.
,
Weber
 
G.
, and
Anand
 
L.
,
1989
, “
An Implicit Time-Integration Procedure for a Set of Internal Variable Constitutive Equations for Isotropic Elasto-Viscoplasticity
,”
International Journal of Plasticity
, Vol.
5
, pp.
521
549
.
13.
Marin, E. B., 1994, “A Critical Study of Finite Strain Porous Inelasticity,” Ph.D. thesis, School of Mechanical Engineering, Georgia Institute of Technology, Atlanta, GA.
14.
McDowell
 
D. L.
,
1992
, “
A Nonlinear Kinematic Hardening Theory for Cyclic Thermoplasticity and Thermoviscoplasticity
,”
International Journal of Plasticity
, Vol.
8
, pp.
695
728
.
15.
McDowell, D. L., Miller, M. P., and Brooks, D. C., 1994, “A Unified Creep-Plasticity Theory for Pb-Sn Solder,” ASTM STP 153, S. A. Schroeder and M.R. Mitchell, eds., ASTM, Philadelphia, pp. 42–59.
16.
Narasimhan
 
R.
,
Rosakis
 
A. J.
, and
Moran
 
B.
,
1992
, “
A Three-Dimensional Numerical Investigation of Fracture Initiation by Ductile Failure Mechanisms in a 4340 Steel
,”
Advances in Fracture/Damage Models for the Analysis of Engineering Problems
, ASME AMD, Vol.
137
, pp.
13
53
.
17.
Odqvist, F. K. G., 1974, “Mathematical Theory of Creep and Creep Rupture,” 2nd ed., Clarendon Press, Oxford, London, pp. 20–35.
18.
Ortiz
 
M.
, and
Popov
 
E. P.
,
1985
, “
Accuracy and Stability of Integration Algorithms for Elasto-Plastic Constitutive Relations
,”
International Journal for Numerical Methods in Engineering
, Vol.
21
, pp.
1561
1576
.
19.
Pierce
 
D.
,
Shih
 
C. F.
, and
Needleman
 
A.
,
1984
, “
A Tangent Modulus Method for Rate Dependent Solids
,”
Computers & Structures
, Vol.
18
, pp.
875
887
.
20.
Sotolongo
 
W.
, and
McDowell
 
D. L.
,
1986
, “
On The Numerical Integration Of Elasto-Plastic Constitutive Model Structures for Nonproportional Cyclic Loading
,”
Computers & Structures
, Vol.
24
, No.
4
, pp.
595
605
.
21.
Tanaka
 
T. G.
, and
Miller
 
A. K.
,
1988
, “
Development of a Method for Integrating Time-Dependent Constitutive Equations With Large, Small, or Negative Strain Rate Sensitivity
,”
International Journal for Numerical Methods in Engineering
, Vol.
26
, pp.
2457
2485
.
This content is only available via PDF.
You do not currently have access to this content.