• 한국지구물리탐사학회:학술대회논문집
• /
• 한국지구물리탐사학회 2003년도 Proceedings of the international symposium on the fusion technology
• /
• pp.184-191
• /
• 2003

In this paper, a mathematical framework based on a homogenisation technique to simulate 'umbrella arch reinforcement system' (UARS) and its implementation into a 3D Finite Element program that can consider stage construction situations are presented. The constitutive model developed allows considering the main design parameters of the problem and only requires geometrical and mechanical properties of the constituents. Additionally, the use of a homogenisation approach implies that the generation of the Finite Element mesh can be easily produced and that re-meshing is not required as basic geometrical parameters such as the orientation of the pipes are changed. The model developed is used to simulate tunnelling with the UARS. From the analyses, the effects of the main design parameters on the elastic and the elastoplastic analyses are considered.

• 한국지반공학회:학술대회논문집
• /
• 한국지반공학회 2000년도 가을 학술발표회 논문집
• /
• pp.327-334
• /
• 2000

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

• 한국콘크리트학회:학술대회논문집
• /
• 한국콘크리트학회 2006년도 춘계학술발표회 논문집(I)
• /
• pp.386-389
• /
• 2006

Prestressed composite girder bridges with concrete infilled steel tube at negative flexural moment region takes the advantages of enhancing local buckling and flexural resistances resulting from the lateral confining effect of concrete due to the interactive reaction in the interface layer of steel tube and concrete. The interface behavior in concrete infilled tube of the test composite girder is analyzed by 8-node zero thickness interface finite element combined with 3-D. elastoplastic concrete constitutive model and 3-D. elastoplastic Mindlin shell element. The interface effects between infillled concrete and steel tube are investigated through the comparision of the experimental and numerical results.

• Geomechanics and Engineering
• /
• 제9권5호
• /
• pp.669-685
• /
• 2015

A three-dimensional elastoplastic constitutive model, also known as a UH model (Yao et al. 2009), was developed to describe the stress-strain relationship for normally consolidated and over-consolidated soils. In this paper, an acoustic tensor and discriminator of bifurcation for the UH model are derived for the strain localization of saturated clays under undrained and fully and partially drained conditions. Analytical analysis is performed to illustrate the points of bifurcation for the UH model with different three-dimensional stress paths. Numerical analyses of cubic specimens for the bifurcation of saturated clays under undrained and fully and partially drained conditions are conducted using ABAQUS with the UH model. Analytical and numerical analyses show the similar bifurcation behaviour of overconsolidated clays in three-dimensional stress states and various drainage conditions. The results of analytical and numerical analyses show that (1) the occurrence of bifurcation is dependent on the stress path and drainage condition; and (2) bifurcation can appear in either a strain-hardening or strain-softening regime.

• 대한기계학회논문집A
• /
• 제28권8호
• /
• pp.1140-1151
• /
• 2004

This paper reports on the finite elements analysis for die compaction process of cemented carbide tool parts. Experimental data were obtained under die compaction and triaxial compression with various loading conditions. The elastoplastic constitutive equations based on the yield function of Shima and Oyane were implemented into an explicit finite element program (ABAQUS/Explicit) and implicit finite element program (PMsolver/Compaction-3D) to simulate compaction response of cemented carbide powder during die compaction. For simulation of die compaction, the material parameters for Shima and Oyane model were obtained by uniaxial die compaction test. Explicit finite element results were compared with implicit results for cemented carbide powder.

• Smart Structures and Systems
• /
• 제30권3호
• /
• pp.245-253
• /
• 2022

A two-layer beam consisting of an elastoplastic layer and a functional layer made of shape memory alloy (SMA) TiNi is considered. Constitutive relations for SMA are set by a microstructural model capable to calculate strain increment produced by arbitrary increments of stress and temperature. This model exploits the approximation of small strains. The equations to calculate the variations of the strain and the internal variables are based on the experimentally registered temperature kinetics of the martensitic transformations with an account of the crystallographic features of the transformation and the laws of equilibrium thermodynamics. Stress and phase distributions over the beam height are calculated by steps, by solving on each step the boundary-value problem for given increments of the bending moment (or curvature) and the tensile force (or relative elongation). Simplifying Bernoulli's hypotheses are applied. The temperature is considered homogeneous. The first stage of the numerical experiment is modeling of preliminary deformation of the beam by bending or stretching at a temperature corresponding to the martensitic state of the SMA layer. The second stage simulates heating and subsequent cooling across the temperature interval of the martensitic transformation. The curvature variation depends both on the total thickness of the beam and on the ratio of the layer's thicknesses.

• Structural Engineering and Mechanics
• /
• 제68권3호
• /
• pp.299-312
• /
• 2018

In this paper, an extended finite element method is proposed to analyze both geometric and material non-linear behavior of general Functionally Graded Material (FGM) plate-shell type structures. A user defined subroutine (UMAT) is developed and implemented in Abaqus/Standard to study the elastoplastic behavior of the ceramic particle-reinforced metal-matrix FGM plates-shells. The standard quadrilateral 4-nodes shell element with three rotational and three translational degrees of freedom per node, S4, is extended in the present study, to deal with elasto-plastic analysis of geometrically non-linear FGM plate-shell structures. The elastoplastic material properties are assumed to vary smoothly through the thickness of the plate-shell type structures. The nonlinear approach is based on Mori-Tanaka model to underline micromechanics and locally determine the effective FGM properties and self-consistent method of Suquet for the homogenization of the stress-field. The elasto-plastic behavior of the ceramic/metal FGM is assumed to follow Ludwik hardening law. An incremental formulation of the elasto-plastic constitutive relation is developed to predict the tangent operator. In order to to highlight the effectiveness and the accuracy of the present finite element procedure, numerical examples of geometrically non-linear elastoplastic functionally graded plates and shells are presented. The effects of the geometrical parameters and the volume fraction index on nonlinear responses are performed.

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

• Computers and Concrete
• /
• 제10권3호
• /
• pp.307-330
• /
• 2012

This work intends to contribute for the improvement of the procedure suggested by Brazilian Technical Code that takes into account the cracked concrete stiffness in the estimative of the displacement of reinforced concrete beams submitted to service loads. A damage constitutive model accounting for induced anisotropy, plastic deformations and bimodular elastic response is used in order to simulate the concrete behaviour, while an elastoplastic behaviour is admitted for the reinforcement. The constitutive models were implemented in a program for bars structures analysis with layered finite elements. Initially, the damage model is briefly presented as well as the parametric identification of the materials that have been used in the reinforced concrete beams. After that, beams with different geometries and reinforcement area are analyzed and a statistical method (ANOVA) is employed in order to identify the main variables in the problem. Soon after, the same procedure is used with another resistance of concrete, where the compression strength is changed. The numerical responses are compared with the ones obtained by Brazilian Technical Code and experimental tests in order to validate the use of the damage model. Finally, some remarks are discussed based on responses presented in this work.

• Geomechanics and Engineering
• /
• 제36권4호
• /
• pp.355-366
• /
• 2024

The original elastoplastic Hardening Soil model is formulated actually partly under hexagonal pyramidal Mohr-Coulomb failure criterion, and can be only used in specific stress paths. It must be completely generalized under Mohr-Coulomb criterion before its usage in engineering practice. A set of generalized constitutive equations under this criterion, including shear and volumetric yield surfaces and hardening laws, is proposed for Hardening Soil model in principal stress space. On the other hand, a Mohr-Coulumb type yield surface in principal stress space comprises six corners and an apex that make singularity for the normal integration approach of constitutive equations. With respect to the isotropic nature of the material, a technique for processing these singularities by means of Koiter's rule, along with a transforming approach between both stress spaces for both stress tensor and consistent stiffness matrix based on spectral decomposition method, is introduced to provide such an approach for developing generalized Hardening Soil model in finite element analysis code ABAQUS. The implemented model is verified in comparison with the results after the original simulations of oedometer and triaxial tests by means of this model, for volumetric and shear hardenings respectively. Results from the simulation of oedometer test show similar shape of primary loading curve to the original one, while maximum vertical strain is a little overestimated for about 0.5% probably due to the selection of relationships for cap parameters. In simulation of triaxial test, the stress-strain and dilation curves are both in very good agreement with the original curves as well as test data.