• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2003.03a
• /
• pp.16-22
• /
• 2003

The focal mechanisms of the earthquakes occurred in 2001 and 2002 are analyzed to understand the regional stress and tectonics in and around Korean Peninsula. The forty -three fault plane solutions are derived using the polarities of first arrival P-waves recorded at KIGAM, Bmh and KEPRI stations. The result suggest that thrust motion with significant amount of strike slip component is dominant mode of faulting. The larger population of mechanism is characterized by WNW-ESE striking nodal planes. NE-SW direction is interpreted as dominant compressional axis orientation of stress field. These solutions are similar to those of medium size earthquakes studied previously, which is known as typical regional tectonic stress orientation in and around Korean Peninsula.

• Structural Engineering and Mechanics
• /
• v.13 no.2
• /
• pp.155-170
• /
• 2002

In this study, a new functional is obtained for folded plates with geometric (kinematic) and dynamic (natural) boundary conditions. This functional is the combination of two different functionals. Both functionals are obtained for thick plates which carry in-plane and lateral forces. A new mixed finite element is developed with $4{\times}13$ nodal parameters for folded plates (REC52). Forces and moments which are the necessary unknowns in engineering problems are obtained directly using the technique suggested here. The use of the global co-ordinate system causes time consuming operations and therefore the Lagrange multiplier method is used to relate the components of the parameters on the fold line. Numerical results are presented for folded plates and compared with experimental results.

• Structural Engineering and Mechanics
• /
• v.51 no.4
• /
• pp.601-625
• /
• 2014

Based on continuum mechanics and the principle of virtual displacements, incremental total Lagrangian formulation (T.L.) and incremental updated Lagrangian formulation (U.L.) were presented. Both T.L. and U.L. considered the large displacement stiffness matrix, which was modified to be symmetrical matrix. According to the incremental updated Lagrangian formulation, small strain, large displacement, finite rotation of three dimensional Timoshenko fiber beam element tangent stiffness matrix was developed. Considering large displacement and finite rotation, a new type of tangent stiffness matrix of the beam element was developed. According to the basic assumption of plane section, the displacement field of an arbitrary fiber was presented in terms of nodal displacement of centroid of cross-area. In addition, shear deformation effect was taken account. Furthermore, a nonlinear finite element method program has been developed and several examples were tested to demonstrate the accuracy and generality of the three dimensional beam element.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.2
• /
• pp.189-196
• /
• 2009

This paper presents the subparametric finite element model formulated by partial-linear layerwise theory for the analysis of laminate composites. The proposed model is based on refined approximations of two dimensional plane for orthotropic thick laminate plate as well as thin case. Three dimensional problem can be reduced to two dimensional case by assuming piecewise linear variation of in-plane displacement and a constant value of out-of-plane displacement across the thickness. The integrals of Legendre polynomials are chosen to define displacement fields and Gauss-Lobatto numerical integration is implemented in order to directly obtain maximum values occurred at the nodal points of each layer without other extrapolation techniques. The validity and characteristics of the proposed model have been tested by using orthotropic multilayered plate problem as compared to the values available in the published references. In this study, the convergence test has been carried out to determine the optimal layer model in terms of central deflection and stresses. Also, the distribution of displacements and stresses across the thickness has been investigated as the number of layer is increased.

• 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.

• Computational Structural Engineering
• /
• v.7 no.4
• /
• pp.103-111
• /
• 1994

Based on a variational principle of the consistent shear deformable discrete laminate theory derived in the companion paper Part I, a finite element procedure for the vibration analysis of laminated composite plates is presented. The present formulation takes the in-plane displacements of an arbitrary layer, the rotations of the cross section of each layer and transverse displacement of the plate as the state variables at a nodal point of finite element, resulting in total nodal degree of freedom of 2(n+l) +1 for the n-layered laminate. Thus, it allows to specify displacement boundary conditions of layer stretching and/or rotation of layer cross sections around the plate edge and/or lateral displacement. The developed procedure is applied to the free vibration problem for sandwich-type hybrid laminates composed of layers with drastically different material properties whose elasticity solutions are known. Comparison of analysis results with other FEM solutions showed that the present formulation yields better accuracy.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.13 no.1
• /
• pp.379-395
• /
• 2021

For the complicated structural details in ships and offshore structures, the traditional hotspot stress approaches are known to be sensitive to the element variables of element topologies, sizes, and integration schemes. This motivated to develop a new approach for predicting reasonable hotspot stresses, which is less sensitive to the element variables and easy to be implemented the real marine structures. The three-point bending tests were conducted for the longitudinal attachments with the round and rectangular weld toes. The tests were reproduced in the numerical simulations using the solid and shell element models, and the simulation technique was validated by comparing the experimental stresses with the simulated ones. This paper considered three hotspot stress approaches: the ESM method based on surface stress extrapolation, the Dong's method based on nodal forces along a weld toe, and the proposed method based on nodal forces perpendicular to an imaginary vertical plane at a weld toe. In order to study the element sensitivities of each method, 16 solid element models and 8 shell element models were generated under the bending and tension loads, respectively. The element sensitivity was analyzed in terms of Stress Concentration Factors (SCFs) in viewpoints of two statistical quantities of mean and bias with respect to the reference SCFs. The average SCFs predicted by the proposed method were remarkably in good agreement with the reference SCFs based on the experiments and the ship rules. Negligibly small Coefficients of Variation (CVs) of the SCFs, which is measure of statistical bias, were drawn by the proposed method.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.37 no.3
• /
• pp.531-539
• /
• 2017

The strut-tie model approach has been recognized as an efficient methodology for the design of all types of concrete members with D-regions, and the approach has been accepted in design codes globally. However, the design of concrete members with the approach requires many iterative numerical structural analyses, numerous graphical calculations, enormous times and efforts, and designer's subjective decisions in terms of the development of appropriate strut-tie model, determination of required areas of struts and ties, and verification of strength conditions of struts and nodal zones. In this study, a computer graphics program, that enables the design of concrete members efficiently and professionally by overcoming the forementioned limitations of the strut-tie model approach, is developed. In the computer graphics program, the numerical programs that are essential in the strut-tie model analysis and design of concrete members including finite element analysis programs for the plane truss and solid problems with all kinds of boundary conditions, a program for automatic determination of effective strengths of struts and nodal zones, and a program for graphical verification of developed strut-tie model's appropriateness by displaying various geometrical shapes of struts and nodal zones, are loaded. Great efficiency and convenience during the application of the strut-tie model approach may be provided by the various graphics environment-based functions of the proposed program.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.13 no.1
• /
• pp.25-35
• /
• 2000

The superstructure of general bridge like slab bridge and slab on girder bridge is composed of elastically supported isotropic plate. The objective of this study is to develop the new analysis method for elastically supported plate with general edge beam or girder(boundaries) under arbitrary out of plane loading. The displacement solutions for the macro-element of plate and beam are obtained by solving for the unknown interactive forces and moments at the beam or nodal line locations after satisfying equilibrium equation along the nodal line. The displacement functions for macro-elements ate proposed in single Fourier series using harmonic analysis, and the equilibrium equations of nodal line are composed by using slope-deflection method. The proposed analysis method is programmed by MS-Fortran and can be applied to all types of isotropic decks with bridge-type boundaries. Numerical examples involving elastically supported plates with various aspect ratio, loading cases, and bridge-type boundary conditions are presented to demonstrate the accuracy of this program. The major advantage of this new analysis method is the development of a simple solution algorithm, leads to obtain rapidly responses of bridge deck system. This proposed method can be used in parametric study of behavior of bridge decks.

• Journal of Korean Society of Steel Construction
• /
• v.12 no.5 s.48
• /
• pp.515-525
• /
• 2000

The purpose of this paper was to develop size & shape optimization programming algorithm of plane trusses. The optimum techniques applied in this study were extended penalty method of Sequential Unconstrained Minimization Techniques(SUMT) and direct search method with multi-variables proposed by Hooke & Jeeves. Upper mentioned two methods were used iteratively at each level of size and shape optimization routines. The design variables of size optimization were circular steel tube(structural member) diameter and thickness, those of shape optimization were joint coordinates, and the objective function was represented as total weight of truss. During the optimum design, two level procedures of size and shape optimization were interacted iteratively until the final optimum values were attained. At the previous studies about shape optimization of truss, the member sectional areas and coordinates were applied as design variables. So that they could not apply the buckling effect of compression member. In this paper, actual sizes of member and nodal coordinates are used as design variables to consider the buckling effect of compression member properly.