• 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.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.50 no.10
• /
• pp.691-699
• /
• 2022

Large aperture antennas with long focal lengths in space have important application for telecommunications, Earth observation and science missions. This paper aims to understand the dynamics of deployment of large mesh antennas and to provide a multibody model for determining the driving forces for the design of reflectors and booms. The modular deployable reflector and boom are designed based on the deployment unit cell. A multibody dynamic model is formulated with Kane's equation and simulated using the pseudo upper triangular decomposition (PUTD) method for solving the constrained problem. Based on the multibody dynamic model, the kinetics of the deployment, the motor driving forces, and the structural dynamic deformation are investigated.

• Journal of the Korea Institute of Information and Communication Engineering
• /
• v.26 no.11
• /
• pp.1586-1591
• /
• 2022

Kirchhoff-Love plate (KLP) equation is a well established theory for a description of a deformation of a thin plate under certain outer source. Meanwhile, analysis of a vibrating plate in a frequency domain is important in terms of obtaining the main frequency/eigenfunctions and predicting the vibration of plate. Among various modal analysis methods, dynamic mode decomposition (DMD) is one of the efficient data-driven methods. In this work, we carry out DMD based modal analysis for KLP where thin plate is under effects of sine-type outer force. We first construct discrete time series of KLP solutions based on a finite difference method (FDM). Over 720,000 number of FDM-generated solutions, we select only 500 number of solutions for the DMD implementation. We report the resulting DMD-modes for KLP. Also, we show how DMD can be used to predict KLP solutions in an efficient way.

• Coupled systems mechanics
• /
• v.7 no.5
• /
• pp.583-610
• /
• 2018

This paper discusses the issues associated with modeling frictional contact between solid bodies undergoing large deformations. The most common model for friction on contact interfaces in solid mechanics is the Coulomb friction model, in which two distinct responses are possible: stick and slip. Handling the transition between these two phases computationally has been a source of algorithmic instability, lack of convergence and non-unique solutions, particularly in the presence of large deformations. Most computational models for frictional contact have used penalty or updated Lagrangian approaches to enforce frictional contact conditions. These two approaches, however, present some computational challenges due to conditioning issues in penalty-type implementations and the iterative nature of the updated Lagrangian formulation, which, particularly in large simulations, may lead to relatively slow convergence. Alternatively, a plasticity-inspired implementation of frictional contact has been shown to handle the stick-slip conditions in a local, algorithmically efficient manner that substantially reduces computational cost and successfully avoids the issues of instability and lack of convergence often reported with other methods (Laursen and Simo 1993). The formulation of this approach, however, has been limited to the small deformations realm, a fact that severely limited its application to contact problems where large deformations are expected. In this paper, we present an algorithmically consistent formulation of this method that preserves its key advantages, while extending its application to the realm of large-deformation contact problems. We show that the method produces results similar to the augmented Lagrangian formulation at a reduced computational cost.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2009.05a
• /
• pp.334-336
• /
• 2009

An elastic-plasticity model during the austenitic decomposition was derived and implemented to incorporate the two important deformation behaviors observed during the phase transformations: the volumetric strain and transformation induced plasticity due to the temperature change and phase transformation. To obtain transformed phase volume fractions during cooling, the fourth order Runge-Kutta method was used to solve the Kirkaldy's phase kinetics model which is function of temperature, austenitic grain size and chemical composition. The volumetric strain was calculated by considering the densities of constituent phases, while the transformation induced plasticity was based on the micro-plasticity due to the volume mismatch between soft austenitic phase and other harder phases. The constitutive equations were implemented into the implicit finite element software and a simple boundary value problem was chosen as a model problem to validate the effect of transformation plasticity on the deformation behavior of steel under cooling from high temperature. It was preliminary concluded that the transformation plasticity plays a critical role in relaxing the developed stress during forming and thus reducing the magnitude of springback.

• Structural Engineering and Mechanics
• /
• v.86 no.1
• /
• pp.69-83
• /
• 2023

Analytical solutions to problems are crucial because they provide high-quality comparison data for assessing the accuracy of numerical solutions. Benchmark analytical solutions for the vibrations of cracked functionally graded material (FGM) plates are not available in the literature because of the high level of complexity of such solutions. On the basis of first-order shear deformation plate theory (FSDT), this study proposes analytical series solutions for the vibrations of FGM rectangular plates with side or internal cracks parallel to an edge of the plates by using Fourier cosine series and the domain decomposition technique. The distributions of FGM properties along the thickness direction are assumed to follow a simple power law. The proposed analytical series solutions are validated by performing comprehensive convergence studies on the vibration frequencies of cracked square plates with various crack lengths and under various boundary condition combinations and by performing comparisons with published results based on various plate theories and the theory of three-dimensional elasticity. The results reveal that the proposed solutions are in excellent agreement with literature results obtained using the Ritz method on the basis of FSDT. The paper also presents tabulations of the first six nondimensional frequencies of cracked rectangular Al/Al2O3 FGM plates with various aspect ratios, thickness-to-width ratios, crack lengths, and FGM power law indices under six boundary condition combinations, the tabulated frequencies can serve as benchmark data for assessing the accuracy of numerical approaches based on FSDT.

• Journal of Powder Materials
• /
• v.28 no.3
• /
• pp.246-252
• /
• 2021

The effect of sintering conditions on the austenite stability and strain-induced martensitic transformation of nanocrystalline FeCrC alloy is investigated. Nanocrystalline FeCrC alloys are successfully fabricated by spark plasma sintering with an extremely short densification time to obtain the theoretical density value and prevent grain growth. The nanocrystallite size in the sintered alloys contributes to increased austenite stability. The phase fraction of the FeCrC sintered alloy before and after deformation according to the sintering holding time is measured using X-ray diffraction and electron backscatter diffraction analysis. During compressive deformation, the volume fraction of strain-induced martensite resulting from austenite decomposition is increased. The transformation kinetics of the strain-induced martensite is evaluated using an empirical equation considering the austenite stability factor. The hardness of the S0W and S10W samples increase to 62.4-67.5 and 58.9-63.4 HRC before and after deformation. The hardness results confirmed that the mechanical properties are improved owing to the effects of grain refinement and strain-induced martensitic transformation in the nanocrystalline FeCrC alloy.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.1154-1158
• /
• 2003

Finite element analysis of welding processes, which entail phase evolution, heat transfer and deformations, is considered in this paper. Attention focuses on numerical implementation of the thermo-elastic-plastic constitutive equation proposed by Leblond in consideration of the transformation plasticity. Based upon the multiplicative decomposition of deformation gradient, hyperelastic formulation is employed for efficient numerical integration, and the algorithmic consistent moduli for elastic-plastic deformations including transformation plasticity are obtained in the closed form. The convergence behavior of the present implementation is demonstrated via a couple of numerical examples.

• Steel and Composite Structures
• /
• v.7 no.5
• /
• pp.355-376
• /
• 2007

The purpose of this study is to predict and investigate the time-dependent creep behavior of composite materials. For this, firstly the evaluation method for the modulus of elasticity of whole fiber and matrix is presented from the limited information on fiber volume fraction using the singular value decomposition method. Then, the effects of fiber volume fraction on modulus of elasticity of GFRP are verified. Also, as a creep model, the nonlinear curve fitting method based on the Marquardt algorithm is proposed. Using the existing Findley's power creep model and the proposed creep model, the effect of fiber volume fraction on the nonlinear creep behavior of composite materials is verified. Then, for the time-dependent analysis of a composite material subjected to uniaxial tension and simple shear loadings, a user-provided subroutine UMAT is developed to run within ABAQUS. Finally, the creep behavior of center loaded beam structure is investigated using the Hermitian beam elements with shear deformation effect and with time-dependent elastic and shear moduli.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2010.04a
• /
• pp.542-545
• /
• 2010

재료의 마이크로 스케일 해석에서 결정의 geometrically necessary dislocation (GND) 효과에 의한 소성구배(plastic gradient)를 고려하는 것은 재료의 소성 거동을 분석하는데 영향을 미친다. 본 연구에서는 먼거리(long range)에서 전위(dislocation)의 영향을 고려하는 GND의 효과를 적용하여 소성 구배의 영향을 받는 다결정(polycrystal) 고체의 거동을 유한요소해석을 이용하여 살펴보았다. 재료의 거동을 분석하기 위해 탄성(elastic)과 소성(plastic) 변형에 먼 거리 변형률(long range strain)을 고려한 항(term)이 포함된 변형 구배(deformation gradient)의 multiplicative decomposition 모델을 사용하였다. 먼 거리 변형률에 의한 영향을 고려하기 위해 구배 경화 계수(gradient hardness coefficient)와 먼 거리 변형률 길이에 대한 재료변수(parameter)가 사용되었다. 각각의 계수들이 다결정 고체의 거동에 미치는 영향을 확인하기 위해 두 변수의 적용에 따른 다결정 고체의 거동을 분석하였다. 다결정 재료의 GND 효과에 의한 소성 구배 효과를 고려해서, 고려하지 않은 경우와 비교하여 발생하는 경화(hardening)의 차이를 분석함으로서 GND에 의한 다결정 고체 거동의 영향을 확인하였다.