• Structural Engineering and Mechanics
• /
• v.6 no.8
• /
• pp.857-870
• /
• 1998

This paper deals with a special class of inverse problems in discrete structural plasticity involving the identification of elastic limits and hardening moduli on the basis of information on displacements. The governing equations lead naturally to a special and challenging optimization problem known as a Mathematical Program with Equilibrium Constraints (MPEC), a key feature of which is the orthogonality of two sign-constrained vectors or so-called "complementarity" condition. We investigate numerically the application of two simple algorithms, both based on the use of the general purpose nonlinear programming code CONOPT accessed via the GAMS modeling language, for solving the suitably reformulated problem. Application is illustrated by means of two numerical examples.

• Fibers and Polymers
• /
• v.7 no.3
• /
• pp.276-285
• /
• 2006

Development and numerical implementation for an elastoplastic constitutive model for anisotropic and asymmetric materials are presented in this paper. The Coulomb-Mohr yield criterion was modified to consider both the anisotropic and asymmetric properties. The modified yield criterion is an isotropic function of the principal values of a symmetric matrix which is linearly transformed from the Cauchy stress space. In addition to the constitutive equation, the numerical treatment for the singularity in the vertex region of yield surface and stress integration algorithm based on elastoplasticity were presented. In order to assess the accuracy of numerical algorithm, isoerror maps were considered. Also, extension of a strip with a circular hole was simulated and results compared with those obtained using the (smooth) Mises yield criterion to validate stress output for a complex stress state.

• Communications of the Korean Mathematical Society
• /
• v.20 no.3
• /
• pp.589-602
• /
• 2005

We present the numerical algorithm for the model for high-strain rate deformation in hyperelastic-viscoplastic materials based on a fully conservative Eulerian formulation by Plohr and Sharp. We use a hyperelastic equation of state and the modified Steinberg and Lund's rate dependent plasticity model for plasticity. A two-dimensional approximate Riemann solver is constructed in an unsplit manner to resolve the complex wave structure and combined with the second order TVD flux. Numerical results are also presented.

• Journal of the korean Society of Automotive Engineers
• /
• v.13 no.3
• /
• pp.58-67
• /
• 1991

A theory of continuum damage mechanics based on the theory of materials of type N was developed and its nonlinear finite element approximation and numerical simulation was carried out. To solve the finite elastoplasticity problems, reasonable kinematics of large deformed solids was introduced and constitutive relations based on the theory of materials of type-N were derived. These highly nonlinear equations were reduced to the incremental weak formulation and approximated by the theory of nonlinear finite element method. Two types of problems, compression moulding problem and pure bending problem, were solved for aluminum 2024.

• Structural Engineering and Mechanics
• /
• v.30 no.6
• /
• pp.665-678
• /
• 2008

The beams components subjected to the loading such as axial, bending and cyclic thermal loads were studied in this research. The used constitutive equations are those of elasto-plasticity coupled to ductile and/or creep damage. The nonlinear kinematic hardening behavior was considered in elastoplasticity modeling. The unified damage law proposed for ductile failure and fatigue by the author of Sermage et al. (2000) and Kachanov's creep damage model applied to cyclic creep and low cycle fatigue of beams. Based on the results of the analysis, the shakedown limit loads were determined through the calculation of the residual strains developed in the beam analysis. The iterative technique determines the shakedown limit load in an iterative manner by performing a series of full coupled elastic-plastic and continuum damage cyclic loading modeling. The maximum load carrying capacity of the beam can withstand, were determined and imposed on the Bree's interaction diagram. Comparison between the shakedown diagrams generated by or without creep and/or ductile damage for the loading patterns was presented.

• Proceedings of the Korean Geotechical Society Conference
• /
• 1991.10a
• /
• pp.166-178
• /
• 1991

Using FLAC, which is an explicit finite difference code written for analysis of problems in geotechnical engineering, a particular example 2 in the Korean Geothechnical Society News has been analysed. The elastoplasticity formulation in FLAC assumes an elastic, perfectly plastic solid in plane strain which conforms to a Mohr-Coulomb yield condition. During tunnel excavations by stages, stresses and diaplacements in region around the tunnel varies according to distance from the face of tunnel and installation of tunnel supports, and soon. In this analysis, the three dimensional support effect of the rock mass during the process of excavation is simulated by using the stress distribution method, and varying the material constants of shotcrete in each stage also considered. The maximum convergency is occured at the crown of the tunnel and estimated to be about 12mm.

• Structural Engineering and Mechanics
• /
• v.30 no.1
• /
• pp.119-129
• /
• 2008

A small strain and elastoplastic formulation of Polygonal Element Method (PEM) is developed for efficient analysis of elastoplastic solids. In this work, the polygonal elements are constructed based on traditional triangular finite meshes. The construction method of polygonal mesh can directly utilize the sophisticated triangularization algorithm and reduce the difficulty in generating polygonal elements. The Wachspress rational finite element basis function is used to construct the approximations of polygonal elements. The incremental variational form and a von Mises type model are used for non-linear elastoplastic analysis. Several small strain elastoplastic numerical examples are presented to verify the advantages and the accuracy of the numerical formulation.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2000.11a
• /
• pp.327-334
• /
• 2000

응력-변형률 관계의 모델링에 있어서 creep, stress relaxation, strain rate effect 등의 묘사는 중요한 지반거동중의 하나인 시간 의존적 거동의 simulation은 있어서 대단히 중요한 요소라 할 수 있다. 특히 지반은 변형률 속도에 대하여 때로는 매우 다른 거동 특성을 보이기 때문에 지반의 모델링에 있어서 변형율 속도를 고려한 구성방정식의 제시는 큰 비중을 차지한다 하겠다. 본 연구에서는 변형율에 따라 변화하는 지반의 거동특성을 보다 현실에 가깝게 묘사하기 위한 구성모델을 제시하였다. 이를 위하여 Bounding Surface Model의 점탄소성 부분을 Perzyna(1966)와 Adachi and Oka(1982)의 구성방정식 이론을 이용하여 발전시켰다. 제안된 구성모델은 기존의 모델에 비하여 다양한 변형율 속도에 적용할 수 있는 모델 정수를 비교적 간단히 결정할 수 있다는 장점이 있으며, 변형율 속도의 영향뿐 아니라 creep, stress relaxation등의 현상도 잘 simulation 할 수 있다. 본 모델은 후에 엄격히 실시되는 실내시험을 통하여 검증될 예정이다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.14 no.1
• /
• pp.11-19
• /
• 2001

기능경사재(FGM)는 구성 물질의 체적분율(volume fraction)이 복합재 전체에 걸쳐 연속적 그리고 기능적으로 분포되어 있어, 기존의 이종물질 접합식(bi-material-type) 복합재보다 현저히 우수한 열기계적 특성을 가진다. 하지만, 기능경사 내열복합재의 열-탄소성 거동은 체적분율의 분포형태와 경사층이 차지하는 상대두께비에 따라 절대적으로 좌우된다. 본 연구는 기능경사 내열복합재의 열-탄소성 특성의 이들 두 설계인자에 대한 파라메트릭 FEM해석을 다룬 것이다. 열-탄소성 이론과 유한요소 근사화에 따라 연구용 2차원 FEM 프로그램을 개발하고, 대표적인 3층 구조의 2차원 기능경사 내열복합재의 열-탄소성 특성을 설계변수의 다양한 조합에 따라 분석하였다.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2001.03a
• /
• pp.207-214
• /
• 2001

응력-변형률 관계의 모델링에 있어서 creep, stress relaxation, strain rate effect 등의 묘사는 중요한 지반거동중의 하나인 시간 의존적 거동에 대한 simulation은 있어서 대단히 중요한 요소라 할 수 있다. 특히 지반은 변형률 속도에 대하여 때로는 매우 다른 거동 특성을 보이기 때문에 지반의 모델링에 있어서 변형율 속도를 고려한 구성방정식의 제시는 큰 비중을 차지한다 하겠다. 본 연구에서는 변형율에 따라 변화하는 지반의 거동특성을 보다 현실에 가갈게 묘사하기 위한 시간 의존적 구성모델을 제시하는데 있다. Bounding Surface Model의 Stress Invariant 부분을 Perzyna(1966)와 Adachi and Oka(1982)의 변형율 속도 의존적인 구성관계 이론을 이용하여 발전시켰다 제안된 구성모델은 다양한 변형율 속도에 적용에 있어서 기존의 방식보다 간단히 모델 정수들을 결정 할 수 있다. 지반거동의 수치적인 해석을 위하여 기존의 Bounding Surface Model에 사용되었던 Program Code를 발전 시켜 사용하였으며, 엄격히 시행된 실내시험의 결과와 비교/검증하였다.