• Geomechanics and Engineering
• /
• v.11 no.1
• /
• pp.77-94
• /
• 2016

A forked tunnel, as a special complicated underground structure, is composed of big-arch tunnel, multi-arch tunnel, neighborhood tunnels and separate tunnels according to the different distances between two separate tunnels. Due to the complicated process of design and construction, surrounding jointed rock mass stability of the big-arch tunnel which belongs to the forked tunnel during excavation is a hot issue that needs special attentions. In this paper, an elasto-plastic damage constitutive model for jointed rock mass is proposed based on the coupling method considering elasto-plastic and damage theories, and the irreversible thermodynamics theory. Based on this elasto-plastic damage constitutive model, a three dimensional elasto-plastic damage finite element code (D-FEM) is implemented using Visual Fortran language, which can numerically simulate the whole excavation process of underground project and perform the structural stability of the surrounding rock mass. Comparing with a popular commercial computer code, three dimensional fast Lagrangian analysis of continua (FLAC3D), this D-FEM has advantages in terms of rapid computing process, element grouping function and providing more material models. After that, FLAC3D and D-FEM are simultaneously used to perform the structural stability analysis of the surrounding rock mass in the forked tunnel considering three different computing schemes. The final numerical results behave almost consistent using both FLAC3D and D-FEM. But from the point of numerically obtained damage softening areas, the numerical results obtained by D-FEM more closely approach the practical behaviors of in-situ surrounding rock mass.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.8 no.2
• /
• pp.19-25
• /
• 2004

The linear analysis for the balance of linear momentum of a structure is relatively easy to perform, but the error becomes large when the structure experiences large deformation. Therefore, the material and geometric nonlinearity need to be considered for the precise calculations in that case. The plastic flow of a ductile steel-like metal mainly transforms its dissipated mechanical energy into heat, which transfers under the first and second law of thermodynamics. This heat increases the temperature of the material and the strength of the material decreases accordingly, which affects mechanical behavior of the given structure. This paper presents a finite-strain thermo-elasto-plastic steel model. This model can handle large deformation and thermal load simultaneously, which is common during earthquake periods. Two 3-dimensional finite element analyses verify this formulation.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.11 no.6
• /
• pp.145-151
• /
• 2012

In this study, finite element analysis is carried out to estimate horizontal forces needed for the required power calculation and vertical forces applied on the structural analysis model for the development of automatic serrated surface at metal gasket core machining system. By considering of elasto-plastic material characteristics, nonlinear contact analysis was conducted to compute these loads according to the change of roll reduction, frictional coefficient and core thickness. As the result, horizontal and vertical reaction force variations are found according to parameters and maximum reaction force is also confirmed to be most affected by roll reduction.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.3
• /
• pp.353-364
• /
• 2007

In this paper, a nonlinear finite element analysis program NUCAS, which has been developed for assessment of ultimate pressure capacity and failure mode for nuclear containment building is described. Degenerated shell element with assumed strain method and low-order solid element with enhanced assumed strain method is adapted to microscopic material and elasto-plastic material model, respectively. Finally, the performance of the developed program is tested and demonstrated with several examples. From the numerical tests, the present results show a good agreement with experimental data or other numerical results.

• Computers and Concrete
• /
• v.15 no.4
• /
• pp.515-545
• /
• 2015

The paper presents results of FE simulations of the strain-rate sensitive concrete behaviour under dynamic loading at the macroscopic level. To take the loading velocity effect into account, viscosity, stress modifications and inertial effects were included into a rate-independent elasto-plastic formulation. In addition, a decrease of the material stiffness was considered for a very high loading velocity to simulate fragmentation. In order to ensure the mesh-independence and to properly reproduce strain localization in the entire range of loading velocities, a constitutive formulation was enhanced by a characteristic length of micro-structure using a non-local theory. Numerical results were compared with corresponding laboratory tests and available analytical formulae.

• Computers and Concrete
• /
• v.20 no.6
• /
• pp.689-704
• /
• 2017

The paper presents experimental and numerical investigations of prefabricated composite structural building reinforced concrete slabs with the insulating material for a residential building construction. The building slabs were composed of concrete and expanded polystyrene. In experiments, the slabs in the full-scale 1:1 were subjected to vertical concentrated loads and failed along a diagonal shear crack. The experiments were numerically evaluated using the finite element method based on two different constitutive continuum models for concrete. First, an elasto-plastic model with the Drucker-Prager criterion defined in compression and with the Rankine criterion defined in tension was used. Second, a coupled elasto-plastic-damage formulation based on the strain equivalence hypothesis was used. In order to describe strain localization in concrete, both models were enhanced in the softening regime by a characteristic length of micro-structure by means of a non-local theory. Attention was paid to the formation of critical diagonal shear crack which was a failure precursor.

• Journal of Welding and Joining
• /
• v.22 no.6
• /
• pp.25-29
• /
• 2004

Some industrial steel structures are composed by components linked by several welding joints to constitute an assembly. The main interest of assembly simulation is to evaluate the global distortion of welded structure. The general method, thermo-elasto-plastic analysis, leads to excessive model size and computation time. In this study, a simplified method called "Local and Global approach" was developed to break down this limit and to provide a accurate solution for distortion. Local and global approach is composed of 3 steps; 1) Local simulation of each welding joint on a dedicated mesh (usually very fine due to high thermal gradients), taking into account for the non linearity of the material properties and the moving heat source. 2) Transfer to the global model of the effects of the welding joints by projection of the plastic strain tensors. 3) Elastic simulation to determine final distortions in global model. The welding deformation test for mock-up structure was performed to verify this approach. The predicted welding distortion by this approach had a good agreement with experiment results.

• Computers and Concrete
• /
• v.2 no.2
• /
• pp.97-110
• /
• 2005

Rheological behaviour of fresh concrete is an important factor in controlling concrete quality. It is recognized that the measurement of the slump is not a sufficient test method to adequately characterize the rheology of fresh concrete. To further understand the slump measurement and its relationship to the rheological properties, an elasto-viscoplastic, 2-D axisymmetric finite element (FE) model is developed. The FE model employs the Bingham material model to simulate the flow of a slump test. An experimental program is carried out using the Slump Rate Machine (SLRM_II) to evaluate the finite element simulation results. The simulated slump-versus-time curves are found to be in good agreement with the measured data. A sensitivity study is performed to evaluate the effects of yield stress, plastic viscosity and cone withdrawal rate on the measured flow curve using the FE model. The results demonstrate that the computed yield stress compares well with reported experimental data. The flow behaviour is shown to be influenced by the yield stress, plastic viscosity and the cone withdrawal rate. Further, it is found that the value of the apparent plastic viscosity is different from the true viscosity, with the difference depending on the cone withdrawal rate. It is also confirmed that the value of the final slump is most influenced by the yield stress.

• Structural Engineering and Mechanics
• /
• v.27 no.2
• /
• pp.135-148
• /
• 2007

A new finite shear wall element model and a method for calculation of 3D multi-storied only shear walled or shear walled - framed structures using finite shear wall elements assumed ideal elasto - plastic material are developed. The collapse load of the system subjected to factored constant gravity loads and proportionally increasing lateral loads is calculated with a method of load increments. The shape functions over the element are determined as a cubic variation along the story height and a linear variation in horizontal direction because of the rigid behavior of the floor slab. In case shear walls are chosen as only one element in every floor, correct solutions are obtained by using this developed element. Because of the rigid behavior of the floor slabs, the number of unknowns are reduced substantially. While in framed structures, classical plastic hinge hypothesis is used, in nodes of shear wall elements when vertical deformation parameter is exceeded ${\varepsilon}_e$, this node is accepted as a plastic node. While the system is calculated with matrix displacement method, for determination of collapse safety, plastic displacements and plastic deformations are taken as additional unknowns. Rows and columns are added to the system stiffness matrix for additional unknowns.

• Journal of Korean Society of Steel Construction
• /
• v.21 no.6
• /
• pp.597-608
• /
• 2009

In this paper, large-deformation elasto-plastic analysis of space frames that considersjoint connection properties is presented. This method is based on the large-deformation formula with finite rotation, which was developed initially for elastic systems, and is extended herein to include the elasto-plastic effect and the member joint connection properties of semi-rigid what?. The analytical method was derived from the Eulerian concept, which takes into consideration the effects of large joint translations and rotations. The localmember force-deformation relationships were obtained from the beam-column approach, and the change caused by the axial strain in the member chord lengths and flexural bowing were taken into account. The effect of the axial force of the member on bending and torsional stiffness, and on the plastic moment capacity, is included in the analysis. The material is assumed to be ideally elasto-plastic, and yielding is considered concentrated at the member ends in the form of plastic hinges. The semi-rigid properties of the member joint connection are considered based on the power or linear model. The arc length method is usedto trace the post-buckling range of the elastic and elasto-plastic problems with the semi-rigid connection. A sample non-linear buckling analysis was carried out with the proposed space frame formulations to demonstrate the potential of the developed method in terms of its accuracy and efficiency.