[논문 리뷰] Fast Fourier Transform computations and build-up of plastic deformation in 2D, elastic-perfectly plastic, pixelwise disordered porous media

François Willot, Yves-Patrick Pellegrini|ArXiv.org|2008. 02. 18.

Composite Material Mechanics참고 문헌 5인용 수 32

한 줄 요약

이 연구는 등방성 또는 비틀림 하중 하에서 2차원의 탄성-완전 플라스틱, 픽셀 수준에서 무질서한 다공성 미디어에서의 플라스틱 변형을 시뮬레이션하기 위해 이산 그린 함수를 활용한 빠른 푸리에 변환(FFT) 방법을 사용한다. 플라스틱 영역의 핵형성-성장-융합 메커니즘이 밝혀지며, 이는 거시적 응력-변형률 거동과 연관된다. 주요 통찰은 플라스틱 클러스터 형성과 밀도에 따른 응력 포화 현상에 기인한다.

ABSTRACT

Stress and strain fields in a two-dimensional pixelwise disordered system are computed by a Fast Fourier Transform method. The system, a model for a ductile damaged medium, consists of an elastic-perfectly matrix containing void pixels. Its behavior is investigated under equibiaxial or shear loading. We monitor the evolution with loading of plastically deformed zones, and we exhibit a nucleation / growth / coalescence scenario of the latter. Identification of plastic ``clusters'' is eased by using a discrete Green function implementing equilibrium and continuity at the level of one pixel. Observed morphological regimes are put into correspondence with some features of the macroscopic stress / strain curves.

연구 동기 및 목표

• 기계적 하중 하에서 무질서하고 다공성이며 탄성-완전 플라스틱인 2차원 미디어에서의 플라스틱 변형의 누적 메커니즘을 규명하는 것.
• FFT 기반 시뮬레이션에서 수렴성과 셀러브룩 패tern(체스보드 패턴)의 오류를 해결하기 위해 이산 그린 함수 수식을 도입하는 것.
• 플라스틱 클러스터 및 갈라진 영역-플라스틱 영역의 형태학적 변화를 거시적 응력-변형률 곡선과 연관짓는 것.
• 효율 매질 모델링의 대체 지표로 사용 가능한 기하학적 지표(예: 클러스터 수, 플라스틱 체적 분율 등)의 잠재력을 탐색하는 것.

제안 방법

• 푸리에 공간에서 리플만-슐레르거 방정식을 기반으로 한 FFT 방법을 채택하여, 역푸리에 변환을 통해 반복적으로 해를 구한다.
• 픽셀 수준에서 평형 및 일致 조건을 확보하기 위해 정방향/역방향 유한차분 기반의 이산 그린 함수를 사용하여 임의의 체스보드 패턴을 방지한다.
• 대칭성을 복원하고 수치적 안정성을 향상시키기 위해 이산 그린 함수(DG1)의 대칭형을 구현한다.
• 반복 해법의 수렴을 가속하기 위해 증강 라그랑주법과 우자와 알고리즘을 적용한다.
• 미세한 평형 조건을 기반으로 한 클러스터 식별을 통해, 응력의 등가 미제스 응력이 유동 응력 Y에 도달할 때 플라스틱 영역을 추적한다.
• 실증 공식을 활용하여 응력-변형률 거동을 클러스터 수 및 플라스틱 체적 분율과 같은 기하학적 지표와 연결한다.

실험 결과

연구 질문

• RQ1 등방성 또는 비틀림 하중 하에서, 무질서하고 다공성이며 탄성-완전 플라스틱인 2차원 미디어에서 플라스틱 변형은 어떻게 핵형성하고 진화하는가?
• RQ2 연속형 대비 이산형 그린 함수의 선택이 FFT 기반 플라스틱성 시뮬레이션의 정확도와 수렴성에 미치는 영향은 무엇인가?
• RQ3 플라스틱 클러스터의 형태학적 특징(예: 수, 크기, 융합)은 거시적 응력-변형률 곡선과 어떻게 관련되는가?
• RQ4 클러스터 수나 플라스틱 체적 분율과 같은 기하학적 지표가 다공성 미디어에서의 효과적 응력 반응을 예측할 수 있는 대체 지표로 기능할 수 있는가?
• RQ5 플라스틱 영역의 투과성(percolation)이 응력 포화를 유도하는 변형률 수준은 언제이며, 이는 응력-변형률 곡선의 이阶도 도함수와 어떻게 관련되는가?

주요 결과

• 플라스틱 변형의 핵형성-성장-융합 메커니즘이 관찰되었으며, 응력-변형률 곡선 상에서 명확한 영역이 식별되었다.
• 융합 속도는 고립된 플라스틱 영역에서 연결된 클러스터로의 전이 기간에 최고에 도달하며, 이는 거시적 응력-변형률 곡선의 이阶도 도함수와 강하게 상관된다.
• 플라스틱 클러스터 형성과 함께 강한 비틀림 변형 국소화가 발생하며, 변형장에서 프랙탈 유사 패턴이 나타난다.
• 거시적 유동 응력은 플라스틱 영역의 투과성 임계점에 도달할 무렵 최종적인 주요 크기로 수렴하며, 이는 응력 포화의 시작을 나타낸다.
• 실증 공식을 통해 기하학적 지표(클러스터 수는 곡률, 플라스틱 체적 분율은 포화도를 반영)를 사용해 응력-변형률 거동을 성공적으로 재현하였다.
• 단순 비틀림 하중의 경우 안정된 융합 영역이 사라지며, 낮은 다공성에서 직선적 비틀림 밴드가 형성되어 국소화 메커니즘이 다를 수 있음을 시사한다.