• 대한기계학회논문집A
• /
• 제34권8호
• /
• pp.1007-1013
• /
• 2010

본 연구에서는 NURBS 의 국부 세분화 방법 중 하나인 계층적 B-스플라인을 이용해 스플라인 유한요소법의 국부 세분화를 수행하는 방법을 제안한다. 세분화가 필요한 영역에 전역 격자로부터 계층적으로 생성된 국소 격자를 중첩시켜 국부 세분화를 수행한다. 국소 격자의 매듭 벡터와 제어점은 전역 격자로부터 추출된 후 세분화 되는 과정을 거친다. 생성된 국소 격자에 적절한 연속성 조건을 부여 함으로써 전역 격자와 국소 격자의 연속성을 유지 한다. 제안된 방법을 이용해 수치 예제의 해석을 수행하였다. 이를 통해 기존 NURBS 기반 스플라인 유한요소법에 비해 제안된 방법의 효율성을 검증하였다.

• 전산구조공학
• /
• 제10권1호
• /
• pp.201-211
• /
• 1997

공간뼈대의 구조에 대하여 기하학적 비선형성이 고려될 수 있는 유한요소이론 및 해석법을 제시한다. 이를 위하여 가상일의 원리를 이용하여 대변형효과를 고려한 3차원 연소체의 평형방정식으로부터, 구속된(restrained warping)효과를 무시하고 유한한 회전각의 2차항의 효과를 포함하는 변위장을 도입하여 초기응력을 받는 공간뼈대요소의 증분평형방정식을 유도한다. 공간뼈대구조를 유한요소로 나누어 요소의 변위장을 요소변위 벡터에 관한 Hermitian다항식으로 나타내고 이를 평형방정식에 대입함으로써 탄성 및 가하학적인 강도행렬을 유도한다. 또한 updated Lagrangian co-rotational formulation에 근거하여, 증분변위로부터 강체회전변위와 순수변형성분을 분리시켜서 강체회전은 요소의 방향변화를 결정하고, 순수변형은 부재력증분을 산정하는 불평형하중 산정법을 제시한다. 공간뼈대구조의 횡-비틂좌굴 및 후좌굴 거동에 대한 예제들을 통하여 본 연구에 대한 해석결과와 문헌의 결과를 비교 검토함으로써 본 연구에서 제시된 이론 및 해석방법의 정당성을 입증한다.

• 한국군사과학기술학회지
• /
• 제1권1호
• /
• pp.128-144
• /
• 1998

일반적인 선형 탄성해석 경계 요소법이 초 고속 회전과 정상 열전도에 의한 열 탄성 효과가 고려된 문제에 적용되었다. 축대칭 경계 요소법 구성이 요약되었고, 등가 경계 적분 방정식의 물체력 핵 함수의 체적 적분 전환방법에 일반화된 내적과 벡터 연산법 개념이 도입되었다. 고차 경계 요소 적용을 위한 이산화 수치 해석법이 요약되었고, 소형 젯트 엔진(ADD 500)의 터어빈 로우터 디스크의 해석 결과가 유한 요소해와 비교되었다.

• 한국지진공학회논문집
• /
• 제3권1호
• /
• pp.41-54
• /
• 1999

비대칭 박벽단면을 갖는 곡선보의 자유진동해석을 수행할 수 있는 유한요소 이론 및 엄밀해를 제시하기 위하여 가상일의 원리를 이용한 3차원 연속체의 운동방정식을 제시한다 박벽단면의 구속된 비틂효과를 고려하는 박벽 곡선보의 변위장을 도입하고 이를 연소체의 운동방정식에 대입하여 단면에 대해 적분함으로써 박벽 곡선보의 운동방정식을 유도한다. 단순지지되고 일축대칭단면을 갖는 박벽 곡선보의 면내 자유진동 모드에 대응하는엄밀해를 산정하였으며 곡선보를 유한요소로 분할하여 요소의 변위장을 요소 변위벡터에 관한 3차의 Hermitian 다항식으로 나타내고 이를 운동방정식에 대입함으로써 탄성강도행렬과 질량 행렬을 유도한다 또한 본 연구에서 얻어진 엄밀해와 곡선보요소를 이용한 유한요소 해석결과를 직선보요소 및 ABAQUS의 쉘요소를 이용하여 얻어진 결과와 비교 검토를 함으로써 본 연구의 타당성을 입증한다.

• 한국전산구조공학회논문집
• /
• 제24권3호
• /
• pp.329-335
• /
• 2011

유한요소법은 수학과 공학을 비롯한 다양한 분야에서 활용되고 있으나 해석대상을 유한 개의 다각형 요소로 분할하여 모델링하기 때문에 기하학적인 형상을 정확하게 기술하지 못하는 어려움이 있다. 그러나 최근에는 NURBS(Non-Uniform Rational B-Spline)를 기저함수로 사용하는 아이소-지오메트릭 해석법(Isogeometric analysis)이 개발되었는데 NURBS는 기하학적 모델을 정확하게 표현할 수 있을 뿐만 아니라 해석의 기저함수로서 응답해석에 사용될 수 있다. 그러나 NURBS 기저함수를 해석에서 따로 구성하는 일은 유한요소해석에서 요소망을 구성하는 만큼 시간과 노력이 많이 요구된다. 아이소-지오메트릭 해석법은 CAD(Computer-Aided Design)와 기하학적 정보를 공유할 수 있기 때문에 CAD 코드로 부터 해석모델의 정보를 직접 얻는 것이 가능하다. 본 논문에서는 상용 CAD 코드인 Rhinoceros 3D를 이용하여 CAD 모델을 작성하고 이를 STEP 파일로 출력하여 NURBS의 노트벡터와 조정점 등의 정보를 아이소-지오메트릭 해석법에 활용하는 기법을 소개한다. 몇몇 수치예제를 통하여 아이소-지오메트릭 해석법의 정확도를 유한요소해석 결과와 비교하여 검증하고, 상용 CAD와 CAE(Computer-Aided Engineering)가 결합된 아이소-지오메트릭 해석법의 효율성을 입증한다.

• 한국강구조학회 논문집
• /
• 제10권3호통권36호
• /
• pp.509-523
• /
• 1998

본 연구에서는 이축 대칭단면을 갖는 박벽 원형아치의 안정성해석을 수행할 수 있는 유한요소 이론 및 엄밀해를 제시하기 위하여, 가상일의 원리를 이용한 3차원 연속체의 운동방정식을 제시한다. 박벽단면의 구속된 비틂(restrained warping)효과를 고려하는 박벽 곡선보의 변위장을 도입하고 이를 연속체의 운동방정식에 대입하여 단면에 대해 적분함으로써 박벽 곡선보의 운동방정식을 유도한다. 단순지지되고 이축 대칭단면을 갖는 박벽 곡선보의 면외좌굴에 대한 엄밀해를 제시하고 박벽 곡선보를 유한요소로 분할하여 요소의 변위장을 요소 변위벡터에 관한 3차의 Hermitian 다항식으로 나타내어 운동방정식에 대입함으로써 탄성 강도행렬과 기하학적 강도행렬을 유도한다. 또한 본 연구에서 얻어진 엄밀해와 박벽 곡선보요소를 이용한 유한요소해석결과를 다른 연구자들의 결과 및 직선 박벽보 요소를 이용한 해석결과와 비교 검토를 함으로써 분 연구의 타당성을 입증한다.

• 대한기계학회논문집B
• /
• 제40권9호
• /
• pp.597-604
• /
• 2016

본 연구에서는 GPU를 이용한 비압축성 유동장의 병렬연산을 위하여, P2P1 유한요소를 이용한 분리 알고리즘 내의 행렬 해법인 이중공액구배법(Bi-Conjugate Gradient)의 CUDA 기반 알고리즘을 개발하였다. 개발된 알고리즘을 이용해 비대칭 협착관 유동을 해석하고, 단일 CPU와의 계산시간을 비교하여 GPU 병렬 연산의 성능 향상을 측정하였다. 또한, 비대칭 협착관 유동 문제와 다른 행렬 패턴을 가지는 유체구조 상호작용 문제에 대하여 이중공액구배법 내의 희소 행렬과 벡터의 곱에 대한 GPU의 병렬성능을 확인하였다. 개발된 코드는 희소 행렬의 1개의 행과 벡터의 내적을 병렬 연산하는 커널(Kernel)로 구성되며, 최적화는 병렬 감소 연산(Parallel Reduction), 메모리 코얼레싱(Coalescing) 효과를 이용하여 구현하였다. 또한, 커널 생성 시 워프(Warp)의 크기에 따른 성능 차이를 확인하였다. 표준예제들에 대한 GPU 병렬연산속도는 CPU 대비 약 7배 이상 향상됨을 확인하였다.

• 대한기계학회논문집
• /
• 제11권1호
• /
• pp.82-87
• /
• 1987

본 논문에서는 열전달문제 역시 변분형으로 전환될 수 있음에 착안하여 전미 분 개념을 도입해서 전도와 대류가 있는 열전달모델에서 주어진 면적 제한조건을 만족 시키며 지정된 경계에서의 온도가 주어진 온도에 가장 근접할 수 있는 모델의 형상을 찾는 방법을 연구하였다. 어떤 물질의 열전달 상태를 바꾸어 경계에서의 온도를 원 하는 바대로 조정하는 문제는 실제 공정에서 중요한 경우가 많다.해석시 열전달 상 태 방정식과 adjoint식은 6절점 삼각형 등계수 요소의 유한 요소법을 이용하여 해석하 였다.설계민감도의 정확한 계산을 위해서는 임의의 형상변화에 따른 경계에서의 수 치적분이 정확해야 하므로 경계를 곡선으로 표시할 수 있는 등계수 요소가 필요하다. 설계 민감도 해석이 진행된 후에는 최적화 기법의 하나인 미분벡터 투영법(Gradient Projection Method)을 사용하여 최적화를 시도했다. 최적설계 과정중 매번 계산결과 에 의해 형상의 변화가 진행되므로 그때마다 유한 요소 모델을 적절히 변화시켜 주어 야 한다. 모델의 경계는 3차함수로 근사화하여 형상이 부드러운 곡선이 되도록 했 으며 설계변수는 근사화한 3차함수를 결정할 수 있도록 정하면 되나 본 연구에서는 모 델의 변화에 따른 y좌표의 변화는 없다고 가정하여 모델경계의 세점을 취해 그 점들의 x좌표를 설계변수로 했다.

• 한국항공우주학회지
• /
• 제33권9호
• /
• pp.56-65
• /
• 2005

시간유한요소법은 시간영역을 고정시키고 행렬 미분방정식 형태의 공간전파 관계식을 풂으로써 시간과 공간에 대한 동적 해석을 수행하는 방법이다. 이 방법은 공간이산화 유한요소법이나 시/공간 동시이산화 유한요소법에 비해 공간에 관한 자유도가 발생하는 것이 두드러진 특징으로, 이를 이용하여 분포형 구동기의 공간에 따른 특성을 최적화하는 데에 효율적으로 사용될 수 있다. 본 논문에서는 임의의 초기조건을 반영할 수 있도록 구성된 상태변수 벡터를 이용하여 구조물을 시간영역에서 이산화하고, 공간영역에서 전파관계식 및 경계조건을 이용하여 공간전파 관계식을 형성하였다. 이 때 구동기의 공간에 따른 형상 분포는 설계되어야 할 변수의 함수이고, 시간반응은 형상함수를 이용하여 이산화 하였다. 포텐셜 에너지 및 운동에너지를 구조물의 변위제어에 적절한 최적의 성능지수로 설정하고, 이를 최소화하도록 미지의 함수인 구동기의 분포형상을 구하였다. 일반적으로 구조물은 임의의 초기조건에서 외란을 받게 되나, 본 연구에서는 구현가능한 제어법칙을 이용하여 최종시간에서 안정화(rest) 조건을 만족한다고 가정하였다. 구동기 분포형상 최적화를 위해 상태/준상태 방정식을 유도하였다. 서브행렬 재형상화와 시/공간 경계조건을 통해 상태변수와 준상태변수에 대한 Ricatti 미분방정식을 유도하였다. 이를 통해 구동기 분포형상 최적화를 구현하였으며, 수치 시뮬레이션을 통해 적절한 구동기의 분포형상 최적화를 수행할 수 있음을 보였다.

• 전자공학회논문지
• /
• 제50권11호
• /
• pp.28-35
• /
• 2013

Krylov-Schur 반복법을 활용하여 2-차원 사각 도파관에서 나타나는 고유특성을 밝혔다. 고유 행렬 방정식은 삼각형 그물 요소의 접선을 기저벡터로 사용한 FEM(유한요소법)으로 구성하였다. 우선 Arnoldi 분해법을 이용하여 이 방정식에 대한 상위 Hessenberg 행렬을 구하였다. 그리고 QR 알골리즘을 통하여 이것을 삼각형 대각 행렬인 Shur 형태로 변형하였다. 수렴 조건에 부합된 몇몇 고유 값들이 삼각형 대각 행렬의 대각 요소에 나타났다. 이들에 대응하는 고유 모드들을 역-반복법으로 구하였다. 수렴조건에 부합되는 고유 값들은 Shur 행렬의 대각선 선두 부분으로 재배열시켰다. 이들은 나머지 고유값 및 고유모드의 쌍을 구하는 반복 과정에서 변형되지 않도록 배제되었다. 이 과정이 연속하여 서너 번 반복되었는데, 그 결과 충분한 신뢰도를 갖는 주요한 몇 개의 TM 및 TE 고유 쌍들이 구하여졌다.