• Journal of Korean Society of Steel Construction
• /
• v.24 no.3
• /
• pp.337-348
• /
• 2012

This study investigated the characteristics of bifurcation and the instability due to the initial imperfection of the space truss, which is sensitive to the initial conditions, and the calculated buckling load by the analysis of Eigen-values and the determinant of tangential stiffness. A two-free nodes model, a star dome, and a three-ring dome model were selected as case studies in order to examine the unstable phenomenon due to the sensitivity to Eigen mode, and the influence of the rise-span ratio and the load parameter on the buckling load were analyzed. The sensitivity to the imperfection of the two-free nodes model changed the critical path after reaching the limit point through the bifurcation mode, and the buckling load level was reduced by the increase in the amount of imperfection. The two sensitive buckling patterns for the model can be explained by investigating the displaced position of the free node, and the asymmetric Eigen mode was a major influence on the unstable behavior due to the initial imperfection. The sensitive mode was similar to the in-extensional mechanism basis of the simplified model. Since the rise-span ratio was higher, the effect of local buckling is more prominent than the global buckling in the star dome, and bifurcation on the equilibrium path occurring as the value of the load parameter was higher. Additionally, the buckling load levels of the star dome and the three-ring model were about 50-70% and 80-90% of the limit point, respectively.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.7 no.1
• /
• pp.65-73
• /
• 1987

An efficient plate bending finite element has been developed using higher-order inplane displacement profiles of the plate. The 6-noded, 21-d.o.f. triangular element including shear deformation effect has been derived from the plate-like continuum by the Galerkin's weighted residual method. Square plate examples were tested with selected element meshes and several aspect ratios for their static behavior under uniformly distributed load. The result of the example tests indicated consistently good performance of the present higher-order plate bending element in comparison with the thin and thick plate solution and other existing finite element solutions.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.35 no.5
• /
• pp.404-411
• /
• 2007

new domain/boundary decomposition method is suggested to perform efficient finite element analyses of contact problems. A penalty method is used for connecting an interface or contact interfaces with neighboring subdomains that satisfy continuity conditions. As a result, the derived effective stiffness matrices are always positive definite, and computational efficiency can be improved to a considerable degree. Moreover, any complex-shaped domain can be divided into independently modeled subdomains without considering the conformity of meshes along the interface. Using a computer code based on the present method, these advantageous features are confirmed through a set of numerical examples.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.10
• /
• pp.2640-2652
• /
• 1994

The stiffness of 6-node isotropic element is stiffer than that of 8-node isotropic element of same configuration. This phenomenon was called 'Relative Stiffness Stiffening Phenomenon'. In this paper, an equation of sampling point modification which correct this phenomenon was derived for the composite plate, as well as an equation for an isotropic plate. The relative stiffness stiffening phenomena of an isotropic plate element could be corrected by modifying Gauss sampling points in the numerical integration of stiffness matrix. This technique could also be successfully applied to the static analyses of composite plate modeled by the 3-dimensional 16-node elements. We predicted theoretical errors of stiffness versus the number of layers that result from the reduction of numerical integration order. These errors coincide very well with the actual errors of stiffness. Therefore, we can choose full integration of reduced integration based upon the permissible error criterion and the number of layers by using the thoretically predicted error.

• Computational Structural Engineering
• /
• v.9 no.4
• /
• pp.117-126
• /
• 1996

Damage detection methods using structural tests can be divided into two methods, i.e., static and dynamic. The static methods which use the stiffness properties of the structure are simpler than the dynamic methods. However, static approaches are very sensitive to the displacement measurement noises and modeling errors. The dynamic methods also have limitations in acquiring the natural frequencies and mode shapes of the high frequencies. In this study, a method for the structural damage assessment using sensitivity matrices is developed, in which the drawbacks of the static and dynamic methods can be compensated. Based on the measurement data for the static displacements and dynamic modal properties, the damage locations and the degree of damage are determined using the presented sensitivity matrix method. The efficiency of the proposed method has been examined through numerical simulation studies on truss type structures.

• Geophysics and Geophysical Exploration
• /
• v.26 no.4
• /
• pp.268-274
• /
• 2023

Finite source inversion is performed with a Green's function matrix and geodetic coseismic displacement. Conventionally, the Green's function matrix is constructed using the Okada model (Okada, 1985). However, for more realistic earthquake simulations, recent research has widely adopted the physics-based model, which can consider various material properties such as elasticity, viscoelasticity, and elastoplasticity. We used the physics-based software PyLith, which is suitable for earthquake modeling. However, the PyLith does not provide a mesh generator, which makes it difficult to perform finite source inversions that require numerous subfaults and observation points within the model. Therefore, in this study, we developed CPInterface (COMSOL-PyLith Interface) to improve the convenience of finite source inversion by combining the processes of creating a numerical model including sub-faults and observation points, simulating earthquake modeling, and constructing a Green's function matrix. CPInterface combines the grid generator of COMSOL with PyLith to generate the Green's function matrix automatically. CPInterface controls model and fault information with simple parameters. In addition, elastic subsurface anomalies and GPS observations can be placed flexibly in the model. CPInterface is expected to enhance the accessibility of physics-based finite source inversions by automatically generating the Green's function matrix.

• Computational Structural Engineering
• /
• v.4 no.2
• /
• pp.95-102
• /
• 1991

Some sections of concrete highway pavements have been viewed with great concern by highway officials and engineers due to the severe cracking and failure problems. This is mainly due to the traffic loads in addition to temperature variations between top and bottom of concrete slab, which cause the concrete slab to curl up and down depending on the thermal gradient, respectively. Subsequently, a major consideration was given to the derivation of stiffness matrix and equivalent nodal loads due to the uniform gravity load, temperature and shrinkage of concrete. And the structural behavior of concrete highway pavement due to the temperature variations throughout the nations has been emphasized.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2006.05a
• /
• pp.304-308
• /
• 2006

The focus of this study is modeling technique for a bellows in vehicle exhaust system. Bellows was developed using tile finite element model by replacing with the equivalent beam. The equivalent beam model were studied in detail. Non-structural node in the cross section of original model is given to expressing their motion. Equivalent mass matrix and stiffness matrix calculated using Guyan reduction method. Material Properties of beam was obtained from the direct comparison between equivalent model and that of Timoshenko beam model. The calculated natural frequencies and mode shape are compared with the reference results and coincided well. The results were compared with the confirmed results, which were in good agreement.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.5
• /
• pp.411-420
• /
• 2009

This study presents an extended finite difference method based on moving least squares(MLS) method for solving potential problems with interfacial boundary. The approximation constructed from the MLS Taylor polynomial is modified by inserting of wedge functions for the interface modeling. Governing equations are node-wisely discretized without involving element or grid; immersion of interfacial condition into the approximation circumvents numerical difficulties owing to geometrical modeling of interface. Interface modeling introduces no additional unknowns in the system of equations but makes the system overdetermined. So, the numbers of unknowns and equations are equalized by the symmetrization of the stiffness matrix. Increase in computational effort is the trade-off for ease of interface modeling. Numerical results clearly show that the developed numerical scheme sharply describes the wedge behavior as well as jumps and efficiently and accurately solves potential problems with interface.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.12 no.4_1
• /
• pp.67-76
• /
• 1992

A hierarchical formulation based on p-version of the finite element method for linear elastic axisymmetric stress analysis is presented. This is accomplished by introducing additional nodal variables in the element displacement approximation on the basis of integrals of Legendre polynomials. Since the displacement approximation is hierarchical, the resulting element stiffness matrix and equivalent nodal load vectors are hierarchical also. The merits of the propoosed element are as follow: i) improved conditioning, ii) ease of joining finite elements of different polynomial order, and iii) utilizing previous solutions and computation when attempting a refinement. Numerical examples are presented to demonstrate the accuracy, efficiency, modeling convenience, robustness and overall superiority of the present formulation. The results obtained from the present formulation are also compared with those available in the literature as well as with the analytical solutions.