• Geomechanics and Engineering
• /
• v.17 no.6
• /
• pp.573-581
• /
• 2019

Excavation of underground cavern and shaft was proposed for the construction of a ventilation facility in an urban area. A shaft connects the street-level air plenum to an underground cavern, which extends down approximately 46 m below the street surface. At the project site, the rock mass was relatively strong and well-defined joint sets were present. A kinematic block stability analysis was first performed to estimate the required reinforcement system. Then a 3-D discontinuum numerical analysis was conducted to evaluate the capacity of the initial support and the overall stability of the required excavation, followed by a 3-D continuum numerical analysis to complement the calculated result. This paper illustrates the application of detailed numerical analyses to the design of the required initial support system for the stability of underground hard rock mining at a relatively shallow depth.

• Tunnel and Underground Space
• /
• v.20 no.1
• /
• pp.28-38
• /
• 2010

In order to investigate the tendency of general displacements and behaviors with respect to each construction process as well as the applicability of numerical analysis schemes, this research has focused on not only analyzing a variety of field observations made in a NATM tunnel, such as displacement of top and side, stress of shotcrete and axial strength of rock bolt, but also carrying out a series of numerical analyses. It was established from the investigation that the 2-dimensional continuum numerical analysis was the one which could more accurately predict displacement of crown and side in the area of one step excavation (patten, P1-P3), while the 2-dimensional discontinuum analysis was the most suitable scheme to study that of two step excavation (patten, P4-P6). In addition, the 2-dimensional continuum analysis enabled to appropriately predict the axial strength of rock bolt and stress of shotcrete in all the area of the tunnel. Finally, it has been possible to conclude from the study that the 3-dimensional continuum analysis should be applied to inspect the behavior and tendency with respect to each stage of the construction as well as in the case of joints, such as large turnouts where relaxation loads in both of horizontal and vertical direction are piled up.

• Journal of the Korean Geotechnical Society
• /
• v.19 no.3
• /
• pp.65-74
• /
• 2003

The dragload on pile groups in consolidating ground was investigated based on a numerical analysis. The case of a single pile and subsequently the response of groups were analyzed by 2D and 3D finite element studies. Conventional continuum elements and special slip elements were used in the analyses for comparison. Based on a limited parametric study, it is shown that dragload for a single pile and group effect are normally overestimated by continuum analyses, compared with the predictions by the slip analyses. The group effect was examined from the slip analysis by considering various factors such as pile configurations, surface loading, interface friction coefficient, and axial loading on piles. An examplary analysis and one previous experimental observation of dragload and group effects were back-analysed. The case histories demonstrated that the slip analysis might predict a better estimate of dragload and group effect compared to the no-slip continuum analysis.

• Steel and Composite Structures
• /
• v.38 no.6
• /
• pp.717-737
• /
• 2021

A theoretical energy-based model to capture the mechanical response of thick woven composite laminates, which are used in such applications as maritime or aerospace, to high-velocity impact was developed. The dependences of the impact phenomenon on material and geometrical parameters were analysed making use of the Vaschy-Buckingham Theorem to provide a non-dimensional framework. The model was divided in three different stages splitting the physical interpretation of the perforation process: a first where different dissipative mechanisms such as compression or shear plugging were considered, a second where a transference of linear momentum was assumed and a third where only friction took place. The model was validated against experimental data along with a 3D finite element model. The numerical simulations were used to validate some of the new hypotheses assumed in the theoretical model to provide a more accurate explanation of the phenomena taking place during a high-velocity impact.

• Structural Engineering and Mechanics
• /
• v.40 no.2
• /
• pp.257-277
• /
• 2011

In this paper the geometrically nonlinear continuum plate finite element model, hitherto not reported in the literature, is developed using the total Lagrange formulation. With the layerwise displacement field of Reddy, nonlinear Green-Lagrange small strain large displacements relations (in the von Karman sense) and linear elastic orthotropic material properties for each lamina, the 3D elasticity equations are reduced to 2D problem and the nonlinear equilibrium integral form is obtained. By performing the linearization on nonlinear integral form and then the discretization on linearized integral form, tangent stiffness matrix is obtained with less manipulation and in more consistent form, compared to the one obtained using laminated element approach. Symmetric tangent stiffness matrixes, together with internal force vector are then utilized in Newton Raphson's method for the numerical solution of nonlinear incremental finite element equilibrium equations. Despite of its complex layer dependent numerical nature, the present model has no shear locking problems, compared to ESL (Equivalent Single Layer) models, or aspect ratio problems, as the 3D finite element may have when analyzing thin plate behavior. The originally coded MATLAB computer program for the finite element solution is used to verify the accuracy of the numerical model, by calculating nonlinear response of plates with different mechanical properties, which are isotropic, orthotropic and anisotropic (cross ply and angle ply), different plate thickness, different boundary conditions and different load direction (unloading/loading). The obtained results are compared with available results from the literature and the linear solutions from the author's previous papers.

• Journal of Welding and Joining
• /
• v.25 no.1
• /
• pp.63-69
• /
• 2007

Once assessment of material failure characteristics is captured precisely in a unified way, it can bedirectly incorporated into the structural failure assessment under various loading environments, based on the theoretical backgrounds so called Local Approach to Fracture. The aim of this study is to develop a numerical fatigue test method by continuum damage mechanics applicable for the assessment of structural integrity throughout crack initiation and structural failure based on the Local Approach to Fracture. The generalized elasto-visco-plastic constitutive equation, which can consider the internal damage evolution behavior, is developed and employed in the 3-D FEA code in order to numerically evaluate the material and/or structural responses. Explicit information of the relationships between the mechanical properties and material constants, which are required for the mechanical constitutive and damage evolution equations for each material, are implemented in numerical fatigue test method. The material constants selected from constitutive equations are used directly in the failure assessment of material and/or structures. The performance of the developed system has been evaluated with assessing the S-N diagram of stainless steel materials.

• Structural Engineering and Mechanics
• /
• v.62 no.1
• /
• pp.33-42
• /
• 2017

We present in this paper a finite element formulation for nonlinear torsional analysis of 3D beams with arbitrary composite cross-sections. Since the proposed formulation employs a continuum mechanics based beam element with kinematics enriched by the extended St. Venant solutions, it can precisely account higher order warping effect and its 3D couplings. We propose a numerical procedure to calculate the extended St. Venant equation and the twisting center of an arbitrary composite cross-section simultaneously. The accuracy and efficiency of the proposed formulation are thoroughly investigated through representative numerical examples.

• Structural Engineering and Mechanics
• /
• v.33 no.6
• /
• pp.657-674
• /
• 2009

In this paper a 3-D continuum damage mechanics formulation for composite laminates and its implementation into a finite element model that is based on the layer-wise laminate plate theory are described. In the damage formulation, each composite ply is treated as a homogeneous orthotropic material exhibiting orthotropic damage in the form of distributed microscopic cracks that are normal to the three principal material directions. The progressive damage of different angle ply composite laminates under quasi-static loading that exhibit the free edge effects are investigated. The effects of various numerical modeling parameters on the progressive damage response are investigated. It will be shown that the dominant damage mechanism in the lay-ups of [+30/-30]s and [+45/-45]s is matrix cracking. However, the lay-up of [+15/-15] may be delaminated in the vicinity of the edges and at $+{\theta}/-{\theta}$ layers interfaces.

• Journal of the Korean Society of Manufacturing Technology Engineers
• /
• v.7 no.3
• /
• pp.110-117
• /
• 1998

A continuum-based configuration design sensitivity analysis method is developed for kinematically driven mechanical systems. The configuration design variable for mechanical systems is defined. The 3-1-3 Euler angle is employed as the orientation design variable. Kinematic admissibility conditions of configuration design change. Direct differentiation method is used to derive the governing equations of the design sensitivity. Numerical examples are presented to demonstrate the validity and effectiveness of the proposed method.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.1
• /
• pp.39-46
• /
• 1993

The need to develop method which can improve the shape design efficiency using high quality approximation is being brought up. In this study, to perform shape optimal design of three-dimensional continuum structures an efficient approximation method for stress constraints is proposed, based on expanding the nodal forces in Taylor series with respect to shape variables. Numerical examples are performed using the 3-D cantilever beam and fixed-fixed beam and compared with other method to demonstrate the efficiency and convergence rate of the Force Approximation method. It is shown that by taking advantage of this high quality approximation, the total number of finite element analysis required for shape optimization of 3-D continuum structures can be reduced significantly, resulting to the same level of efficiency achieved previously in sizing optimization problems. Also, shape representation by super curve technique applied to obtain optimal shape finds useful method.