등매개변수 요소

편집 역사 토론
구조해석 수치해석 마지막 수정: 2026-07-11 04:17:22

1. 개요[편집]

등매개변수 요소(Isoparametric element)는 유한요소법에서 요소의 기하학적 형상과 그 안의 물리량(변위 등)을 똑같은 형상함수(shape function)로 보간하는 요소를 말한다. “등(等)“은 같다는 뜻, “매개변수(parametric)“는 매개변수 좌표를 쓴다는 뜻이니, 이름 그대로 “형상과 해를 같은 함수로 매개한다”는 아이디어를 통째로 담고 있다.

이 개념 하나 덕분에 유한요소법은 직사각형·정삼각형 같은 반듯한 요소를 넘어, 곡선 경계와 찌그러진 형상까지 자유롭게 다룰 수 있게 되었다. 실무에서 쓰는 상용 유한요소법 코드의 요소는 사실상 전부 등매개변수 정식화(isoparametric formulation) 위에 서 있다. 1968년 이르온스(Bruce Irons)가 정립한 이 아이디어는1 오늘날 CAE 소프트웨어의 보이지 않는 뼈대다.

2. 핵심 아이디어: 하나의 함수, 두 가지 역할[편집]

찌그러진 사각형 요소를 생각해 보자. 실제 좌표계(물리 좌표 x,yx, y)에서 이 요소는 삐뚤빼뚤하지만, 매개변수 좌표계(자연 좌표 ξ,η\xi, \eta, 각각 −1-1에서 +1+1)에서는 반듯한 정사각형이다. 등매개변수 정식화의 핵심은 이 반듯한 기준 요소(reference element) 위에서만 계산하고, 실제 형상으로는 사상(mapping)을 통해 넘어가는 것이다.

기하 형상의 보간:

x=∑i=1nNi(ξ,η) xi,y=∑i=1nNi(ξ,η) yix = \sum_{i=1}^{n} N_i(\xi, \eta)\, x_i, \qquad y = \sum_{i=1}^{n} N_i(\xi, \eta)\, y_i

변위장의 보간:

u=∑i=1nNi(ξ,η) ui,v=∑i=1nNi(ξ,η) viu = \sum_{i=1}^{n} N_i(\xi, \eta)\, u_i, \qquad v = \sum_{i=1}^{n} N_i(\xi, \eta)\, v_i

두 식에서 형상함수 NiN_i가 완전히 같다는 점이 등매개변수의 정의 그 자체다. 절점 좌표 xi,yix_i, y_i를 넣으면 형상이 되고, 절점 변위 ui,viu_i, v_i를 넣으면 해가 된다. 만약 기하 보간의 차수가 해 보간보다 낮으면 준매개변수(subparametric), 높으면 초매개변수(superparametric) 요소라 부르지만, 실무의 절대다수는 등매개변수다.

3. 형상함수[편집]

형상함수는 기저함수의 일종으로, 각 절점에서 1이고 나머지 절점에서 0이 되는 성질(크로네커 델타 성질)을 갖는다. 4절점 사각형(Q4) 요소의 형상함수는 다음과 같이 깔끔하다.

Ni(ξ,η)=14(1+ξiξ)(1+ηiη)N_i(\xi, \eta) = \frac{1}{4}(1 + \xi_i \xi)(1 + \eta_i \eta)

여기서 (ξi,ηi)(\xi_i, \eta_i)는 ii번째 절점의 자연 좌표다. 8절점(Q8)이나 9절점(Q9) 같은 2차 요소로 가면 형상함수가 테일러 급수의 고차항을 담아 곡선 경계를 표현할 수 있게 된다. 2차 요소의 변이 곡선을 그릴 수 있다는 것이 등매개변수의 진짜 위력이다. 원통이나 날개 단면 같은 곡면을 직선 요소로 근사하려면 요소를 무진장 잘게 쪼개야 하지만, 2차 등매개변수 요소는 하나로 곡선 한 조각을 흉내 낸다.

형상함수가 만족해야 할 조건이 몇 가지 있다. 절점에서 값이 1/0이어야 하고(보간성), 어느 점에서든 모든 형상함수의 합이 1이어야 하며(분할의 통일성, partition of unity), 강체 운동과 일정 변형률 상태를 정확히 재현할 수 있어야 한다(완전성). 이 조건들이 깨지면 해가 수렴하지 않는다.

4. 자코비안: 두 좌표계를 잇는 다리[편집]

문제는 유한요소 정식화에서 필요한 것이 변위의 물리 좌표 미분(∂u/∂x\partial u / \partial x)인데, 형상함수는 자연 좌표 ξ,η\xi, \eta의 함수라는 점이다. 이 간극을 메우는 것이 자코비안 행렬이다.

J=[∂x∂ξ∂y∂ξ∂x∂η∂y∂η]\mathbf{J} = \begin{bmatrix} \dfrac{\partial x}{\partial \xi} & \dfrac{\partial y}{\partial \xi} \\[2mm] \dfrac{\partial x}{\partial \eta} & \dfrac{\partial y}{\partial \eta} \end{bmatrix}

연쇄 법칙으로 자연 좌표 미분을 물리 좌표 미분으로 바꾸려면 자코비안의 역행렬 J−1\mathbf{J}^{-1}이 필요하다. 그래서 자코비안 행렬식 det⁡J\det \mathbf{J}가 0이 되거나 음수가 되면 사상이 뒤집히면서 계산이 통째로 깨진다. 이것이 실무에서 그렇게 잔소리를 듣는 “요소 뒤틀림(distortion)” 경고의 정체다.2 요소가 심하게 찌그러지거나 오목해지면 det⁡J\det \mathbf{J}가 위험해지고, 해의 정확도가 급락한다. 격자 품질 지표 중 자코비안 비(Jacobian ratio)가 바로 이걸 감시하는 숫자다.

5. 가우스 적분점에서의 수치적분[편집]

등매개변수 요소의 강성행렬은 다음과 같은 적분으로 정의된다.

Ke=∫−11∫−11BTD B det⁡J  dξ dη\mathbf{K}^e = \int_{-1}^{1} \int_{-1}^{1} \mathbf{B}^{T} \mathbf{D}\, \mathbf{B}\, \det\mathbf{J}\; d\xi\, d\eta

여기서 B\mathbf{B}는 변형률-변위 행렬, D\mathbf{D}는 재료 물성 행렬이다. 이 적분을 손으로 정확히 풀기는 거의 불가능하다. B\mathbf{B}와 det⁡J\det\mathbf{J}가 자연 좌표의 복잡한 유리함수라서다. 그래서 수치적분, 특히 가우스-르장드르 구적법(Gauss-Legendre quadrature)을 쓴다.

Ke≈∑i∑jwiwj BT(ξi,ηj) D B(ξi,ηj) det⁡J(ξi,ηj)\mathbf{K}^e \approx \sum_{i} \sum_{j} w_i w_j\, \mathbf{B}^{T}(\xi_i, \eta_j)\, \mathbf{D}\, \mathbf{B}(\xi_i, \eta_j)\, \det\mathbf{J}(\xi_i, \eta_j)

기준 요소가 반듯한 정사각형이라 적분 구간이 항상 [−1,1][-1, 1]로 고정된다는 점이 여기서 빛을 발한다. 어떤 찌그러진 요소든 가우스 적분점의 위치와 가중치 ww가 미리 정해진 표로 통일된다.3 가우스 적분점은 함수를 실제로 평가하는 지점인데, 재미있게도 응력·변형률 같은 물리량이 가장 정확한 곳도 절점이 아니라 바로 이 적분점이다(Barlow 점). 그래서 후처리에서 응력을 절점으로 옮길 때는 적분점 값을 외삽한다.

적분점을 정식보다 하나 적게 쓰는 저감적분(reduced integration)은 계산을 빠르게 하고 과강성(shear locking)을 완화하지만, 대신 모래시계 모드(hourglass mode)라는 가짜 0-에너지 변형이 생길 수 있어 주의가 필요하다. 이 트레이드오프는 명시적 동해석 코드에서 특히 자주 등장한다.

6. 왜 이 개념이 중요한가[편집]

등매개변수 요소가 없었다면 유한요소법은 반듯한 격자에 갇힌 학술적 장난감에 머물렀을 것이다. 이 개념 덕분에 임의의 곡선 형상을 소수의 요소로 표현하고, 모든 요소를 하나의 기준 요소 계산 루틴으로 통일할 수 있게 되었다. 프로그래밍 관점에서도 축복인데, 삼각형이든 사각형이든 육면체든 형상함수와 적분점 표만 바꿔 끼우면 같은 조립 코드가 돌아간다.

정리하면 등매개변수 요소는 세 가지를 한꺼번에 해낸다. 형상과 해를 같은 함수로 다루는 우아함, 자코비안 사상을 통한 임의 형상 대응, 그리고 가우스 적분점 기반 수치적분의 표준화다. 이 삼박자가 구조해석뿐 아니라 전자기해석, 열해석 등 모든 유한요소 분야의 공통 언어가 되었다.

7. 관련 문서[편집]

8. Footnotes[편집]

  1. Irons, B. M. (1968)의 등매개변수 사각형 요소 논문이 출발점으로 꼽힌다. “요소를 반듯한 기준 공간에서 계산하고 실제로는 구부린다”는 발상은, 지금 보면 당연하지만 당시엔 유한요소법을 실무로 끌어올린 결정적 한 방이었다. ↩

  2. 격자 생성기가 “Jacobian check failed”를 뱉으면 십중팔구 어딘가에 오목하거나 뒤집힌 요소가 숨어 있다. 이걸 무시하고 돌리면 해석이 발산하거나, 더 무섭게는 그럴싸한 쓰레기 결과를 내놓는다. Garbage in, garbage out의 요소 버전. ↩

  3. 2×22 \times 2 가우스 적분이면 4개 점, 각 방향 좌표가 ±1/3\pm 1/\sqrt{3}이다. 이 네 자리 숫자를 외우고 다니는 유한요소 개발자가 실제로 있다. 그게 바로 등매개변수의 표준화가 준 선물이자 저주다. ↩