• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.11
• /
• pp.2922-2931
• /
• 1994

We propose a modified 6-node element, where the sampling point of Gauss quadrature moved in the thickness direction. The modified 6-node element has been applied to static problems and forced motion analyses. In this study, this method is extended to the finite element analysis of the natural frequencies of two dimensional problems. We also propose a modified 16-node element for three dimensional problems, which behaves much like a 20-node element with smaller degree of freedom. The modified 6-node and 16-node elements have been applied to the modal analyses of beams and plates, respectively. The results agree well with the results of the 8-node or 20-node element models.

• Computational Structural Engineering
• /
• v.8 no.4
• /
• pp.87-97
• /
• 1995

For the same configuration of two-dimensional finite element models, 6-node element exhibits stiffer bending stiffness than 8-node element. This is true in the relation between 16-node element and 20-node element for three-dimensional model. This stiffening phenomenon comes from the elimination of several mid nodes from full-node elements. Therefore, this may be called 'relative stiffness stiffening phenomenon'. It seems that there are a couple of ways to correct the stiffening effect, however, we could find only one effective method-the method of modification of Gauss sampling points-which passes the patch test and does not alter other kinds of stiffness, such as extensional stiffness. The quantity of modification is a function of Poisson's ratios of the constituent materials. We could obtain two modification equations, one for plane stress case and the other for plane strain case. This method can be extended to 3-dimensional solid elements. Except the exact plane strain cases, most 3-dimensional plates could be modeled successfully with 16-node element modified by the equation for the plane stress case. The effectiveness of the modification method is checked by applying it to several examples with excellent improvements. In numerical examples, beams with various boundary conditions are subjected to static and time-dependent loads. Free and forced motion analyses of beams and plates are also tested. The beam and plate may be composed of isotropic multilayers as well as a single layer.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2009.04a
• /
• pp.304-307
• /
• 2009

본 논문에서는 3차원 모델링을 이용하여 원공을 갖는 알루미늄 판의 패치 보강효과에 대해 알아보고자 한다. 구조물의 노후화로 인해 높은 응력을 받는 부재의 응력 특이점에서 내구력이 급격하게 저하되거나 때로는 부재의 정적파괴를 유발시키는 원인을 제공한다. 이로 인해 과거에는 손상된 모재에 보강 재료를 연결시키기 위하여 리벳 또는 볼트와 같은 기계적 연결을 통해 보강하였으나 최근에는 접착패치보강 기법이 그 주류를 이루고 있다. 패치 보강시 일면 패치 보강으로 인하여 면외 휨 효과가 발생된다. 판의 두께 방향에 따른 응력집중계수를 별도로 분석하였다. 기존의 3차원 솔리드 요소는 해의 정확성은 뛰어난 반면에 상당한 컴퓨터 시간을 요구하는 단점을 가지고 있다. 이러한 문제를 극복하기 위해서, 본 논문에서는 각 층의 변위장을 2차원 형상함수와 1차원 형상함수의 조합으로 구성하여, 면내거동에 대한 p-세분화와 면외거동에 대한 p-세분화를 분리시키는 방식을 취한다. 또한, 에너지 함수의 적분시 Gauss-Lobatto 적분법을 사용하여 절점의 위치에서의 응력점을 구하는 경우, 외삽과정을 계산하는 단계를 생략하면서도, 해의 정확성 측면에서는 거의 차이가 없기 때문에 좀 더 효율적인 수치적분이 될 수 있다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.21 no.3
• /
• pp.295-304
• /
• 2008

This paper comprises two specific objectives. The first is to examine the applicability of Ordinary Kriging interpolation(OK) to finite element method that is based on variogram modeling in conjunction with different allowable limits of separation distance. The second is to investigate the accuracy according to theoretical variograms such as polynomial, Gauss, and spherical models. For this purpose, the weighted least square method is applied to obtain the estimated new stress field from the stress data at the Gauss points. The weight factor is determined by experimental and theoretical variograms for interpolation of stress data apart from the conventional interpolation methods that use an equal weight factor. The validity of the proposed approach has been tested by analyzing two numerical examples. It is noted that the numerical results by Gauss model using 25% allowable limit of separation distance show an excellent agreement with theoretical solutions in literature.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.12 no.1
• /
• pp.92-98
• /
• 2001

When analyzing the scatting problem of multi-layered planar structures using closed-form Green's function, one of the main difficulties is that the numerical integrations for the evaluation of diagonal matrix elements converge slowly and are not so stable. Accordingly, even when the integration fur the singularity of type $e^{-jkr}/{\gamma}$, corresponding to the source dipole itself, is performed using such a method, this difficulty persists in the integration corresponding to the finite number of complex images. In order to resolve this difficulty, a new technique based upon the Gaussian quadrature in polar coordinates for the evaluation of the two-dimensional generalized exponential integral is presented. Stability of the algorithm and convergence is discussed. Performance is demonstrated for the example of a microstrip patch antenna.

• Computational Structural Engineering
• /
• v.10 no.4
• /
• pp.117-130
• /
• 1997

In this paper, a modified 10-node equivalent solid element(MQM10 element), which has smallest degrees of freedom among 3-dimensional solid elements accounting bending deformation as well as extensional and shear deformations of isotropic plates, is proposed. The proposed MQM10 element exhibits stiffer bending stiffness due to the reduction of degrees of freedom from 20-node element or Q11 element. As an effective way to correct the relative stiffness stiffening phenomenon, the modification equation of Gauss sampling points is proposed. The quantity of modification is a function of Poisson's ratio. The effectiveness of MQM10 element is tested by applying it to several examples. It is noted that the results of static and free vibration analysis of isotropic plates using MQM10 elements show a good agreement with those using 20-node element.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.4
• /
• pp.1320-1332
• /
• 1996

In the present study, we consider modified 5-node equivalent solid element which has smallest degree of freedom among 2-dimensional solid elements accounting bending deformation as well as extensional and shear deformations, We shall investigate static and dynamic characteristics of this element, which is very effective in thin beam, thick beam, large displacement problems, beam of variable thickness, and asymmetrically stepped beam, etc., as well as relatively simple problems of beam. The degree of freedom of this element is 10, which is smaller than 18 of 9-node element, 16 of 8-node elemtns, 12 of modified 6-node element and Q6 element. Therefore, this element is expected to broaden the effective range of application of the solid elements in the bending problems further.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.15 no.3
• /
• pp.491-499
• /
• 2002

A p-version finite element model based on degenerate shell element is proposed tot the analysis of orthotropic laminated plates. In the nonlinear formulation of the model, the total Lagrangian formulation is adopted with large deflection and moderate rotation being accounted tot in the sense of yon Karman hypothesis. The material model is based on the Huber-Mises yield criterion and Prandtl-Reuss flow rule in accordance with the theory of strain hardening yield function, which is generalized lot anisotropic materials by introducing the parameters of anisotropy. The model is also based on extension of equivalent-single layer laminate theory(ESL theory) with shear deformation, leading to continuous shear strain at the interface of two layers. The integrals of Legendre polynomials are used for shape functions with p-level varying from 1 to 10. Gauss-Lobatto numerical quadrature is used to calculate the stresses at the nodal points instead of Gauss points. The validity of the proposed P-version finite element model is demonstrated through several comparative points of iew in terms of ultimate load, convergence characteristics, nonlinear effect, and shape of plastic tone.

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

• Journal of Ocean Engineering and Technology
• /
• v.12 no.2 s.28
• /
• pp.1-21
• /
• 1998

In this paper, a non-conforming 3-D 8-node solid element(MQM10) has beets applied to the analyses of static and dynamic bending problems of laminated composite plates The QM10 element exhibits stiffer bending stiffness which is caused by the reduction of degree of freedom from Q11 element. As an effective way to correct the relative stiffness stiffening phenomenon the modification of Gauss sampling points for composite plates is proposed. The quantity of modification is a function of material properties. Also, another two modified equations are obtained, one is modification for stress, and the other is modification of coefficient of shear modulus in free vibration. It is noted that MQM10 element can analyse the static and free vibration problems of various 3-dimensional composite plates composed of unidirectional laminae, woven laminae or braided laminae. The results of MQM10 element are in good agreement with those of 20-node element.