A Finite Element Procedure of a Cyclic Thermoviscoplasticity Model for Solder and Copper Interconnects

1998 ◽  
Vol 120 (1) ◽  
pp. 24-34 ◽  
Author(s):  
C. Fu ◽  
D. L. McDowell ◽  
I. C. Ume
Keyword(s):  
Finite Element ◽  
Electronic Packaging ◽  
Copper Interconnects ◽  
Implicit Time ◽  
Integration Scheme ◽  
Single Element ◽  
Ofhc Copper ◽  
Computation Efficiency ◽  
Backward Euler ◽  
Finite Element Procedure

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.

scholarly journals Investigation of Air Pocket Behavior in Pipelines Using Rigid Column Model and Contributions of Time Integration Schemes

Water ◽  
2021 ◽  
Vol 13 (6) ◽  
pp. 785
Author(s):  
Arman Rokhzadi ◽  
Musandji Fuamba
Keyword(s):  
Transient Flow ◽  
Time Integration ◽  
Driving Pressure ◽  
Implicit Time ◽  
Numerical Dissipation ◽  
Integration Scheme ◽  
Time Step ◽  
Column Model ◽  
Integration Schemes ◽  
Backward Euler

This paper studies the air pressurization problem caused by a partially pressurized transient flow in a reservoir-pipe system. The purpose of this study is to analyze the performance of the rigid column model in predicting the attenuation of the air pressure distribution. In this regard, an analytic formula for the amplitude and frequency will be derived, in which the influential parameters, particularly, the driving pressure and the air and water lengths, on the damping can be seen. The direct effect of the driving pressure and inverse effect of the product of the air and water lengths on the damping will be numerically examined. In addition, these numerical observations will be examined by solving different test cases and by comparing to available experimental data to show that the rigid column model is able to predict the damping. However, due to simplified assumptions associated with the rigid column model, the energy dissipation, as well as the damping, is underestimated. In this regard, using the backward Euler implicit time integration scheme, instead of the classical fourth order explicit Runge–Kutta scheme, will be proposed so that the numerical dissipation of the backward Euler implicit scheme represents the physical dissipation. In addition, a formula will be derived to calculate the appropriate time step size, by which the dissipation of the heat transfer can be compensated.

Multiscale Modeling of Large Deformation Processes in Polycrystalline Metals

2003 ◽  
Author(s):  
Antoinette Maniatty ◽  
Karel Matous ◽  
Jing Lu
Keyword(s):  
Finite Element ◽  
Three Dimensional ◽  
Grain Structure ◽  
Current Work ◽  
Integration Scheme ◽  
Two Dimensional ◽  
Scale Bridging ◽  
Grain Structures ◽  
Backward Euler ◽  
Bulk Response

A mesoscale model for predicting the evolution of the grain structure and the mechanical response of polycrystalline aggregates subject to large deformations, such as arise in bulk metal forming processes, is presented. The gain structures modeled are either experimentally observed or are computer generated and statistically similar to experimentally observed grain structures. In order to capture the inhomogeneous deformations and the resulting grain structure characteristics, a discretized model at the mesoscale is used. This work focuses on Al-Mg-Si alloys. Scale bridging is used to link to the macroscale. Examples involving two-dimensional grain structures and current work on three-dimensional grain structures are presented. The present work provides a framework to model the mesoscopic behavior and interactions between grains during finite strains. The mesoscale is characterized by a statistically representative voluem element (RVE), which contains the grains of a polycrystal. Experimentally observed grain structures are used both as models directly (for two-dimensional cases) and to define statistical characteristics to verify the similarity of computer generated grain structures (for three-dimensional cases). A Monte Carlo method based on the Potts model is used to define three-dimensional grain structures. In order to make the representative grain structure appropriate for scale-bridging, we design them with periodicity. A three-field, updated Lagrangian finite element formulation with a kinematic split of the deformation gradient into volume preserving and volumetric parts is used to create a stable finite element method in the context of nearly incompressible behavior. A fully implicit two-level backward Euler integration scheme is derived for integrating the constitutive equations, and consistent linearization is used in Newton’s method to solve the resulting equations. In addition, the average of the boundary conditions and bulk response must match the macroscopically measured bulk response. To illustrate and verify the proposed model, we analyze examples involving two-dimensional grain structures and compare with results from a Taylor model. Current work on three-dimensional grain structures ara also presented.

Failure and Reliability Analysis Influenced by Various Factors Using Accelerated Disturbed State Concept

2003 ◽  
Author(s):  
Russell D. Whitenack ◽  
Chandra S. Desai
Keyword(s):  
Finite Element ◽  
Reliability Analysis ◽  
Constitutive Modeling ◽  
Electronic Packaging ◽  
Analysis Tool ◽  
Disturbed State Concept ◽  
Powerful Approach ◽  
Finite Element Procedure ◽  
Cycles To Failure

The disturbed state concept (DSC) presented herein represents a unified and powerful approach for constitutive modeling of materials and interfaces in electronic packaging. Together with the computer finite element procedure it provides an analysis tool for calculation of stresses, strains, disturbance and cycles to failure. The accelerated procedure allows economical approximation of cycles to failure and distribution of disturbance at different cycles for design and reliability.

Modeling of the Long-Term Behavior of Glassy Polymers

2012 ◽  
Vol 135 (1) ◽  
Author(s):  
Sami Holopainen ◽  
Mathias Wallin
Keyword(s):  
Experimental Data ◽  
Finite Element ◽  
Glassy Polymers ◽  
Integration Scheme ◽  
Element Formulation ◽  
Algebraic Equations ◽  
Plastic Spin ◽  
Backward Euler ◽  
Long Term Behavior

The constitutive model for glassy polymers proposed by Arruda and Boyce (BPA model) is reviewed and compared to experimental data for long-term loading. The BPA model has previously been shown to capture monotonic loading accurately, but for unloading and long-term behavior, the response of the BPA model is found to deviate from experimental data. In the present paper, we suggest an efficient extension that significantly improves the predictive capability of the BPA model during unloading and long-term recovery. The new, extended BPA model (EBPA model) is calibrated to experimental data of polycarbonate (PC) in various loading–unloading situations and deformation states. The numerical treatment of the BPA model associated with the finite element analysis is also discussed. As a consequence of the anisotropic hardening, the plastic spin enters the model. In order to handle the plastic spin in a finite element formulation, an algorithmic plastic spin is introduced. In conjunction with the backward Euler integration scheme use of the algorithmic plastic spin leads to a set of algebraic equations that provides the updated state. Numerical examples reveal that the proposed numerical algorithm is robust and well suited for finite element simulations.

scholarly journals Development of a continuum plasticity model for the commercial finite element code ABAQUS

2011 ◽  
Vol 2 (2) ◽  
pp. 275-283
Author(s):  
M. Safaei ◽  
W. De Waele
Keyword(s):  
Finite Element ◽  
Yield Surface ◽  
Euler Method ◽  
Integration Scheme ◽  
Isotropic Hardening ◽  
Implicit Integration ◽  
Backward Euler ◽  
User Material Subroutine ◽  
Von Mises ◽  
Future Work

The present work relates to the development of computational material models for sheet metalforming simulations. In this specific study, an implicit scheme with consistent Jacobian is used forintegration of large deformation formulation and plane stress elements. As a privilege to the explicitscheme, the implicit integration scheme is unconditionally stable. The backward Euler method is used toupdate trial stress values lying outside the yield surface by correcting them back to the yield surface atevery time increment. In this study, the implicit integration of isotropic hardening with the von Mises yieldcriterion is discussed in detail. In future work it will be implemented into the commercial finite element codeABAQUS by means of a user material subroutine.

Thermomechanical Finite Element Analysis of Problems in Electronic Packaging Using the Disturbed State Concept: Part 1—Theory and Formulation

1998 ◽  
Vol 120 (1) ◽  
pp. 41-47 ◽  
Author(s):  
C. Basaran ◽  
C. S. Desai ◽  
T. Kundu
Keyword(s):  
Finite Element ◽  
Constitutive Model ◽  
Electronic Packaging ◽  
Constitutive Models ◽  
Cyclic Behavior ◽  
Element Analysis ◽  
Theoretical Formulation ◽  
Finite Element Procedure ◽  
Reliable Design ◽  
Number Of Cycles

Accurate prediction of the thermomechanical cyclic behavior of joints and interfaces in semiconductor devices is essential for their reliable design. In order to understand and predict the behavior of such interfaces there is a need for improved and unified constitutive models that can include elastic, inelastic, viscous, and temperature dependent microstructural behavior. Furthermore, such unified material models should be implemented in finite element procedures so as to yield accurate and reliable predictions of stresses, strains, deformations, microcracking, damage, and number of cycles to failure due to thermomechanical loading. The main objective of this paper is to present implementation of such an unified constitutive model in a finite element procedure and its application to typical problems in electronic packaging; details of the constitutive model are given by Desai et al. (1995). Details of the theoretical formulation is presented in this Part 1, while its applications and validations are presented in Part 2, Basaran et al. (1998).

scholarly journals Nonlinear finite element system simulation of piezoelectric vibration-based energy harvesters

2021 ◽  
pp. 1045389X2110482
Author(s):  
Andreas Hegendörfer ◽  
Paul Steinmann ◽  
Julia Mergheim
Keyword(s):  
Finite Element ◽  
System Simulation ◽  
Electric Circuit ◽  
Software Tool ◽  
Simulation Method ◽  
Implicit Time ◽  
Integration Scheme ◽  
Power Circuit ◽  
Energy Harvesters ◽  
Piezoelectric Vibration

Piezoelectric vibration-based energy harvesters consist of an electromechanical structure and an electric circuitry, influencing each other. We propose a novel approach that allows a finite element based system simulation of nonlinear electromechanical structures coupled to nonlinear electric circuitries. In the finite element simulation the influence of the electric circuit on the electromechanical structure is considered via the vector of external forces, using an implicit time integration scheme. To demonstrate the applicability of the new simulation method an active power circuit is considered. Several examples of piezoelectric vibration-based energy harvesters, connected to standard or synchronized switch harvesting on inductor (SSHI) circuits, showing linear or nonlinear mechanical behavior, are studied to validate the proposed simulation method against numerical results reported in the literature. The advocated method allows for consistent and efficient simulations of complete nonlinear energy harvesters using only one software tool.

Modeling of 3D Magnetostrictive Systems with Application to Galfenol and Terfenol-D Actuators

2012 ◽  
Vol 77 ◽  
pp. 11-28
Author(s):  
Marcelo J. Dapino ◽  
Suryarghya Chakrabarti
Keyword(s):  
Finite Element ◽  
Time Integration ◽  
Three Dimensional ◽  
Nonlinear Dynamic System ◽  
Element Model ◽  
Implicit Time ◽  
Integration Scheme ◽  
Modeling Framework ◽  
Time Integration Scheme ◽  
Parametric Analyses

This work presents a unified approach to model three dimensional magnetostrictive transducers. Generalized procedures are developed for incorporating nonlinear coupled constitutive behavior of magnetostrictive materials into an electro-magneto-mechanical finite element modeling framework. The finite element model is based on weak forms of Maxwell's equations for electromagnetics and Navier's equations for mechanical systems. An implicit time integration scheme is implemented to obtain nonlinear dynamic system responses. The model is implemented into a finite element (FE) solver and applied to two case studies, a Galfenol unimorph actuator and a magnetohydraulic Terfenol-D actuator for active engine mounts. Model results are compared with experiments, and parametric analyses are conducted which provide guidelines for optimization of actuator design.

scholarly journals Two Scales, Hybrid Model for Soils, Involving Artificial Neural Network and Finite Element Procedure

2015 ◽  
Vol 36 (2) ◽  
pp. 29-36
Author(s):  
Marcin Krasiński ◽  
Marek Lefik
Keyword(s):  
Finite Element ◽  
Engineering Properties ◽  
Single Element ◽  
Laboratory Sample ◽  
Engineering Problems ◽  
Geotechnical Problem ◽  
Finite Element Procedure ◽  
Correct Representation ◽  
Constitutive Parameters ◽  
Level Analysis

Abstract A hybrid ANN-FE solution is presented as a result of two level analysis of soils: a level of a laboratory sample and a level of engineering geotechnical problem. Engineering properties of soils (sands) are represented directly in the form of ANN (this is in contrast with our former paper where ANN approximated constitutive relationships). Initially the ANN is trained with Duncan formula (Duncan and Chang [2]), then it is re-trained (calibrated) with some available experimental data, specific for the soil considered. The obtained approximation of the constitutive parameters is used directly in finite element method at the level of a single element at the scale of the laboratory sample to check the correct representation of the laboratory test. Then, the finite element that was successfully tested at the level of laboratory sample is used at the macro level to solve engineering problems involving the soil for which it was calibrated.

Sign in / Sign up

Export Citation Format

Share Document