본 논문은 8절점 고체요소를 이용하여 항공기 충돌에 의한 원전 격납건물의 동적 거동을 분석하고 그 결과를 기술하였다. 콘크리트의 재료적 특성을 표현하기 위하여 Drucker-Prager항복기준을 바탕으로 항복면과 파괴면을 형성하였다. 이때 항복면과 파괴면은 콘크리트의 소성변형이 누적되면 가변하는 것으로 가정하였다. 철근의 재료특성은 변형도에 의존적인 탄성/점소성모델을 이용하여 표현하였다. 표준고체요소의 성능저하를 방지하기 위하여 Hughes가 제시한 B bar법을 바탕으로 변형도-변위관계 행렬을 형성하였다. 동적 시간이력해석을 수행하기 위하여 안정적인 수렴성을 가지는 암시적인 Newmark법을 도입하였다. 마지막으로 시간이력해석을 통하여 콘크리트 균열변형도의 수준과 충돌하는 항공기의 종류에 따른 격납건물의 동적거동변화를 조사하고 이를 정량적으로 기술하였다.

Research about spherical shells been applying most usually is achieved by many investigators already and generalized equation has been derived. But, existent research is limited in case that spherical shell's roof rigidity is isotropy or orthotropy, and research that consider periodicity of rigidity-distribution that can happen by doing spherical shell's roof system by lattice system is not gone entirely. The purpose of this paper is applying Galerkin method to spherical shell that model periodicity of roof rigidity distribution that appear by roof lattice form of large space structure and develop structural analysis program that formularize. Rigidity-model of this research selects that of spherical shell which has 2-way grid. In this paper, buckling-strength and deformation distribution of isotopic spherical shell and 2-way grid spherical shell obtained by developed program could confirm the reliability by comparison with result of existent research.

속이 빈 축대칭 회전체인 두꺼운 쉘의 정확한 고유진동수와 모우드형상을 결정하기 위해서 3차원적인 해석방법이 사용되었다. 이 축대칭 회전쉘의 모선을 직선으로 한정하지 않았으며, 쉘의 두께 또한 일정한 것으로 제한하지 않았다. 이 쉘의 중앙면은 임의의 곡율을 가지며, 쉘의 두께도 임의적으로 변한다. 자오선방향, 두께방향, 원주방향으로의 변위 성분인는 시간반응의 정현성(sinusoidal)과방향으로의 주기성을 지니며,와 z 방향으로는 대수다항 식의 형태로 가정되었다. 이 쉘의 변형률에너지와 운동에너지를 공식화하였으며, 진동수의 최소화를 통해 상위경계치의 진동수를 구하고 다항식의 차수를 증가시켜 엄밀해에 수렴된 진동수를 구할 수 있다. 선형적으로 두께가 변하는 두꺼운 원추형쉘과 구형쉘에 대한 예를 통하여 하위 다섯 개의 진동수에 대해서 유효 숫자 4자리까지의 정확한 수렴성연구가 이루어졌다. 이 해석 방법은 두께가 매우 두꺼운 쉘 뿐만이 아니라 얇은 쉘에도 적용이 가능하다

It is often hard to obtain analytical solutions of boundary value problems of shells. Introducing some approximations into the governing equations may allow us to get analytical solutions of boundary value problems. Instead of an analytical procedure, we can apply a numerical method to the governing equations. Since the governing equations of shells of revolution under symmetric load are expressed by ordinary differential equations, a numerical solution of ordinary differential equations is applicable to solve the equations. In this paper, the governing equations of orthotropic spherical shells under symmetric load are derived from the classical theory based on differential geometry, and the analysis is numerically carried out by computer program of Runge-Kutta methods. The numerical results are compared to the solutions of a commercial analysis program, SAP2000, and show good agreement.

본 연구에서는 복합재료 원통쉘의 정적 구조해석 결과를 신경회로망에 적용하여 원통쉘에 가해진 하중특성을 추론하였다. 신경회로망 알고리즘은 역전파 학습법의 학습율이 가변적으로 조정될 수 있도록 프로그램을 개발하였으며, 입력패턴은 원통쉘에 하중이 가해졌을 때, 원통쉘의 측면에서 발생하는 9지점의 변형률을 이용하였다 출력층은 가해진 하중특성으로 설정하였으며, 학습결과 원통쉘의 하중특성 추론 학습에 성공하였다. 은닉층의 층수를 1층에서 3층까지 학습결과를 비교분석하였으며, 하중특성은 0.5% 이내로 추론이 가능해졌다. 본 연구 결과 신경회로망을 이용한 복합재료 원통쉘의 역공학이 가능해졌다.

본 논문은 철근 콘크리트 구조물의 비탄성 해석을 수행하기 위하여 개발된 9절점 퇴화 쉘 요소에 대하여 기술하였다. 개발된 쉘 유한요소는 퇴화고체기법과 함께 구조물에서 발생하는 횡 전단 변형효과를 고려하기 위하여 Reissner-Mindlin (RM)가정을 도입하였다. RM가정을 바탕으로 한 퇴화 쉘 요소는 쉘의 두께가 얇거나, 즉 종횡비가 작거나, 균일하지 않은 유한요소망을 사용할 경우 구조물의 강성이 과대하게 계산되는 강성과대현상(Locking phenomenon)이 나타나게 된다. 강성과대현상은 선택적 감차 적분, 비 적합 모드, 가변형도 등을 사용하여 개선하는데 특히 가변형도법에 바탕을 둔 대체변형도는 많은 유한요소에 성공적으로 적용되어 왔다. 그러나 이와는 대조적으로 콘크리트의 비탄성 해석에 가변형도법을 도입하고 그 성능을 조사한 사례는 매우 적다. 따라서 본 연구에서는 가변형도법과 미시적 재료모델을 바탕으로 RM 쉘 요소를 정식화하고 미를 유한요소프로그램으로 개발하였다. 개발된 철근 콘크리트 쉘 요소의 성능은 수치예제를 통하여 검증하였다. 수치예제로부터 개발된 쉘요소를 이용하여 구해진 해석결과가 실험결과 또는 다른 해석결과에 근접함을 알 수 있다.

이 논문에서는 철근콘크리트 쉘구조의 동적해석을 위한 비선형 유한요소 해법을 제시하였다. 사용된 프로그램은 철근콘크리트 구조물의 해석을 위한 RCAHEST이다. 유한요소로서는 면내회전강성도를 갖는 4절점 평면 쉘요소가 사용되었다. 두께방향에 대한 철근과 콘크리트의 재료성질을 고려하기 위하여 층상화기법이 도입되었다. 재료적 비선형성에 대해서는 균열콘크리트에 대한 인장, 압축, 전단모델과 콘크리트 속에 있는 철근모델을 조합하여 고려하였다. 이에 대한 콘크리트의 균열모델로서는 분산균열모델을 사용하였다. 동적 평형방정식의 해는 HHT법에 의한 수치적분으로 구하였다. 신뢰성 있는 해석결과와 비교를 통하여 이 논문의 제안방법이 철근 콘트리트 쉘구조의 비선형 동적해석에 적합한 방법임을 입증하고자 한다.

보강된 판 및 쉘구조의 동적 비선형해석을 수행하기 위하여, 유한회전을 고려한 변형된 쉘유한요소를 이용하여 total Lagrangian formulation이 제시된다. 전단구속 (shear locking) 현상과 가상의 제로에너지 모우드를 동시에 제거하기 위하여 가정변형도 개념을 채용한다. 탄소성해석에서는 return mapping 미해rithm이 쉘구조의 붕괴 해석에 적용된다. Newmark 직접적분법을 사용하여 동하중 및 지진하중을 받는 쉘구조의 동적 비선형해석 결과를 제시한다.

기존의 Mindlin이론을 바탕으로 한 감절점 쉘 요소는 두께가 두꺼운 쉘 구조를 해석할 경우 매우 좋은 결과가 얻어지는 것으로 알려져 있다. 그러나 두께가 얇은 쉘의 경우 전단구속 및 막구속 현상으로 인해 요소의 강성이 과대하게 평가되는 것으로 나타났다. 이러한 문제를 해결하기 위해 본 연구에서는 서로 다른 두 가지 막 및 전단 구속 해결방법을 조합하여 쉘 요소에 적용하였다. 첫 번째 요소는 전단 변형률에 가정된 전단변형률을 적용하였고, 막 변형률에는 감차적분법을 적용하였다. 두 번째 요소의 경우, 막 변형률과 전단 변형률 모두에 가정된 전단 변형률을 적용하여 쉘 구조해석을 수행하였다. 여러 가지 수치 해석을 수행한 결과, 두 가지 쉘 요소가 빠른 수렴성을 가지고 있는 것으로 나타났으며, 막 구속과 전단 구속 현상이 제거된 것을 나타났다.

스위프 기하학적 모델은 곡선, 면 또는 입방체를 주어진 경로를 따라 이동시킴으로써 기하학적 모델을 생성하는 기법이다. 따라서 스위핑을 사용하면 프리즘 쉘의 곡면을 쉽게 정의할 수 있다. 본 논문은 스위프 기하학적 모델을 프리즘 쉘의 최적화에 적용하는 절차에 대하여 기술하였다. 제시한 스위프 기하학적 모델을 유한요소법과 융합하였고 프리즘 쉘의 반응을 계산하기 위해 9절점 퇴화쉘요소를 채용하였다. 본 연구에서 제시한 최적화과정을 증명하기 위하여 수치예제를 풀어 보았다. 수치예제를 통하여 제시한 스위프 기하학적 모델이 많은 종류의 프리즘 쉘을 최적화하는데 효율적이고 신뢰적인 방법인 것으로 나타났다.

보강된 판 및 쉘구조의 기하학적 비선형해석을 수행하기 위하여, total lagrangian formulation에 근거한 증분 평형방정식을 적용하고, 강도행렬 산정시 회전각의 2차항을 포함시켜 기하학적 비선형 해석시 해의 수렴성을 향상시켰으며, 보강된 쉘 구조의 해석시 보강재를 쉘 요소로 모델링하고 주부재와 보강재의 연결점에서 일반적인 변환관계를 이용하였다. 등매개 쉘 유한요소의 단점인 locking 현상을 극복하기 위하여 가정 변형률장을 적용하여 감차적분 또는 선택적분시 나타날 수 있는 제로 에너지 모드를 제거하였다. 수치해석 예제를 통하여 가정 변형률장에 근거한 쉘유한요소에 대한 효율성 및 적용성을 확인하였다.

최 등1)은 total lagrangian formulation에 근거한 증분 평형방정식을 적용하고, 강도행렬 산정시 회전각의 2차항을 포함시켜 기하학적 비선형 해석시 해의 수렴성을 향상시켰다. 또한 등매개 쉘 유한요소의 단점인 전단구속 현상과 제로 에너지 모드가 발생하는 문제를 극복하기 위하여 가정 변형률장을 적용하여 보강된 판 및 쉘 구조의 비선형 해석법을 개발하였다. 본 연구에서는 잔류응력을 고려한 쉘구조의 극한강도 해석을 수행하기 위하여, 대변형거동과 함께 소성붕괴거동을 추적할 수 있는 알고리즘을 제시한다. 잔류응력을 고려한 증분평형방정식에 return mapping algorithm을 이용한 탄소성 해석법을 결합시켜서 보강된 판 및 쉘구조의 극한거동을 파악한다. 수치해석 예제를 통하여 본 연구에서 제시된 유한요소 및 비선형 해석 알고리즘에 대한 효율성 및 적용성을 확인하였다.

비축대칭 하중을 받는 축대칭 쉘의 해석시 구조의 축대칭성을 고려하면 시간과 노력을 절약할 수 있다. 하중과 변위에 대하여 원주방향으로 Fourier 급수전개를 고려함으로서 비축대칭하중을 받는 축대칭 쉘의 해석은 뼈대요소처럼 취급할 수 있다. 본 논문에서는 Fourier 급수전개를 이용한 통상의 유한요소법에 의하여 비축대칭 하중을 받는 원형탱크의 강성행렬을 유도하고, 이 강성행렬을 행렬의 조작에 의해 전달행렬로 전환하여 전달행렬법을 적용하였다. 이 논문은 비축대칭하중을 받는 축대칭 쉘의 해석을 위한 연립방정식의 수를 최소화하는데 그 목적이 있다. 제안된 방법에 의한 풍하중과 물하중을 받는 원형탱크의 해석결과는 타 방법에 의한 해석결과와 잘 일치하고 있다.

자오방향 및 주변방향으로 피르스트레스트 하중이 작용된 축대칭 쉘 구조물을 기하학적으로 축대칭인 구조물의 특성을 최대한으로 이용할 수 있도록 회전 링요소로 모델화하였다 보강링 요소의 모델은 축대칭 쉘요소를 이용하였으며 본체 구조물과 절점에서 부착되있는 것으로 가정하여 이의 편심을 고려하였다 유체-구조물의 상호관계는 접촉면에서 구조물의 가속도에 비례한 부가질량으로 표현하였으며 부가질량은 유체를 비점성 비압축 및 비회전을 가정하여 유한요소법에 의해 구하였다 이에 대한 수치해석을 통하여 고유진동해석 및 지진하중을 주하중으로 한 동적해석을 실시하였다 프로그램을 통하여 해석한 결과를 프리스트레스 하중 하에서 고유진동수에 대한 정해와 비교한 결과 20개의 요소로 모델링한 경우에서도 정해와 근접한 해를 얻을 수 있었다 또한 내부유체가 있는 경우와 링보강을 한 경우에 대한 고유진동수를 문헌과 비교한 결과 근접한 해를 얻을 수 있었다.

본 논문에서는 반복하중을 받는 철근 콘크리트 쉘구조물의 해석을 위한 비선형 유한요소 해법을 제시하였다 유한 요소로서는 충상화기법을 이용한 부재회전강성도를 갖는 4절점 평면 쉘요소가 개발되었다 두께 방향에 대한 철근과 콘크리트의 재료성질을 고려하기 위하여 충상화기법이 도입되었다. 재료적 비선형성에 대해서는 균열콘크리트에 대한 인장, 압축, 전단모델과 콘크리트중에 있는 철근모델을 조합하여 고려하였다. 이에 대한 콘크리트의 균열모델로서는 분산균열모델을 사용하였으며 철근에 대해서는 1축 응력상태로가정하여 등가의 분산 분포된 철근량으로 모델화하였다 구성모델은 재하, 제하 그리고 재재하과정을 포함하여 요소는 반복하중하에서 철근콘크리트 쉘의 거동을 파악할 수 있다 신뢰성 있는 실험결과와 비교를 통하여 본 논문의 해석방법이 반복하중을 받는 철근콘크리트 쉘구조의 비선형 해석에 적합한 방법임을 입증하고자 한다.

알루미노실리케이트 계 졸이 지르콘 쉘 몰드의 뮬라이트 층 생성에 미치는 영향을 고찰하였다. 알루미노실리케이트 졸은 콜로이달 실리카와 수용성 질산알루미늄을 혼합하여 제조하였으며 50˚C에서 48시간 조건에서 겔화 하였다. 이러한 겔은 깁사이트 및 알루미노실리케이트 복합 겔로 구성되어 있었으며, Si 이온과 결합하는 모든 Al 이온의 배위수가 4임을 확인하였다. 뮬라이트 상은 1300˚C 이상에서 소결하였을 때 관찰되었으며 뮬라이트의 XRD 피크는 소결 온도와 질산알루미늄의 농도가 증가할수록 노실리케이트 졸 슬러리를 2차 층에 코팅하였다. 그 결과 1차 및 3차 층의 분리가 일어났으며, 이는 소결시 졸의 1차 및 3차 층으로 침윤 발생과 잔류 실리카와 뮬라이트 간의 열팽창계수 차이에 기인한 것으로 판단된다.

본 논문에서는 쉘 최적화에 대한 연구 결과를 기술하였다. 본 연구의 주목적은 쉘 구조물의 최적형상과 두께 분포를 찾는데 있다. 쉘의 변형에너지를 목적함수로 사용하고 초기 쉘의 부피를 제약조건을 고려하였다. 본 연구에서는 Computer-Aided Geometric Design (CAGD) 기법을 이용하여 쉘의 형상과 그 두께 분포를 표현하였고 쉘의 변형에너지를 측정하기 위해서 가변형 도를 채용한 퇴화 쉘 요소(Degenerated Shell Element)를 도입하였다. 최적 값을 구하기 위해서 세 가지 수학적 프로그래밍 기법을 제공하는 프로그램 DOT를 사용하였다. 마지막으로 새로이 개발된 쉘 최적화시스템의 효율성을 최적화예제로써 증명하였다.

지르콘 쉘 몰드를 제조하여 소결 온도 및 시간에 따른 소결거동이 기계적특성에 미치는 영향을 관찰하였다. 1000˚C에서 1.5시간동안 소결 한 쉘 몰드의 1차 코팅 표면층에서 미세균열의 크기 및 수가 극대화되었으며 소결온도 및 시간이 증가할수록 상온강도와 기공의 비표면은 감소하였다. 1차 소결 후 1500˚C에서 4시간동안 소결 처리한 주형의 고온변형거동은 백업층과 하중의 역방향으로 스터코와 지르콘 슬러리의 경계면을 따라 변형이 발생하였다. 1차 코팅층에서 알루미나 스터코와 지르콘 슬러리의 열팽창계수 차이와 백업 코팅층에서 지르콘 슬러리간의 입도 차이가 고온변형시험 중 역변형을 일으킨 것으로 판단된다.

압력절점은 요소의 균등한 압력증분을 1개의 자유도로 갖는 절점이며, 유한요소의 하중-변위 평형방정식에 체적과 압력의 관계를 추가하여 한계압력 이후에서도 체적변화에 따른 압력증분을 직접적으로 제저할 수 있는 절점이다. 본 연구에서는 철근콘크리트의 평면 구성 방정식과 적층정식화에 적용한 쉘 요소에 압력절점을 추가하고 해석시 체적을 제어함으로써 철근콘크리트 원통형 구조에 대해 파괴까지의 극한내압 능력을 해석할 수 있는 체적제어 비선형 해석기법을 개발하였다. 본 논문에서 제안한 해석기법을 이용하여 철근콘크리트 원통형 구조물에 대하여 비선형 해석을 수행하여 한계압력과 한계압력 이후의 구조물의 거동을 예측하였으며 실험결과와 비교 검증하였다.

본 연구에서는 대변형 쉘 구조물에 효과적인 적응적 유한요소 자동생성 기법을 제안한다. 사후 오차평가에 기초하여 기하학적 비선형 해석시 각 하중 단계에서의 요소 재생성에 초점을 맞추고 있다. 응력오차로부터 얻어진 요소크기 함수로 등고선을 구성하고, 요소 재생성 기법으로 advancing front method의 일종인 패이빙법(paving method)을 이용하여 적응적 요소 자동생성을 수행한 결과, 그 유용성을 확인하였다.