• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.9
• /
• pp.1555-1561
• /
• 2003

Finite element analysis using modern constitutive equation is one of the most general tools to simulate the deformation behavior and to predict the life of the structure. Constitutive equation becomes complicated so as to predict the material behavior more accurately than the classical models. Because of the complexity of constitutive model, numerical treatment becomes so difficult that the calculation should be verified carefully. One-element tests, simple tension or simple shear, are usually used to verify the accuracy of finite element analysis using complicated constitutive model. Since this test is mainly focused on the time integration scheme, it is also necessary to verify the equilibrium iteration using material stiffness matrix and to compare FE results with solution of structures. In this investigation, viscoplastic solution of thick walled cylinder was derived considering axial constraints and was compared with the finite element analysis. All the numerical solutions showed a good coincidence with FE results. This numerical solution can be used as a verification tool for newly developed FE code with complicated constitutive model.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.12 no.3
• /
• pp.139-144
• /
• 2004

Dynamic crack initiation in ductile steel is investigated by means of impact loaded 3 point bend(PB) specimens. Results from non-viscoplastic and viscoplastic materials are compared. Their materials are applied with various impact velocities and static strain rates. The specimen has the size 320${\times}$750 mm with a thickness of 10 mm. A modified 3PB specimen design with reduced width at the ends has been developed in order to avoid the initial compressive load of the crack tip and also to avoid the uncertain boundary conditions at the impact heads. Numerical simulations are made by using the FEM code ABAQUS. Therefore, their results are plotted by shapes of the von Mises plastic stress and equivalent plastic strain of the specimens applied by various impact velocities.

• Transactions on Electrical and Electronic Materials
• /
• v.5 no.6
• /
• pp.215-218
• /
• 2004

This paper presents a three-dimensional numerical simulation for thermal oxidation process. A new elasto-viscoplastic model for robust numerical oxidation simulation is proposed. The three-dimensional effects of oxidation process such as mask lifting effect and corner effects are analyzed. In nano-scale process, the oxidant diffusion is punched through to the other side of the mask. The mask is lifted so the thickness of oxide region is greatly enhanced. The compressive pressure during the oxidation is largest in the mask corner of the island structure. This is because the masked area near the corner is surrounded by an area larger than the others in the island structure. This stress induces the retardation of the oxide growth, especially at the masked corner in the island structure.

• Journal of Theoretical and Applied Mechanics
• /
• v.3 no.1
• /
• pp.45-65
• /
• 2002

A constitutive model on oorthotropic thermo-elasto-viscoplasticity for fiber-reinforced composite materials Is illustrated, and their thermomechanical responses are predicted with the fully-coupled finite element formulation. The unmixing-mixing scheme can be adopted with the multipartite matrix method as the constitutive model. Basic assumptions based upon the composite micromechanics are postulated, and the strain components of thermal expansion due to temperature change are included In the formulation. Also. more than two sets of mechanical variables, which represent the deformation states of multipartite matrix can be introduced arbitrarily. In particular, the unmixing-mixing scheme can be used with any well-known isotropic viscoplastic theory of the matrix material. The scheme unnecessitates the complex processes for developing an orthotropic viscoplastic theory. The governing equations based on fully-coupled thermomechanics are derived with constitutive arrangement by the unmixing-mixing concept. By considering some auxiliary conditions, the Initial-boundary value problem Is completely set up. As a tool of numerical analyses, the finite element method Is used with isoparametric Interpolation fer the displacement and the temperature fields. The equation of mutton and the energy conservation equation are spatially discretized, and then the time marching techniques such as the Newmark method and the Crank-Nicolson technique are applied. To solve the ultimate nonlinear simultaneous equations, a successive iteration algorithm is constructed with subincrementing technique. As a numerical study, a series of analyses are performed with the main focus on the thermomechanical coupling effect in composite materials. The progress of viscoplastic deformation, the stress-strain relation, and the temperature History are careful1y examined when composite laminates are subjected to repeated cyclic loading.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.2
• /
• pp.488-499
• /
• 1996

SMC(Sheet molding compound) is a thermosetting material reinforced with chopped fiberglass. The compression molding of SMC was analyzed based on a rigid thermo-viscoplastic approach using a three dimensional finite element program coupled with temperatures. Only the temperature analysis part was tested in this paper by solving one-dimensional heat transfer problem and comparing with the exact solutions available in the literature. Based on this comparison the program was proved to be valid and was further applied in solving compression molding of SMC between flat dies. To investigate the usefulness of a rigid thermo-viscoplastic approach in the compression molding analysis of SMC charge, compression of rectangular shaped SMC charge at plane strain and three dimensionalde formation condition was analyzed under the same condition as given in the literature. From this comparison it was found out that the rigid thermo-viscoplastic approach was useful in analyzing SMC compression molding between flat dies.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.24 no.1
• /
• pp.97-106
• /
• 2011

In this study, viscoplastic material properties of solder alloy which is used in the electronics packages are inversely estimated. A specimen is fabricated to this end, and an experiment is conducted to examine deformation by Moir$\acute{e}$ interferometry. As a result of the experiment, bending displacement of the specimen and shear strain of the solder are obtained. A viscoplastic finite element analysis procedure is established, and the material parameters are determined to match closely with the experiments. The uncertainties which include inherent experimental error and insufficient data of experiments are addressed by using the method of computer model calibration. As a result, material parameters are identified in the form of confidence interval, and the displacements and strains using these parameters are predicted in the form the prediction interval.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.11
• /
• pp.3516-3524
• /
• 1996

Service conditions for structures at elevated temperatures in nuclear power plant involve transient thermal and mechanical load levels that are severe enough to caeuse inelastic deformations due to creep and plasticity. Therefore, a systematic mehtod of inelastic analysis is needed for the design of structural components in nuclear poser plants subjected to such loading conditions. In the present investigation, the Chabodhe model, one of the unified viscoplastic constitutive equations, was selected for systematic inelastic analysis. The material response was integrated based on GMR ( generallized mid-point rule) time integral scheme and provided to ABAQUS as a material subroutine, UMAT program. By comparing results obtaned from uniaxial analysis using the developed UMAT program with those from Runge-Kutta solutions and experimentaiton, the validity of the adopted Chaboche model and the numerical stability and accuracy of the developed UMAT program were verified. In addition, the developed material subroutine was applied for uniaxial creep and tension analyses for the plate with a hole in the center. The application further demonstrates usefulness of the developed program.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.33 no.4
• /
• pp.380-387
• /
• 2009

The reliability concerns of solder interconnections in flip chip PBGA packages are produced mainly by the mismatch of coefficient of thermal expansion(CTE) between the module and PCB. Finite element analysis has been employed extensively to simulate thermal loading for solder joint reliability and deformation of packages in electronic packages. The objective of this paper is to study the thermo-mechanical behavior of FC-PBGA package assemblies subjected to temperature change, with an emphasis on the effect of the finite element model, material models and temperature conditions. Numerical results are compared with the experimental results by using $moir{\acute{e}}$ interferometry. Result shows that the bending displacements of the chip calculated by the finite element analysis with viscoplastic material model is in good agreement with those by $moir{\acute{e}}$ inteferometry.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.5 s.176
• /
• pp.1075-1083
• /
• 2000

The prediction of the inelastic behavior of the structure is an essential part of reliability assessment procedure, because most of the failures are induced by the inelastic deformation, such as creep and plastic deformation. During decades, there has been much progress in understanding of the inelastic behavior of the materials and a lot of inelastic constitutive equations have been developed. These equations consist of the definition of inelastic strain and the evolution of the state variables introduced to quantify the irreversible processes occurred in the material. With respect to the definition of the inelastic strain, the inelastic constitutive models can be categorized into elastoplastic model, unified viscoplastic model and separated viscoplastic model and the different integration methods have been applied to each category. In the present investigation, the generalized integration method applicable for various types of constitutive equations is developed and implemented into ABAQUS by means of UMAT subroutine. The solution of the non-linear system of algebraic equations arising from time discretization with the generalized midpoint rule is determined using line-search technique in combination with Newton method. The strategy to control the time increment for the improvement of the accuracy of the numerical integration is proposed. Several numerical examples are considered to demonstrate the efficiency and applicability of the present method. The prediction of the inelastic behavior of the structure is an essential part of reliability assessment procedure, because most of the failures are induced by the inelastic deformation, such as creep and plastic deformation. During decades, there has been much progress in understanding of the inelastic behavior of the materials and a lot of inelastic constitutive equations have been developed. These equations consist of the definition of inelastic strain and the evolution of the state variables introduced to quantify the irreversible processes occurred in the material. With respect to the definition of the inelastic strain, the inelastic constitutive models can be categorized into elastoplastic model, unified viscoplastic model and separated viscoplastic model and the different integration methods have been applied to each category. In the present investigation, the generalized integration method applicable for various types of constitutive equations is developed and implemented into ABAQUS by means of UMAT subroutine. The solution of the non-linear system of algebraic equations arising from time discretization with the generalized midpoint rule is determined using line-search technique in combination with Newton method. The strategy to control the time increment for the improvement of the accuracy of the numerical integration is proposed. Several numerical examples are considered to demonstrate the efficiency and applicability of the present method.

• The Journal of Engineering Geology
• /
• v.3 no.3
• /
• pp.215-230
• /
• 1993

In-situ rock mass demonstrates the variety of structural features, and especially the mechanical and spatial characteristics of joint (or joint system) greatly affect the deformation and fallure strength of the rock mass. In this study finite element model capable of analyzing the viscoplastic behavior of reinforced jointed rock mass has been developed based on equivalent material approach. Accuracy and reliability of the numerical model have verified by simuiating the behavior of simplified block model and comparing the results with analytic solutions. Practical applicability was also demonstrated by analyzing the time-dependent behavior of underground oil storage tunnel and assessing the reinforcement effect of rockbolt.