쉘형 구조물의 동적 불안정 문제를 다룬 연구결과는 다소 발표되고 있으나 위상곡면을 이용하여 카오스 생성에 관한 기본적 현상을 다룬 연구는 거의 없는 실정이다. 동적 비선형 문제에서 여러 가지 초기조건에 의해 불안정 현상이 민감하게 발생하는 이유를 파악하기 위해서는 위상곡면에서 끌개의 특성을 조사하여 동적 불안정 생성과정을 검토하는 일은 매우 중요하다. 본 연구에서는 얕은 EP 쉘이 스텝하중을 받을 때, 직접 좌굴과 간접 좌굴의 발생 경로를 파악하기 위하여 Galerkin 법에 의해 전개된 이산화 방정식을 구한다. 이를 수치해석 기법으로 위상곡선과 연속응답스펙트럼을 구해 동적 불안정 특성을 규명한다.

쉘형 구조물은 외력에 대해 효과적으로 저항할 수 있어, 두께를 얇게하면서 대공간 구조물을 만들 수 있는 장점이 있다. 이러한 장점은 구조형태에 크게 의존한다. 그러므로 많은 설계자들은 최적 형태를 설계에 반영하고자 하지만, 이는 간단치 않다. 지금까지 보다 최적의 형태를 얻기 위한 많은 기법들이 발표되어 왔고, 이들은 각각 장단점을 나타낸다. 본 논문에서는 최적의 곡면을 얻기 위해 기하학적 비선형을 고려한 유한요소법을 이용한 비교적 간단한 방법을 제안한다. 이러한 방법을 이용하여 다양한 쉘형 곡면에 적용하여 최적곡면을 얻고, 이를 비교한다.

PSC 박스 교량은 콘크리트, 철근과 텐던으로 구성된 구조물로서 콘크리트의 인장 균열, 철근의 비선형 거동 등 재료의 비선형성 거동 특성 및 콘크리트의 시간 의존적 특성을 가지고 있는 복합 구조물이다. PSC 박스 교량의 시공 중 거동 특성을 고려하기 위하여 뼈대 요소(프레임 요소)를 이용한 시공단계의 설계가 수행되고 있다. 그러나 PSC 박스 교량 중 곡선램프교 등의 경우는 교량의 외측 및 내측의 변위 및 응력 값이 현저히 다르다. 따라서 PSC 박스 교량의 텐던량 및 시공 중 긴장력이 외측 및 내측에서 다르게 산정되어야 함에도 불구하고 현실적으로는 계산이 불가능하여 같은 양의 텐던과 부적절한 긴장력을 사용하고 있어 시공 중 항상 안전사고에 노출되고 있다. 이러한 단점을 해결하기 위하여 3차원 해석이 필수적으로 요구되고 있으며 본 연구에서는 PSC 박스 교량의 해석 기법에 필요한 준 적합 쉘 요소를 제안하고자 한다.

본 논문에서는 적응적 h-유한요소 세분화에 의한 박스형 절판 구조물의 선형좌굴 유한요소해석법을 제안한다. 면내회전 자유도를 갖는 변절점 평판쉘유한요소를 사용하여 유한요소의 거동을 개선하고 6자유도를 갖는 다른 유한요소와의 자유도의 연결을 용이하게 한다. 이와 같이 개발된 평판쉘유한요소에 의하여 박스형 절판구조물의 정확한 구조해석이 가능한데, 변절점유한요소를 정식화함으로써 적응적 h-유한요소 세분화시에 발생하는 다른 패턴의 사각형 유한요소 세분화망의 연결을 용이하게 해결한다. 오차평가에 대한 개선된 응력장을 얻기 위하여 상위수렴 조각회복법을 적용한다. 이와 같이 상위수렴 조각회복법에 의한 개선된 응력장에 의하여 구성된 유한요소 세분화망을 이용하여 좌굴하중과 좌굴모드를 자동적으로 구할 수 있도록 한다.

본 연구에서는 NURBS곡면식을 바탕으로 하는 곡면 모델링과 쉘 유한요소해석의 효율적인 연동체계를 개발하고자 한다. 기하학적으로 정확한 쉘 유한요소해석에서 정확한 기하량의 계산은 필수적이며, 따라서 곡면을 표현하는 일반적인 방법인 NURBS곡면식으로 부터 필요한 기하량을 직접 계산한다면 보간에 의해 발생할 수 있는 기하학적 오차를 줄임으로써 해의 수렴성을 높일 수 있다. 아울러 기하학적으로 정확한 쉘 유한요소를 일반적인 곡면에 적용하기 힘들었던 한계점을 극복하여 수학적으로 표현 가능한 단순한 곡면들뿐만 아니라, NURB곡면식으로 표현 가능한 일반적인 곡면의 해석이 가능하게 되어 적용범위를 확장할 수 있다. 본 연구에서는 곡면을 생성함에 있어 주어진 데이터 점들을 보간하여 NURBS곡면을 생성하는데, 이러한 데이터들은 일반적으로 곡면의 스캐닝을 통해 얻을 수 있다. 곡면을 보간하여 NURBS곡면을 생성하는 과정에서 사용되는 매개변수 정의방식에 때라 생성된 곡면의 정확성이 차이를 보이므로 곡면의 형상에 따라 적합한 방식을 사용하여 곡면을 보간 할 필요가 있다. 몇 가지 잘 알려진 수치예제를 통하여 개발된 연동체계의 성능과 정확성을 검증하고 그 결과를 비교 분석하였다.

물에 잠긴 다공 원통형 쉘의 경우 물에 잠긴 상태로 유한요소해석을 하기에는 거의 불가능하므로 등가물성치를 사용하여야 한다. 다공 평판의 경우 이에 대한 등가물성치를 ASME 코드에서 제시하고 있지만, 다공 원통형 쉘의 등가물성치에 대한 연구는 아직까지 수행된 적이 없다. 따라서 본 연구에서는 유한요소해석을 이용하여 다공 원통형 쉘의 동적 해석에 이용할 수 있는 등가물성치를 제안하였고 그 타당성을 검증하였다.

격납건물은 원자로 사고발생시 방사능물질의 외부 유출을 막는 최후의 방벽이므로 가동 중 원전의 격납건물에 대한 안전성평가는 반드시 수행되어야 된다. 이러한 맥락에서 이 논문은 원전 격납건물의 비선형해석을 위해 탄소성 모델을 바탕으로 개발된 8절점 가변형도 쉘 요소와 이를 이용한 구조물의 비선형해석에 대하여 기술하였다. 비선형해석을 위해 콘크리트의 압축거동에 Drucker-Prager 파괴기준을 적용하였고 파괴포락선의 형상을 결정짓는 재료매개변수는 이축응력 실험으로부터 도출하였다. 개발된 쉘 유한요소는 퇴화 고체기법과 횡 전단변형도를 고려하기 위하여 Reissner-Mindlin(RM)가정을 도입하였고 쉘의 두께가 얇거나, 즉 종횡비가 작거나, 균일하지 않은 유한요소망을 사용할 경우 구조물의 강성이 과대하게 평가되는 묶임현상(locking phenomenon)을 제거하기 위해 본 논문에서는 가변형도법을 도입하였다. 개발된 철근콘크리트 쉘 요소의 성능검증을 위해서 벤치마크 테스트를 수행하였고 그 결과 이 논문에서 도출한 유한요소해석 결과는 실험결과와 잘 일치 하였다

Researches on spherical shell which is most usually applied have been completed by many investigators already and generalized numerical formula was derived. But the existent researches are limited to those on spherical shell with isotropic or orthotropic roof stiffness, periodic distribution of roof stiffness that can be caused by spherical and latticed roof system is not considered. Therefore, the object of this study is to develop a structural analysis program to analyze spherical shells that have periodicity of roof stiffness distribution caused by latticed roof of large space structure, grasp buckling characteristics and behavior of structure.

등방성 혹은 비등방성 적층복합판 및 쉘의 선형 정적 문제와 자유진동 해석이 새로운 변형률 변위 관계가 도입된 개선된 9절점 쉘 요소에 의하여 수행되었다. 그 관계에서 새롭게 추가된 휨 변형률과 변위사이의 관계 항들에 의한 효과는 비틀어진 보 문제에서 검토되었다. 정식화의 전 과정을 통해, 식들의 모든 항들은 자연 좌표계에 기초하고 있다. 가정 자연 변형률 방법이 막 잠김과 전단 잠김 거동을 제거하기 위하여 사용하였다. 적층 복합판 및 쉘의 고유치의 계산을 위해 Lanczos방법을 사용하였고 질량행렬을 구성하기 위하여 Gauss적분법을 사용하였다. 정식화의 유효성을 평가하기 위해 수치 예제를 해석적 해와 비교하였으며, 제시된 결과는 자유진동 조건하에서 적층체의 거동을 이해하는데 유용할 것이다.

본 연구의 목적은 축대칭 하중을 받는 원통형 펄의 엄밀해를 구하는데 있어서, 간략하면서도 엄밀한 해를 구하는 방법을 제시하고자 하는데 있다. 이는 임의 형상의 구조해석을 위한 강력한 도구이긴 하지만 여전히 근사해석인 유한요소법에 대체될 수 있을 것이다. 이를 위해 본 논문은 반복법의 일종인 인도행렬법을 이용한 절점역계의 분배방식을 사용하였다. 원통형 쉘의 분배와 전달인자는 한성지반상의 보에 대한 미분방정식으로부터 구해진 것이다. 이러한 방법을 축대칭 집중하중과 정수압을 받는 원통형 쉘에 각각 적용해 보았고, 그 결과는 BEF 이론해와 비교할 때 만족할 만 하였다.

This paper is to study the energy absorption characteristics of CF/Epoxy(Carbon Fiber/Epoxy Resin) laminated shell with the various curvatures subjected to transverse impact loadings under the low impact velocity in consideration of design of structural members for use of transportation machine, which are consisted of the characteristics of high stiffness, strength and lightweight. The curvature radius are associated with the energy absorption characteristics of CF/Epoxy laminated shell which is brittleness material. In all tests, maximum load of CF/Epoxy laminated plate is higher than that of laminated shell with curvature, but maximum deflection is lower. And then absorbed energy of laminated shell with curvature is higher than laminated plate(curvature radius is unlimited), As curvature radius is increased, the absorbed energy is increased in laminated shell with curvature.

두꺼운 축대칭 쌍곡형 쉘의 고유진동수를 결정하는 3차원 해석법이 제시되었다. 수학적으로 2차원적인 전통적인 쉘 이론과는 달리, 본 연구의 해석법은 3차원적인 동탄성방정식을 근간으로 하였다. 반경방향, 원주방향, 축방향으로의 변위성분인 u, u, uz를 시간에 대해서는 정현적으로, 에 대해서는 주기적으로, r과 z방향으로는 대수 다항식으로 표현하였다. 쌍곡형 쉘의 위치(변형률)에너지와 운동에너지를 정식화하고 리츠법을 사용하여 고유치문제를 해결하였으며, 진동수의 최소화과정을 통해 고유진동수를 엄밀해의 상위경계치로 구하였다. 대수 다항식의 차수가 증가하면 진동수는 엄밀해에 수렴하게 된다. 축대칭 쌍곡형 쉘의 하위 5개의 진동수에 대해서 유효숫자 4자리까지의 수렴성 연구가 이루어졌다. 쌍곡형 쉘의 서로 다른 2개의 두께 비, 3개 의 축비(axis ratio), 3개의 shv이 비를 가진 총 18개의 형상을 지닌 자유 경계의 축대칭 쌍곡형 쉘의 수치결과를 도표화하였다. 프와송 비는 0.3으로 고정하였다. 본 연구의 해석법은 매우 두꺼운 쉘 뿐만 아니라 얇은 쉘에도 적용이 가능하다.

In this work, a finite element model is presented for geometrically non-linear analysis of shell structures. Finite element by using a three-node flat triangular shell element is formulated. The non-linear incremental equilibrium equations are formulated by using an updated Lagrangian formulation and the solutions are obtained with the incremental/iterative Newton-Raphson method and arc length method. Some of results are presented for shell structures. The obtained results are in good agreement with the results available in existing literature.

임의의 경계조건과 변두께를 갖는 축대칭 반구형 쉘과 반구형체의 진동수와 모우드형상을 결정하는 3차원적 해석법이 소개되었다. 수학적으로 2차원적인 전통적인 쉘이론과는 달리 본 연구의 해석법은 3차원 동적 탄성방정식을 사용하였다 자오선방향 (Φ), 법선방향(z), 원주방향()으로의 변위성분인 Φ, z, 는 시간에 대해서는 정현적으로, 에 대해서는 주기적으로, 와 z 방향에 대해서는 대수다항식으로 표현될 수 있다. 축대칭 반구형 쉘의 변형률 에너지와 운동 에너지를 정식화하고, 리츠법으로 고유치문제를 계산하였다. 진동수의 최소화과정을 통해 엄밀해의 상위 경계치 진동수를 구하였으며, 이 때, 다항식의 차수를 증가시키면 진동수는 엄밀해에 수렴하게 된다. 자오선방향으로 선형적으로 꿩 두께가 변하는 반구형 쉘과 반구형체치 3차원적 진동수를 최초로 계산하였으며, 축방향으로 난 조그만 원추형 구멍이 진동수에 미치는 영향도 분석하였다. 상두께와 자유경계조건을 갖는 두꺼운 축대칭 반구형 쉘에 대한 3차원적 리츠해와 3차원적 유한요소법에 의한 진동수를 서로 비교하였다.

본 연구에서는 축대칭 쉘구조물의 정동적해석을 효과적으로 수행할 수 있는 새로운 유한요소를 제안하였다. 본 유한요소는 축대칭 쉘의 전단변형률을 고려하였으며, 쉘의 경계에서 기술할 수 있는 변수들만으로 표현되는 효율적인 형태의 수정된 혼합 변분이론에 바탕하여 유한요소정식화를 수행하였다. 또한, 변위장에 대해 무절점 자유도를 추가적으로 도입하여 요소의 수치적 성능을 크게 향상시켰다 계산의 효율성을 위해, 요소정식화의 최종단계에서 정치조건으로부터 응력매개변수들을 제거하고, 동적축약을 통하여 무절점 자유도 성분들 또한 최종적인 유한요소방정식에서 제거되어짐으로써, 일반적인 변위기저 요소와 같은 크기의 유한요소방정식을 얻을 수 있다. 몇 가지 수치예제들에 대한 해석을 통하여, 무절점 자유도와 변위장에 일치하는 적절한 응력매개변수를 가지는 제안된 혼합 축대칭 쉘요소가 정동적해석에서 대단히 정확하고 효율적임을 확인할 수 있었다

본 논문은 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.