판과 셸 이론

구조해석 고체역학 마지막 수정: 2026-07-13 04:21:52

1. 개요[편집]

판과 셸 이론(plate and shell theory)은 두께가 다른 두 치수에 비해 얇은 판상(板狀) 구조를, 두께 방향을 적분으로 날려버리고 중립면(mid-surface)이라는 2차원 곡면 하나의 거동으로 축약해 해석하는 고체역학 이론이다. 보 이론이 3차원 물체를 1차원 선으로 줄였다면, 판·셸 이론은 3차원을 2차원 면으로 줄이는 한 차원 위의 형제뻘이다. 평평하면 판(plate), 곡률이 있으면 셸(shell)이라 부른다. 자동차 차체 패널, 항공기 동체와 날개 외피, 압력용기, 선박 외판, 심지어 콜라 캔까지 — 얇은 껍데기로 하중을 버티는 것들이 전부 이 이론의 관할이다.¹

핵심은 얇음이 주는 축복이다. 두께 방향으로는 응력이 거의 선형이라 가정하고 적분해 치우면, 미지수가 두께 좌표에서 해방되어 계산량이 극적으로 줄어든다. 물론 그 대가로 얇음의 정도에 따라 서로 다른 근사 모델을 골라 써야 하는 번거로움이 따라온다.

2. 키르히호프-러브 판[편집]

얇은 판의 고전 모델은 키르히호프-러브 판(Kirchhoff–Love plate) 이론이다. 보 이론의 오일러-베르누이 가정을 그대로 2차원으로 확장한 것으로, 두께 방향 법선이 변형 후에도 곧고 중립면에 수직을 유지한다고 본다. 즉 횡전단 변형을 무시한다. 이 가정 아래 판의 처짐 $w(x, y)$ 는 다음 4계 편미분방정식, 이른바 판 방정식(또는 소피 제르맹-라그랑주 방정식)을 만족한다.

D \nabla^4 w = q

여기서 $\nabla^4 = \nabla^2\nabla^2$ 은 중조화(biharmonic) 연산자, $q$ 는 단위 면적당 하중이다. $D$ 는 판의 굽힘 강성(flexural rigidity)으로, 다음처럼 정의된다.

D = \frac{E h^3}{12 (1 - \nu^2)}

$h$ 는 판 두께, $\nu$ 는 포아송비다. 보의 $EI$ 에 대응하는 자리인데, 두께의 세제곱 의존성은 여기서도 그대로다. 판이 얇아질수록 강성이 세제곱으로 곤두박질치니, 얇은 철판이 손으로도 휘청거리는 이유가 여기 있다. 분모의 $(1-\nu^2)$ 은 판이 한 방향으로 휘면 포아송 효과로 옆 방향 변형이 구속되어 실제보다 뻣뻣해지는 2차원 특유의 보정이다.²

3. 민들린-라이스너 판[편집]

판이 두꺼워지거나(두께/변 길이 비가 대략 1/10을 넘거나), 복합재 해석처럼 층간 전단이 중요해지면 키르히호프 가정이 처짐을 과소평가한다. 이때 등장하는 것이 민들린-라이스너 판(Mindlin–Reissner plate) 이론이다. 보 이론의 티모셴코 보에 정확히 대응하는 모델로, 법선이 변형 후 중립면에 수직일 필요가 없게 풀어주어 횡전단 변형을 허용한다.

그 결과 처짐 $w$ 와 두 개의 단면 회전각 $\theta_x, \theta_y$ 가 서로 독립된 미지수가 된다. 미지수가 늘어 계산은 무거워지지만, 두꺼운 판과 모드 해석에서의 고차 진동에 훨씬 정확하다. 실무 유한요소법 코드가 판·셸에 대해 민들린 계열을 기본으로 채택하는 이유는, 회전각을 독립 절점 자유도로 두는 편이 요소 정식화가 깔끔하고 저차 보간으로도 처리되기 때문이다.

4. 셸: 막 거동과 굽힘 거동[편집]

셸(shell)은 곡률을 가진 판이다. 곡면이라는 사실 하나가 거동을 근본적으로 바꾼다. 평평한 판은 횡하중을 오직 굽힘으로만 버티지만, 곡면 셸은 하중을 면 내부의 인장·압축, 즉 막(membrane) 응력으로도 나눠 진다. 달걀 껍데기가 그 얇기로도 손아귀 압력을 버티는 것, 돔 지붕이 기둥 없이 넓은 공간을 덮는 것이 전부 막 거동이 하중의 대부분을 인장·압축으로 흘려보내기 때문이다.³

그래서 셸의 총 거동은 대략 막 거동 + 굽힘 거동의 중첩으로 이해된다. 곡률이 클수록 막 거동이 지배적이라 대단히 효율적이지만, 국소적으로 집중하중이나 경계 조건이 걸리면 그 부근에서 굽힘이 급격히 되살아나는 경계층(boundary layer)이 생긴다. 얇은 셸일수록 이 경계층이 얇고 응력이 뾰족해서, 격자를 촘촘히 깔지 않으면 응력집중을 놓친다. 셸 좌굴은 초기 형상 결함에 극도로 민감한 것으로 악명 높아, 좌굴 해석에서 판·셸은 늘 골칫거리 단골이다.

5. 전단 잠김과 MITC[편집]

민들린 판·셸 요소를 순진하게 짜서 얇은 판에 적용하면 처참한 일이 벌어진다. 두께가 0으로 가는 극한에서 전단 변형률이 0이 되어야 하는데, 저차 요소는 이 조건을 정확히 만족시키지 못해 존재하지도 않는 전단 에너지를 잔뜩 만들어낸다. 그 결과 요소가 실제보다 수십·수백 배 뻣뻣해져 처짐이 거의 안 나오는 현상이 전단 잠김(shear locking)이다.⁴

이 잠김을 잡으려 여러 처방이 개발되었다. 전단 항만 낮은 차수로 적분하는 축소 적분(reduced integration)은 간단하지만 가짜 변형 모드(hourglass mode)를 부르는 부작용이 있다. 오늘날 상용 코드의 주력은 전단 변형률을 요소 내부 특정 점(tying point)에서 따로 보간해 잠김을 원천 차단하는 MITC(Mixed Interpolation of Tensorial Components) 계열 요소다. MITC4(4절점), MITC9 같은 이름이 Abaqus, NASTRAN 등의 셸 요소 목록에서 흔히 보이는 이유가 이것이다. 잠김·수렴성·강건성 사이의 줄타기는 판·셸 유한요소법 연구가 지금도 논문을 쏟아내는 분야다.

6. 관련 문서[편집]

“판이냐 셸이냐”의 경계는 곡률뿐이다. 곡률이 0이면 판, 0이 아니면 셸. 그래서 유한요소 코드에서는 아예 둘을 합쳐 “셸 요소”로 부르고 평판은 곡률 0인 특수 경우로 취급하는 경우가 많다. 판을 별도 이름으로 부르는 건 다분히 역사적 관성이다. ↩
소피 제르맹(Sophie Germain)은 여성이 과학원에 발도 못 들이던 19세기 초에 판의 진동 방정식으로 파리 과학원 상을 받았다. 다만 초기 유도에 오류가 있어 라그랑주가 손을 봤고, 그래서 방정식에 두 사람 이름이 함께 붙는다. ↩
이 원리를 극한까지 밀어붙인 게 얇은 셸 구조물이다. 두께 대비 반경 비가 1000을 넘는 냉각탑이나 얇은 콘크리트 셸 지붕이 서 있을 수 있는 것은, 하중이 거의 전부 막 응력으로 흐르도록 형상을 잡았기 때문이다. ↩
잠김은 판·셸만의 문제가 아니다. 비압축성 재료에서의 체적 잠김(volumetric locking) 등 사촌이 여럿 있다. 공통점은 “요소가 특정 물리적 극한을 표현할 자유도가 부족해 가짜 강성을 만들어낸다”는 것. 유한요소 연구의 절반은 이 잠김들과의 전쟁이라 해도 과언이 아니다. ↩