• Tunnel and Underground Space
• /
• v.10 no.3
• /
• pp.355-365
• /
• 2000

As large underground projects such as radioactive waste disposal, hot water and heat storage, and geothermal energy become influential, the study, which consider all aspects of thermics, hydraulics and mechanics would be needed. Thermo-Hydro-Mechanical coupling analysis is one of the most complex numerical technique because it should be implemented with the combined three governing equations to analyze the behavior of rock mass. In this study, finite element code, which is based on Biot's consolidation theory, was developed to analyze the thermo-hydro-mechanical coupling in continuum rock mass. To verify the implemented program, one-dimensional consolidation model under the isothermal and non-isothermal conditions was analyzed and was compared with the analytic solution. The parametric study on two-dimensional consolidation was also performed and the effects of several factors such as poisson's ratio and hydraulic anisotropy on rock mass behavior were investigated. In the future, this program would be revised to be used for analysis of general discontinuous media with incorporating discrete joint model.

• Structural Engineering and Mechanics
• /
• v.48 no.1
• /
• pp.41-56
• /
• 2013

Initiation and growth of delamination is a great concern of designers of composite structures. Interface elements with de-cohesive constitutive law in the content of continuum damage mechanics can be used to predict initiation and growth of delamination in single and mixed mode conditions. In this paper, an interface element based on the cohesive zone method has been developed to simulate delaminatoin growth of post-buckled laminate under fatigue loading. The model was programmed as the user element and user material by the "User Programmable Features" in ANSYS finite element software. The interface element is a three-dimensional 20 node brick with small thickness. Because of mixed-mode condition of stress field at the delamination-front of post-buckled laminates, a mixed-mode bilinear constitutive law has been used as user material in this model. The constitutive law of interface element has been verified by modelling of a single element. A composite laminate with initial delamination under quasi-static compressive Loading available from literature has been remodeled with the present approach. Moreover, it will be shown that, the closer the delamination to the free surface of laminate, the slower the delamination growth under compressive fatigue loading. The effects of laminate configuration on delamination growth are also investigated.

• Computers and Concrete
• /
• v.7 no.6
• /
• pp.483-496
• /
• 2010

This work deals with the determination of crack openings in 2D reinforced concrete structures using the Finite Element Method with a smeared rotating crack model or an embedded crack model. In the smeared crack model, the strong discontinuity associated with the crack is spread throughout the finite element. As is well known, the continuity of the displacement field assumed for these models is incompatible with the actual discontinuity. However, this type of model has been used extensively due to the relative computational simplicity it provides by treating cracks in a continuum framework, as well as the reportedly good predictions of reinforced concrete members' structural behavior. On the other hand, by enriching the displacement field within each finite element crossed by the crack path, the embedded crack model is able to describe the effects of actual discontinuities (cracks). This paper presents a comparative study of the abilities of these 2D models in predicting the mechanical behavior of reinforced concrete structures. Structural responses are compared with experimental results from the literature, including crack patterns, crack openings and rebar stresses predicted by both models.

• Transactions of Materials Processing
• /
• v.31 no.5
• /
• pp.261-266
• /
• 2022

To reduce the weight of the vehicle, the application of the high strength steel sheets to chassis parts is increased. High forming load is induced during the shearing process of steel chassis parts made of high strength steel, and the possibility of an eccentric load is increased depending on the product seating condition on the die, which decreases the stability and lifespan of the die. In this paper, a three-dimensional finite element analysis with the continuum element was conducted using the damage theory for the cam-trimming process of the front lower arm. The structural analysis of the trimming die was performed with the forming load result obtained from the analysis, and the amount of deflection and the stress distribution of the die during the shearing process were evaluated for the confirmation of the tool stability. The shape of the weak region of the die was modified according to structural analysis and then the stability was confirmed with the finite element analysis. The analysis result showed that the possibility of tool failure during cam-trimming process was remarkably reduced, and the reliability of the proposed modified design was validated.

• Transactions of Materials Processing
• /
• v.26 no.5
• /
• pp.286-292
• /
• 2017

Magnesium alloy shows strong anisotropy and asymmetric behavior in tension and compression curve, especially at room temperature. These characteristics limit the application of finite element method (FEM) which is based on conventional continuum mechanics. To accurately predict the material behavior of magnesium alloy at microstructural level, a methodology of fully coupled multiscale simulation is presented and a crystal plasticity model as a constitutive equation in the simulation of metal forming process is introduced in this study. The existing constitutive equation for rigid plastic FEM is modified to accommodate deviatoric stress component and its derivatives with respect to strain rate components. Viscoplastic self-consistent (VPSC) polycrystal model was selected as a constitutive model because it was regarded as the most robust model compared to Taylor model or Sachs model. Stiffness matrix and load vector were derived based on the new approach and implemented into $DEFORM^{TM}-3D$ via a user subroutine handling stiffness matrix at an elemental level. The application to extrusion and rolling process of pure magnesium is presented in this study to assess the validity of the proposed multiscale process.

• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.1
• /
• pp.145-156
• /
• 1993

An improved analysis model for material nonlinearity induced by elasto-plastic deformation and damage including large strain response was proposed. The elasto-plastic-damage constitutive model based on the continuum damage mechanics approach was adopted to overcome limitations of the conventional plastic theory, which can manage the anisotropic tonsorial damages evolved during time-independent plastic deformation process of materials. Updated Lagrangian finite element formulation for elasto-plastic damage coupling problem including large deformation, large rotation and large strain problems was completed to develop a numerical model which can predict all kinds of structural nonlinearities and damage rationally. Finally, a finite element analysis code for the 2-dimensional plane problem was developed and the applicability and validity of the numerical model was investigated through some numerial examples. Calculations showed reasonable results in both geometrical nonlinear problem due to large deformation and material nonlinearity including the damage effect.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.8 no.3
• /
• pp.1-10
• /
• 1988

A higher order plate bending finite element using cubic in-plane displacement profiles is proposed for geometrically nonlinear analysis of thin and thick plates. The higher order plate bending element has been derived from the three dimensional plate-like continuum by discretization of the equations of motion by Galerkin weighted residual method, together with enforcing higher order plate assumptions. Total Lagrangian formulation has been used for geometrically nonlinear analysis of plates and consistent linearization by Newton-Raphson method has been performed to solve the nonlinear equations. The element characteristics have been computed by, selective reduced integration technique using Gauss quadrature to avoid shear locking phenomenon in case of extremely thin plates. Several numerical examples were solved with FEAP macro program to demonstrate versatility and accuracy of the present higher order plate bending element.

• Journal of the Korean Society for Heat Treatment
• /
• v.7 no.4
• /
• pp.233-243
• /
• 1994

Constitutive relation of thermoelasto-plastic material undergoing phase transformation during quenching process were developed on the basic of continuum thermodynamics. The metallic structure, temperature and residual stresses distributions were numerically calculated by the finite element technique. The metallic structure were defined by transformation from austenite to pearlite and characterized as a fuction of thermal history and mixture rule of phase. On the distribution of thermal stress along the radial direction, axial and tangential stresses are compressive in the surface, and tential in the inner part. Radial stress is tensile in the whole body. The reversion of residual stress takes plase at 11.5~15.5mm from the center.

• Proceedings of the KSME Conference
• /
• 2007.05a
• /
• pp.780-785
• /
• 2007

Accurate and fast haptic simulations of deformable objects are desired in many applications such as medical virtual reality. In haptic interactions with a coarse model, the number of nodes near the haptic interaction region is too few to generate detailed deformation. Thus, local refinement techniques need to be developed. Many approaches have employed purely geometric subdivision schemes, but they are not proper in describing the deformation behavior of deformable objects. This paper presents a continuum mechanics-based finite element adaptive method to perform haptic interaction with a deformable object. This method superimposes a local fine mesh upon a global coarse model, which consists of the entire deformable object. The local mesh and the global mesh are coupled by the s-version finite element method (s-FEM), which is generally used to enhance accurate solutions near the target points even more. The s-FEM can demonstrate a reliable deformation to users in real-time.

• Steel and Composite Structures
• /
• v.27 no.5
• /
• pp.567-576
• /
• 2018

The objective of this work is to analyze large deflections of a fiber reinforced composite cantilever beam under point loads. In the solution of the problem, finite element method is used in conjunction with two dimensional (2-D) continuum model. It is known that large deflection problems are geometrically nonlinear problems. The considered non-linear problem is solved considering the total Lagrangian approach with Newton-Raphson iteration method. In the numerical results, the effects of the volume fraction and orientation angles of the fibre on the large deflections of the composite beam are examined and discussed. Also, the difference between the geometrically linear and nonlinear analysis of fiber reinforced composite beam is investigated in detail.