• Journal of the KSME
• /
• v.23 no.3
• /
• pp.200-206
• /
• 1983

FAM의 기본적인 구상은 해석 하고자하는 선형 또는 비선형 편미분 방정식을 국부적으로 해석 적인 해를 구하여 이용하자는 것이다. 그러기 위하여 유한차분법(FDM)과 유한변분법(FEM)에 서와 같이 전체유동장을 작은 요소로 나누고 그 요소 내에서 국부해를 구한 다음 이들 요소를 중첩시킴으로써 각 요소의 미지수에 대한 대수식을 얻어서 수치해를 구하자는 것이다. 그러나 FDM에서와 같이 국부요소에서 미분항을 구하지 않고, FEM 에서와 같이 요소에서 형상함수를 도입하지 않는 상태에서 해석적인 해를 구하고 있기 때문에 수치해석에서 얻어지는 미분양들은 비교적 정확하게 구해진다. 따라서 Navier-Stokes 방정식이나 에너지 방정식에서 최고차항이 작은 파라메타, 즉 레이놀즈수나 피크리수의 역수로 곱하여서 있는 경우에도 안정된 해를 구할 수 있다고 알려져 있다. 요소자체의 계수를 구하는 데는 계산시간이 많이 소요되지만 수치해석 상의 안정성이나 수렴성이 좋기 때문에 전체계산시간은 오히려 적게 걸리는 경우도 있다고 한다.

• 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 Korea Institute of Military Science and Technology
• /
• v.1 no.1
• /
• pp.128-144
• /
• 1998

The BEM for linear elastic stress analysis is applied to the highly rotating axisymmetric body problem which also involves the thermoelastic effects due to steady-state thermal conduction. The axisymmetric BEM formulation is briefly summarized and an alternative approach for transforming the volume integrals associated with such body force kernels into equivalent boundary integrals is described in a way of using the concept of inner product and vector identity. A discretization scheme for higher order BE is outlined for numerical treatment of the resulting boundary integral equations, and it is consequently illustrated by determining the stress distributions of the turbine rotor disk of the small turbojet engine(ADD 500) for which a FEM stress solution has been furnished by author.

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

• Journal of Ocean Engineering and Technology
• /
• v.12 no.1
• /
• pp.99-106
• /
• 1998

3차원 파수조에서 완전 비선형파를 시뮬레이션하기 위하여 우선 랜킨 소스를 기저로한 적분방정식을 고차경계요소법을 이용하여 이산화하였다. 그리고 방사경계조건은 파흡수 비치와 포텐셜 스트레칭 기법을 이용하여 모델링하였으며, 비선형 자유표면과 경계조건식은 고차 예측 및 보정 기법을 이용하여 시간 적분하였다. 파흡수 비치는 파의 진행방향 특성에 따라 수조내에 다양하게 배치할 수 있으며 비칭서 흡수가 덜된 파는 수조의 길이 방향 끝단에서 포텐셜 스트레칭 기법에 의하여 반사없이 진행하도록 하였다. 수치실험 결과 일-에너지 보존법칙과 모멘텀-임펄스 보존 법칙이 만족됨으로써 본 수치기법의 효용성이 검증되었다.

• Journal of IKEEE
• /
• v.19 no.4
• /
• pp.520-525
• /
• 2015

This paper investigates the consensus problem for high-order integrators with an arbitrary large communication delay. In order to solve this problem, new consensus controller with an additional design parameter that can eliminate the effect of a communication delay on the consensus problem is proposed. Also, it is proved that the proposed consensus controller can always solve the consensus problem of high-order integrators even in the presence of an arbitrarily large communication delay. Finally, an illustrative example is given in order to show the effectiveness of our design method.

• Journal of Korean Society on Water Environment
• /
• v.23 no.5
• /
• pp.683-688
• /
• 2007

Water flow into unsaturated soils is most often modeled by Richards equation consisting of the mass conservation law and Darcy's law. Three standard forms of Richards equation are presented as the head (${\Psi}$)-based form, the moisture content (${\theta}$) based form, and the mixed form. Numerical solutions of these partial differential equations with highly nonlinear terms can cause poor results along with significant mass balance errors. The numerical solution based on the mixed form of Richards equation is known that the mass is perfectly conserved without any additional computational efforts. The aim of this study is to develop fully implicit numerical scheme of Richards equation for one-dimensional vertical unsaturated flow in homogeneous soils using the finite difference approximation, and then to perform sensitivity analysis of infiltration to the variations in the unsaturated soil properties and to different soil types.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.36 no.9
• /
• pp.907-913
• /
• 2012

High-order perturbation solutions have been obtained for the simple physical model describing the liquid pool spreading with a continuous spill, and these are shown to improve over first-order perturbation solutions. The non-dimensional governing equations for the model are derived to obtain more general solutions. Non-dimensional parameters are sought as the governing parameters for the non-dimensional equations, and the non-dimensional evaporation rate is used as the perturbation parameter. The results show that the high-order solutions exhibit an improvement over the first-order solutions with respect to the pool volume as well as the spreading radius. In addition, as the order of the perturbation solutions increases, the difference between the numerical solutions and the perturbation solutions is significantly reduced. Finally, it is revealed that the third-order solutions are reasonable because they almost agree with the numerical solutions.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.2 no.3
• /
• pp.75-86
• /
• 1982

The degeneration of two classes of thick beam elements has been conducted, one (DB6) based on the conventional Timoshenko beam assumptions whereas the other (DB7) based on the assumed cubic axial displacement profile. While an adjustable shear correction factor is required for the DB 6 element to compensate for the unrealistic distribution of shear strain across the thickness, the DB 7 element assumes the more realistic quadratic profile of shear strain at the outset. With the plane-stress continuum solution as reference, solutions obtained by these two element models are compared with the analytical Timoshenko solution, the analytical thin beam solution and several available solutions of other existing beam elements. The result indicates that the performance of the higher order beam element DB 7 is consistently superior to any others. This is true for the whole range of aspect ratios of beam, in both static and free vibration analyses.

• Journal of Korean Society of Steel Construction
• /
• v.23 no.4
• /
• pp.475-481
• /
• 2011

In this study, a simplified method of improving not only transverse shear stress but also shear strain based on the first-order shear deformation theory was developed. Unlike many established methods, such as the higher-order shear deformation and layerwise theories, this method can easily apply to finite elements as only $C^0$ continuity is necessary and the formulation of equations is very simple. The basic concept in this method, however, must be corrected:the distribution of the transverse shear stresses and shear strains through the thickness from the formulation based on the higher-order shear deformation theory. Therefore, the shear correction factors are no longer required, based on the first-order shear deformation theory. Numerical analyses were conducted to verify the validity of the proposed formulations. The solutions based on the simplified method were in very good agreement with the results considering the higher-order shear deformation theory.