• 전산구조공학
• /
• 제11권1호
• /
• pp.217-226
• /
• 1998

본 논문은 변위제약모드를 갖는 트러스구조물의 형태해석을 목적으로 하였으며, 이를 위하여 해의 존재조건과 무어-펜로즈(Moore-Penrose) 일반역행렬을 이용하였다. 또한, 수치해석과정에서의 변위제약모드로는 호몰로지변형(homologous deformation)을 고려하여 해석하였고, 다음으로 다양한 변위제약모드와 절점에 작용하는 하중비를 만족하는 구조물의 형태를 구하였다. 본 논문에서의 형태해석문제는 지정된 변위를 만족하는 구조물의 형태를 찾는 일종의 역문제(inverse problem)로서 일반적인 구조해석과정과는 반대되는 입장에서 접근하였다. 또한, 본 논문에서는 수치해석과정에서 근사해의 정도를 향상시키기 위하여 뉴튼-랩슨법을 사용하였고, 수치해석예제로서 부재의 배열형태에 따라 3가지모델을 선택하였으며, 이들 모델을 통하여 적용한 해석기법의 정확성과 효율성을 검증하였다.

• 터널과지하공간
• /
• 제9권4호
• /
• pp.296-305
• /
• 1999

많은 터널 현장에서 계측을 수행하고 있지만, 일반적인 역해석 방법으로 이 계측자료를 통해 라이닝에 작용하는 변형율과 응력 등의 역학량을 추정하기는 어렵다. 그 이유는 현장에서 계측이 가능한 적은 변위수로부터 라이닝의 응력상태를 추정할 수 없기 때문이다. 계측의 시간과 비용을 고려해 볼 때, 계측자료를 이용할 수 있는 기법의 개발은 시급한 연구과제라 할 수 있다. 본 연구에서는 터널계측에서 계측 가능한 소수의 변위로 전테 라이닝의 변위와 응력 등의 역학량을 계산할 수 있는 역해석 모델을 개발하였다. 역해석 방법으로서 복공 외면에 작용하는 외력의 자승이 최소화될 때 해를 구하는 방법을 사용하였다. 라이닝 외면에 작용하는 하중을 가정하여 유한요소법으로 변위와 응력을 계산하고, 여기서 얻은 몇개의 변위를 계측 변위로 가정하여 역해석 모델에 입력함으로써 절점 변위, 요소의 변형율 및 응력 등을 다시 구하였다. 처음 계산된 변위와 역계산으로 구한 변위를 비교함으로써 모델의 신뢰성을 평가하였다. 자동계측으로 라이닝의 비교적 정확한 계측을 수행하고 있는 지하철 한강 하저터널 현장에서 얻은 계측변위를 개발한 역해석 모델을 적용하여 실제 계측한 응력과 유사한 결과를 얻었다.

• 한국전산구조공학회논문집
• /
• 제12권4호
• /
• pp.591-597
• /
• 1999

본 연구에서는 비 균질 퇴적층으로 인한 지진파의 증폭에 대한 경계/유한요소 해석을 수행하였다. 수치해석을 위해, 비 균질 퇴적층은 8절점 등 매개 변수 유한요소 사용하여 모델링하였고, 그 주위의 균질 반무한 지반은 3절점 등매개 변수 경계요소를 사용하여 모델링하였다. 경계요소와 유한요소의 접촉면에서, 표면 력의 평형조건과 변위의 적합 조건에 의해 두 개의 요소를 결합하는 알고리듬을 개발하였다. 수치해석의 영향인자로서 SH파, P파와 SV파의 입사각, 무 차원 진동수 그리고 반무한 지반과 퇴적층사이의 전단 파 속도 비와 질량밀도 비를 고려하였다.

• 수산해양기술연구
• /
• 제38권2호
• /
• pp.110-118
• /
• 2002

본 논문에서는 무요소 이론을 정식화하였고 이를 이용한 1차원 및 2차원 EFG 프로그램을 Visual Basic과 C언어를 이용하여 작성해 보았다. 그리고 각각의 EFG 수치해석의 예를 작성된 프로그램을 이용하여 해를 구하였다. 해석결과는 다른 문헌의 결과와 일치하였으며 해석결과에서 나타나듯이 무요소 해의 정도는 영향영역의 비례축소인자 dmax와 가중함수의 종류, 절점 배치형태에 의해 좌우된다는 사실을 알 수 있었다. 특히 1, 2차원 EFG 해석결과에서 가장 최적의 해를 보이며 정해(exact solution)에 가장 근접한 조건은 dmax = 2 이고 가중함수가 3차 Spline형일 때로 나타났으며 유한요소법과 마찬가지로 절점의 수가 많을 수록 그리고 절점을 균일하게 배치할수록 높은 정도를 나타내는 것을 알 수 있었다. 특히 2차원의 경우 3차 Spline형 이외의 다른 가중함수를 사용할 경우에 상당히 큰 오차를 나타내는 점은 1차원 EFG 해석의 결과와는 다른 점이었지만 그 외 대부분 같은 결과를 나타내었다. 1차원에서 절점을 임의로 배치한 경우는 비교적 균일하게 배치한 경우가 해에 근접하는 형태를 나타내었으며 절점 간격이 상대적으로 적은 곳에서 큰 오차를 나타내었다. 그리고 절점을 임의로 선택할 때 변위가 모두 ‘0’의 값을 가지는 경우를 볼 수 있는데, 이는 화면상의 좌표계산에서 생긴 미소한 오차가 절점들에 의해 반복됨으로서 발생하는 것으로 보인다. 또한 탄성계수 값이 클 경우 dmax 에서 계산이 제대로 수행되지 못하는 경우가 있는데, 이는 수치가 double형의 크기를 초과하기 때문인 것으로 보인다. 결과에서 나타나듯이 무요소법에서 적당한 가중함수와 비례축소 인자를 사용하면 정해에 가까운 우수한 해를 얻을 수 있다는 것을 알 수 있다. 비록 프로그래밍 과정이나 이론의 정식화가 유한요소법에 비해 상당히 어려운 점은 있으나 무요소법은 요소의 정보를 필요치 않으므로 사용자 입장에서는 매우 편리할 것이다. 앞으로 경계조건을 효과적으로 만족시키는 문제를 해결하고 효과적인 알고리즘이 개발된다면 실용적으로 유한요소법을 대신할 수 있는 좋은 대안이 될 수 있을 것이라 생각된다.

• 한국강구조학회 논문집
• /
• 제16권3호통권70호
• /
• pp.295-303
• /
• 2004

본 연구에서는 등방성 및 복합적층판 해석을 위해 추가변형률을 갖는 개선된 4절점 Reissner-Mindlin 평판휨요소를 제안하였다. 전단잠김현상과 가상적인 제로에너지모드를 제거하기 위해 비적합 변위모드와 Bubble 함수식에 근거한 새로운 형태의 전단변형률을 추가함으로써 횡방향 전단거동을 개선하였다. Andelfinger와 Ramm(1993)이 제시한 기본적인 추가변형률은 면내거동을 개선시키고자 그대로 적용하였다. 1차전단변형이론에 근거한 새로 개발된 4절점 평판요소를 '14EASP'라 명하였다. 14EASP 유한요소의 특징과 성능을 평가하고자 몇가지의 수치해석예제를 적용하였으며, 다른 유한요소 및 해석적인 해와 비교하였다. 그 결과 본 연구에서 제안한 14EASP는 보다 안정적이고, 수렴성이 빠르며, 특히 요소형상이 왜곡된 경우에도 정확한 결과를 도출하였다.

• 한국지반신소재학회논문집
• /
• 제20권3호
• /
• pp.11-20
• /
• 2021

본 논문에서는 유한요소법과 경계요소법을 결합하여 기하학적으로 급변 부위가 있는 이차원 탄성 정적 문제에 대하여 효율적이고 정확한 해석 결과를 얻기 위한 유한요소법과 경계요소법의 근사 결합 방법을 제시한다. 이차원 문제의 유한요소로서는 3절점, 4절점 평면응력 요소를 적용하고, 이차원 문제의 경계요소로는 3절점 경계요소를 적용한다. 모델링 단계에서는 우선 전체 해석 대상을 유한요소로 모델링한 후에 기학학적 급변 부위를 경계요소로 모델링 하는데, 유한요소의 모델링을 위하여 정의된 절점을 이용하여 경계요소를 정의한다. 해석 단계에서는 전체 해석 대상에 대하여 유한요소 해석을 우선적으로 수행하고, 이후에 경계요소 해석을 자동으로 수행하는데, 경계부에서의 경계조건은 유한요소 해석 결과인 변위 조건과 응력 조건을 적용한다. 수치예제로서 이차원 탄성 정적 문제인 균열이 있는 평판에 대한 해석 결과를 제시하고 고찰한다.

• 대한기계학회논문집
• /
• 제15권3호
• /
• pp.751-756
• /
• 1991

본 연구에서는 Ashwell이 제시한 변형률요소를 전단효과를 고려한 두꺼운 곡 선보 요소에 적용 하였다. 막 변형률, 곡률, 전단변형률 각각에 독립된 변형률 함수 를 가정하여 미분 방정식의 일반해를 구하면 정확한 강체변위의 표현은 물론, 강성과 잉현상을 피할 수 있고 얇은 곡선보에서 두꺼운 곡선보에 이르기까지 보의 해석에 있 어서, 2절점으로 구성되는 적은 자유도수에서 높은 정확도를 보여주는 간편하고도 효 율적인 요소를 개발하고자 하였다.

• 한국원자력학회:학술대회논문집
• /
• 한국원자력학회 1997년도 춘계학술발표회논문집(2)
• /
• pp.30-35
• /
• 1997

물성치가 상이한 계에서 접합면을 포함한 전 영역내의 연속 응력장 개선방안을 제안하였다. 재료의 물성치가 100배 차이가 나는 직선보 예제에 대해 변위형 유한요소해에서는 접착면 상ㆍ하측에서 응력들이 불연속이며 상당한 차이를 보이고 있는데 반하여 본 연구에서의 응력장은 연속이며 접합면의 절점응력들이 이론해에 근접하고 있다 또한 본 연구에서의 연속 응력장으로 계산한 변형률 에너지는 수회 이내의 반복계산에서 이론해에 수렴하고 있다.

• EDISON SW 활용 경진대회 논문집
• /
• 제4회(2015년)
• /
• pp.269-273
• /
• 2015

처음 에디슨 경진대회를 준비할 때 업로드 되어 있는 다양한 프로그램을 보면서 어떤 주제를 선정할 지 고민을 하던 중, 눈에 들어온 프로그램은 서울대학교에서 제작한 Co-rotational Plane beam-Dynamic tip load이란 프로그램이었다. 위 프로그램은 한쪽 단이 고정된 외팔보(cantilever beam)의 끝단에 하중을 시간에 관련된 함수(sin, cos)로 주어 각 절점에서의 변위(X, Y), 속도 그리고 가속도를 알려주는 프로그램이다. 우리는 이 프로그램을 이용해서 외팔보 끝단의 진동 형태와 주기가 어떤 요인(단면의 모양, 재료의 성질, 가하는 하중의 진동수 등)에 의해서 변하는지 분석할 수 있었다. 거기서 더 나아가 하중을 실제 지진상황에서의 크기와 유사하게 주고 위 프로그램에서의 보의 거동(진동 주기, 진동 변위 등)을 건물의 거동이라 가정했을 때 면진 주기를 얼마로 설정해야하는지도 알 수 있었다.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2011년도 정기 학술대회
• /
• pp.764-767
• /
• 2011

본 논문에서는 유한요소법을 이용하여 임의의 단면을 가지는 빔의 비틀림 문제를 해석 할 수 있는 방법론을 제시하였다. 빔 유한요소에서 연속적인 뒤틀림함수를 얻기 위해 각 절점에서 뒤틀림자유도를 정의한 후 빔의 길이 방향으로 보간하였다. 이러한 방법의 사용은 뒤틀림구속효과와 비선형문제에 쉽게 접근 할 수 있게 한다. 또한, 임의의 단면에 대한 뒤틀림함수는 각 단면에서 St.Venant 방정식을 유한요소법을 통해 수치적으로 계산된다. 단면에서 계산된 해는 3차원 일반 빔 요소의 변위장에 매핑된다. 위와 같은 절차를 통해 개발된 빔 유한요소를 사용하면 임의의 단면을 가진 빔 구조물을 자유/구속 뒤틀림조건에서 비틀림, 굽힘, 신축 변형이 복합적으로 고려하여 해석해 낼 수 있다. 이렇게 해석된 결과를 검증하기 위하여 사각단면과 L단면에서의 결과 값을 고찰하였다.