용어 목록 > 응집영역모델 (Cohesive Zone Model)

응집영역모델 (Cohesive Zone Model)

쉽게 풀면

종이를 천천히 찢을 때 완전히 두 조각으로 갈라지기 전에도 찢어지는 부분 근처는 서로 잡아당기며 힘을 어느 정도 버티는 구간이 있습니다. 응집영역모델은 이렇게 균열이 완전히 벌어지기 전, 균열 끝부분에서 아직 서로 붙잡고 있는 힘과 벌어짐 정도의 관계를 수학적으로 표현해서 파괴 과정을 모사하는 방법입니다.

왜 중요한가

복합재료의 층간분리나 접착 계면의 파괴처럼 균열이 점진적으로 진전하는 현상을 시뮬레이션하려면 이 모델이 자주 사용됩니다. 유한요소해석과 결합하여 실제 구조물에서 균열이 언제, 어떻게 진전할지 예측하는 데 널리 활용됩니다. 균열 선단의 특이 응력장을 직접 다루는 대신 견인력-변위 관계로 손상 과정을 표현하므로, 유한요소해석에서 초기 균열을 미리 넣어 두지 않고도 균열의 발생과 성장 거동을 함께 모사할 수 있습니다. 다만 일반적인 계면 응집요소는 균열이 진행할 경로에 미리 배치해야 하므로, 경로를 알 수 없는 문제에는 다른 기법이 필요합니다.

어떤 절차인가

1. 계면(또는 예상 균열 경로)에 응집요소를 둡니다. 이 요소는 두께가 거의 없고, 벌어짐(변위 차)과 그 사이에 작용하는 견인력(traction)의 관계로 거동을 정의합니다.
2. 견인력–변위 관계(응집법칙)를 고릅니다. 삼각형(이선형)이 가장 흔하고, 사다리꼴이나 다항식·지수형도 씁니다(Needleman, 1987). 처음에는 벌어짐에 비례해 견인력이 커지다가 최대 견인력에서 손상이 시작되고, 이후 견인력이 줄어 0이 되면 완전히 분리됩니다.
3. 곡선 아래 면적이 파괴에너지입니다. 따라서 최대 견인력과 파괴에너지, 두 값으로 삼각형 법칙이 정해집니다.
4. 모드 혼합(열림과 전단이 섞인 경우)이면 혼합 모드 기준을 더합니다(Camanho et al., 2003).
5. 유한요소해석에서 하중을 늘리며 손상의 시작과 진전을 함께 계산합니다.

이 모델의 뿌리에는 균열 선단 앞쪽에 응력이 작용하는 영역을 두는 생각이 있습니다(Dugdale, 1960; Barenblatt, 1962). 콘크리트 균열에 적용한 가상균열모델도 같은 계열입니다(Hillerborg et al., 1976).

언제 쓰고 무엇을 얻나

균열 선단의 무한대 응력(특이점)을 직접 다루지 않고 균열의 발생과 진전을 하나의 모델로 표현할 수 있고, 접착 조인트의 박리나 복합재료의 층간분리처럼 균열이 퍼지는 경로가 계면으로 정해져 있는 문제에 잘 맞습니다. 이와 달리 전통적인 파괴역학은 이미 있는 균열이 자랄지 판정하는 데 J-적분 같은 양과 파괴인성을 씁니다.

한계도 분명합니다. 일반적인 계면 응집요소는 균열이 진행할 경로를 미리 알고 그 경로에 배치해야 합니다. 또 요소망이 응집영역 길이에 비해 너무 거칠면 결과가 부정확해집니다(Turon et al., 2007).

아래는 가상의 접착층 값으로 삼각형 법칙을 계산한 예입니다.

• 최대 견인력 σ_max = 50 MPa, 파괴에너지 G_c = 0.5 N/mm, 초기 강성 K = 100,000 N/mm³
• 손상 시작 변위 δ₀ = σ_max / K = 0.0005 mm
• 완전 분리 변위 δ_f = 2G_c / σ_max = 0.02 mm (삼각형 면적 ½·σ_max·δ_f = G_c)
• 특성 길이 l = E·G_c / σ_max² = 0.6 mm (E = 3000 MPa 가정). 응집영역 길이는 이 길이에 모델에 따른 상수(대략 0.2~1)를 곱한 정도이고, 그 길이 안에 요소가 여러 개 들어가야 합니다.

논문에는 이렇게 적는다

• 응집법칙의 형태(삼각형, 사다리꼴 등)와 최대 견인력, 파괴에너지, 초기 강성을 모두 적습니다. 이 값들에 결과가 민감하므로 얻은 방법(실험 보정 여부)도 밝혀야 합니다.
• 응집요소를 어디에 배치했는지(미리 정한 경로인지), 요소 종류와 응집영역 안의 요소 수, 요소망 민감도 확인을 적습니다.
• 하중–변위 곡선, 파괴 하중, 균열 길이 등을 실험과 비교해 모델을 보정·검증한 결과를 보고합니다.
• 소프트웨어에서 쓰는 손상 안정화(점성 정규화) 설정이 있으면 그 값도 적습니다.

논문에서는 이렇게 쓰입니다

실제 논문에서

이 용어를 직접 쓴 국내 논문 문장을 아직 수집하지 못했습니다. 수집되는 대로 출처와 함께 싣겠습니다.

예시 문장

아래는 용어의 쓰임을 보여 주려고 만든 예시로, 실제 논문에서 가져온 문장이 아닙니다.

컴퓨터 시뮬레이션에서 응집영역모델을 사용해 층간분리가 퍼져나가는 과정을 재현했다는 뜻입니다.

모델의 세부 값을 실험 결과에 맞춰 조정해서 접착 부위가 언제 파괴되는지 예측했다는 의미입니다.

조금 더 깊게 보면

응집영역모델은 흔히 응력-변위 관계를 삼각형이나 사다리꼴 형태의 곡선으로 정의하며, 이 곡선 아래 면적이 해당 계면의 파괴에너지에 해당합니다. 유한요소해석에서는 균열이 진전할 것으로 예상되는 경로에 응집요소를 배치하여 손상의 시작과 진전을 함께 계산합니다.

균열 선단 앞쪽에 응집영역을 설정해 견인력과 균열면 변위 사이의 관계를 정의함으로써, 특이점 없이도 균열의 발생과 전파를 하나의 손상 모델 안에서 함께 다룰 수 있게 해줍니다. 삼각형 형태의 견인력-변위 관계가 가장 널리 쓰이지만, 재료의 연화 거동에 따라 사다리꼴이나 지수함수 형태의 관계식을 대신 사용하기도 합니다.

주의할 점

응집영역모델의 결과는 모델에 입력하는 응력-변위 곡선의 형태와 파라미터 값에 민감하게 반응하므로, 이를 실험으로 신중하게 보정하지 않으면 예측이 실제와 크게 어긋날 수 있습니다. 모델의 예측 결과는 견인력-변위 곡선의 형태와 최대 견인력, 파괴에너지 등 입력 파라미터 설정에 민감하므로, 실험으로 이 값들을 신중하게 보정해야 신뢰할 수 있는 해석 결과를 얻을 수 있습니다.

헷갈리기 쉬운 용어

파괴를 다루는 틀은 균열의 시작을 어떻게 보느냐가 다릅니다.

응집영역모델의 파괴에너지는 파괴역학의 임계 에너지해방률과 대응시켜 입력합니다.

참고 문헌

1. Dugdale, D. S. (1960). Yielding of steel sheets containing slits. Journal of the Mechanics and Physics of Solids, 8(2), 100–104.
2. Barenblatt, G. I. (1962). The mathematical theory of equilibrium cracks in brittle fracture. Advances in Applied Mechanics, 7, 55–129.
3. Hillerborg, A., Modéer, M., & Petersson, P.-E. (1976). Analysis of crack formation and crack growth in concrete by means of fracture mechanics and finite elements. Cement and Concrete Research, 6(6), 773–781.
4. Needleman, A. (1987). A continuum model for void nucleation by inclusion debonding. Journal of Applied Mechanics, 54(3), 525–531.
5. Camanho, P. P., Dávila, C. G., & de Moura, M. F. (2003). Numerical simulation of mixed-mode progressive delamination in composite materials. Journal of Composite Materials, 37(16), 1415–1438.
6. Turon, A., Dávila, C. G., Camanho, P. P., & Costa, J. (2007). An engineering solution for mesh size effects in the simulation of delamination using cohesive zone models. Engineering Fracture Mechanics, 74(10), 1665–1682.

표의 기준값은 위 문헌에 근거한 관례이며, 분야·학술지마다 다른 기준을 쓸 수 있습니다.

관련 용어