![The Korean Society for Power System Engineering](/images/basic/pubreader_back.gif)
유한요소-전달강성계수법을 이용한 다양한 형상을 갖는 축대칭 셸의 정적해석
초록
축대칭 셸의 정적해석은 공학 분야에서 중요한 주제 중의 하나이다. 저자들은 유한요소법의 모델링 기법과 전달강성계수법의 전달 기법을 조합한 유한요소-전달강성계수법을 이용하여 다양한 형상을 갖는 축대칭 셸의 정적해석을 수행하기 위한 전산 알고리즘을 제안한다. 저자들은 유한요소-전달강성계수법으로 다양한 형상을 갖는 축대칭 셸의 정적해석을 수행하기위한 전산알고리즘을 정식화한 후, 전산 프로그램을 개발하였다. 저자들은 유한요소-전달강성계수법과 유한요소법으로 3가지 계산모델(원형 돔, 물탱크, 압력용기)에 대하여 정적해석을 수행하였다. 양 방법의 계산 결과의 비교를 통해, 유한요소-전달강성계수법은 다양한 형상을 갖는 축대칭 셸의 정적해석에 이용될 수 있는 신뢰성 있는 해석 기법임을 확인하였다.
Abstract
Static analysis of axisymmetric shells is one of the important topics in engineering. The authors suggest a computational algorithm for carrying out the static analysis of axisymmetric shells with various shapes by using the finite element-transfer stiffness coefficient method (FE-TSCM) which combines both the modeling technique of the finite element method (FEM) and the transfer technique of the transfer stiffness coefficient method. In this study, the authors developed the computational program after computational algorithm was formulated by the FE-TSCM for the static analysis of the axisymmetric shell with various shapes. The authors carried out the static analysis about three computational models, which were a spherical shell cap, a water tank, and a pressure vessel by using the FE-TSCM and the FEM. From comparing the computational results of both methods, it was confirmed that the FE-TSCM is a reliable analysis method for static analysis of axisymmetric shells with various shapes.
Keywords:
Static Analysis, Axisymetric Shell, Finite Element, Transfer Stiffness Coefficient Method키워드:
정적해석, 축대칭 셸, 유한요소, 전달강성계수법1. 서 론
두께가 얇고 굽은 판인 셸(shell)로 구성된 구조물은 가벼우면서 튼튼하여 기계, 항공, 해양, 토목, 원자력 등 다양한 산업분야에서 널리 사용되고 있어 공학적으로 중요한 구조물 중 하나이다. 셸 구조물의 설계에서 가장 기본이 되는 것이 하중에 대한 셸 구조물의 변위나 응력을 계산하는 정적해석(static analysis)이다.1,2)
셸 구조물의 정적해석을 수행하기 위해서는 모델링 과정이 선행되어야 하는데, 여기에 솔리드 요소, 셸 요소, 축대칭 요소 등이 사용될 수 있다.(3-6) 솔리드 요소는 3차원 캐드 소프트웨어로 만든 형상에 사용자가 쉽게 적용할 수는 있지만 계산결과의 신뢰성을 확보하기 위해서는 엄청나게 많은 요소로 모델링해야 하므로 계산 시간이나 컴퓨터 메모리 측면에서 취약점이 있다. 이러한 문제점을 해결하기 위하여 쉘 요소를 사용하면 해석에 필요한 전체 자유도 수를 크게 줄일 수 있다. 그러나 쉘 구조물이 축대칭 형상과 하중을 갖는 경우에는 축대칭 요소를 사용하는 것이 자유도를 대폭 줄일 수 있어 아주 효과적이다.
저자들은 축대칭 셸 구조물의 정적해석 알고리즘의 개발을 위해 오랜 기간 동안 연구를 해 왔고, 그 결과 단순한 원통형 셸과 원추형 셸의 정적해석을 위한 전산 알고리즘과 프로그램을 개발하여 제안한 바 있다.7,8) 그러나 과거의 연구는 특정 원통형 셸과 원추형 셸에만 적용이 제한되어 원형 돔, 구형 셸, 반구형 셸, 원통-원추-반구형 조립 셸 등 다양한 형상을 갖는 축대칭 셸 구조물에는 적용이 불가능하였다. 따라서 이 연구의 목적은 단순한 원통형 또는 원추형 셸과 같은 직선형 셸뿐만 아니라 원형 돔, 반구형 셸과 같은 곡선형 셸을 포함한 다양한 형상을 갖는 축대칭 셸의 정적해석을 수행할 수 있는 전산 알고리즘 및 프로그램을 개발하는 것이다.
2. 정적해석 알고리즘
2.1 모델링
이 연구에서는 Fig. 1과 같은 원통, 원추, 반구형 셸이 조립된 축대칭 셸 구조물을 Fig. 2와 같은 절두 원추요소(truncated cone element)의 직렬 결합으로 고려한다. Fig. 1과 같이 축대칭 셸 구조물 전체를 n개의 절두 원추요소로 분할할 경우, 원추요소의 좌우측에는 분할 원이 존재하는데, 이 원을 링 절점(ring node) 혹은 간략히 절점(node)이라 하고, 좌단에서 우단으로 총 n+1개의 절점을 순차적으로 절점 1, 절점 2, ⋯, 절점 n+1로 부른다. 각 절점은 원추요소의 자오선(xi축)방향과 이에 수직인 법선(yi축)방향 그리고 회전(θ)방향의 3자유도를 갖는다. 한편, 원추요소는 셸의 자오선방향으로 직선이므로 곡선 셸을 모델링할 경우 분할수를 증가할수록 해석 결과가 정확해진다.
구조물이 탄성지지 된 경우, 기초와 절점 사이에 3개의 기초지지 스프링을 삽입하고, 스프링상수 값을 통해 모델링한다. 구조물 좌우단의 경계조건은 절점 1과 절점 n+1의 기초지지 스프링을 통해 모델링한다. 절점에는 하중 또는 모멘트가 직접 작용할 수 있고, 원추요소에 작용하는 압력을 통해 외력을 정의할 수 있다.
유한요소해석을 통해 Fig. 2와 같은 i번째 축대칭 절두 원추요소의 강성행렬은
(1) |
와 같고, i번째 원추요소에 작용하는 압력의 등가절점하중은 다음 식과 같다(부록 참고).3)
(2) |
2.2 강성계수행렬 및 힘보정벡터의 전달
절점 i의 좌측과 우측에서 힘벡터와 변위벡터사이의 관계를 식 (3)과 (4)로 각각 정의한다.
(3) |
(4) |
Fig. 3처럼 원추요소 i가 원추요소 i-1보다 βi만큼 반시계 방향으로 기울어져 있으면
(5) |
가 되고, 여기서 좌표변환행렬
(6) |
이고, 다음과 같은 성질을 갖는다.
(7) |
식 (4), (5), (7)로부터 다음 식과 같은 강성계수행렬 및 힘보정벡터의 좌표 변환식을 유도할 수 있다.
(8) |
기초지지 스프링을 갖는 절점 i에 집중하중이나 압력에 의한 등가집중하중이 작용할 경우, 절점 i에서 힘의 평형식은
(9) |
가 되고, 여기서
(10) |
식 (3), (4), (9)로부터 절점 i 좌우측의 강성계수행렬 및 힘보정벡터 사이의 관계(포인트 전달식)를 구하면 다음과 같다.
(11) |
일반적으로 유한요소법에서 사용하는 힘벡터의 정방향과 전달강성계수법에서 사용하는 힘벡터의 정방향이 다르므로, 식 (1)을 변형하면 다음과 같다.
(12) |
식 (3), (4), (12)로부터 원추요소 i 좌우측의 강성계수행렬 및 힘보정벡터 사이의 관계(필드 전달식)를 구하면 다음과 같다.
(13) |
여기서
(14) |
2.3 절점 변위 및 힘벡터 계산
구조물의 좌단에 해당하는 절점 1에는 기초지지 스프링이 존재할 수 있으나 절점 1의 좌측에 원추 요소는 존재하지 않는다. 따라서
(15) |
을 구하고, 필드 전달식인 식 (13)으로부터 와 를 구하고, 좌표 변환식인 식 (8)로부터 와 를 구한다. 그리고 포인트 전달식인 식 (11)로부터 와 를 구할 수 있다. 반복적으로 상기 좌표 변환식과 포인트 전달식 그리고 필드 전달식을 적용하면 최종적으로 마지막 절점인 절점 n+1 우측의 와 도 유도할 수 있다.
우단의 경계조건도 스프링으로 모델링하였으므로 마지막 절점인 절점 n+1의 우측은 자유단으로 고려할 수 있고, 그 결과 은 영벡터가 된다. 따라서 식 (4)로부터
(16) |
가 된다. 식 (3), (4), (12)로부터 식 (17)을 유도할 수 있는데, 이 식은 원추요소 좌우측의 변위벡터 사이의 관계식이다.
(17) |
따라서 식 (16)으로 마지막 절점의 변위벡터를 구한 후, 식 (17)과 식 (5)를 반복적으로 적용하면 모든 절점의 변위벡터를 구할 수 있다. 절점 좌우측의 힘벡터는 절점의 변위벡터를 이용하여 식 (3), (4)로부터 구할 수 있다.
3. 해석 결과
3.1 균일한 분포하중을 받는 원형 돔
Fig. 4에 나타낸 첫 번째 계산모델은 중심각이 70°인 원형 돔(dome)이고, 상부에 균일한 분포하중을 받고 있으며 경계조건은 고정이다. 돔의 반경은 2.286 m, 두께는 7.62 cm, 탄성계수는 20.68 GPa, 푸아송비는 1/6, 분포하중은 6.894 kPa이다.
계산모델은 돔의 중앙부를 기준으로 좌우 대칭이므로 10개의 축대칭 원추요소로 모델링하였고, 돔의 중앙부를 절점 1, 고정점을 절점 11로 한다. 계산모델의 대칭성을 고려하여 절점 1의 변위 와 각변위 θ1을 추가로 구속하여 계산하였다.
유한요소-전달강성계수법으로 각 절점의 변위를 계산한 결과가 Table 1이다. 유한요소법에 근거하여 작성된 Ross3)의 전산 프로그램으로 동일 모델에 대하여 계산한 결과가 Table 2이다. 양 방법에 의한 계산 결과는 아주 잘 일치하였다.
Table 3은 돔의 9번째와 10번째 요소 내외부의 자오선방향 응력을 유한요소-전달강성계수법과 유한요소법으로 계산한 결과이고, Table 4는 원주방향 응력(hoop stress)을 계산한 결과이다. 양 방법의 결과가 아주 잘 일치하였다.
3.2 원통-원추형 결합 셸 형태의 물탱크
Fig. 5에 나타낸 두 번째 계산모델은 내경 5 m, 높이 10 m의 원통형 셸과 높이 5 m의 원추형 셸이 결합된 구조물이고, 내부에는 10 m 높이의 물이 저장된 물탱크이다. 경계조건은 구조물의 하부가 고정이고, 셸의 두께는 1 cm, 탄성계수는 200 GPa, 푸아송비는 0.3이다.
저장된 물의 정수압을 하중으로 고려하고, 원통형 셸은 10개, 원추형 셸은 5개의 축대칭 원추요소로 모델링한다. 절점 1은 물탱크의 가장 하부에 해당되고, 절점 11은 원통형 셸과 원추형 셸의 결합부 그리고 절점 16은 가장 상부에 해당된다. 대칭성을 고려하여 절점 16의 변위 와 각변위 θ16을 추가로 구속하여 계산하였다.
유한요소-전달강성계수법으로 홀수 절점의 변위를 계산한 결과를 Table 5에 나타내었다. 유한요소법에 근거하여 작성된 Ross3)의 전산 프로그램으로 동일 모델에 대하여 계산한 결과를 Table 6에 나타내었다. 양 방법의 계산 결과를 비교해 보면 대체로 잘 일치하고 있음을 알 수 있었다.
Fig. 6은 계산 결과인 각 절점의 변위와 구조물의 대칭성을 이용하여 3차원 형태로 나타낸 물탱크의 변형된 형상이다. 변형 결과를 확연히 보여주기 위해 각 절점의 변위를 1,000배 확대하였고, 색상과 격자는 3차원 형상을 명확히 표현하기 위해 제공된 것 외에는 특별한 의미가 없다.
3.3 내압을 받는 압력 용기
Fig. 7에 나타낸 세 번째 계산모델은 좌우에 2개의 반구형 셸과 2개의 원통형 셸 그리고 하나의 원추형 셸이 결합된 구조물로서 내부에서 일정한 압력을 받는 압력용기이다. 이 계산모델의 전체 길이는 5.6 m, 두께는 1 cm, 탄성계수는 210 GPa, 푸아송비는 0.3이다. 경계조건은 좌단이 고정이고, 내부 게이지 압력은 1 MPa이다.
좌측부터 반구형 셸은 10개, 첫 번째 원통형 셸은 20개, 원추형 셸은 10개, 두 번째 원통형 셸은 10개, 마지막 반구형 셸은 6개, 따라서 총 56개의 축대칭 원추요소로 모델링한다. 절점 1은 압력용기의 좌단에 해당되고, 절점 57은 우단에 해당된다. 대칭성을 고려하여 절점 57의 변위 와 각변위 θ57을 추가로 구속하여 계산하였다.
유한요소-전달강성계수법으로 주요 절점의 변위를 계산한 결과를 Table 7에 나타내었다. 유한요소법에 근거하여 작성된 Ross3)의 전산 프로그램으로 동일 모델에 대하여 계산한 결과를 Table 8에 나타내었다. 양 방법의 계산 결과를 비교해 보면 아주 잘 일치하고 있음을 알 수 있다.
Fig. 8은 계산 결과로부터 나타낸 압력용기의 변형된 형상이다. 변형 결과를 확연히 보여 주기 위해 각 절점의 변위를 500배 확대하였고, 색상과 격자는 3차원 형상을 명확히 표현하기 위해 제공된 것 외에는 특별한 의미가 없다.
5. 결 론
축대칭 셸의 정적해석은 공학분야에서 중요한 주제 중의 하나이다. 이 논문에서 저자들은 유한요소법의 모델링 기법과 전달강성계수법의 전달 기법을 조합한 유한요소-전달강성계수법을 이용하여 다양한 형상을 갖는 축대칭 셸의 정적해석을 수행하기 위한 전산 알고리즘을 제안한다.
이 연구에서는 유한요소-전달강성계수법과 유한요소법으로 다양한 형상을 갖는 축대칭 셸의 정적해석을 수행하기 위한 전산알고리즘을 정식화한 후, 유한요소-전달강성계수법의 계산 정확성을 확인하기 위하여 3개의 계산모델(균일한 압력을 받는 원형 돔, 정수압을 받는 물탱크, 내부 압력을 받는 압력용기)에 대해 정적해석을 수행하였다. 그리고 동일한 모델에 대하여 유한요소법으로 해석을 수행한 후, 양 방법의 계산 결과를 비교한 결과 이 연구에서 제안한 유한요소-전달강성계수법이 다양한 형상을 갖는 축대칭 셸의 정적해석에 신뢰성 있는 해석 기법임을 확인할 수 있었다.
― 기 호 설 명 ―
kx, ky : | x축 및 y축방향 스프링의 스프링상수 |
kθ : | 회전방향 스프링의 스프링상수 |
l : | 원추요소의 자오선방향 길이 |
px, py : | x축 및 y축방향의 압력 |
qx, qy : | x축 및 y축방향의 집중 절점하중 |
qθ : | 절점에 작용하는 모멘트 하중 |
r : | 원추요소 내부의 반경 |
R : | 절점의 반경 |
t : | 원추요소의 두께 |
u, w : | x축 및 y축방향의 변위 |
d : | 절점의 변위벡터(3×1) =[u, w, θ]T |
: | 절점 좌측과 우측의 힘보정벡터(3×1) |
: | 절점 좌측과 우측의 힘벡터(3×1) |
p, r : | 원추요소 좌측 및 우측 절점에 부여되는 압력의 등가절점 하중벡터(3×1) |
q : | 절점의 하중벡터(3×1) =[qx, qy, qθ]T |
A, B, C, D : | 요소강성행렬의 부분행렬(3×3) |
H : | 기초지지 강성행렬(3×3) |
K : | 원추요소의 강성행렬(6×6) |
R : | 좌표변환행렬(3×3) |
: | 절점 좌측, 우측의 강성계수행렬(3×3) |
그리스 문자
α : | 원추요소의 반정각 |
β : | 인접 원추요소 사이의 경사각 |
θ : | 각변위 |
ξ : | 요소의 정규화 좌표 (=x/l) |
상첨자
i, j : | 국부좌표계 번호 |
T : | 전치행렬 |
O : | 전체좌표계 |
-1 : | 역행렬 |
하첨자
i : | 절점 번호 |
(i) : | 요소 번호 |
Author contributions
M. S. Choi; Conceptualization, Investigation, Methodology, Software and Writing-original draft. K. U. Yang; Project administration, Validation, Visualization, Writing-review & editing and Supervision.
References
- H. Kraus, 1967, “Thin Elastic Shells”, John Wiley & Sons, New York.
-
J. N. Reddy, 2007, “Theory and Analysis of Elastic Plates and Shells (Second Edition)”, CRC Press, Boca Raton, 403-478.
[https://doi.org/10.1201/9780849384165]
- C. T. F. Ross, 1982, “Computational Methods in Structural and Continuum Mechanics”, John Wiley & Sons, New York, 79-89.
-
W. Jr. Weaver and P. R. Hohnston, 1984, “Finite elements for Structural Analysis”, Prentice-Hall, New Jersey.
[https://doi.org/10.1115/1.3167704]
-
H. S. Ryu and H. C. Sin, 1997, “Two Node Meridional Strain-based Axisymmetric Shell Elements”, Transactions of the KSME A, 21(6), 925-932.
[https://doi.org/10.22634/KSME-A.1997.21.6.925]
- Y. W. Kwon and H. C. Bang, 2000, “The Finite element method using MATLAB (2nd Edition)”, CRC Press, Boca Raton.
-
M. S. Choi and D. J. Yeo, 2012, “Stress Analysis of Axisymmetric Cylindrical Shell”, Journal of The Korean Society for Power System Engineering, 16(6), 45-51.
[https://doi.org/10.9726/kspse.2012.16.6.045]
-
M. S. Choi, J. H. Byun and D. J. Yeo, 2015, “Structural Analysis of Axisymmetric Conical Shells Using Finite Element-Transfer Stiffness Coefficient Method”, Journal of The Korean Society for Power System Engineering, 19(1), 38-44.
[https://doi.org/10.9726/kspse.2015.19.1.038]
Appendix
부 록
축대칭 셸의 해석에 사용되는 절두 원추요소의 강성행렬과 압력의 등가절점하중을 구하는 방법은 여러 가지가 있지만 이 연구에서 채택한 Ross3)가 사용한 방법은 다음에 기초한다.
i번째 절두 원추요소 내부의 임의의 위치 ξ = xi/li (0 ≤ ξ ≤ 1)에서 자오선방향 및 법선방향의 변위 ui, wi를 형상함수행렬과 절점변위벡터를 이용하여 다음과 같이 가정한다.
여기서
Novozhilov의 식으로 변형률벡터를 다음 식과 같이 유도할 수 있다.
여기서
탄성이론으로부터 합응력(stress resultant)을 표시하면 다음 식과 같다.
여기서 E와 ν는 탄성계수와 푸아송비이다.