• 대한기계학회논문집
• /
• 제15권1호
• /
• pp.190-200
• /
• 1991

본 연구에서는 고온에 노출되는 열-점소성 거동의 해석을 위해 소성유동, 크 리프, 응력풀림(stress relaxation)등의 거동을 동시에 다룰 수 있는 통일구성방정식 모델에 대해 논하고 적절한 모델을 선정한다. 이 모델은 미소변형이론에 근거한 것 이므로 구조물의 거동을 소변형률(small strain)과 소회전(small rotation)의 범위내 로 가정하여 해석한다.선정된 모델에 대해서 시간변화율 형태의 방정식으로부터 유 한요소법을 통한 수치화와 사용된 구성방정식을 효율적으로 처리할 수 있는 수치해석 법상의 알고리듬을 제안한다. 제안된 알고리즘을 사용하여 유한요소법 전산코드를 적상하고, 작성된 코드를 이용하여 고온에서 하중을 받는 단순보와 국부적으로 심한 가열을 받는 구조물에 적용하여 고전적인 구성방정식으로 복합적인 해석이 어려웠던 열-점소성 거동을 효과적으로 해석할 수 있음을 보인다. 본 논문은 응력해석에 주안 점을 두었으므로 열해석에 관한 상세한 논의는 가급적 생략하기로 한다.

• Geomechanics and Engineering
• /
• 제37권2호
• /
• pp.179-187
• /
• 2024

In this paper, a unified time-dependent constitutive model of Darcy flow and non-Darcy flow is proposed. The influencing factors of flow velocity are discussed, which demonstrates that permeability coefficient is the most significant factor. Based on this, the dynamic evolution characteristics of water inrush during tunneling through fault fracture zone is analyzed under the constant permeability coefficient condition (CPCC). It indicates that the curves of flow velocity and hydrostatic pressure can be divided into typical three stages: approximate high-velocity zone inside the fault fracture zone, velocity-rising zone near the tunnel excavation face and attenuation-low velocity zone in the tunnel. Furthermore, given the variation of permeability coefficient of the fault fracture zone with depth and time, the dynamic evolution of water flow in the fault fracture zone under the variable permeability coefficient condition (VPCC) is also studied. The results show that the time-related factor (α) affects the dynamic evolution distribution of flow velocity with time, the depth-related factor (A) is the key factor to the dynamic evolution of hydrostatic pressure.

• 한국농공학회지
• /
• 제30권4호
• /
• pp.134-153
• /
• 1988

This paper was carried out to investigate the existence of a unique stress- strain behavior by obtaining some factors influencing the time dependent stress- strain behavior of clay. The results obtained from this study were summarized as follows ; 1. The relationship between stress ratro and strain in normally consolidated clay was in- dependent on pre-shear consolidation pressure. Therefore, shear strain could be expressed as a function with stress ratio. 2. The constitutive equation of shear strain on Modified Carn Clay Model coincided better with the observed value than Cam Clay Model. 3. The relationships between deviator stress and shear strain, between pore water pressure and shear strain were unified by the mean equivalent pressure. 4. The shear strain contour in norrnally consolidated clay was increased linearly through origin, but that in overconsolidated clay was not in accordance with the result of the former. 5. Because the effective stress path of normally consolidated clay was unified by the mean equivalent pressure, state boundary surface in (e,p,q) space was transformed into two dimensional surface. But it was considered to be suitable that the unified stress- strain in overconsolidated clay be expressed by a function with overconsolidation ratio. 6. The deviator for constant strain was increased linearly with increment of strain rate ($\varepsilon$) on semi-log scale, but pore water pressure was decreased. 7. The behavior of stress relaxation was transformed from linear to curvilinear with inc - rement of strain rate before stress relaxation test, and pore water pressure was increased in total range. 8. The strain of creep was increased linearly with increment of time on semi-log scale. The greater the strain rate before creep test became, the greater the increment of strain of creep became. And the pore water pressure during creep test was increased generally with increment of time on semi-log scale.

• 대한기계학회논문집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.

• Geomechanics and Engineering
• /
• 제25권3호
• /
• pp.195-205
• /
• 2021

Time effect on the deformation and strength characteristics of geogrid reinforced sand retaining wall has become an important issue in geotechnical and transportation engineering. Three physical model tests on geogrid reinforced sand retaining walls performed under various loading conditions were simulated to study their rate-dependent behaviors, using the presented nonlinear finite element method (FEM) analysis procedure. This FEM was based on the dynamic relaxation method and return mapping scheme, in which the combined effects of the rate-dependent behaviors of both the backfill soil and the geosynthetic reinforcement have been included. The rate-dependent behaviors of sands and geogrids should be attributed to the viscous property of materials, which can be described by the unified three-component elasto-viscoplastic constitutive model. By comparing the FEM simulations and the test results, it can be found that the present FEM was able to be successfully extended to the boundary value problems of geosynthetic reinforced soil retaining walls. The deformation and strength characteristics of the geogrid reinforced sand retaining walls can be well reproduced. Loading rate effect, the trends of jump in footing pressure upon the step-changes in the loading rate, occurred not only on sands and geogrids but also on geogrid reinforced sands retaining walls. The lateral earth pressure distributions against the back of retaining wall, the local tensile force in the geogrid arranged in the retaining wall and the local stresses beneath the footing under various loading conditions can also be predicted well in the FEM simulations.