• 대한기계학회논문집A
• /
• 제29권2호
• /
• pp.253-261
• /
• 2005

In decades, a substantial body of work on a unified viscoplastic model which considers the mechanism of plastic deformation and creep deformation has developed. The systematic scheme for numerical analysis of unified model is necessary because the dominant failure mechanism is the defect growth and coalescence in materials. In the present study, the unified viscoplastic model for materials with defects suggested by Suquet and Michel was employed for numerical analysis. The constitutive equations are integrated based on the generalized mid-point rule and implemented into a finite element program (ABAQUS) by means of user-defined subroutine (UMAT). To evaluate the validity of the developed UMAT code and the assessment of the adopted viscoplastic model, the results obtained from the UMAT code was compared with the numerical reference solution and experimental data. The unit cell analysis also has been investigated to study the effect of strain rate, temperature, stress triaxiality and initial defect volume fraction on the growth and coalescence of the defect.

• Geomechanics and Engineering
• /
• 제17권2호
• /
• pp.119-131
• /
• 2019

Time dependent creep settlements are one of the most important causes of material deteriorations for the huge water structures such as concrete faced rockfill dams (CFRDs). For this reason, performing creep analyses of CFRDs is vital important for monitoring and evaluating of the future and safety of such dams. In this study, it is observed how changes viscoplastic behaviour of a CFR dam depending the time. Ilısu dam that is the longest concrete faced rockfill dam (1775 m) in the world is selected for the three dimensional (3D) analyses. 3D finite difference model of Ilısu dam is modelled using FLAC3D software based on the finite difference method. Two different special creep material models are considered in the numerical analyses. Wipp-creep viscoplastic material model and burger-creep viscoplastic material model were rarely used for the creep analyses of CFRDs in the last are taken into account for the concrete slab and rockfill materials-foundation, respectively. Moreover, interface elements are defined between the concrete slab-rockfill materials and rockfill materials-foundation to provide interaction condition for 3D model. Firstly, dam and foundation are collapsed under its self-weight and static behaviour of the dam is evaluated for the empty reservoir conditions. Then, reservoir water is modelled considering maximum water level of the dam and time-dependent creep analyses are performed for maximum reservoir condition. In this paper, maximum principal stresses, vertical-horizontal displacements and pore pressures that may occur on the dam body surface during 30 years (from 2017 to 2047) are evaluated in detail. According to numerical analyses, empty and maximum reservoir conditions of Ilısu dam are compared with each other in detail. 4 various nodal points are selected under the concrete slab to better seen viscoplastic behaviour changes of the dam and viscoplastic behaviour differences of these points during 30 years are graphically presented. It is clearly seen that horizontal-vertical displacements and principal stresses for maximum reservoir condition are more than the empty reservoir condition of the dam and significant pore pressures are observed during 30 years for maximum reservoir condition. In addition, horizontal-vertical displacements, principal stresses and pore pressures for 4 nodal points obviously increased until a certain time and changes decreased after this time.

• Computers and Concrete
• /
• 제8권4호
• /
• pp.371-388
• /
• 2011

In this paper, time-continuous constitutive equations for strain rate-dependent materials are presented first, among which those for the overstress and the consistency viscoplastic models are considered. By allowing the stress states to be outside the yield surface, the overstress viscoplastic model directly defines the flow rule for viscoplastic strain rate. In comparison, a rate-dependent yield surface is defined in the consistency viscoplastic model, so that the standard Kuhn-Tucker loading/unloading condition still remains true for rate-dependent plasticity. Based on the formulation of the consistency viscoplasticity, a computational elasto-viscoplastic constitutive model is proposed for the short fiber-reinforced fresh cementitious paste for extrusion purpose. The proposed constitutive model adopts the von-Mises yield criterion, the associated flow rule and nonlinear strain rate-hardening law. It is found that the predicted flow stresses of the extrudable fresh cementitious paste agree well with experimental results. The rate-form constitutive equations are then integrated into an incremental formulation, which is implemented into a numerical framework based on ANSYS/LS-DYNA finite element code. Then, a series of upsetting and ram extrusion processes are simulated. It is found that the predicted forming load-time data are in good agreement with experimental results, suggesting that the proposed constitutive model could describe the elasto-viscoplastic behavior of the short fiber-reinforced extrudable fresh cementitious paste.

• 한국지반공학회논문집
• /
• 제18권4호
• /
• pp.55-64
• /
• 2002

본 논문에서는 점탄-점소성 구성모델을 다층지반에서의 대형지진 발생시의 조건에 적용하기 위해 일본 고베 포트아일랜드에서 발생한 1995 Hyogoken Nanbu 지진에 대한 지진응답 해석을 수행하였다. 지진응답해석 결과 점성토의탄-점소성 모델과 점탄-점소성 모델로 계산된 가속도 기록은 포트아일랜드에서 계측된 가속도 기록과 거의 일치함을 알 수 있었으며, 점성토 지반 부근에서 점탄-점소성 모델과 탄-점소성 모델은 미세하게 다른 거동 특성을 나타내어 점소성 모델의 타당성을 확인하였다. 따라서 동적 점탄-점소성 구성모델은 대변형률 영역에서 점성토의 소성변형을 유발하는 대형 지진 등의 발생시 점성토의 증폭 및 감쇠특성의 파악을 위해 적용가능한 모델임이 입증되었다.

• Journal of Theoretical and Applied Mechanics
• /
• 제3권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.

• 소성∙가공
• /
• 제4권3호
• /
• pp.204-213
• /
• 1995

A finite element program is developed to analyze the flow phenomena in SMC compression molding as a viscoplastic model. The calculation of temperature distribution is also carried out by uncoupling the thermal analysis from the flow analysis. SMC molding processes with a flat plate substructure and the one with a T-shaped rib are considered in numerical simulation. The numerical results provide deformed shapes, temperature distribution in a SMC charge, and the forming load. The simulation of compression molding of a flat plate with a T-shaped rib requires a remeshing technique for the whole process.

• Computers and Concrete
• /
• 제8권3호
• /
• pp.243-278
• /
• 2011

The paper presents simulations of the Brazilian test using two numerical models. Both models are regularized in order to obtain results independent of discretization. The first one, called gradient damage, is refined by additional averaging equation which contains gradient terms and an internal length scale as localization limiter. In the second one, called viscoplastic consistency model, the yield function depends on the viscoplastic strain rate. In this model regularization properties are governed by the assumed strain rate. The two models are implemented in the FEAP finite element package and compared in this paper. Parameter studies of the split test are performed in order to point out the features of each model.

• 한국농공학회지
• /
• 제44권4호
• /
• pp.114-122
• /
• 2002

Up to now, many constitutive models for clay have been proposed and studied based on the elasto-plastic or elasto-viscoplastic theory and it has been recognized that the effect of time on the loading process is a salient feature. In the present study, cyclic behavior characteristics of clay was studied with a viscoelastic-viscoplastic constitutive model for clay based on the non-linear kinematic hardening rule. In order to examine the behavior of clay several cyclic untrained triaxial tests and also their numerical simulations were performed. As results of that, it was found that the proposed model can well describe cyclic behaviors of clay such as frequency dependent characteristics, and have the high feasibility of numerical simulation for dynamic analysis.

• Transactions on Electrical and Electronic Materials
• /
• 제5권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.

• 대한기계학회논문집A
• /
• 제27권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.