감소 적분과 모래시계 제어

구조해석 유한요소법 마지막 수정: 2026-07-20 04:12:47

1. 개요[편집]

감소 적분(reduced integration)은 유한요소법에서 요소 강성행렬을 만들 때, 완전적분(full integration)에 필요한 수치적분 점 수보다 하나 적은 적분점을 쓰는 기법이다. 얼핏 “정확도를 깎아 먹는 게으른 짓” 같지만, 실은 요소 잠김을 완화하고 계산 비용을 줄이는 강력한 처방이다. 대신 공짜는 아니어서, 모래시계 모드(hourglass mode)라는 제로 에너지 모드(zero-energy mode)가 딸려 나오고, 이걸 억제하는 것이 모래시계 제어(hourglass control)다.¹

이 조합은 오늘날 명시적 동해석 코드의 사실상 표준이다. LS-DYNA의 기본 솔리드 요소, Abaqus/Explicit의 C3D8R, 그리고 대부분의 충돌·성형·낙하 해석이 8절점 육면체 1점 감소적분 요소 위에서 돌아간다. 이유는 명확하다. 완전적분 대비 강성 계산이 8배 빠르니까.

2. 완전적분 대비 장점[편집]

8절점 육면체(HEX8) 요소를 완전적분하려면 $2\times2\times2 = 8$ 개의 가우스 적분점이 필요하다. 감소적분은 이걸 요소 중심 단 1점으로 줄인다. 얻는 이득은 두 가지다.

비용 절감: 적분점마다 변형률·응력을 계산하고 재료 모델을 호출해야 하는데, 명시적 해석에서는 이 물성 갱신이 전체 연산의 큰 비중을 차지한다. 적분점이 1/8이면 강성/내력 계산도 대략 1/8로 준다. 수백만 요소를 수십만 스텝 돌리는 충돌 해석에서 이 차이는 며칠과 몇 시간을 가른다.
잠김 완화: 완전적분은 요소 내부 모든 적분점에서 굽힘·비압축 제약을 촘촘히 평가해 전단 잠김·체적 잠김을 유발한다. 감소적분은 그 제약이 평가되던 지점을 제거하므로 요소가 굽힘·비압축 변형을 훨씬 자연스럽게 표현한다. 얇은 셸이나 준비압축성 재료에서 변위가 정답에 가깝게 나온다.

3. 부작용: 제로 에너지 모드[편집]

문제는 적분점 하나로는 요소가 표현할 수 있는 모든 변형 모드를 “감지”하지 못한다는 것이다. 1점 적분 HEX8 요소에는 요소 중심에서 변형률이 정확히 0으로 계산되는데도 실제로는 절점이 움직이는 변형 패턴이 존재한다. 즉 변형 에너지 $U = \frac{1}{2}\mathbf{d}^\top \mathbf{K} \mathbf{d} = 0$ 인데도 $\mathbf{d} \neq \mathbf{0}$ 인 모드다. 강성 행렬 관점에서는 강체 운동이 아닌데도 고유값이 0인 고유벡터에 해당한다.

이 모드가 만드는 변형 형상이 마름모가 번갈아 찌그러지는 모래시계(모래통) 모양이라 hourglass mode라 부른다. 아무 저항 없이 이 모드로 마음껏 흔들릴 수 있으므로, 제어하지 않으면 해가 격자 전체에 걸쳐 지그재그로 폭주하며 물리적 의미를 완전히 잃는다. 해석 결과 컨투어에 체스판·바둑판 무늬가 보인다면 십중팔구 모래시계 모드가 새고 있는 것이다.²

제로 에너지 모드의 개수는 요소마다 정해져 있다. 1점 적분 HEX8은 강체 운동 6개를 제외하고도 12개의 모래시계 모드를 가지며, 4절점 사각형 평면 요소(QUAD4)는 2개를 가진다. 완전적분이었다면 이 모드들이 모두 적분점에서 감지되어 0이 아닌 강성을 가졌을 텐데, 적분점을 줄이는 대가로 이들이 강성 없는 “구멍”으로 남은 것이다. 결국 감소적분은 잠김을 유발하던 과잉 강성을 덜어내면서, 그 과정에서 필요한 강성 일부까지 함께 날려버린 셈이다.

4. 모래시계 제어[편집]

모래시계 모드를 잡는 방법은 그 제로 에너지 모드에만 선택적으로 인공 저항을 부여하는 것이다. 정상적인 물리 모드(굽힘·인장 등)에는 손을 대지 않고, 에너지가 0으로 계산되는 방향으로만 복원력을 더해 준다. 크게 두 계열이 있다.

강성 기반(stiffness-type): 모래시계 모드에 인공 강성을 부여한다. 준정적·저속 문제에서 정확하지만, 계수가 크면 오히려 인위적으로 뻣뻣해진다.
점성 기반(viscosity-type): 모래시계 변형 속도에 비례하는 감쇠력을 준다. 고속 충돌처럼 변형률 속도가 높은 문제에 적합하다.

두 경우 모두 사용자가 모래시계 제어 계수(hourglass coefficient)를 지정하는데, 보통 0.03~0.1 범위의 작은 값을 쓴다.³ 너무 작으면 모래시계가 새고, 너무 크면 요소가 잠김처럼 뻣뻣해져 감소적분의 이점을 스스로 갉아먹는다. 이 계수 튜닝은 감쇠 튜닝과 함께 명시적 해석 엔지니어의 대표적인 손맛 영역이다.

정교한 모래시계 제어는 단순히 인공력을 더하는 데 그치지 않고, 그 힘이 강체 회전이나 정상적인 변형 모드와 직교하도록 설계한다. Flanagan-Belytschko가 제안한 방식은 모래시계 모드 형상 벡터에 절점 속도를 투영해 순수 모래시계 성분만 추출한 뒤 그 성분에만 저항을 주므로, 물리 응답을 오염시키지 않는다. 이런 정식화 덕분에 오늘날 감소적분 요소는 정확도를 크게 희생하지 않고도 완전적분의 8배 속도를 누린다.

4.1. 완전적분의 재부상[편집]

그렇다고 완전적분이 죽은 건 아니다. 준정적 성형이나 정밀 응력 평가가 필요한 곳에서는, 모래시계 자체가 아예 없는 완전적분 요소를 선호하기도 한다. 다만 순수 완전적분은 잠김에 취약하므로, 실무 요소는 대개 완전적분에 비적합 모드나 향상 변형률(enhanced assumed strain)을 얹은 개선판이다. Abaqus의 C3D8I(비적합 모드 요소)가 대표적인 절충으로, 모래시계 걱정 없이 굽힘도 잘 잡는다. 결국 요소 선택은 “속도와 잠김 회피(감소적분) 대 견고함과 정밀도(개선된 완전적분)” 사이의 판단이며, 정답은 문제 성격에 달렸다.

5. 실무 진단: 모래시계 에너지[편집]

건전성을 판정하는 표준 지표는 모래시계 에너지 비율이다. 해석이 끝난 뒤 전체 내부 에너지(internal energy) 대비 모래시계 에너지(hourglass energy)의 비를 확인한다.

\frac{E_\text{hourglass}}{E_\text{internal}} < 0.1

업계 관례상 이 비율이 대략 5~10%를 넘으면 결과를 신뢰하지 않고 재해석한다.⁴ 비율이 높을 때의 처방은 (1) 격자를 세분화해 모래시계 모드의 영향을 줄이거나, (2) 완전적분 요소나 개선된 정식화 요소로 바꾸거나, (3) 하중을 한 요소에 점하중으로 때리지 않고 분산시키는 것이다. 점하중이 감소적분 요소 하나에 직접 걸리면 모래시계가 거의 반드시 터진다.

한 가지 덧붙이면, 감소적분과 모래시계 제어의 균형은 요소 잠김과 정확히 반대 방향의 줄다리기다. 완전적분은 잠기고, 순수 감소적분은 흐물댄다. 실무 요소는 그 사이에서 “잠기지도 새지도 않는” 최적점을 노린다. 감소적분+모래시계 제어는 그 절충의 가장 성공적이고 널리 쓰이는 형태다.

6. 관련 문서[편집]

“모래시계 제어”는 영어 hourglass control의 직역인데, 현장에서는 그냥 “아워글래스 잡는다”고 한다. 신입이 “결과에 바둑판 무늬가 떠요” 하면 십중팔구 아워글래스다. ↩
모래시계 모드는 요소 하나에서는 조용하다가 격자 전체에 걸쳐 인접 요소끼리 위상이 맞물리면 거대한 지그재그 변형으로 증폭된다. 그래서 “요소 하나만 보면 멀쩡한데 전체 그림이 미쳐 있는” 특유의 증상이 나온다. ↩
코드마다 기본값과 알고리즘이 다르다. LS-DYNA는 *HOURGLASS 카드로 여러 타입(표준 점성, Flanagan-Belytschko, 강성 등)을 고를 수 있고, 성형 해석에서는 타입 6(강성 기반) 조합이 국룰처럼 쓰인다. ↩
“모래시계 에너지 10% 룰”은 논문에 못 박힌 절대 기준이라기보다 오랜 관행이다. 정밀한 응력값이 필요한 국소 부위라면 훨씬 더 엄격하게, 전체 거동만 보는 개념 설계 단계라면 다소 느슨하게 본다. ↩