• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.34 no.3
• /
• pp.895-906
• /
• 2014

A fundamental study of the dynamically penetrating anchor (DPA - colloquially known as torpedo anchor) embedded into deep seabed was conducted using measurement data and numerical approaches. Numerical simulation of such a structure penetration was often suffered by severe mesh distortion arising from very large soil deformation, complex contact condition and nonlinear soil behavior. In recent years, a Coupled Eulerian-Lagrangian method (CEL) has been used to solve geomechanical boundary value problems involving large deformations. In this study, 3D finite element analyses using the CEL formulation are carried out to simulate the construction process of dynamic anchors. Through comparisons with results of field measurements, the CEL method in the present study is in good agreement with the general trend observed by in-situ measurements and thus, predicts a realistic large deformation movement for the dynamic anchors by free-fall dropping, which the conventional FE method cannot. Additionally, the appropriate parametric studies needed for verifying the characteristic of dynamic anchor are also discussed.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.42 no.5
• /
• pp.151-159
• /
• 2000

Dynamic loading of structures often causes excursions of stresses will into the inelastic range and the influence of geometry changes on the response is also significant in may cases. In general , the shell structures designed according to quasi-Static analysis may collapse under condition of dynamic loading. Therefore, for a more realistic prediction on the lad carrying capacity of these shell. both material and geometric nonlinear effects should be considered. In this study , the material nonlinearity effect on the dynamic response is formulated by the elasto-viscoplastic model highly corresponding to the real behavior of the material. Also, the geometrically nonlinear behavior is taken into account using a Total Lagrangian formulation. the reinforcing bars are modeled by the equivalent steel layer at the location of reinforcements, and Von Mises yield criteria is adopted for the steel layer behavior. Also, Drucker-Prager yield criteria is applied for the behavior of concrete. the shape imperfection of dome is assumed as 'dimple type' which can be expressed Wd1=Wd0(1-(r-a)m)n while the shape imperfection of wall is assumed as sinusoidal curve which is Wwi =Wwo sin(n $\pi$y/l). In numerical test, three cases of shape imperfection of 0.0 -5.0cm(opposite direction to loading ; inner shape imperfection)and 5cm (direction to loading : outward shape imperfection) and thickness of steel layer determined by steel ratio of 0,3, and 5% were analyzed. The effect of shape imperfection and steel ratio and behavior characteristics of perfect shape shell and imperfect shape shell are identified through analysis of above mentioned numerical test. Dynamic behaviors of dome and wall according toe combination of shape imperfection and steel ratio are also discussed in this paper.

• Journal of the Korean GEO-environmental Society
• /
• v.7 no.6
• /
• pp.133-143
• /
• 2006

A fully coupled effective stress dynamic analysis procedure for modeling seismic liquefaction on slope is presented. An elasto-plastic formulation is used for the constitutive model UBCSAND in which the yield loci are radial lines of constant stress ratio and the flow rule is non-associated. This is incorporated into the 2D version of Fast Lagrangian Analysis of Continua (FLAC) by modifying the existing Mohr-Coulomb model. This numerical procedure is used to simulate centrifuge test data from the Rensselaer Polytechnic Institute (RPI). UBCSAND is first calibrated to cyclic direct simple shear tests performed on Nevada sand. Both pre- and post-liquefaction behaviour is captured. The centrifuge test is then modeled and the predicted accelerations, excess porewater pressures, and displacements are compared with the measurements. The results are shown to be in general agreement. The procedure is currently being used in the design of liquefaction remediation measures for a number of dam, bridge, tunnel, and pipeline projects in Western Canada.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.26 no.1
• /
• pp.24-38
• /
• 1989

In this study, an attempt for formulating a new and simplified rectangular finite element having only four corner nodal points is made to analyze the elastic-plastic large deformation behaviour up to the ultimate limit state of plates with initial imperfections. The present finite element contains the geometric nonlinearity caused by both in-plane and out-of-plane large deformation because for very thin plates the influence of the former may not be negligible. Treatment of expanded plastic zone in the plate thickness direction of the element is simplified based upon the concept of plastic node method so that the elastic-plastic stiffness matrix of the element is derived by the simple matrix operation without performing complicated numerical integration. Thus, a considerable saving of the computational efforts is expected. A computer program is also completed based on the present formulation and numerical calculation for some examples is performed so as to verify the accuracy and validity of the program.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.13 no.4
• /
• pp.437-448
• /
• 2000

In this study, a numerical analysis method for soil-pile interaction in frequency domain problem is presented. The total soil-pile interaction system is divided into two parts so called near field and far field. In the near field, beam elements are used lot a pile and plain strain finite elements for soil. In the far field, dynamic fundamental solution for multi-layered half planes based on boundary element formulation is adopted lot soil. These two fields are coupled using FE-BE coupling technique. In order to verify the proposed soil-pile interaction analysis, the dynamic responses of pile on multi-layered half planes are simulated and the results are compared with the experimental results. Also, various numerical analyses of piles considering different conditions of soil-pile interaction system are performed to examine the dynamic behavior of the system. It has been found that the developed method which satisfies the radiation conditions of multi-layered half planes can be applied to various structure systems effectively in frequency domain.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.13 no.4
• /
• pp.449-459
• /
• 2000

The original Mindlin-type degenerated shell element perform reasonably well for moderately thick shell structures. However, when full integration for analysis of thin shell is used to evaluate the stiffness matrix, the stiffness of shell element is often over-estimated due to shear or membrane locking phenomena. To correct this problem, the formulation of the new degenerated shell element is derived by the combination of two different techniques. The first type of elements(TypeⅠ) has used assumed shear strains in the natural coordinate system to overcome the shear locking problem, the reduced integration technique in in-plane strains(membrane strains) to avoid membrane locking behaviour. Another element(TypeⅡ) has applied the assumed strains to both of membrane strain and transverse shear strains. The improved degenerated shell element has been tested by several numerical problems of shell structures. Numerical results indicate that this shell element shows fast convergence and reliable solutions.

• Journal of Korean Association for Spatial Structures
• /
• v.11 no.4
• /
• pp.49-59
• /
• 2011

This paper presents a new numerical method by using Singular Value Decomposition force method for form-finding of tensegrity structures with constraints. In order to maintain the tensegrity structures stable, state of self-equilibrium stress should be introduced. The existing force method has some advantages about form-finding of tensegrity structures. However, this method has complex formulation. In force density method, dummy members are required for form -finding of tensegrity structures with constraints. Therefore this study proposes new force method using Singular value decomposition. The proposed method is both having easy basic consept and simple computation than existing force method. A numerical example is presented to demonstrate the efficiency and robustness in searching initial single integral feasible self-stress mode for tensegrity structures with constraints.

• Journal of the Korean Geotechnical Society
• /
• v.21 no.8
• /
• pp.37-43
• /
• 2005

Hoek-Brown model, which was developed in order to predict the behavior of rock mass, has widely been utilized and revised by many researchers to solve various problems encountered in tunnelling and slope stability analysis. However, there is no schematic investigation on the application of the Hoek-Brown model to numerical analysis including finite element simulations. In this paper the Hoek-Brown model was formulated as a constitutive model according to the procedure of generalized plasticity theory, and a Rounded Hoek-Brown model, which could overcome the numerical difficulties by modifying the edge part of the yield surface as a curve shape, was newly proposed. The new model could satisfy the requirements as an elasto-plastic constitutive soil model and follow the yield surface of the original Hoek-Brown model in the compression mode. The constitutive equation for the proposed model herein was established and presented to be applicable to the generalized nonlinear finite element analysis.

• Steel and Composite Structures
• /
• v.10 no.3
• /
• pp.207-222
• /
• 2010

Fire incident in buildings is common, so the fire safety design of the framed structure is imperative, especially for the unprotected or partly protected bare steel frames. However, software for structural fire analysis is not widely available. As a result, the performance-based structural fire design is urged on the basis of using user-friendly and conventional nonlinear computer analysis programs so that engineers do not need to acquire new structural analysis software for structural fire analysis and design. The tool is desired to have the capacity of simulating the different fire scenarios and associated detrimental effects efficiently, which includes second-order P-D and P-d effects and material yielding. Also the nonlinear behaviour of large-scale structure becomes complicated when under fire, and thus its simulation relies on an efficient and effective numerical analysis to cope with intricate nonlinear effects due to fire. To this end, the present fire study utilizes a second-order elastic/plastic analysis software NIDA to predict structural behaviour of bare steel framed structures at elevated temperatures. This fire study considers thermal expansion and material degradation due to heating. Degradation of material strength with increasing temperature is included by a set of temperature-stress-strain curves according to BS5950 Part 8 mainly, which implicitly allows for creep deformation. This finite element stiffness formulation of beam-column elements is derived from the fifth-order PEP element which facilitates the computer modeling by one member per element. The Newton-Raphson method is used in the nonlinear solution procedure in order to trace the nonlinear equilibrium path at specified elevated temperatures. Several numerical and experimental verifications of framed structures are presented and compared against solutions in literature. The proposed method permits engineers to adopt the performance-based structural fire analysis and design using typical second-order nonlinear structural analysis software.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.1
• /
• pp.139-149
• /
• 1994

A numerical simulation of unsteady axisymmetric turbulent flow was performed for a reciprocating engine including port/valve assembly. The governing equations based on a nonorthogonal coordinate formulation with Cartesian velocity components were used and discretised by the finite volume method with non-staggered variable arrangements. The modified $\kappa-\xi$. turbulence model which included the effect of compressibility was used. The results of twodimensional transient calculation for the axisymmetric configuration were compared with the experimental data. Although slightly low rms velocity was predicted compared to the experimental data, predicted velocity distributions at the valve exit and in-cylinder region showed good agreements with the experimental data. The flow at the valve exit was separated at the same valve lift position with the experimental data. Two vortices incylinder region were generated during the initial intake process. The clockwise main vortex became strong and moved upward to the top wall. The counter-clockwise second vortex became weak and stick to the upper left corner of the cylinder. After middle intake process, new vortex adjacent to upper cylinder wall appeared by the piston motion and therefore, the in-cylinder flow was formed into three vortices. The cylinder pressure just before bottom dead center of piston was higher than inlet pressure and then the reverse flow occured at the valve exit. The in-cylinder flow characteristics were strongly dependent on piston motion, but insensitive to valve motion.