• Journal of the Earthquake Engineering Society of Korea
• /
• v.15 no.5
• /
• pp.11-24
• /
• 2011

A real-time hybrid test of a 3 story-3 bay reinforced concrete frame which is divided into numerical and physical substructure models under uniaxial earthquake excitation was run using a fixed iteration implicit HHT time integration method. The first story inner non-ductile column was selected as the physical substructure model, and uniaxial earthquake excitation was applied to the numerical model until the specimen failed due to severe damage. A finite-element analysis program, Mercury, was newly developed and optimized for a real-time hybrid test. The drift ratio based on the top horizontal displacement of the physical substructure model was compared with the result of a numerical simulation by OpenSees and the result of a shaking table test. The experiment in this paper is one of the most complex real-time hybrid tests, and the description of the hardware, algorithm and models is presented in detail. If there is an improvement in the numerical model, the evaluation of the tangent stiffness matrix of the physical substructure model in the finite element analysis program and better software to reduce the computational time of the element state determination for the force-based beam-column element, then the comparison with the results of the real-time hybrid test and the shaking table test deserves to make a recommendation. In addition, for the goal of a "Numerical simulation of the complex structures under dynamic loading", the real time hybrid test has enough merit as an alternative to dynamic experiments of large and complex structures.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.3
• /
• pp.25-34
• /
• 1993

The long-term behavior, such as in excavation problems of weak medium, can be dealt with by the elasto-viscoplasticity models. In this paper, a combined formulation of elasto-viscoplasticity using boundary elements and finite elements without using internal cells is presented. The domain integral introduced due to the viscoplastic stresses is transformed into a boundary integral applying direct integration in cylindrical coordinates. The results of the developed boundary element analysis are compared with those from the explicit solution and from the finite element analysis. It is observed that the boundary element analysis without internal cells results in some error because of its deficiency in handling the nonlinearity in local stress concentration. Therefore, a coupled analysis of boundary elements and finite elements, in which finite elements are used in the area of stress concentration, is developed. The coupled method is applied to a time dependent inelastic problem with semi-infinite boundaries. It results in reasonable solution compared with other methods where relatively higher degree of freedoms are employed. Thus, it is concluded that the combined analysis may be used for such problems in the effective manner.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2000.04a
• /
• pp.16-16
• /
• 2000

충격파를 포함하는 초음속 유동을 해석하는 수치해법 중에서 많이 사용되어진 것은 엄밀 및 근사 리만 해법과 플럭스 분할 기법들로서 이들은 Euler 방정식에 기반을 두고 선형 또는 비선형파의 상호작용을 풍상 차분법으로 기술하는 방법들이다. 이러한 수치기법들은 과거 광범위하게 사용되어 왔으나 최근 여러 가지 단점이 발견되었다. 이와 같은 문제점을 극복하고자 입자의 통계적인 운동을 기술하는 기체 운동론에 근거하여 BGK 수치기법이 제시되었다. 이는 비충돌 볼츠만 방정식으로부터 입자의 수준에서 플럭스 분할 기법 형태의 풍상차분법을 구현하는 것으로 볼츠만 방정식의 충돌항을 BGK 모델로 대치하고 이것의 적분해로부터 수치 플럭스를 구한다. 이 수치기법은 기존의 리만해법에 비하여 수치적으로나 물리적으로 매우 타당한 성질인 강건성, 정확성, 엔트로피 조건, 양수보존성 등을 가지고 있음이 밝혀졌다. 이와 같은 수치기법을 사용하여 로켓 노즐 내의 아음속, 천이음속, 초음속에서의 유동장 해석을 위한 프로그램을 작성하였다. 시간 적분에 대하여는 정상 상태의 계산을 위하여 내재적 시간 적분 방법을 사용하였으며, 공간 이산화 방법으로는 임의의 제어체적에 대하여 적분형 보존 방정식을 적용하는 유한 체적법을 사용하였다. 초음속 입구 유동과 출구에서 초음속과 저음속 유동의 두가지 경우를 고려하여 얻은 결과를 기존의 연구 결과와 비교하여 본 결과 잘 일치하였다. 입구 유동이 저음속이고 출구 유동이 초음속인 경우에 대하여도 해석결과가 실험결과와 잘 일치하였다. 상대적으로 낮은 온도, 압력 조건과 높은 온도, 압력 조건을 가지는 고체 로켓 모터 노즐 내의 유동을 해석하였다. 이들 해석 결과를 전압, 전온도로 표준화시킨 결과 서로 일치하였으며, 파라서 저온, 저압에서 얻은 결과도 표준화시킬 경우, 고온, 고압에서도 사용될 수 있음을 알 수 있었다.의 영향에 초점을 맞추었다.다고 판단되며 배기 가스 자체에 대기 공기중에 함유되어 있던 습기가 얼어붙는(Icing화) 문제가 발생하기 때문에 배기가스의 Icing을 방지하기 위하여 압축기 끝단에서 공기를 추출하여 배기부분에 송출할 필요성이 있는 것으로 판단되었다. 출구가스의 기체 유동속도가 매우 빠르므로 (100-l10m.sec) 이를 완화하기 위한 디퓨저의 설계가 요구된다고 판단된다. 또 연소기 후방에 물을 주입하는 경우 열교환기 및 기타 부분품에 발생할 수 있는 부식 및 열교환 효율 저하도 간과할 수 없는 문제로 파악되었다. 이러한 기술적 문제가 적절히 해결되는 경우 비활성 가스 제너레이터는 민수용으로는 대형 빌딩, 산림, 유조선 등의 화재에 매우 적절히 사용되어 질 수 있을 뿐 아니라 군사적으로도 군사작전 중 및 공군 기지의 화재 그리고 지하벙커에 설치되어 있는 고급 첨단 군사 장비 등의 화재 뿐 아니라 대간첩작전 등에 효과적으로 활용될 수 있을 것으로 판단된다.가 작으며, 본 연소관에 충전된 RDX/AP계 추진제의 경우 추진제의 습기투과에 의한 추진제 물성 변화는 미미한 것으로 나타났다.의 향상으로, 음성개선에 효과적이라고 사료되었으며, 이 방법이 편측 성대마비 환자의 효과적인 음성개선의 치료방법의 하나로 응용될 수 있으리라 생각된다..7%), 혈액투석, 식도부분절제술 및 위루술·위회장문합술을 시행한 경우가 각 1례(2.9%)씩이었다. 13) 심각한 합병증은 9례(26.5%)에서 보였는데 그중 식도협착증이 6례(17.6%), 급성신부전증 1례(2.9%), 종격동기흉과 폐염이 병발한 경우와 폐염이 각 1례(2.9%)였다. 14) 식도경 시행회수는 1회가 17례(54.8%), 2회가 9례(29.0%), 3회 이상이 5례(16.1%)였다.EX>$IC_{50}$/ 값이 210 $\mu\textrm{g}$<

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.26 no.5A
• /
• pp.855-862
• /
• 2001

본 논문에서는 백색잡음과 페이딩 환경하의 통신시스템의 성능평가시 백색잡음에 대해서만 모멘트 기법을 이용조건부 오류확률을 구하고 페이딩에 대해서는 해석적 접근방식을 취함으로써 전체적인 컴퓨터 수행시간을 획기적으로 개선시킬 수 있는 기법을 제안하고, 제안된 기법의 효율성을 컴퓨터 시뮬레이션 수행시간 비교를 통해 입증하였다. 현재 널리 사용되는 Monte Carlo 시뮬레이션 기법은 IS-95 순방향 채널과 같이 많은 샘플수를 필요로 하는 시스템에서 상당한 시뮬레이션 수행시간을 요구하게 된다. 이러한 문제점을 해결하기 위하여, 제안된 방식에서는 백색잡음하의 수신신호의 N차 모멘트를 측정하여 이산확률밀도함수를 구함으로써 수신신호의 확률적 특성을 모사하고, 이로부터 내삽법과 외삽법을 적용하여 조건부 누적확률분포함수를 산출함으로써 채널오류율을 구하게 된다. 그리고 페이딩 성능평가시 조건부 채널오류율을 이용하여 산출된 조건부 누적확률분포함수를 환경조건에 맞춘 후에 페이딩의 확률밀도함수와 수치 해석적으로 간단히 적분하여 성능평가를 함으로써 수행시간의 줄임과 동시에 정확한 채널오류율을 구하게 된다.

• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.586-592
• /
• 2000

터보머신 태부에 존재하는 정익 - 동익의 상호작용 유동현상을 수치모사 하는 코드를 병렬화 하였다 정익 - 동익의 상호작용을 해석하는 데에 편리하도륵 Multi-Block Grid System을 도입하여 계산영역을 형성하였고, 동익의 움직임으로 인해 발생하는 Sliding Interface부분은 Patched 알고리즘을 적용하여 해석하였다. 정익과 동익의 수를 1대 1로 단순화시켜 수치모사한 결과와 정익과 동익의 수를 실제 조건과 더 비슷하게 설정한 3대 4의 비율로 맞추어 수치모사한 결과를 비교하였다. 또한, 병렬컴퓨팅으로 인해 단축된 계산시간을 다른 연구에서의 계산시간들과 서로 비교하였다. 2차원 비정상 압축성 Navier-Stokes 방정식이 이용되었고, 난류모델링에는 K-w SST 모델링이 적응되었다. Roe의 FDS 기법을 사용하여 플럭스를 계산하였고, MUSCL 기법을 적용하여 3차의 공간정확도를 갖도록 하였다. 시간적분에는 이보성의 DP-SGS를 사용하였다. 해석결과의 분석에는 Time-averaged pressure distribution과 Pressure amplitude distribution 데이터를 사용했다.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1996.05d
• /
• pp.241-246
• /
• 1996

지진시 미끄럼 등과 같이 전형적인 비선형거동을 하는 구조물에 대한 동적해석을 수행하는 경우 먼저 입력지진의 가속도 시간이력에 대한 기준선 조정이 필요할 때가 있다. 인공적으로 작성된 지진기록의 경우 때로 가속도 시간이력을 적분하여 속도 및 변위 시간이력을 얻었을 때 증가하는 형태로 나타나 이로 인하여 비선형응답이 비정상적으로 커질 수 있기 때문이다. 본 논문에서는 바닥이 마찰거동을 하는 구조물에 대해 간단히 모델하여 이러한 예를 보였으며 또한 주로 사용되는 기준선 조정방법들의 응답영향을 비교하였다. 그 결과 입력지진의 기준선 조정을 하지 않는 것이 항상 보수적인 결과만을 보여 주지 않는다는 점과 기준선 조정의 방범에도 표준화가 필요하다는 점을 파악할 수 있었다.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.11 no.3
• /
• pp.419-428
• /
• 2000

Multiconductor transmission line(MTL) equations are solved by FDTD(Finite-Difference Time-Domain) method to predict crosstalk and fields to transmission line coupling on lossy multiconductor transmission lines terminated in arbitrary linear loads. Skin effect losses as well as dc losses are included in the analysis. In order to increase computational efficiency, the convolution integral of internal impedance of conductors and the line currents is computed by using Prony method. For boundary conditions of MTLs terminated in linear loads, state-variable formulation is adopted. The simulated results by FDTD method are compared with the measured ones obtained by using TEM cell. The predictions are in good agreement with the measurements. In addition, it has been found that skin effect losses as well as dc losses of the conductors should be included for accurate predictions on relatively high loss transmission lines such as PCB. It has also been found that dc losses and skin-effect losses affect late-time responses and early-time responses, respectively.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.25 no.2
• /
• pp.139-148
• /
• 2012

The MLS(Moving Least Squares) Difference Method is a numerical scheme that combines the MLS method of Meshfree method and Taylor expansion involving not numerical quadrature or mesh structure but only nodes. This paper presents an dynamic algorithm of MLS difference method for solving transient solid mechanics problems. The developed algorithm performs time integration by using Newmark method and directly discretizes strong forms. It is very convenient to increase the order of Taylor polynomial because derivative approximations are obtained by the Taylor series expanded by MLS method without real differentiation. The accuracy and efficiency of the dynamic algorithm are verified through numerical experiments. Numerical results converge very well to the closed-form solutions and show less oscillation and periodic error than FEM(Finite Element Method).

• Geophysics and Geophysical Exploration
• /
• v.10 no.2
• /
• pp.154-160
• /
• 2007

This paper compares the characteristics of three different algorithms for three-dimensional (3D) magnetotelluric (MT) modeling. These methods are developed by Mackie et al. (1994), Sasaki (1999) and Nam et al. (2007). The first and second methods are based on the finite difference method (FDM), while the last one the finite-element method (FEM). MT responses, apparent resistivities and phases, for a COMMEMI 3D-2 model show a good agreement with integral equation solutions and only minor discrepancies are found over the anomalous bodies in the 3D model. The computation time of the two methods based on FDM is short and the static divergence correction contributes to speed up. The FEM modeling scheme is accurate but slow.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.11 no.1
• /
• pp.7-19
• /
• 1999

The energy transfer between sea-wave components by way of nonlinear wave-wave interactions plays a central role in spectral evolution. Since huge calculation time is required to exact computation of the resulting Boltzmann integral, however, the exact nonlinear energy transfer has not been directly introduced into operational wave models. Thus, effective calculation methods were examined in the present study which exploit the scale property of a scattering coefficient and the detailed balance of interactions. The improved Webb's method (IWM) has inherent stability because singularities degenerate into a negligible point. The improved Masuda's method (IMM) makes a quasi-analytical treatment of the inherent singularities and requires only 1.3 seconds of computer time via Pentium 300MHz processor. The IMM is, therefore, projected to be very useful for theoretical researches in spectral evolution with fetch- or duration-limited situations.