[논문 리뷰] Rupture processes in fiber bundle models

P. C. Hemmer, Alex Hansen외 1명|arXiv (Cornell University)|2006. 02. 16.

Optical Network Technologies참고 문헌 5인용 수 4

한 줄 요약

이 논문은 다양한 응력 재분배 메커니즘을 가진 섬유 번들 모델에서의 파손 과정을 조사하며, 붕괴 폭발 크기 통계에 초점을 맞춘다. 분석적 및 수치적 방법을 사용하여, 등가하중분배 모델은 열역학적 극한에서 지수 ξ = 5/2인 거듭제곱 법칙 폭발 분포를 보이며, 국소하중분배 모델은 교차 행동을 보이고 새로운 임계성 클래스를 나타내며, 임계점에서 지수 −1.6인 거듭제곱 법칙 군집 크기 분포를 보인다.

ABSTRACT

Fiber bundles with statistically distributed thresholds for breakdown of individual fibers are interesting models of the static and dynamics of failures in materials under stress. They can be analyzed to an extent that is not possible for more complex materials. During the rupture process in a fiber bundle avalanches, in which several fibers fail simultaneously, occur. We study by analytic and numerical methods the statistics of such avalanches, and the breakdown process for several models of fiber bundles. The models differ primarily in the way the extra stress caused by a fiber failure is redistributed among the surviving fibers.

연구 동기 및 목표

• 다양한 응력 재분배 규칙 하에서 섬유 번들 모델의 고장 폭발 통계적 성질을 이해하기 위해.
• 다양한 하중분배 메커니즘에서 폭발 크기 분포가 보편적인 거듭제곱 법칙을 따르는지 확인하기 위해.
• 완전한 붕괴에 가까운 임계 행동을 조사하며, 특히 자가조직 임계성의 발생을 다루기 위해.
• 폭발 통계 및 고장 역학 측면에서 등가하중분배(GLOBAL) 모델과 국소하중분배 모델을 비교하기 위해.
• 탄성 결합과 공간 상관관계가 치명적 고장을 유도하는 데 미치는 영향을 탐구하기 위해.

제안 방법

• 증가하는 하중 하에서 섬유 번들 파손을 시뮬레이션하기 위해 등가하중분배(GLOBAL) 및 국소하중분배 모델을 사용한다.
• 각 하중 단계에서 번들의 평균 강도를 계산하기 위해 반복 고정점 알고리즘을 적용한다.
• 행렬의 역행렬 (I + K G)⁻¹의 테일러 전개를 기반으로 한 조건부 반복 해법을 사용하여 효율적인 계산을 수행한다.
• 그린 함수 G를 포함하는 행렬 곱셈을 가속하기 위해 빠른 푸리에 변환(FFTs)을 사용한다.
• 국소화가 시작될 때 고장 났던 결합의 군집을 식별하고 분석하기 위해 호센-코플라만 알고리즘을 구현한다.
• 고장을 다른 불안정성에 따라 제어하기 위해 적용력 F 또는 변위 D를 사용한다.

실험 결과

연구 질문

• RQ1 등가하중분배 섬유 번들에서 폭발 크기 분포가 열역학적 극한에서 거듭제곱 법칙을 따르는가?
• RQ2 시스템이 임계점 근처이거나 국소하중재분배일 경우 폭발 지수 ξ는 어떻게 변화하는가?
• RQ3 국소하중분배 모델에서 치명적 고장의 시작 시 고장 군집 크기 분포의 성격은 무엇인가?
• RQ4 국소하중분배 시스템의 고장 과정은 자가조직화되어 있으며, 새로운 임계성 클래스에 속하는가?
• RQ5 탄성 상수와 공간 상관관계는 치명적 고장으로의 전이에 어떻게 영향을 미치는가?

주요 결과

• 등가하중분배 모델에서는 열역학적 극한에서 폭발 크기 분포가 지수 ξ = 5/2인 거듭제곱 법칙을 따른다.
• 시스템이 임계점 근처일 경우 폭발 지수는 ξ = 3/2로 교차하여 스케일링 행동의 변화를 나타낸다.
• 유한한 탄성 상수를 가진 국소하중분배 모델에서는 여전히 지수 ξ ≈ 2.6인 거듭제곱 법칙 폭발 분포를 보이며, 이는 글로벌 모델의 ξ = 2.5와 일치한다.
• 변위 제어 하에서 국소화가 시작될 때 군집 크기 분포는 지수 −1.6인 거듭제곱 법칙을 따르며, 이는 퍼콜레이션(−2.05)과는 다름없는 새로운 임계성 클래스를 나타낸다.
• 탄성 상수 e가 작을 경우 공간 상관관계와 국소화로 인해 치명적 고장이 발생하며, 이 불안정성은 큰 e에서 힘 제어 하에서는 존재하지 않는다.
• 행렬 (I + K G)⁻¹의 테일러 전개를 기반으로 한 조건부 반복 해법은 안정적이고 효율적인 시뮬레이션을 가능하게 하며, 계산 복잡도는 nL²ln(L)로 표현되며 여기서 n은 전개의 항 수이다.