• Geotechnical Engineering
• /
• v.14 no.4
• /
• pp.129-140
• /
• 1998

Elasto-plastic finite element analysis of geotechnical boundary value problems necessitate the stress integration for the known strain increments. For the elasto-plastic constitutive model, the stress integration is generally achieved by numerical schemes, because analytical integration is impossible for general strain path. In this case, the accuracy of numerical stress integration has an important role on the overall accuracy of nonlinear finite element solution. In this study, the Sloan's substepping method which is one of explicit integration methods has been adopted and iris applicability has been checked. The unstability and inaccuracy of ifs results initiated from initial stress level were revealed. So. a new modified numerical integration method which employs the basic concept of modified Euler scheme for error control is proposed and accuracy and stability of the solutions are confirmed by triaxial test simulation.

• Fire Science and Engineering
• /
• v.19 no.4 s.60
• /
• pp.18-25
• /
• 2005

In this study, numerical analysis of two-dimensional unsteady natural convection of air in a square enclosure heated from below, was performed as a basic research of fire science. SIMPLE algorithm was used to the pressure term of momentum equations in the numerical analysis. The numerical analysis were studied for the two model cases and two heat conditions, respectively, which are different with insulation of enclosures and position of heat applied. Also, the ceiling temperatures of enclosure were measured to compare the accuracy of numerical analysis, and it is found that the temperature predicted by numerical analysis were agreed well with the measurements. Streamline and isotherm of the each model case were acquired for each time step.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.15 no.1
• /
• pp.1-8
• /
• 1990

In this paper, the transient electromagneti wave scattering at dielectric cylinder is studied by new numerical analysis method. Basic formulation of boundary integral equation (BIE) for numerical method is started weighted residual technique. BIE is made to two simultaneous equation at surface inner and outside point of dielectric cylinder in extended boundary condition (EBC) and surface boundary condition (SBC). Numerical method is used Boundary element method (BEM) that is two form, one is direct method and the other is indirect method, so that this method that transformes operator inversion martics is used numerical analysis. A good agreement of this numerical solution and the other results is obtained.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.11 no.2
• /
• pp.186-196
• /
• 2000

In this paper, modified square loop frequency selective surface for dualband communication antenna systems is proposed, and the scattering characteristics is discussed. The analysis for the problem of scattering by periodic structures with a dielectric slab is formulated using the spectral-domain immittance approach and Floquet's theorem. The method of moments which uses rooftop subdomain basis function is employed to solve the equations. Numerical results include the comparison between the transmission characteristic of general square loop and that of modified square loop. Also, the transmission characteristics of modified square loop for arbitrary incident angle and polarization is presented. To verify analysis results, modified square loop frequency selective surface was fabricated and the calculated results were compared with the measured results. The measured results showed good agreement with the calculated results.

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

본 논문은 기하학적 비선형성을 가진 보존적 단일 하중 매개변수의 탄성 상태 공간구조의 분기이론에 관한 수치 해석적 기본 방법 및 경로 추적, pin-pointing, 경로 전환을 기술하고 있다. 비선형 탄성 불안정 상태는 극한점과 분기점으로 분류될 수 있으며, 평형경로상의 평형점의 계산 및 평형경로상의 특이점을 찾기 위한 pin-pointing 반복계산을 수행하는 일반적인 비선형 수치해석법으로 극한점을 계산할 수 있다. 그러나 분기좌굴 해석을 위해서는 좌굴 후 분기경로의 추적을 위한 분기경로 전환 알고리즘이 추가적으로 필요하다. 본문에서는 에너지이론에 기초한 일반 탄성안정이론을 소개하고, 평형경로 추적, 분기 좌굴점을 찾기 위한 직접법과 분기경로 전환에 관한 이론을 전개한다. 분기좌굴 해석예제로 트러스로 이루어진 스타돔, 핀지지의 평면아치, 평면프레임, 3차원 공간프레임의 분기좌굴 해석을 수행하여 본문에서 제시한 수치해석법의 정확성 및 실용성을 검증한다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.5C
• /
• pp.239-250
• /
• 2009

Tunnelling-induced settlements and lateral ground movements have been investigated by numerical parametric studies considering ground condition, excavation characteristics, and ground layers. Before the numerical study the existing methods of ground movement estimation have been collected and analysed to have some information of ground movements and to improve them providing a fundamental material for the numerical study. Numerical model simulation has been performed of a physical model test of tunnel excavation in which the ground movements were measured reliably and the results have been used to determine the numerical approach and the appropriate soil constitutive mode. With this procedure done, the results of numerical parametric studies have been put together to analyze and understand tunnelling-induced settlements and lateral ground movements.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.10 no.1
• /
• pp.91-104
• /
• 2008

In tunnels, safety factor concept has been suggested to estimate their stability quantitatively. It is merely limited in the framework of mechanical analysis. However safety factor concept has not been applied in hydro-mechanical coupled analyses due to their modelling complexity. Recently studies on this topic are being actively made. In this study, induced drainage modelling methods for hydro-mechanical coupled analyses are compared and analyzed to estimate safety factor of a subsea tunnel exactly. To this end, methods both controlling hydraulic characteristic of shotcrete and using a drainage well are considered. Sensitivity analysis were carried out on rock class, thickness of shotcrete, and hydraulic properties of rock mass. As the results of this study, it turned out that the induced drainage modelling using a drainage well would give more reliable results than that of controlling hydraulic characteristic of shotcrete in estimating tunnel stability in hydro-mechanical coupled analyses.

• Proceedings of the Korea Society for Energy Engineering kosee Conference
• /
• 1994.05a
• /
• pp.28-38
• /
• 1994

분류층 가스화기 설계를 위한 일차연구로서 가스화기 이차공기 주입방법에 따른 비반응 난류장에 대한 수치해석을 검사체적에 기초한 유한차분방법을 이용하여 수행하였다. 압력과 속도의 연계문제는 SIMPLEC 알고리즘을, 레이놀즈 전단력은 k-$\varepsilon$ 난류모델을 사용하였다. 입자궤적 계산은 공기역학적 항력만을 고려하였으며 비선형적인 공기저항력에 의한 난류변동상관모델은 고려치 않았다. 이차공기 주입방법(parallel injection과 nonparallel 3$0^{\circ}C$ injection)에 따른 수치해석을 수행하여 Ar tracer의 질량분율에 대한 실험자료와 비교하여 만족할 만한 결과를 얻었으며 이차공기의 주입각 및 기타 제반변수에 따른 유동장 변화를 분석하였다.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1992.10a
• /
• pp.130-135
• /
• 1992

복잡한 구조물의 동특성 해석은 수치 해석적인 방법과 실험적 방법 모두 단 순 구조물의 동특성 해석에 비해 정확도가 떨어지는 반면 계산시간과 노력 은 크게 증가하게 된다. 이 경우 구조물 전체를 여러개의 간단한 부분 구조 계로 부할하고, 각 부분 구조계에 대해서 해석후 그 결과들을 적절한 결합 조건하에서 다시 조합하여 전체 구조계에 대해 동특성 해석을 수행하거나, 고감쇠 처리된 구조물과 같이 고전적 이론 해석기법의 적용이 어려운 경우 실험적인 해석방법과 이론적인 해서방법을 혼합 사용할 수 있는 부분 구조 합성법(Substructural Synthesis Method)을 사용하는 것이 효과적이다. 부분 구조합성법은 1) 응답특성을 이용한 방법 2) 모우드 특성을 이용한 방법 3) 특성행렬을 이용한 방법 등이 있으며 본 연구에서는 부분 구조계 응답함수 로부터 직접 특성행렬을 산출하는 방법을 이용하여 전체 구조계의 동특성을 해석할 수 있는 부분 구조 합성법을 제시하였다.

• Proceedings of the Korea Information Processing Society Conference
• /
• 2004.05a
• /
• pp.989-992
• /
• 2004

수학적 모델을 컴퓨터 상에 실현시키는데 있어 보다 효율적인 알고리즘을 구현하고 개발하는 것이 수치해석 연구의 궁극적인 목표이다. 일반적으로 수치해석을 사용할 때 컴퓨터 상에서 구한 계산 결과, 즉 근사 값과 수학적으로 구한 값인 참값은 정확하게 같지 않다. 따라서 근사 값이 얼마나 참값에 가까운가에 따라 알고리즘의 효율성을 평가하는 오차평가는 수치해석의 가장 중요한 과제라 할 수 있다. 대부분의 경우 오차평가에 있어 오차의 한계를 이용하지만 주어진 문제의 참값을 모르기 때문에 정확한 오차평가를 할 수 없다. 본 논문에서는 수치등각사상을 구하기 위한 해법중 하나인 Wegmann해법에 있어 몇 가지 수학적 이론에 근거하여 참값을 모르더라도 오차평가를 할 수 있는 방법을 제안하고 수치실험을 통해 그 유효성을 입증한다.