• 제목/요약/키워드: MLS difference method

검색결과 23건 처리시간 0.025초

1차 미분 근사를 이용한 MLS차분법의 동적해석 (Dynamic Analysis of MLS Difference Method using First Order Differential Approximation)

  • 김경환;윤영철;이상호
    • 한국전산구조공학회논문집
    • /
    • 제31권6호
    • /
    • pp.331-337
    • /
    • 2018
  • 본 논문은 MLS(moving least squares) 차분법의 1차 미분 근사함수를 바탕으로 시간에 따른 수치해석이 가능한 해석기법을 제시한다. 오직 1차 미분 근사함수로만 지배방정식을 이산화했으며, 근사함수를 조립하는 형태로 전체 시스템 방정식을 구성하여 차분법으로 이산화된 운동방정식이 유한요소법(finite element method)과 유사한 모습을 갖게 되었다. 운동방정식을 시간적분하기 위해서 중앙차분법(central difference method)을 사용하였다. 유한요소 알고리즘을 통해서 MLS 차분법과 유한요소법의 고유진동 해석을 수행하였으며, 두 해석결과를 비교하였다. 또한, 동적해석 결과를 기존의 2차 미분 근사함수를 활용한 해석결과와 함께 도시함으로써 제안된 수치기법의 정확성을 검증하였다. 1차 미분 근사함수를 조립하는 과정에서 해석결과의 떨림현상이 억제되었으며 상대적으로 균일한 응력분포를 구할 수 있었다.

MLS차분법을 이용한 재료비선형 문제 해석 (Development of MLS Difference Method for Material Nonlinear Problem)

  • 윤영철
    • 한국전산구조공학회논문집
    • /
    • 제29권3호
    • /
    • pp.237-244
    • /
    • 2016
  • 본 연구는 재료비선형 문제를 다루기 위한 비선형 MLS 차분법의 정식화 과정을 제시한다. MLS 차분법은 절점모델을 기반으로 고속 미분근사식을 활용하여 지배 미분방정식을 직접 이산화 하는데, 변수를 변위로 일원화한 Navier 방정식을 사용하여 탄성재료 문제를 다룬 기존의 MLS 차분법은 재료의 구성방정식을 별도로 고려할 수 없다. 본 연구에서는 비선형 재료의 구성방정식을 반영할 수 있는 강정식화를 위해 1차 미분근사를 반복 사용하는 겹미분근사를 고안했다. 응력의 발산으로 표현되는 평형방정식을 그대로 이산화하고 Newton 방법을 적용하여 반복계산을 통해 수렴해를 찾는 비선형 알고리즘을 제시했다. 응력 계산과 내부변수의 갱신은 return mapping 알고리즘을 활용하였고, 알고리즘 접선계수(algorithmic tangent modulus)의 적용을 통해 빠르고 안정적인 반복계산이 가능하도록 하였다. 재생성 시험을 통해 겹미분근사의 정당성을 검증했고, 비선형재료에 대한 인장문제의 해석을 통해 개발된 비선형 MLS 차분 알고리즘의 정확성과 안정성을 확인하였다.

MLS 차분법을 이용한 고체역학 문제의 동적해석 (Dynamic Algorithm for Solid Problems using MLS Difference Method)

  • 윤영철;김경환;이상호
    • 한국전산구조공학회논문집
    • /
    • 제25권2호
    • /
    • pp.139-148
    • /
    • 2012
  • MLS(Moving Least Squares) 차분법은 무요소법의 이동최소제곱법과 Taylor 전개를 이용하여 요소망의 제약 및 수치 적분이 없이 절점만을 이용하여 미분방정식을 수치해석할 수 있는 방법이다. 본 연구에서는 고체역학 문제의 동적해석을 위하여 MLS 차분법의 시간이력해석 알고리즘을 제시한다. 개발된 알고리즘은 Newmark 방법으로 시간적분을 하였으며, 강형식을 그대로 이산화하여 해석을 수행했다. 이동최소제곱법을 이용해 Taylor 전개식을 근사하여 실제 미분계산없이 미분근사식을 얻기 때문에 고차까지 Taylor 다항식의 차수를 증가하는 것이 용이하다. 1차원과 2차원 수치예제들을 통하여 동적해석을 위한 MLS 차분법의 정확성과 효율성을 검증하였다. 수치결과들이 정확해에 잘 수렴하였으며, 유한요소법(FEM)의 해석결과와 비교하여 떨림현상(oscillation) 및 주기성(periodicity) 오차에 대해 보다 안정적인 모습을 보였다.

MLS 차분법을 이용한 동적균열전파 해석 (Analysis of Dynamic Crack Propagation using MLS Difference Method)

  • 윤영철;김경환;이상호
    • 한국전산구조공학회논문집
    • /
    • 제27권1호
    • /
    • pp.17-26
    • /
    • 2014
  • 본 논문은 MLS(Moving Least Squares) 차분법을 바탕으로 동적균열전파 해석을 수행하기 위한 알고리즘을 제시한다. MLS 차분법은 절점만으로 이루어진 수치모델을 사용하며, 이동최소제곱법을 이용하여 전개한 Taylor 다항식을 기초로 미분근사식을 유도하기 때문에, 요소망의 제약에서 완벽하게 벗어난 절점해석이 가능하다. 시간항을 포함하는 동적 평형방정식은 Newmark 방법으로 시간적분 하였다. 동적하중을 받는 균열이 전파할 때, 매 시간단계마다 절점모델을 재구성하지 않고 균열선단 주변에서 국부적인 수정을 통해 해석이 가능하다. 동적균열을 묘사하기 위해 가시한계법(visibility criterion)을 적용하였고, 동적 에너지해방률을 산정하여 균열의 진전유무와 그에 상응하는 진전방향을 결정하였다. 모드 I 상태와 혼합모드 상태에서 균열이 진전하는 현상을 모사하였고, 이론해와 Element-Free Galerkin법으로 계산한 결과와의 비교를 통해 개발된 알고리즘의 정확성과 안정성을 검증하였다.

복합 불연속면을 갖는 포텐셜 문제 해석을 위한 확장된 MLS 차분법 (Extended MLS Difference Method for Potential Problem with Weak and Strong Discontinuities)

  • 윤영철;노혁천
    • 한국전산구조공학회논문집
    • /
    • 제24권5호
    • /
    • pp.577-588
    • /
    • 2011
  • 본 논문은 복합 불연속면을 갖는 포텐셜 문제의 해석을 위해 확장된 MLS(Moving Least Squares) 차분법을 제시한다. 계면경계를 따라 해(solution)와 수직방향, 접선방향 미분들이 모두 불연속 특이성을 나타내는 복합 불연속면을 묘사하기 위해 계단함수, 쐐기함수, 가위함수와 같은 불연속 특이함수를 추가하여 기존의 MLS 차분법을 개선했다. 계면경계조건은 기지의 조건으로서 지배방정식의 이산화과정에서 추가의 미지계수를 발생시키지 않는다. 포아송 방정식 형태의 지배미분 방정식을 풀기 위해 내부영역과 경계에 절점을 배치하고 차분식을 구성한다. 차분식을 조립한 계 방정식을 직접 풀기 때문에 계산효율성이 매우 우수하다. 수치예제는 제시된 해석기법의 우수성을 잘 보여주며, 균열전파, 이동경계, 상호작용 문제 등 다양한 불연속 문제로의 확장이 기대된다.

MLS 차분법의 결정 변수에 따른 정확도 분석 및 혼합변분이론을 통한 미분근사 성능향상 (On the Improvement of the Accuracy of Higher Order Derivatives in the MLS(Moving Least Square) Difference Method via Mixed Formulation)

  • 김현영;김준식
    • 한국전산구조공학회논문집
    • /
    • 제33권5호
    • /
    • pp.279-286
    • /
    • 2020
  • 본 연구에서는 점근해석 및 논로컬 이론에서 요구하는 4차 이상의 고차 미분근사를 수행하기 위하여 계방정식에 혼합변분이론을 적용하여 MLS 차분법으로부터 구해지는 고차 미분근사의 정확도를 큰 폭으로 향상시킨다. 또한, MLS 차분법에 존재하는 세 가지 조건변수에 따른 고차미분근사의 정확도를 비교·분석한다. 혼합변분이론의 합응력을 후처리하여 변위의 미분을 근사할 경우 기존의 변위장 기반 계방정식의 차분 결과에 비해 미분 차수가 2차 낮아진다. 해석 범위내 절점 수가 과도하게 많거나 기저 차수가 클 경우 MLS 차분법의 영향영역 내에서 과적합(overfitting)이 발생한다. 또한 영향영역이 최적 범위 이상으로 넓어질 경우 근사의 정확도가 떨어진다. 위 내용을 사인 하중을 받는 단순지지보 수치예제로부터 확인하였다.

MLS 차분법을 위한 Explicit 동적해석 알고리즘 개발 (Development of Explicit Dynamic Algorithm for MLS Difference scheme)

  • 김경환;윤영철;이상호
    • 한국전산구조공학회:학술대회논문집
    • /
    • 한국전산구조공학회 2010년도 정기 학술대회
    • /
    • pp.179-182
    • /
    • 2010
  • 본 연구에서는 MLS 차분법을 이용하여 동역학 문제를 해석하기 위한 explicit 동적해석 알고리즘을 제시한다. 격자망이 없는 장점을 부각시키기 위해 이동최소제곱법에 근거한 Taylor 전개로부터 미분근사를 얻고 차분식을 구성했다. 지배 미분방정식의 시간항을 CDM(Central difference Method) 차분하여 빠른 속도로 동적해석을 수행하였다. 수치결과를 통해 본 연구에서 제시한 알고리즘의 정확성과 안정성을 확인할 수 있었다.

  • PDF

MLS 차분법을 활용한 동적 균열전파해석의 Rayleigh 감쇠영향 분석 (A Study of Rayleigh Damping Effect on Dynamic Crack Propagation Analysis using MLS Difference Method)

  • 김경환;이상호;윤영철
    • 한국전산구조공학회논문집
    • /
    • 제29권6호
    • /
    • pp.583-590
    • /
    • 2016
  • 본 논문은 강형식 기반의 MLS 차분법에 Rayleigh 감쇠효과를 적용한 동적균열진전 해석기법을 제시한다. Rayleigh 감쇠 효과가 반영된 동적 평형방정식과 구성방정식을 도출하고, MLS 미분근사식을 이용하여 지배방정식들을 이산화하였다. 평형방정식뿐만 아니라 구성방정식에서도 감쇠효과를 적절하게 고려하여 기존의 무요소 강정식화 기법에서 고려하지 못했던 비례감쇠 알고리즘을 구현하였다. 시간관련 항을 포함한 동적 평형방정식은 중앙차분법(central difference method)을 이용하여 시간적분 하였고, 속도에 대한 차분식을 lagging시켜 이산화 방정식을 간소화시켰다. 균열의 기하학적 특성은 표면력 '0'인 자연경계 조건을 균열면에 놓인 절점들에 부과하여 묘사하였으며, 균열성장으로 인해 해석단계마다 변하는 절점의 생성 및 이동 효과를 계방정식 구성에 반영하였다. 단일균열과 다중균열을 갖는 수치예제를 통해서 제안된 수치기법의 정확성을 검증하였으며, 비례감쇠 효과의 고려가 동적균열진전 해석결과에 미치는 영향을 보였다.

MLS 유한차분법을 이용한 복합재료의 열전달문제 해석 (Heat Transfer Analysis of Composite Materials Using MLS Finite Difference Method)

  • 윤영철
    • 한국전산구조공학회:학술대회논문집
    • /
    • 한국전산구조공학회 2008년도 정기 학술대회
    • /
    • pp.2-7
    • /
    • 2008
  • A highly efficient moving least squares finite difference method (MLS FDM) for heat transfer analysis of composite material with interface. In the MLS FDM, governing differential equations are directly discretized at each node. No grid structure is required in the solution procedure. The discretization of governing equations are done by Taylor expansion based on moving least squares method. A wedge function is designed for the modeling of the derivative jump across the interface. Numerical examples showed that the numerical scheme shows very good computational efficiency together with high aocuracy so that the scheme for heat transfer problem with different heat conductivities was successfully verified.

  • PDF

머리전달함수 측정법의 실험적 비교 (Comparison of Measurement Methods for Head-related Transfer Function(HRTF))

  • 안태수;이두호
    • 한국소음진동공학회논문집
    • /
    • 제19권12호
    • /
    • pp.1260-1268
    • /
    • 2009
  • Three methods(the stepped sine method, the statistical method(random excitation method) and the maximum-length sequence(MLS) method) for head-related transfer functions(HRTFs) are experimentally compared in view point of accuracy and efficiency. First, the stepped sine method has high signal-to-noise ratio, but low efficiency. Second, the statistical method is fast measurement speed, but weak to noise than the other methods. Finally, the MLS method shows both good efficiency and high signal-to-noise ratio, but it needs additional software or equipment such as MLS signal generator. For comparison of measurement accuracy, HRTFs of KEMAR dummy are measured for various azimuths and elevations. Error norms for magnitude and phase of HRTFs are defined and calculated for the measured HRTFs. The calculated error norms show that the methods give similar results in magnitude and phase except a little phase difference in the MLS method.