[논문 리뷰] The nature of quasistatic deformation in granular materials

Jean-Noël Roux|ArXiv.org|2009. 01. 15.

Granular flow and fluidized beds참고 문헌 4인용 수 7

한 줄 요약

이 논문은 구형 입자 집합의 이산 요소 시뮬레이션을 통해 준정적 입자 변형에서 두 가지 별개의 유변학적 영역—제1형(접촉 스프링 강성에 의한 탄성 변형)과 제2형(접촉 네트워크의 비탄성 재배열)—을 규명한다. 주요 기여는 제1형이 작은 변형 및 하중 반전 조건에서 지배적이며, 이 경우 변형은 접촉 강성에 반비례함을 보여주는 바, 제2형은 더 큰 변형에서 나타나며 힘의 결정성 부족과 구조적 이질성이 뚜렷하다는 점이다.

ABSTRACT

Strain in granular materials in quasistatic conditions under varying stress originate in (I) contact deformation and (II) rearrangements of the contact network. Depending on sample history and applied load, either mechanism might dominate. One may thus define rheological regimes I and II accordingly. Their properties are presented and illustrated here with discrete numerical simulation results on sphere packings. Understanding the microscopic physical origin of strain enables one to clarify such issues as the existence of macroscopic elasticity, the approach to stress-strain relations in the large system limit and the sensitivity to noise.

연구 동기 및 목표

• 준정적 입자 변형에서 변형의 미세구조적 기원을 규명하여, 접촉 변형과 네트워크 재배열을 구분한다.
• 두 가지 별개의 유변학적 영역—제1형(탄성 접촉 변형이 지배함)과 제2형(접촉 네트워크 재배열이 지배함)을 정의하고 특성화한다.
• macroscopic 기계적 거동(탄성 및 플라스틱성 포함)이 접촉 네트워크 구조와 하중 역사를 어떻게 좌우하는지 조사한다.
• 힘의 결정성 부족과 협업 수준이 영역 간 전이를 결정하는 데 어떤 역할을 하는지 고찰한다.
• 하중 반전과 외란이 변형의 비가역성과 영역 지배도에 미치는 영향을 평가한다.

제안 방법

• Hertz-Mindlin 접촉 법칙을 적용한 단일 크기의 구형 입자 집합에 대한 이산 요소 방법(DEM) 시뮬레이션, 삼축 압축 조건 하에서 수행.
• 준정적 조건에서의 이탈을 정량화하기 위해 무차원 관성 수치 $ I = \dot{\epsilon}_{a} \sqrt{m/aP} $ 를 사용.
• 수치적 안정성을 확보하기 위해 정상 감쇠의 비율인 $ \zeta $ 를 사용한 법선 점성 감쇠 도입.
• 탄성 변형을 정량화하기 위해 무차원 접촉 강성 매개변수 $ \kappa = (E/P)^{2/3} $ 를 정의.
• 새로운 접촉이 형성되지 않는 고정된 접촉 목록을 사용한 전체 시뮬레이션과의 비교를 통해 제1형 거동를 분리 분석.
• 다양한 샘플 준비 방법(A 대비 C)에 따른 응력-변형 곡선, 탄성 모odulus, 힘 네트워크의 결정성 부족 정도를 분석.

실험 결과

연구 질문

• RQ1 입자 물질에서 변형이 주로 탄성 접촉 변형에 기인하는지, 아니면 접촉 네트워크 재배열에 기인하는지가 무엇에 의해 결정되는가?
• RQ2 샘플 준비 방식, 협업 수준, 하중 경로에 따라 제1형 또는 제2형의 지배 정도가 어떻게 달라지는가?
• RQ3 제1형에서 변형은 접촉 강성에 대해 어떻게 스케일링되는가? 제2형과의 차이는 무엇인가?
• RQ4 하중 반전은 입자 물질에서 비가역(플라스틱) 변형의 범위와 성격에 어떻게 영향을 미치는가?
• RQ5 힘의 결정성 부족과 협업 수준 등의 미세구조적 특성이 영역 간 전이를 얼마나 잘 예측할 수 있는가?

주요 결과

• 제1형은 탄성 접촉 변형이 지배하는 영역으로, 변형은 접촉 강성 $ E_{m} $ 에 반비례하며, 초도 하중 및 하중 반전 시 관찰된다.
• 제1형에서는 응력-변형 곡선의 초도 기울기가 매크로스코픽의 영률 $ E_m $ 와 일치하며, 이는 $ \kappa $ 와 비례한다. 또한 변형은 $ 1/E_m $ 에 비례한다.
• 하중 반전 시, 초도 해제 기울기는 $ E_m $ 와 동일하며, 비가역 변형은 약 $ \sim 0.2P $ 의 간격에서 발생한다. 이는 제1형 변형이 순수 탄성은 아님을 확인한다.
• 제2형에서는 변형의 크기가 접촉 강성에 민감하지 않으며, 오히려 입자 간 간극에 의존한다. 이 영역은 더 큰 변형에서 지속적인 네트워크 고장과 복구를 동반하며 지배한다.
• 영역 간 전이는 힘의 결정성 부족에 의해 규정된다: 제1형은 힘 네트워크가 고도로 초과정적(큰 힘 해 공간)일 때 발생하며, 제2형은 마찰성 슬립으로 인해 이 공간이 자유도의 약 10%로 감소할 때 나타난다.
• 새로운 접촉이 없는 고정된 접촉 목록을 사용한 시뮬레이션은 제1형에서 전체 시뮬레이션과 구분되지 않으며, 이는 이 영역에서 접촉 네트워크 변화가 무시할 만큼 미미하다는 것을 확인한다.