• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.20 no.12
• /
• pp.614-621
• /
• 2019

A generalized laminated composite beam element is presented for the flexural and buckling analysis of laminated composite beams with double and single symmetric cross-sections. Based on shear-deformable beam theory, the present beam model accounts for transverse shear and warping deformations, as well as all coupling terms caused by material anisotropy. The plane stress and plane strain assumptions were used along with the cross-sectional stiffness coefficients obtained from the analytical technique for different cross-sections. Two types of one-dimensional beam elements with seven degrees-of-freedom per node, including warping deformation, i.e., three-node and four-node elements, are proposed to predict the flexural behavior of symmetric or anti-symmetric laminated beams. To alleviate the shear-locking problem, a reduced integration scheme was employed in this study. The buckling load of laminated composite beams under axial compression was then calculated using the derived geometric block stiffness. To demonstrate the accuracy and efficiency of the proposed beam elements, the results based on three-node beam element were compared with those of other researchers and ABAQUS finite elements. The effects of coupling and shear deformation, support conditions, load forms, span-to-height ratio, lamination architecture on the flexural response, and buckling load of composite beams were investigated. The convergence of two different beam elements was also performed.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.30 no.12 s.255
• /
• pp.1626-1633
• /
• 2006

In the present paper, in the context of the meshless interpolation of a moving least squares (MLS) type, a novel method which uses primary and secondary nodes in the domain and on the global boundary is introduced, in order to improve the accuracy of solution. The secondary nodes can be placed at any location where one needs to obtain a better resolution. The support domains for the shape functions in the MLS approximation are defined from the primary nodes, and the secondary nodes use the same support domains. The shape functions based on the MLS approximation, in an integration domain, have a single type of a rational function, which reduces the difficulty of numerical integration to evaluate the weak form. The present method is very useful in an adaptive calculation, because the secondary nodes can be easily added and moved without an additional mesh. Several numerical examples are presented to illustrate the effectiveness of the present method.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.7
• /
• pp.1590-1596
• /
• 1995

The energy release rates for a dynamically growing crack in orthotropic materials are expressed explicitly in terms of dynamic stress intensity factors. The stress functions suitable for the problem are found and the evaluation of the J-integral for the theoretical singular crack tip fields yields energy release rates. The present results are simpler than the existing ones and can be reduced to the well known solutions in special cases. Examples of extracting stress intensity factors from the finite element solution using the present results are given for the dynamically growing crack problem of orthotropic materials.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.8 no.2
• /
• pp.25-32
• /
• 1988

A formulation of an isoparametric quadrilateral higher-order plate bending finite element is presented. The 8-noded 28-d.o.f. plate element has been degenerated from the three-dimensional continuum by introducing the plate assumptions and considering higher-order in-plane displacement profile. The element characteristics have been derived by the Galerkin's weighted residual method and computed by using the selective reduced integration technique to avoid shear-locking phenomenon. Several numerical examples are given to demonstrate the accuracy and versatility of the proposed quadrilateral higher-order plate bending element over the other existing plate finite elements in both static and dynamic analyses.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2001.04a
• /
• pp.53-60
• /
• 2001

Meshfree methods have been attracting issue as computational methods during past a few years. Nowadays, various meshfree methods such as EFGM, RKPM h-p cloud method and etc. were developed and applied in engineering problems. But, most of them were not truly meshless method because background mesh of cell was required for the spatial integration of a weak form. A nodal integration is required for truly meshless methods but it is known that this method gives a little unstable and incorrect solutions. In this paper, an improvement scheme of the existed nodal integration which the weak form can be simply integrated without any stabilization term is proposed. Numerical tests show that the proposed method is more convenient and gives more correct solutions than the previous method.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.12
• /
• pp.3007-3019
• /
• 1993

For the same configuration, the stiffness of 6-node two dimensional isoparametric is stiffer than that of 8-node two dimensional isoparametric element. This phenomenon may be called 'Relative Stiffness Stiffening Phenomenon.' In this paper, the relative stiffness stiffening phenomenon was studied, and could be corrected by modifying the position of Gauss integral points used in the numerical integration of the stiffness matrix. For the same deformation (bending) energy of 6-node and 8-node two dimensional isoparametric elements, Gauss integral points of 6-node element have to move closer, in comparison with those of 8-node element, in the case of numerical integration along the thickness direction.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.14 no.4
• /
• pp.481-494
• /
• 2001

An enhanced degenerated shell finite element (FE), which has been developed for inelastic analysis of reinforced concrete structures is described in this paper. Generally, Reissner-Mindlin (RM) assumptions are adopted to develop the degenerated shell FE so that transverse shear deformation effects is considered. However, it is found that there are serious defects such as locking phenomena in RM degenerated shell FE since the stiffness matrix has been overestimated in some situations. As remedies of locking phenomena, reduced integration, incompatible mode and assumed strain method have been used. Especially, the assumed strain method has been successfully used in many FEs. But contrarily, there is a few investigation on the performance of the assumed strains in the inelastic analysis of concrete structures. Therefore, shell formulation is provided in this paper with emphasis on the terms related to the stiffness matrix based on assumed strain method and microscopic concrete material model. Finally, the performance of the present shell element is tested and demonstrated with several numerical examples. From the numerical tests, the present result shows a good agreement with experimental data or other numerical results.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2006.05a
• /
• pp.588-593
• /
• 2006

자연하천은 일반적으로 만곡수로나 사행수로 형태를 보이고 있으며, 직선수로에서와 달리 원심력에 기인한 이차류 영향을 받게 된다. 이차류에 의해서 수면에서는 만곡부 바깥쪽으로, 하상에서는 만곡부 안쪽으로의 흐름특성을 보이게 된다. 만곡부 안쪽으로 가해지는 하상 전단응력에 기인하여 하상에서의 입자가 만곡부 안쪽으로 이송되며, 결과적으로 만곡부 안쪽에는 점사주가, 바깥쪽에는 소(pool)가 생성된다. 또한 지형경사의 생성으로 입자에 가해지는 중력효과도 변화된다. 따라서 이와 같은 자연하천의 흐름과 하상변동을 수치모의 하기 위해서는 만곡부 이차류 특성을 고려한 모형이 필요하다. 본 연구에서는 수심 적분된 흐름방정식과 하상토 보존방정식 (Exner equation)을 이용한 하상변동을 위한 비연계 수치모형을 위해서 하상토 보존방정식의 유한요소 알고리즘을 개발하였다. 하상토 보존방정식은 흐름 특성에 따른 평형 유사량의 공간변화율을 이용하여 일정 기간 동안의 하상 변화량을 계산한다. 이 때 이차류에 의한 하상 전단응력의 편각 및 지형경사 변화에 따른 실제 입자의 이송방향을 보정하여 평형 유사량이 계산된다. 이러한 보정식을 적용시키기 위해서는 유속성분의 공간변화량 및 지형경사의 공간성분이 필요하다. 유한요소법은 연속성 변수를 이산화시켜 근사해를 구하는 수치기법의 일종이기 때문에, 요소망이 불규칙적으로 구성되었을 경우 임의의 절점에서 연속성을 지닌 변수의 공간변화율을 계산하는데 어려움이 있다. 따라서 본 연구에서는 평형 유사량 계산 시에 절점이 아닌 요소 내부에서 평형 유사량을 계산하는, 하상토 보존방정식의 새로운 유한요소 알고리즘을 개발하고, 새로운 알고리즘을 적용시킨 수치모형의 검증을 행하였다. 경계조건 알고리즘의 검증으로 위해서 Soni 등 (1980)이 행한 상류 유입 유사량에 따른 하상변동을 수치 모의하고 실험치와 비교하였으며, Sutmuller와 Glerum (1980)이 수행한 만곡수로에서의 하상변동을 모의하고 실험과 비교하였다. 새로운 알고리즘을 적용시킨 하상토 보존방정식의 유한요소 수치모형의 결과는 매우 안정적이며, 실험과 매우 유사한 결과를 얻을 수 있었다. 본 수치모델은 현재 균일한 입자의 하상토만을 고려하므로, 입자분급이나 하상 장갑화 현상 등은 무시한다.

• Journal of Korean Association for Spatial Structures
• /
• v.8 no.4
• /
• pp.81-90
• /
• 2008

This paper presents a finite element formulation for the three-dimensional hierarchical solid element using Integrals of Legendre polynomials. The proposed hexahedral solid element is composed of four different modes including vertex, edge, face, and internal mode, respectively. The eigenvalue and patch test have been carried out to confirm the zero-energy mode and constant strain condition. In addition to these, a posteriori error estimation has been studied for the p-adaptive finite element analysis that is based on a smoothing technique to compute a post-processed solution from the finite element solution. The uniform p-refinement and non-uniform p-refinement are compared in terms of convergence rate as the number of degree of freedom is increased. The simple cantilever beam is tested to show the performance of the proposed solid element.

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