환형 전단 실험은 입자 재료(예: 토양, 모래)가 전단력 하에서 어떻게 거동하는지 분석하는 핵심 방법입니다. PFC3D5.0을 활용해 이 실험을 자동화하는 완전한 코드 프레임워크를 제공하며, 모델 구축부터 데이터 추출까지 전 과정을 다룹니다. 아래 내용은 실제 시뮬레이션에서 바로 적용 가능한 기술적 구현에 초점을 맞춥니다.

모델 구축: 3단계 접근법

모델링은 크게 세 단계로 구성됩니다: 벽 생성, 입자 배치, 경계 조건 설정. 원통형 전단 상자를 사용하며, 하단 벽은 고정하고 상단 벽에 회전 운동을 적용합니다. 초기에는 입자의 마찰 계수를 낮게 설정하여 추후 압밀 과정을 용이하게 만듭니다.

; 벽 생성
wall generate id 1 cylinder ...    ; 하단 고정 벽
wall generate id 2 cylinder ...    ; 상단 회전 벽
wall attribute damp 0.7            ; 벽 감쇠 계수

; 입자 분포
ball distribute ...
ball attribute density 2650         ; 예: 석영 밀도

입자 분포 후 목표 간극비까지 압밀이 필요합니다. 이 과정은 서보 제어(servo-control)를 사용하여 가장 안정적으로 수행할 수 있습니다. 서보 제어 중에는 접촉력을 지속적으로 모니터링하여 입자가 과도하게 변형되지 않도록 해야 합니다. 압밀이 완료되면 마찰 계수를 실제 실험 조건에 맞게 조정합니다.

주요 제어 매개변수

전단 속도는 수치적 불안정성을 방지하기 위해 매우 중요합니다. 준정적(quasi-static) 조건을 모방하기 위해 전단 속도는 1e-5 ~ 1e-4 스케일로 설정합니다. 응력 제어 모드를 사용하여 수직 응력을 가하는 것이 변위 제어보다 실제 실험에 더 가깝습니다.

; 응력 제어(서보) 함수 예시
def apply_confining_stress
    loop local i (1, wall.num)
        if wall.id(i) == 2
            wall.force.app = -target_stress * area
        endif
    endloop
end

핵심 Fish 함수: 데이터 수집 및 분석

전단 중 다섯 가지 주요 지표를 실시간으로 기록해야 합니다. 전단 응력, 수직 변위, 부피 변화, 운동 에너지, 겉보기 마찰 계수. 벽의 회전 각도를 활용하여 전단 변형률을 계산하는 것이 효율적인 방법입니다.

; 전단 응력 계산 함수
def calc_shear_stress
    shear_force = 0.0
    loop foreach cnt wall.contact.list
        if wall.isin(cnt.end1)   ; 상단 벽과의 접촉만 고려
            shear_force += cnt.force.shear
        endif
    endloop
    current_shear = shear_force / area
end

부피 변화 계산 시 벽의 움직임으로 인한 영향을 배제해야 합니다. 입자의 중심(centroid) 좌표를 사용하여 유효 영역 내 부피를 통계적으로 계산하는 것을 권장합니다. 운동 에너지 모니터링은 시스템이 안정적인 상태에 도달했는지 판단하는 데 도움이 되며, 이는 데이터 수집 시점을 결정하는 중요한 기준입니다.

결과 해석: 전문가적 접근

겉보기 마찰 계수(전단 응력 ÷ 수직 응력) 곡선이 안정적인 평탄 구간(plateau)에 도달하면 재료가 임계 상태(critical state)에 진입했음을 의미합니다. 이 시점에서 부피 변화 곡선도 안정화되어야 합니다. 특히 밀집된 모래는 팽창(dilatancy)을, 느슨한 모래는 수축을 보일 수 있으므로 주의 깊게 관찰해야 합니다.

수직 응력이 일정하게 유지되는 것은 서보 제어의 성공 여부를 판단하는 중요한 지표입니다. 수직 응력 변동 폭이 5%를 초과할 경우, 벽의 강성 및 감쇠 매개변수를 재검토해야 합니다. 전체 에너지 중 운동 에너지 비율은 항상 1% 미만이어야 하며, 그렇지 않다면 전단 속도가 너무 빠른 것은 아닌지 의심해야 합니다.

코드 패키지 활용 팁

제공되는 코드에는 세 가지 기본 조건(밀집 모래, 느슨한 모래, 과도 상태)이 사전 설정되어 있습니다. 다음 세 가지 변수만 조정하여 케이스를 전환할 수 있습니다.

• 초기 간극비 (0.5 ~ 0.8)
• 입도 분포 계수 (0.1 ~ 0.3)
• 구속압 (50 kPa ~ 1 MPa)

시뮬레이션 완료 시 다음 네 개의 데이터 파일이 자동 생성됩니다.

• 응력-변형률 곡선.dat
• 부피 변화.dat
• 에너지 기록.dat
• 실시간 매개변수.log

시뮬레이션 중 입자가 비산하는 문제가 발생하면, 일반적으로 감쇠 계수가 너무 작은 경우입니다. 전단 단계에서는 국부 감쇠(local damping)를 0.3 정도로 설정하여 수렴 속도와 역학적 응답 사이의 균형을 유지하는 것이 좋습니다. 입자 흐름 시뮬레이션은 마치 요리와 같아서, 적절한 타임스텝(time-step) 설정이 결과의 질을 결정합니다.