• The Korean Journal of Applied Statistics
• /
• v.33 no.6
• /
• pp.665-682
• /
• 2020

This paper introduces numba library in Python, which improves computational efficiency of the provided implemented code written by naive Python language by applying just-in-time (JIT) compilation. To apply just-in-time compilation, the numba only needs to use a decorator on a target Python function. We provide implementation examples with numba for the permutation test and the parameter estimation for Gaussian mixture distribution. We also numerically show the efficiency of numba by comparing the total computation times of the implementation using naive python and the implementation using numba for each application.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.13 no.5
• /
• pp.2403-2407
• /
• 2012

Finite number of pile elements are considered in load transfer method. And section force and movement of each pile element are computed by considering compatibilities between pile displacement and the load transfer along a pile and between displacement and resistance at the tip of the pile. For the conventional load transfer method, large amount of computations due to iterations are needed. Formulation of finite difference equation from the differential equation which depicts pile behavior under axial loading was accomplished in order to simplify the computation for obtaining pile section forces and displacements. By comparing the results between the simplified computation method and the reported data, there was no difference between the two results.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.8
• /
• pp.1513-1518
• /
• 1992

Numerical solution technique is suggested to analyze the vibration of a spring composed of stacked beams fastened together. Bernoulli-Euler beam theory for small deflection is used, and incremental Coulomb friction law is adopted for the interface friction. The validity of the present solution technique is checked for the perfectly bonded case and the perfect sliding case.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.12
• /
• pp.2298-2305
• /
• 1992

To obtain the convenient solution of the multibody dynamic systems composed of flexible beams, linear finite element technique is adopted and the nodal coordinates are interpolated in the global inertia frame. Mass matrix becomes an extremely simple constant matrix and the force vector also becomes extremely simple because Coriolis acceleration and centrifugal force are not required. And the elastic force is also simply computed from the moving frame attached to the material. To solve the global differential algebraic euation. an ODE technique is adopted after Lagrange multiplier is computed by the accelerated iterative technique, and the time demanding procedures such as Newton-Raphson iterations and decomposition of the big matrix are not required. The accuracy of the present solution is checked by a well-known example problem.

• 전기의세계
• /
• v.39 no.3
• /
• pp.47-54
• /
• 1990

유한요소법과 경계요소법의 합성으로 전자계 해석을 하는 기법은 각 방법의 장점을 수용하여 경계가 없는 무한영역의 전자장을 분석하는 기법으로서 어떤 복잡하고 어려운 기하학적 구조의 문제도, 비선형이나 비균질성 재질의 문제도 쉽게 formulation이 가능하여 용이하게 해석할 수 있지만 전체 System matrix방정식이 비대칭이며 부분적인 full matrix를 형성하여 계산시간이 길어 진다는 단점도 있다. 적용예에서 보여 준 것과 같이 합성요소법은 그 해가 실제에 근사한 값을 가질수 있다고 생각되며, 계산시간을 단축시키기 위하여 직접법이나 반복법을 사용한 새로운 해법들이 도입되고 있다. 최근에는 system전체 node의 순서를 고려한 NDRA(Nested Dissection Reordering Algorithm)이 도입되고 있고, System matrix자체를 유한 요소법의 형태로 유지시키며 풀수 있는 방법으로 알려진 Absorbin 경계조건을 사용하여 전자파에 대한 해석을 하고 있다. 유한 및 경계요소 합성법은 초고압 옥외용 전력기기의 전자장 해석과 설계, 레이다나 안테나 등의 전자파 해석문제, 초전도 응용, 전력기기의 전자장해석과 설계, 우주공간에서의 전력전송문제 등을 쉽게 model화하여 적용할 수 있을 것이다.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.1 no.1
• /
• pp.1-10
• /
• 1997

Current seismic design philosophy for reinforced concrete (RC) structures on energy dissipation through large inelastic defomations. A nonlinear dynamic analysis which is used to represent this behavior is time consuming and expensive, particularly if the computations have to be repeated many times. Therefore, the selection of an efficient yet accurate alogorithm becomes important. The main objective of the present study is to propose a new technique herein called the prediction approach with siffness measure (PASM) method in the convetional direct integration methods, the triangular decomposition of matrix is required for solving equations of motion in every time step or every iteration. The PASM method uses a limited number of predetermined decomposed effective matrices obtained from stiffness states of the structure when it is deformed into the nonlinear range by statically applied cyclic loading. The method to be developed herein will reduce the overall numerical effort when compared to approaches which recompute the stiffness in each time step or iteration.

• Journal of the Korean Society of Surveying, Geodesy, Photogrammetry and Cartography
• /
• v.17 no.1
• /
• pp.33-39
• /
• 1999

In order to calculate geoidal height refered to Bessel ellipsoid, methods to translate geoidal heights from a certain coordinate system to an arbitrary system with the corresponding ellipsoid are studied. and geoidal heights refered to Bessel ellipsoid were computed from EGM96 Model refered to GRS80 using iteration method pro-posed in this paper. Transformation parameters between WGS84 and Bessel were calculated using geoidal heights computed from iteration method. The result of coordinate transformation(standard deviation) were 0.009 second in latitude and 0.006 in longitude and 0.393m in orthometric height.

• Computational Structural Engineering
• /
• v.8 no.4
• /
• pp.65-74
• /
• 1995

Direct method developed for the inelastic analysis of planar frames subjected to monotonic loads is extended to cyclic loads. Two frame elements for Direct Method(inelastic truss and inelastic beam) are developed. The accuracy and reliability of the preposed method is verified by comparing the analysis results of example with step-by-step analysis. Direct Method is superior to Step-by-step analysis in view of reliability of solution and analysis cost.

• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 2002.09a
• /
• pp.156-160
• /
• 2002

실험적 베리오그램의 모델링에 SA(Simulated Annealing)기법을 이용하였다. 최소 자승법의 해를 구하기 위하여 기존의 상용 프로그램에서 많이 이용되고 있는 반복법에 근거한 방법에 비해서 SA 기법은 초기 가정값에 크게 영향을 받지 않고 일정한 모델 인자의 값을 제시하였다. 임의의 초기 가정값을 입력하여도 충분한 반복 계산을 통하여 목적함수의 값이 광역적 최소값으로 수렴하는 것을 확인할 수 있었다. 베리오그램 모델이 일반적으로 비선형 모델이기 때문에 목적함수의 지역적 최소값으로의 수렴이 문제가 되고 이로 인하여 구해지는 인자의 값이 정확하지 않을 수 있지만 SA 기법을 이용하여 최소 자승법의 해를 구하게 되면 정확한 인자의 값을 구할 수 있음을 확인하였다.

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

Krylov 부공간 모델차수축소법은 초기 유한요소모델과 축소모델의 전달함수의 계수인 모멘트를 일치시키는 방법을 이용하는 축소기법으로 이미 대형 유한요소모델의 주파수응답 해석의 효율적인 계산에 많이 사용되고 있는 방법 중의 하나이다. 본 논문에서는 Krylov 부공간 축소기법을 이용한 관심 주파수영역에 대한 주파수응답 해석 및 이를 통하여 계산된 주파수응답의 여러 가지 설계변수에 대한 설계민감도 해석방법을 제안하였다. 일반적으로 구조물의 주파수응답을 고려한 최적설계를 위해서는 설계변수에 대한 관심 주파수영역에서의 주파수응답 및 그의 민감도 정보가 요구되므로, 고려하는 유한요소모델이 대형일 경우에 관심 주파수영역에서의 반복적인 해석으로 인한 계산비용의 문제가 대두된다. 본 논문에서는 축소모델을 이용하여 주파수응답과 주파수응답의 설계민감도 해석을 수행하여 계산의 효율성을 극대화하였다. 민감도 계산에는 시간측면과 구현의 용이성 측면에서 장점이 있는 준해석적 방법을 이용하였다. 수치 예제를 통하여 축소기법을 이용한 주파수응답의 설계민감도 해석 결과를 유한차분법에 근거한 민감도 결과와 비교하였다. 본 논문에서 제안된 방법을 이용하는 경우, 주파수응답을 고려한 최적설계를 계산비용 측면에서 매우 효율적으로 수행할 수 있을 것이다.