• 대한기계학회논문집A
• /
• 제24권5호
• /
• pp.1075-1083
• /
• 2000

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. These equations consist of the definition of inelastic strain and the evolution of the state variables introduced to quantify the irreversible processes occurred in the material. With respect to the definition of the inelastic strain, the inelastic constitutive models can be categorized into elastoplastic model, unified viscoplastic model and separated viscoplastic model and the different integration methods have been applied to each category. In the present investigation, the generalized integration method applicable for various types of constitutive equations is developed and implemented 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 line-search technique in combination with Newton method. The strategy to control the time increment for the improvement of the accuracy of the numerical integration is proposed. Several numerical examples are considered to demonstrate the efficiency and applicability of the present method. 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. These equations consist of the definition of inelastic strain and the evolution of the state variables introduced to quantify the irreversible processes occurred in the material. With respect to the definition of the inelastic strain, the inelastic constitutive models can be categorized into elastoplastic model, unified viscoplastic model and separated viscoplastic model and the different integration methods have been applied to each category. In the present investigation, the generalized integration method applicable for various types of constitutive equations is developed and implemented 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 line-search technique in combination with Newton method. The strategy to control the time increment for the improvement of the accuracy of the numerical integration is proposed. Several numerical examples are considered to demonstrate the efficiency and applicability of the present method.

• Structural Engineering and Mechanics
• /
• 제3권2호
• /
• pp.185-192
• /
• 1995

A numerical algorithm for plane stress and shell elasto-plasticity is presented in this paper. The proposed strain decomposition (SD) algorithm is an elastic predictor/plastic corrector algorithm, and in the context of operator splitting, is a return mapping algorithm. However, it differs significantly from other return mapping algorithms in that only the necessary response functions are used without invoking their gradients, and the stress increment is updated only at the end of the time step. This makes the proposed SD algorithm more suitable for materials with complex yield surfaces and will guard against error accumulation during the time step. Comparative analyses of structural systems using the proposed strain decomposition (SD) algorithm and the iterative radial return (IRR) algorithm are presented. The results demonstrate the accuracy and usefulness of the proposed algorithm.

• 대한토목학회논문집
• /
• 제43권4호
• /
• pp.469-476
• /
• 2023

반복하중이나 동적하중에 대한 지반의 소성모델은 지반구조물의 비선형 수치해석에 매우 중요하다. 단일 항복면 모델은 반복하중에 대해 선형적 거동을 보이는 반면, 개발된 탄성영역이 없는 J₂-경계면 소성모델은 동일한 물성치로 효과적으로 지반의 비선형성을 모사할 수 있다. 경계면 내부 항복면의 반경을 0으로 수렴시켜 탄성영역이 사라지도록 수식화하고, 소성경화 계수과 팽창률을 이용하여 소성변형 증분을 정의하였다. 개발된 모델의 응력-변형률 증분식을 제시하고, 쌍곡선 모델에 대한 소성경화 계수를 유도하였다. 삼축압축조건과 얕은기초의 반복하중에 대한 비교해석은 개발된 모델의 안정적인 수렴성, 이론식과의 일치성, 그리고 이력경로을 보여 주었다. 또한, 수정된 쌍곡선함수에 대한 소성경화 계수를 제시하여, 1차원 등가선형모델에 부합하는 모델변수 산정법을 제안하여 지반의 다차원 거동을 모델링할 수 있도록 하였다.

• Earthquakes and Structures
• /
• 제15권2호
• /
• pp.123-132
• /
• 2018

In this paper the strain energy density (SED) model is used to analyze the seismic behavior of suspen-domes and a new criterion is established for judging the seismic failure based on a characteristic point in the SED model. Firstly, a nonlinear time-history response analysis was carried out using the finite-element package ANSYS for typical suspen-domes subjected to different ground motions. The seismic responses including nodal displacements, ratios of yielding members, strain energy density and structural maximum deformation energy were extracted corresponding to the increasing peak ground acceleration (A). Secondly, the SED sum ($I_d$) was calculated which revealed that the $I_d-A$ curve exhibited a relatively large change (called a characteristic point) at a certain value of A with a very small load increment after the structures entered the elastic-plastic state. Thirdly, a SED criterion is proposed to judge the seismic failure load based on the characteristic point. Subsequently, the case study verifies the characteristic point and the proposed SED criterion. Finally, this paper describes the unity and application of the SED criterion. The SED method may open a new way for structural appraisal and the SED criterion might give a unified criterion for predicting the failure loads of various structures subjected to dynamic loads.

• Geomechanics and Engineering
• /
• 제35권6호
• /
• pp.593-601
• /
• 2023

In this paper, an analytical solution of large-strain cylindrical cavity expansion in compaction grouting problem under temperature field is given. Considering the stress increment caused by temperature, the analytical solution of cavity expansion under traditional isothermal conditions is improved by substituting the temperature stress increment into the cavity expansion analysis. Subsequently, combined with the first law of thermodynamics, the energy theory is also introduced into the cylindrical cavity expansion analysis, and the energy dissipation solution of cylindrical cavity expansion is derived. Finally, the validity and reliability of solution are proved by comparing the results of expansion pressure with those in published literatures. The results show that the dimensionless expansion pressure increases with the increase of temperature, and the thermal response increases with the increase of dilation angle. The higher the exothermic temperature of grouting slurry, the greater the plastic deformation energy of the surrounding soil, that is, the greater the influence on the surrounding soil deformation and the surrounding environment. The proposed solution not only enrich the theoretical system of cavity expansion, but also can be used as a theoretical tool for energy geotechnical engineering problems, such as CPT, nuclear waste disposal, energy pile and chemical grouting, etc.

• 한국지반환경공학회 논문집
• /
• 제18권3호
• /
• pp.25-30
• /
• 2017

본 논문에서는 변형률 응답 곡선의 형상을 정량적으로 분석하는 새로운 기법을 소개하였다. 변형률 응답 곡선은 응력 공간에서 단위 응력을 입력하였을 때의 발생한 변형률을 변형률 공간에 작성한 이미지이다. 변형률 응답 곡선의 형상 변화 특성을 퓨리에 기술자를 이용하여 분석하였으며, 탄소성 모델 기반으로 그 의미를 해석하였다. 전산해석결과 응력 경로가 항복면 내부에 머물러 탄성 변형만 존재할 경우, 퓨리에 기술자는 일정하게 그 값을 유지하였으나 응력 경로가 항복면에 도달하여 소성변형이 발생하는 순간 퓨리에 기술자는 변화하였다. 불교란 자연 시료와 재성형 시료에 대해 실시한 응력 경로 시험 결과, 자연 시료에 대해서만 응력 경로가 항복면에 도달하였을 때 명확한 변형률 응답 곡선의 형태 변화가 발견되었으며, 퓨리에 기술자 값도 명확하게 변화하였다. 퓨리에 기술자 중, 비대칭성과 관련 있는 퓨리에 기술자 값이 항복 및 점성토 구조변화의 지표로 활용할 수 있는 것을 확인하였다.

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

• 한국농공학회지
• /
• 제33권2호
• /
• pp.112-119
• /
• 1991

Some experiments were carried out to investigate the effects of freezing and thawing on the strength and strain characteristics of alluvial silty clay under the different temperatures, loading and moisture conditions. The results were as follows; 1. The soil used was proved to be consisted of silty clay with honey-combed structure, and showed higher dilatancy, frost activity and lower stability in natural state. 2. Soil treated with freezing and thawing cycles showed lower compressive strength compared with the non treated, The strength decreased with incement of freezing and thawing cycles. It's shapes of stress-strain curves were flat and did not formulate a peak while the peak strength of higher moisture content soil decreased with the increment of moisture content. It's decrement ratio was most distinctly shown at the first one cycle of freezing and thawing. 3. The cohesion decreased due to freezing and thawing cycles but internal frcition angle was not changed. 4. The liquid limit decreased with increment of freezing and thawing cycles, and became almost constant after three cycles of freezing and thawing. 5. The strength under simple loading at failure mode was appeared to be higher compared with the cyclic loading after freezing and thawing but initial moisture content effect was not observed. 6. Ice lense was not observed within 50% of ice content ratio but observed over 100%. The higher the ice content ratio, the higher the peak strength. As a matter of fact, it seems that an optimum ice content ratio exists for plastic mode and the least compressive strength.

• 한국지반공학회지:지반
• /
• 제12권4호
• /
• pp.21-32
• /
• 1996

흙의 응력-변형률 거동을 더욱 근사적으로 예측하기 위하여,소성 증분이론에 한계상태 토질역학의 개념을 도입한 것이 항복-경계면 소성 모델이다. 이 모텔은 등방 압밀 흙의 거동을 묘사하기 위하여 두개의 타원과 하나의 쌍곡선으로 구성되었다. 따라서 사용된 여러가지 매개변수로 인하여 이 모델은 매우 복잡하다. 그러므로,실제 지반 문제에 이 모델을 적용하기 위하여는 이론의 정확한 이해와 숙련이 요구된다. 본 논문에서는,이 모델에 사용된 여러 매개변수들 중항복-경계면 형상 매개변수 R과 A, 투영 중심 매개변수 C를 변화시켜 그 결과를 수치해석하였다. 최종적으로 단조하중과 주기하중에 대한 각각의 민감도를 분석하였고 각 매개변수 값의 범위를 제안하였다.

• Geomechanics and Engineering
• /
• 제20권5호
• /
• pp.371-384
• /
• 2020

A settlement prediction method based on shear wave velocity measurements and soil nonlinearity was recently developed and verified by means of centrifuge tests. However, the method was only applicable to heavily overconsolidated soil deposits under enlarged yield surfaces. In this study, the settlement evaluation method was refined to consider the stress history of the sublayer, based on an overconsolidation ratio evaluation technique, and thereby incorporate irrecoverable plastic deformation in the settlement calculation. A relationship between the small-strain shear modulus and overconsolidation ratio, which can be determined from laboratory tests, was adopted to describe the stress history of the subsurface. Based on the overconsolidation ratio determined, the value of an empirical coefficient that reflects the effect of plastic deformation over the elastic region is determined by comparing the overconsolidation ratio with the stress increment transmitted by the surface design load. The refined method that incorporate this empirical coefficient was successfully validated by means of centrifuge tests, even under normally consolidated loading conditions.