• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.28 no.2
• /
• pp.133-142
• /
• 2004

In this article, the constitutive relations of three types of piezoelectric benders, which are a unimorph bender, a bimorph bender and a triple-layer bender, are derived based on the beam theory under the quasi-static equilibrium condition. The relation coefficients are described as the geometry and material properties of the benders. More general constitutive relations involving fixed-free, fixed-roll, and fixed-simply supported boundary conditions under the inconsistent length condition between the piezoelectric layer and the nonpiezoelectric one are discussed. The complicated constitutive relations can be easily calculated and checked by using the symbolic function in ‘Mathematica’. The relation coefficients for the benders are plotted in three dimensional graph using the developed program.

• Transactions of Materials Processing
• /
• v.21 no.4
• /
• pp.221-227
• /
• 2012

The automotive and electronic industries have seriously considered the use of magnesium alloys because of their excellent properties such as strength to weight ratio, EMI shielding capability, etc. However, it is difficult to form magnesium alloys at room temperature because of the mechanical deformation related to twinning. Hence, magnesium alloys are normally formed at elevated temperatures. In this study, a temperature dependent constitutive model, the C-H/V model, for the magnesium alloy AZ31B sheet is proposed. A hardening law based on nonlinear kinematic and H/V(Hollomon/Voce) hardening model is used to properly characterize the Bauschinger effect and the stabilization of the flow stress. Material parameters were determined from a series of uni-axial cyclic experiments(C-T-C) with the temperature ranging between 150 and $250^{\circ}C$. The developed models are fit to experimental data and a comparison is made.

• Nuclear Engineering and Technology
• /
• v.55 no.11
• /
• pp.4146-4158
• /
• 2023

This paper proposes a combined plastic and creep constitutive model of A533B1 pressure vessel steel to simulate progressive deformation of nuclear pressure vessels under severe accident conditions. To develop the model, recent tensile test data covering a wide range of temperatures (from RT to 1,100 ℃) and strain rates (from 0.001%/s to 1.0%/s) was used. Comparison with experimental data confirms that the proposed combined plastic and creep model can well reflect effects of temperature and strain rate on tensile behaviour up to failure. In the companion paper (Part II), the proposed model will be used to simulate OECD lower head failure (OLHF) test data.

• Smart Structures and Systems
• /
• v.30 no.3
• /
• pp.221-237
• /
• 2022

This paper is concerned with nonlinear modeling and efficient numerical simulation of electrostrictive materials and structures. Two types of such materials are considered: relaxor ferroelectric ceramics and electrostrictive polymers. For ceramics, a geometrically linear formulation is developed, whereas polymers are studied in a geometrically nonlinear regime. In the paper, we focus on constitutive modeling first. For the reversible constitutive response under consideration, we introduce the augmented Helmholtz free energy, which is composed of a purely elastic part, a dielectric part and an augmentation term. For the elastic part, we involve an additive decomposition of the strain tensor into an elastic strain and an electrostrictive eigenstrain, which depends on the polarization of the material. In the geometrically nonlinear case, a corresponding multiplicative decomposition of the deformation gradient tensor replaces the additive strain decomposition used in the geometrically linear formulation. For the dielectric part, we first introduce the internal energy, to which a Legendre transformation is applied to compute the free energy. The augmentation term accounts for the contribution from vacuum to the energy. In our formulation, the augmented free energy depends not only on the strain and the electric field, but also on the polarization and an internal polarization; the latter two are internal variables. With the constitutive framework established, a Finite Element implementation is briefly discussed. We use high-order elements for the discretization of the independent variables, which include also the internal variables and, in case the material is assumed incompressible, the hydrostatic pressure, which is introduced as a Lagrange multiplier. The elements are implemented in the open source code Netgen/NGSolve. Finally, example problems are solved for both, relaxor ferroelectric ceramics and electrostrictive polymers. We focus on thin plate-type structures to show the efficiency of the numerical scheme and its applicability to thin electrostrictive structures.

• Computers and Concrete
• /
• v.25 no.1
• /
• pp.29-36
• /
• 2020

This work presents a numerical modeling procedure to simulate the refractory concrete lining in fluid catalytic cracking units of oil refineries. The model includes the simulation of the anchors that reinforce the contact between the refractory concrete and the steel casing. For this purpose, the constitutive relations of an interface finite element are set to values that represent the homogenized behavior of the anchored interface. The parameters of this constitutive relation can be obtained by experimental tests. The model includes also multi-surface plasticity, in order to represent the behavior of the refractory concrete lining. Since the complexity of real case applications leads to high computational costs, the models presented here were implemented in a high-performance parallelized finite element platform. A case study representing a riser similar to the ones used by the refinery industry demonstrates the potential of the model.

• Geomechanics and Engineering
• /
• v.28 no.1
• /
• pp.89-105
• /
• 2022

The behavior of the soil subgrade is complex and irregular against loads. When modeling, the soil is often replaced by a more straightforward system called a subgrade model. The Winkler method of linear elastic springs is a popular method of soil modeling in which the spring constant shows the modulus of subgrade reaction. In this research, the factors affecting the distribution of the modulus of subgrade reaction of elastoplastic subgrades are examined. For this purpose, critical theories about the modulus of subgrade reaction were examined. A square raft foundation on a sandy soil subgrade with was analyzed at different internal friction angles and Young's modulus values using ABAQUS software. To accurately model the actual soil behavior, the elastic, perfectly plastic constitutive model was applied to investigate a foundation on discrete springs. In order to increase the accuracy of soil modeling, equations have been proposed for the distribution of the subgrade reaction modulus. The constitutive model of the springs is elastic, perfectly plastic. It was observed that the modulus of subgrade reaction under an elastic load decreased when moving from the corner to the center of the foundation. For the ultimate load, the modulus of subgrade reaction increased as it moved from the corner to the center of the foundation.

• International Journal of Automotive Technology
• /
• v.1 no.1
• /
• pp.35-41
• /
• 2000

The crashworthiness of vehicles with finite element methods depends on the geometry modeling and the material properties. The vehicle body structures are generally composed of various members such as frames, stamped panels and deep-drawn parts from sheet metals. In order to ensure the impact characteristics of auto-body structures, the dynamic behavior of sheet metals must be examined to provide the appropriate constitutive relation. In this paper, high strain-rate tensile tests have been carried out with a tension type split Hopkinson bar apparatus specially designed for sheet metals. Experimental results from both static and dynamic tests with the tension split Hopkinson bar apparatus are interpolated to construct the Johnson-Cook and a modified Johnson-Cook equation as the constitutive relation, that should be applied to simulation of the dynamic behavior of auto-body structures. Simulation of auto-body structures has been carried out with an elasto-plastic finite element method with explicit time integration. The stress integration scheme with the plastic predictor-elastic corrector method is adopted in order to accurately keep track of the stress-strain relation for the rate-dependent model accurately. The crashworthiness of the structure with quasi-static constitutive relation is compared to the one with the rate-dependent constitutive model. Numerical simulation has been carried out for frontal frames and a hood of an automobile. Deformed shapes and the Impact energy absorption of the structure are investigated with the variation of the strain rate.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2004.05a
• /
• pp.233-236
• /
• 2004

The major objective of this paper is to clarify the effect of constitutive laws on bulk forming design based on the ideal flow theory. The latter theory is in general applicable for perfectly/plastic materials. However, its kinematics equations constitute a closed-form system, which are valid for any incompressible materials, therefore enabling us to extend design solutions based on the perfectly/plastic constitutive law to more realistic laws with rate sensitive hardening behavior. In the present paper, several constitutive laws commonly accepted for the modeling of cold and hot metal forming processes are considered and the effect of these laws on one particular plane-strain design is demonstrated. The closed form solution obtained describes a non-trivial nonsteady ideal process. The design solutions based on the ideal flow theory are not unique. To achieve the uniqueness, the criterion that the plastic work required to deform the initial shape of a given class of shapes into a prescribed final shape attains its minimum is adopted. Comparison with a non-ideal process is also made.

• Geomechanics and Engineering
• /
• v.5 no.3
• /
• pp.241-261
• /
• 2013

The bearing mechanism of pile during installation and loading process which controls the deformation and distribution of strain and stress in the soil surrounding pile tip is complex and full of much uncertainty. It is pointed out that particle crushing occurs in significant stress concentrated region such as the area surrounding pile tip. The solution to this problem requires the understanding and modeling of the mechanical behavior of granular soil under high pressures. This study aims to investigate the sand behavior around pile tip considering the characteristics of sand crushing. The numerical analysis of model pile loading test under different surcharge pressure with constitutive model for sand crushing is presented. This constitutive model is capable of predicting the dilatancy of soil from negative to positive under low confining pressure and only negative dilatancy under high confining pressure. The predicted relationships between the normalized bearing stress and normalized displacement are agreeable with the experimental results during the entire loading process. It is estimated from numerical results that the vertical stress beneath pile tip is up to 20 MPa which is large enough to cause sand to be crushed. The predicted distribution area of volumetric strain represents that the distributed area shaped wedge for volumetric contraction is beneath pile tip and distributed area for volumetric expansion is near the pile shaft. It is demonstrated that the finite element formulation incorporating a constitutive model for sand with crushing is capable of producing reasonable results for the pile loading problem.

• Korean Journal of Materials Research
• /
• v.25 no.2
• /
• pp.81-89
• /
• 2015

A strain-gradient crystal plasticity constitutive model was developed in order to predict the Hall-Petch behavior of a Ni-base polycrystalline superalloy. The constitutive model involves statistically stored dislocation and geometrically necessary dislocation densities, which were incorporated into the Bailey-Hirsch type flow stress equation with six strength interaction coefficients. A strain-gradient term (called slip-system lattice incompatibility) developed by Acharya was used to calculate the geometrically necessary dislocation density. The description of Kocks-Argon-Ashby type thermally activated strain rate was also used to represent the shear rate of an individual slip system. The constitutive model was implemented in a user material subroutine for crystal plasticity finite element method simulations. The grain size dependence of the flow stress (viz., the Hall-Petch behavior) was predicted for a Ni-base polycrystalline superalloy NIMONIC PE16. Simulation results showed that the present constitutive model fairly reasonably predicts 0.2%-offset yield stresses in a limited range of the grain size.