• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.25 no.5
• /
• pp.413-420
• /
• 2012

In this study, the Kernel integration scheme for 2D linear elastic direct boundary element method has been discussed on the basis of subparametric element. Usually, the isoparametric based boundary element uses same polynomial order in the both basis function and mapping function. On the other hand, the order of mapping function is lower than the order of basis function to define displacement field when the subparametric concept is used. While the logarithmic numerical integration is generally used to calculate Kernel integration as well as Cauchy principal value approach, new formulation has been derived to improve the accuracy of numerical solution by algebraic modification. The subparametric based direct boundary element has been applied to 2D elliptical partial differential equation, especially for plane stress/strain problems, to demonstrate whether the proposed algebraic expression for integration of singular Kernel function is robust and accurate. The problems including cantilever beam and square plate with a cutout have been tested since those are typical examples of simple connected and multi connected region cases. It is noted that the number of DOFs has been drastically reduced to keep same degree of accuracy in comparison with the conventional isoparametric based BEM. It is expected that the subparametric based BEM associated with singular Kernel function integration scheme may be extended to not only subparametric high order boundary element but also subparametric high order dual boundary element.

• Journal of Korean Society of Steel Construction
• /
• v.16 no.6 s.73
• /
• pp.819-832
• /
• 2004

This study describes an efficient and facile method for configuration optimum design of structures. One of the ways to achieve numerical shape representation and the selection of design variables is using the design element concept. Using this technique, the number of design variables could be drastically reduced. Isoparametric mapping was utilized to automatically generate the finite element mesh during the optimization process, and this made it possible to easily calculate the derivatives of the coordinates of generated finite element nodes w.r.t. the design variables. For the structural analysis, finite element analysis was adopted in the optimization procedure, and two different techniques(the deterministic method, a modified method of feasible direction; and the stochastic method, a genetic algorithms) were applied to obtain the minimum volumes and section areas for an efficient configuration optimization procedure. Futhermore, spline interpolation was introduced to present a realistic optimum configuration that meet the manufacturing requirements. According to the results of several numerical examples(steel structures), the two techniques suggested in this study simplified the process of configuration optimum design of structures, and yielded improved objective function values with a robust convergence rate. This study's applicability and capability have therefore been demonstrated.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.8 no.2
• /
• pp.147-158
• /
• 2004

For stiffened plates composed of composite materials, many researchers have used a finite element method which connected isoparametric plate elements and beam elements. However, the finite element method is difficult to reflect local behavior of stiffener because beam elements are transferred stiffness for nodal point of plate elements, especially the application is limited in case of laminated composite structures. In this paper, for analysis of laminated composite stiffened plates, 3D shell elements for stiffener and plate are employed. Reissner-Mindlin's first order shear deformation theory is considered in this study. But when thickness will be thin, isoparamatric plate bending element based on the theory of Reissner-Mindlin is generated by transverse shear locking. To eliminate the shear locking and virtual zero energy mode, the substitute shear strain field is used. A deflection distribution is investigated for simple supported rectangular and skew stiffened laminated composite plates with arbitrary orientation stiffener as not only variation of slenderness and aspect ratio of the plate but also variation of skew angle of skew stiffened plates.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.19 no.11
• /
• pp.205-212
• /
• 2018

In order to overcome shortcomings of commercial analysis program for solving certain geotechnical problems, finite element method adopting isoparametric quadrilateral element was selected as a tool for analyzing soil behavior and calculating process was programmed. Two examples were considered in order to verify reliability of the developed program. One of the two examples is the case of acting isotropic confining pressure on finite element and the other is the case of acting shear stress on the sides of the finite element. Isoparametric quadrilateral element was considered as the finite element and displacements in the element can be expressed by node displacements and shape functions in the considered element. Calculating process for determining strain which is defined by derivatives using global coordinates was coded using the Jacobian and the natural coordinates. Four point Gauss rule was adopted to convert double integral which defines stiffness of the element into numerical integration. As a result of executing analysis of the finite element under isotropic confining pressure, calculated stress corresponding to four Gauss points and center of the element were equal to the confining pressure. In addition, according to the analyzed results for the element under shear stress, horizontal stresses and vertical stresses were varied with positions in the element and the magnitudes and distribution pattern of the stresses were thought to be rational.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.7 no.2
• /
• pp.21-28
• /
• 1987

The paper describes the analysis of large displacement problems including instability phenomena. The element used in this is a degenerated isoparametric shell element with eight nodes. Total Lagrangian formulation has been adopted in this study using Newton-Raphson iteration method with incremental load. The linear stability analyses performed usually for the initial position can be repeated at several advanced fundamental states on the non-linear buckling path. Thus a current estimate of the failure load is given. The numerical examples of a cylindrical panel under uniform load, simply supported plate under axial load, and clamped plate under uniform load are carried out. The examples applying degenerated isoparametric elements to bifurcation buckling and nonlinear collapse problems are also performed.

• Computational Structural Engineering
• /
• v.10 no.2
• /
• pp.179-188
• /
• 1997

A 4-node element for efficient finite element analysis of plate bending is presented in this paper. This element is formulated based on Mindlin plate theory to take account of shear deformation. To overcome the overestimation of shear stiffness in thin Mindlin plate element, especially in the lower order element, five nonconforming displacement modes are added to the original displacement fields. The proposed nonconforming element does not possess spurious zero-energy mode and does not show shear locking phenomena in very thin plate even for distorted mesh shapes. It was recognized from benchmark numerical tests that the displacement converges to the analytical solutions rapidly and the stress distributions are very smooth. The element also provides good results for the case of high aspect ratio.

• Journal of the Korean Geotechnical Society
• /
• v.36 no.12
• /
• pp.27-34
• /
• 2020

Use of axisymmetric shell element for the structure increases the efficiency and accuracy in finite element analysis of the interaction between the ground and the structure. This paper derived the force balance equation and the moment balance equation for an axisymmetric shell element based on Kirchhoff's theory. The governing equation for the axial deformation used the isoparametric shape function in the Galerkin formulation, and the governing equation for the shell bending used the higher-order shape function. The developed axisymmetric shell element was combined with Geo-COUS, a geotechnical finite element program for the coupled analysis with the ground. The accuracy of the developed element was confirmed through the example analyses of the circular plate and the liquid storage tank. And the energy balance equation for the axisymmetric shell element is presented.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2010.04a
• /
• pp.211-214
• /
• 2010

플랫 플레이트 시스템은 기둥 주위의 국부적인 응력집중 현상으로 인한 뚫림전단 파괴에 대해 취약하다. 따라서 유한요소해석을 통해 이러한 플랫 플레이트 시스템의 뚫림전단 성능을 평가하고자 한다. 슬래브의 전단을 고려하기 위하여 Reissner-Mindlin 가정을 바탕으로 한 등매개변수 감절점 쉘 요소를 적용하였다. 콘크리트의 재료적 비선형 거동을 고려하기 위해 압축거동은 수정압축장 이론을 적용하였으며 인장강성효과 또한 콘크리트 재료모델에 반영하였다. 기존 실험결과와의 비교를 통해 타당성을 검증하고자 하였다. 비교 결과, 약 16%의 오차율을 보였으며 보강비가 낮은 실험체에 비해 보강비가 높은 실험체가 실험결과에 가까운 값을 예측하는 것으로 나타났다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.12 no.4
• /
• pp.591-597
• /
• 1999

본 연구에서는 비 균질 퇴적층으로 인한 지진파의 증폭에 대한 경계/유한요소 해석을 수행하였다. 수치해석을 위해, 비 균질 퇴적층은 8절점 등 매개 변수 유한요소 사용하여 모델링하였고, 그 주위의 균질 반무한 지반은 3절점 등매개 변수 경계요소를 사용하여 모델링하였다. 경계요소와 유한요소의 접촉면에서, 표면 력의 평형조건과 변위의 적합 조건에 의해 두 개의 요소를 결합하는 알고리듬을 개발하였다. 수치해석의 영향인자로서 SH파, P파와 SV파의 입사각, 무 차원 진동수 그리고 반무한 지반과 퇴적층사이의 전단 파 속도 비와 질량밀도 비를 고려하였다.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.32 no.3
• /
• pp.58-65
• /
• 2004

Thermal and structural finite element analyses were performed for the turbopump turbine bladed disk model with shroud of a liquid rocket engine. The only 1/80 part model was analyzed which consists of 3D eight node isoparametric solid elements. The applied loading history consists of a startup condition with a thermal spike and a steady state. Heat transfer coefficient on the blade was predicted using the commercial Navier-Stokes solver, Fluent. Transient thermal responses during startup and steady states were calculated using a 3D finite element code developed. Maximum stress and shroud tip displacement under the influence of centrifugal and thermal loading were also determined.