[논문 리뷰] A threshold-force model for adhesion and mode I fracture

Srivatsan Hulikal, Kaushik Bhattacharya외 1명|arXiv (Cornell University)|2016. 06. 10.

Adhesion, Friction, and Surface Interactions참고 문헌 6인용 수 3

한 줄 요약

이 논문은 인장 임계력 힘을 갖는 이산 스프링 모델을 제안하여 점착 접촉과 모드 I 균열을 시뮬레이션하며, 선형 탄성 균열역학과 일치하는 바에 따라 임계력 힘을 $1/\sqrt{r}$로 스케일링할 경우 매크로스코픽으로 치수에 영향을 받지 않는 거동를 나타낸다. 이 모델은 주요 JKR 점착 스케일링 법칙—최대 인장력이 균열 에너지와 구의 반지름에 선형적으로 의존하는 것—을 재현하여 임계력 힘을 균열 에너지, 엘라스토모듈러스, 푸아송 비율과 같은 물성으로 직접 캘리브레이션할 수 있게 한다.

ABSTRACT

We study the relation between a threshold-force based model at the microscopic scale and mode I fracture at the macroscopic scale in a system of discrete interacting springs. Specifically, we idealize the contact between two surfaces as that between a rigid surface and a collection of springs with long-range interaction and a constant tensile threshold force. We show that a particular scaling similar to that of crack-tip stress in Linear Elastic Fracture Mechanics leads to a macroscopic limit behavior. The model also reproduces the scaling behaviors of the JKR model of adhesive contact. We determine how the threshold force depends on the fracture energy and elastic properties of the material. The model can be used to study rough-surface adhesion.

연구 동기 및 목표

• 이산 스프링에서의 미시적 임계력 힘 거동과 매크로스코픽 모드 I 균열역학 간의 연결 고리를 설정하기 위해.
• 해상도 크기와의 비례 관계에서 매크로스코픽으로 치수에 영향을 받지 않는 반응을 얻기 위해 임계력 힘의 스케일링을 규명하기 위해.
• 특히 최대 인장력이 균열 에너지와 구의 반지름에 선형적으로 의존하는 바를 재현함으로써 모델을 검증하기 위해.
• 균열 에너지, 전단 모odulus, 푸아송 비율과 같은 물성으로 표현되는 임계력 힘의 해석적 표현을 유도하기 위해.
• 거칠은 표면에 특히 적합한 장거리 탄성 상호작용을 수반하는 점착 접촉의 수치 시뮬레이션을 가능하게 하기 위해.

제안 방법

• 반탄성 표면을 장거리 탄성 상호작용을 갖는 정렬된 이산 스프링 집합으로 모델링하며, 반무한체 탄성의 부시앙스크 및 루브 솔루션에 따라 제어된다.
• 각 스프링에 인장 임계력 힘 $F_{\text{th}}$ 를 도입하며, 이에 도달하면 스프링이 분리되고 힘이 0으로 떨어진다.
• 길이 척도 $L^*$ 와 힘 척도 $F^*$ 를 이용한 비차원화를 통해 차원 없는 탄성 계수 커널 $\bar{C}_{ij}$ 를 유도하며, $i \neq j$ 에서는 $1/\bar{r}_{ij}$, $i = j$ 에서는 $3.8/\bar{\Delta}$ 로 표현된다.
• 대규모 시스템에서 장거리 상호작용을 효율적으로 계산하기 위해 빠른 다체법(Fast Multipole Method)을 적용한다.
• 이동 제어 하중 조건 하에서 스프링 변형 및 힘을 업데이트하기 위해 예측-보정 알고리즘을 사용하며, 각 시간 단계에서 동적 접촉 상태 업데이트를 수행한다.
• 해상도 간격 $\Delta$ 에 대해 $F_{\text{th}} \propto \Delta^{\alpha}$ 로 스케일링하여 임계력 힘을 캘리브레이션하며, 매크로스코픽 불변성을 위한 임계값으로 $\alpha = 0.75$ 를 규명한다.

실험 결과

연구 질문

• RQ1 해상도 크기와의 비례 관계에서 어떤 임계력 힘 스케일링이 이산 스프링 모델에서 매크로스코픽으로 치수에 영향을 받지 않는 반응을 유도하는가?
• RQ2 임계력 힘 모델은 JKR 점착 모델의 핵심 스케일링 법칙을 어떻게 재현하는가?
• RQ3 모델의 임계력 힘은 균열 에너지 및 탄성 계수와 같은 기본 물성과 관련이 있는가?
• RQ4 모델의 거동은 특히 균열 끝에서의 $1/\sqrt{r}$ 응력 특이성에 기반한 선형 탄성 균열역학 예측과 일치하는가?
• RQ5 상호작용 커널을 수정하거나 자유도를 추가함으로써 점탄성 거동 또는 혼합 모드 균열을 시뮬레이션하기 위해 모델을 확장할 수 있는가?

주요 결과

• 매크로스코픽 반응이 해상도 크기와 무관해지기 위해서는 임계력 힘이 $\Delta^{0.75}$ 로 스케일링되어야 하며, 이는 선형 탄성 균열역학에서의 $1/\sqrt{r}$ 특이성과 정확히 일치한다.
• 모델은 JKR 스케일링을 재현한다: 최대 인장력은 균열 에너지와 구의 반지름에 선형적으로 증가하며, 점착 접촉 이론과의 일관성을 확인한다.
• 최대 인장력은 효과적 응력 집중 계수 $\bar{K}_{\text{th}}$ 의 제곱에 비례하며, $\bar{F}_{\text{max}} = 7.6 \bar{K}_{\text{th}}^2 \bar{R}$ 로 표현되며, 이는 균열역학 유사성의 검증을 받는다.
• 임계력 힘 매개변수 $K_{\text{th}}$ 는 $K_{\text{th}} = \pi \sqrt{\frac{3G\gamma}{7.6(1-\nu)}}$ 로 물성과 해석적으로 관련되어 있어 측정 가능한 물성 매개변수로부터 직접 캘리브레이션할 수 있다.
• 모델은 다양한 상호작용 커널에 대해 강건하다: 장거리 탄성 상호작용이 제거되거나 국소 상호작용으로 대체되면 치수 불변성이 붕괴되며, 이는 탄성 상호작용 범위의 역할을 확인한다.
• 커널 $C_{ij}(t)$ 를 시간에 따라 변화시키는 방식으로 시간 의존 거동을 확장할 수 있으며, 이는 점탄성 에너지 손실 및 속도 의존 효과의 시뮬레이션을 가능하게 한다.