• Structural Engineering and Mechanics
• /
• v.47 no.1
• /
• pp.27-44
• /
• 2013

A finite element computational procedure for the accurate analysis of quasistatic thermorheological complex structures response is developed. The geometrical nonlinearity, arising from large displacements and rotations (but small strains), is accounted for by the total Lagrangian description of motion. The Schapery's nonlinear single-integral viscoelastic constitutive model is modified for a time-stress-temperature-dependent behavior. The nonlinear thermo-viscoelastic constitutive equations are incrementalized leading to a recursive relationship and thereby the resulting finite element equations necessitate data storage from the previous time step only, and not the entire deformation history. The Newton-Raphson iterative scheme is employed to obtain a converged solution for the non-linear finite element equations. The developed numerical model is verified with the previously published works and a good agreement with them is found. The applicability of the developed model is demonstrated by analyzing two examples with different thermal/mechanical loading histories.

• Transactions of Materials Processing
• /
• v.19 no.8
• /
• pp.468-473
• /
• 2010

Shape memory alloys (SMAs) have the ability to recover their original shape upon thermo-mechanical loading even after large inelastic deformation. The unique feature is known as pseudoelasticity and shape memory effect caused by the crystalline structural transformation between two solid-state phases called austenite and martensite. To support the engineering application, a number of constitutive models, which can be formally classified into either micromechanics-based or phenomenological model, have been developed. Most of the constitutive models include a kinetic law governing the crystallographic transformation. The present work presents a one-dimensional, phenomenological constitutive model for SMAs in the context of the unified viscoplasticity theory. The proposed model does not incorporate the complex mechanisms of phase transformation. Instead, the effects induced by the transformation are depicted through the growth law for the back stress that is an internal state variable of the model.

• Transactions of Materials Processing
• /
• v.13 no.2
• /
• pp.168-173
• /
• 2004

Microforming can be a good application for bulk metallic glasses. It is important to simulate the deformation behaviour of the bulk metallic glasses in a supercooled liquid region for manufacturing micromachine parts. For these purposes, a correct constitutive model which can reproduce viscosity results is essential for good predicting capability. In this paper, we studied deformation behaviour of the bulk metallic glasses using the finite element method in conjunction with the fictive stress constitutive model which can describe non-Newtonian as well as Newtonian behaviour. A combination of kinetic equation which describes the mechanical response of the bulk metallic glasses at a given temperature and evolution equations fur internal variables provide the constitutive equation of the fictive stress model. The internal variables are associated with fictive stress and relation time. The model has a modular structure and can be adjusted to describe a particular type of microforming process. Implementation of the model into the MARC software has shown its versatility and good predictive capability.

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

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. The complexity of these constitutive equations generally requires a stable and accurate numerical method. The radial return mapping is one of the most robust integration scheme currently used. Nonlinear kinematic hardening model of Armstrong-Fredrick type has recovery term and the direction of kinematic hardening increment is not parallel to that of plastic strain increment. In this case, The conventional radial return mapping method cannot be applied directly. In this investigation, we expanded the radial return mapping method to consider the nonlinear kinematic hardening model and implemented this integration scheme 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 Newton method and bisection method. Using dynamic yield condition derived from linearization of flow rule, the integration scheme for elastoplastic and viscoplastic constitutive model was unified. Several numerical examples are considered to demonstrate the efficiency and applicability of the present method.

• Proceedings of the Korean Powder Metallurgy Institute Conference
• /
• 1999.04a
• /
• pp.15-16
• /
• 1999

• Honam Mathematical Journal
• /
• v.36 no.1
• /
• pp.11-27
• /
• 2014

In this paper, we propose closures for multi-phase flow models, which satisfy boundary conditions and conservation constraints. The models governing the evolution of the fluid mixing are derived by applying an ensemble averaging procedure to the microphysical equations characterized by distinct phases. We consider compressible multi species multi-phase flow with surface tension and transport.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.9
• /
• pp.1227-1238
• /
• 2000

This study examines the effects of the discretization schemes on numerical solutions of viscoelastic fluid flows. For this purpose, a temporally evolving mixing layer, a two-dimensional vortex pair interacting with a wall, and a turbulent channel flow are selected as the test cases. We adopt a fourth-order compact scheme (COM4) for polymeric stress derivatives in the momentum equations. For convective derivatives in the constitutive equations, the first-order upwind difference scheme (UD) and artificial diffusion scheme (AD), which are commonly used in the literature, show most stable and smooth solutions even for highly extensional flows. However, the stress fields are smeared too much and the flow fields are quite different from those obtained by higher-order upwind difference schemes for the same flow parameters. Among higher-order upwind difference schemes, a third-order compact upwind difference scheme (CUD3) shows most stable and accurate solutions. Therefore, a combination of CUD3 for the convective derivatives in the constitutive equations and COM4 for the polymeric stress derivatives in the momentum equations is recommended to be used for numerical simulation of highly extensional flows.

• Structural Engineering and Mechanics
• /
• v.64 no.5
• /
• pp.611-620
• /
• 2017

A new theory of weightless sagging planer elasto-flexible cables under point loads is developed earlier by the authors and used for predicting the nonlinear dynamic response of cable-suspended linear elastic beams. However, this theory is not valid for nonlinear elastic cracked concrete beams possessing different positive and negative flexural rigidity. In the present paper, the flexural response of simply supported cracked concrete beams suspended from cables by two hangers is presented. Following a procedure established earlier, rate-type constitutive equations and third order nonlinear differential equations of motion for the structures undergoing small elastic displacements are derived. Upon general quasi-static loading, negative nodal forces, moments and support reactions may be introduced in the cable-suspended concrete beams and linear modal frequencies may abruptly change. Subharmonic resonances are predicted under harmonic loading. Uncoupling of the nodal response is proposed as a more general criterion of crossover phenomenon. Significance of the bilinearity ratio of the concrete beam and elasto-configurational displacements of the cable for the structural response is brought out. The relevance of the proposed theory for the analysis and the design of the cable-suspended bridges is critically evaluated.

• The Journal of the Acoustical Society of Korea
• /
• v.20 no.4
• /
• pp.54-61
• /
• 2001

Magnetostrictive materials have nonlinear elasto-magnetic properties. However the constitutive equations to describe the nonlinear properties are not available, yet. In this study we develope the equation in magnetostrictive materials by use of piezomagnetic constitutive equation which is quasi-linearized. With the wave equation, we determine the propagation velocity inside the magnetostrictive materials when a plane wave propagates along a given magnetic field. Validity of the calculated velocity is verified through comparison with experimental velocity measurement results for the most representative magnetostrictive materials. Terfenol-D.

• Structural Engineering and Mechanics
• /
• v.39 no.2
• /
• pp.169-186
• /
• 2011

This paper is devoted to the development of the equations which describe the elastic response of a curved laminate subjected to in-plane loads and bending moments. Similar to the classic $6{\times}6$ ABD matrix constitutive relation of a flat laminate, a new $6{\times}6$ matrix constitutive relation between force resultants, moment resultants, mid-plane strains and deformed curvatures for a curved laminate is formulated. This curved lamination theory will provide the fundamental basis for the analyses of curved laminated structures. The stress predictions by the present curved lamination theory are compared to those by the curved laminate analysis that neglected the nonlinear terms in the derivation of the constitutive relation. The results show that the curved laminate analysis that neglected the nonlinear terms cannot reflect the effect of curvature and can no longer predict the stresses accurately as the curvature becomes noticeable. In this paper, a curved lamination theory that retains the nonlinear terms and, therefore, accounts for the effect of the non-flat geometry of the structure will be developed.