• 제목/요약/키워드: 적분-미분 정식화

검색결과 8건 처리시간 0.026초

이동최소제곱 유한차분법을 이용한 응력집중문제 해석(I) : 고체문제의 정식화 (Analysis of Stress Concentration Problems Using Moving Least Squares Finite Difference Method(I) : Formulation for Solid Mechanics Problem)

  • 윤영철;김효진;김동조;윙 캠 리우;테드 벨리치코;이상호
    • 한국전산구조공학회논문집
    • /
    • 제20권4호
    • /
    • pp.493-499
    • /
    • 2007
  • 본 연구에서는 미분 가능한 함수가 Taylor 전개로 표현되고 그 계수들은 주어진 함수와 미분에 대한 근사값을 제공할 수 있다는 점에 착안하여 m차 Taylor 다항식을 구성하고 이동최소제곱법을 이용하여 그 계수들을 구했다. 계산된 근사함수와 미분을 콜로케이션 개념을 바탕으로 균열 문제를 포함하는 고체문제에 대한 지배 미분방정식에 적용하여 차분식 형태의 이산화된 계방정식을 구성하였다. 본 연구의 해석기법은 격자망(grid)에 의존적이고 근사함수가 없는 유한차분법과 형상함수의 미분과 약형식의 적분산정, 필수경계조건 처리가 어려운 Galerkin법 기반의 무요소법의 단점을 효과적으로 극복한 새로운 수치기법이다.

안정된 전자파 과도 산란해를 얻기 위한 시간영역 전장 적분방정식 해석 (Time-Domain Electric Field Integral Equation Solving for a Stable Solution of Electromagnetic Transient Scattering)

  • 정백호;김채영
    • 대한전자공학회논문지TC
    • /
    • 제39권4호
    • /
    • pp.201-208
    • /
    • 2002
  • 본 논문에서는 3 차원 임의 형태 도체의 지연 산란 응답을 얻기 위한, 새로운 시간영역 전장 적분방정식(Time-Domain Electric Field Integral Equation: TD-EFIE)을 제안한다. 자기 벡터 전위의 시간 미분항은 중앙 차분으로, 전기 스칼라 전위는 시간에 대한 평균을 취한 두 개의 항으로 근사하였다. 이로부터 도체에 의한 산란 지연 응답 해의 산출시, 기존의 방법보다 정확하고 더욱 안정된 해를 얻을 수 있었다. 제안된 방법의 자세한 정식화 과정을 보였으며, 주파수 영역에서의 이산 푸리에 역변환 (Inverse Discrete Fourier Transform: IDFT) 결과치와 제안된 방법에 의한 수치해를 각각 비교하였다.

와도를 기저로 한 비압축성 점성유동해석 방법 (A Vorticity-Based Method for Incompressible Viscous Flow Analysis)

  • 서정천
    • 한국전산유체공학회지
    • /
    • 제3권1호
    • /
    • pp.11-21
    • /
    • 1998
  • A vorticity-based method for the numerical solution of the two-dimensional incompressible Navier-Stokes equations is presented. The governing equations for vorticity, velocity and pressure variables are expressed in an integro-differential form. The global coupling between the vorticity and the pressure boundary conditions is fully considered in an iterative procedure when numerical schemes are employed. The finite volume method of the second order TVD scheme is implemented to integrate the vorticity transport equation with the dynamic vorticity boundary condition. The velocity field is obtained by using the Biot-Savart integral. The Green's scalar identity is used to solve the total pressure in an integral approach similar to the surface panel methods which have been well established for potential flow analysis. The present formulation is validated by comparison with data from the literature for the two-dimensional cavity flow driven by shear in a square cavity. We take two types of the cavity now: (ⅰ) driven by non-uniform shear on top lid and body forces for which the exact solution exists, and (ⅱ) driven only by uniform shear (of the classical type).

  • PDF

탄성균열 해석을 위한 이동최소제곱 유한차분법의 내적확장 (Intrinsic Enrichment of Moving Least Squares Finite Difference Method for Solving Elastic Crack Problems)

  • 윤영철;이상호
    • 대한토목학회논문집
    • /
    • 제29권5A호
    • /
    • pp.457-465
    • /
    • 2009
  • 본 연구는 균열선단에서 응력특이성을 갖는 탄성균열문제를 해석하기 위한 이동최소제곱 유한차분법을 제시한다. 응력특이성을 유발하는 균열선단 주변장을 모형화하기 위해 근사식에 선단주변함수를 내재적으로 도입하여 이동최소제곱 근사의 틀을 그대로 유지하면서 실제 미분계산을 거의 하지 않고 미분근사를 할 수 있는 이동최소제곱 Taylor 다항식 근사의 장점을 살렸다. 균열문제 정식화시 시간소모적인 적분과정이 필요한 약정식화 대신 해석영역에 배치된 절점에서 지배 미분방정식에 대한 차분식을 직접 구성하는 강정식화를 적용하여 계산 효율성을 향상시켰다. 균열문제 해석을 통해 내적확장된 이동최소제곱 유한차분법이 응력 특이성을 내포한 선단주변 변위장을 정확히 묘사할 수 있을 뿐만 아니라 응력확대계수를 정확히 계산 할 수 있음을 보였다.

경계요소법을 이용한 압축성형다이 가열선의 최적위치 설계 (Optimal Positioning of Heating Lines in a Compression Molding Die Using the Boundary Element Method)

  • 이부윤;조종래
    • 대한기계학회논문집
    • /
    • 제17권6호
    • /
    • pp.1478-1485
    • /
    • 1993
  • 본 연구에서는 압축 성형다이 가열선의 최적위치를 결정하는 문제를 다이의 형상 최적설계 문제로 정식화하고자 한다. 최적화 문제의 목적은 다이 내면의 온도를 지정된 온도로 균일하게 유지하는 것이다. 또한 Lee, Choi와 Kwak의 형상 설계 민감 도 해석을 위한 직접 미분 방법을 응용하여 가열선 위치변화에 관한 민감도 계산을 위한 경계 적분 방정식을 유도하고, 경계요소법으로 온도와 온도의 민감도를 해석하고 자 한다. 수치적 응용의 예로서, 넓은 평판의 성형을 위한 압축 성형다이의 가열선의 최적위치를 결정하는 문제를 다루고, 최적위치에서 가열선의 개수가 온도의 균일성에 미치는 영향을 검토한다.

비압축성 점성유동의 와도와 압력 경계조건 (On the Vorticity and Pressure Boundary Conditions for Viscous Incompressible Flows)

  • 서정천
    • 한국전산유체공학회:학술대회논문집
    • /
    • 한국전산유체공학회 1998년도 춘계 학술대회논문집
    • /
    • pp.15-28
    • /
    • 1998
  • As an alternative for solving the incompressible Navier-Stokes equations, we present a vorticity-based integro-differential formulation for vorticity, velocity and pressure variables. One of the most difficult problems encountered in the vorticity-based methods is the introduction of the proper value-value of vorticity or vorticity flux at the solid surface. A practical computational technique toward solving this problem is presented in connection with the coupling between the vorticity and the pressure boundary conditions. Numerical schemes based on an iterative procedure are employed to solve the governing equations with the boundary conditions for the three variables. A finite volume method is implemented to integrate the vorticity transport equation with the dynamic vorticity boundary condition . The velocity field is obtained by using the Biot-Savart integral derived from the mathematical vector identity. Green's scalar identity is used to solve the total pressure in an integral approach similar to the surface panel methods which have been well-established for potential flow analysis. The calculated results with the present mettled for two test problems are compared with data from the literature in order for its validation. The first test problem is one for the two-dimensional square cavity flow driven by shear on the top lid. Two cases are considered here: (i) one driven both by the specified non-uniform shear on the top lid and by the specified body forces acting through the cavity region, for which we find the exact solution, and (ii) one of the classical type (i.e., driven only by uniform shear). Secondly, the present mettled is applied to deal with the early development of the flow around an impulsively started circular cylinder.

  • PDF

다양한 적층 배열을 갖는 복합 신소재 판 구조물의 유한차분 비선형 해석 (Finite Difference Nonlinear Analysis of Composite Plate Structures with Various Layer Sequences)

  • 이상범;이상열;이래철
    • 한국구조물진단유지관리공학회 논문집
    • /
    • 제9권4호
    • /
    • pp.159-168
    • /
    • 2005
  • 본 연구는 다양한 적층 배열을 갖는 비등방성을 보이는 첨단 복합 신소재 판구조물의 유한 차분 비선형 해석을 수행한다. 복잡한 편미분 방정식으로 표현되는 역학문제들을 수치해석 하는 경우 본 연구에서 사용한 유한차분법은 유한요소법에 비하여 체눈 생성 및 수치적분 과정을 피할 수 장점을 갖는다. 유한차분법을 이용한 많은 연구들은 단지 에너지 방법을 사용한 고정 혹은 단순 경계조건에 대하여 수행되었다. 그러나 이러한 접근방법은 자유경계에 대하여 불가피하게 발생하는 가상점 문제를 충분히 만족시킬 수 없다. 그러므로 본 연구에서는 임의의 경계조건을 갖는 비등방성 복합 적층한의 비선형 거동 문제를 보다 효과적으로 해결할 수 있는 유한차분식을 정식화 하였다. 적층 배열 변화를 비롯한 다양한 매개변수에 대하여 본 연구에서 제안한 접근방법을 사용하여 적층판의 복잡한 비선형 거동을 분석하였다.

등매개변수 사변형요소를 적용한 유한요소해석법 (Finite element method adopting isoparametric formulation of the quadrilateral elements)

  • 이승현;한진태
    • 한국산학기술학회논문지
    • /
    • 제19권11호
    • /
    • pp.205-212
    • /
    • 2018
  • 본 연구에서는 상용 해석프로그램에서 구현하기에 어려움이 있는 지반공학적 문제를 해결하기 위한 쉽고 직관적인 해석 프로그램 개발의 일환으로 계산의 정확도가 상대적으로 높은 요소를 채택한 유한요소법을 정식화 하고 해석과정을 프로그램화 하였다. 개발된 프로그램의 계산과정에 있어서의 신뢰성 확인을 위해 두 가지 예에 대한 해석을 수행하고 결과분석을 해 보았는데 첫 번째 예는 등방구속압이 요소에 작용하는 경우이고 나머지 예는 전단응력이 요소의 측면에 작용하는 경우이다. 유한요소를 구성하는 요소로는 등매개변수 사변형 요소를 채택하였는데 요소내의 변위는 요소의 절점변위와 형상함수로 표현된다. 전체좌표(global coordinate)에 의한 미분계수로 표현되는 변형률을 얻기 위해 자코비언과 자연좌표(natural coordinate)를 이용하는 계산과정을 코딩하였다. 요소의 강성행렬을 정의하는 이중적분식을 수치적분으로 변환시키기 위해 4점 가우스 구적법을 적용하였다. 개발된 프로그램의 계산과정 검증을 위해 등방구속압이 작용하는 요소에 대한 해석을 수행한 결과 요소내의 네 개의 가우스점과 요소 중앙에 대해 계산된 응력값이 등방구속압과 일치됨을 알 수 있었다. 개발된 프로그램의 계산과정 검증을 위해 전단응력이 작용하는 요소에 대한 해석을 수행한 결과 요소내에 발생되는 횡방향응력 및 연직응력이 위치에 따라 변화됨을 알 수 있었으며 외력에 대한 발생응력의 크기 및 분포양상이 합리적임을 알 수 있었다.