점탄소성 유동 모델의 볼록 해석과 통합

초록

볼록 해석 도구를 이용해 선형 점성, 완전 소성, 그리고 비선형 파워‑법칙 점성을 조합한 점탄소성 소산 퍼텐셜을 체계적으로 유도한다. 병렬·직렬 결합의 수학적 구조를 인피멀 컨볼루션으로 표현하고, 기존 경험적 조화 평균 모델과의 차이를 명확히 밝힌다.

상세 분석

본 논문은 점탄소성 재료의 응답을 하나의 볼록 소산 퍼텐셜 ζ_vp(ε) 로 기술하는 체계적 프레임워크를 제시한다. 먼저 완전 소성(1차 동차 퍼텐셜 ζ₁(·)=σ_a|·|)과 선형 점성(2차 동차 퍼텐셜 ζ₂(·)=½D|·|²)을 병렬로 연결하면 σ∈∂ζ₁(ε)+∂ζ₂(ε) 즉 σ∈∂(ζ₁+ζ₂)(ε) 가 된다. 이는 전통적인 Bingham 모델과 동일하며, σ=σ_a ε/|ε|+Dε 로 표현된다. 직렬 결합에서는 전체 변형률을 ε=ε₁+ε₂ 로 분해하고, 각각 ζ₁′(ε₁)=σ, ζ₂′(ε₂)=σ 로 두어 인피멀 컨볼루션 ζ_vp=ζ₁□ζ₂ 로 정의한다. 이때 ζ₁□ζ₂는 Yosida 근사이며, 결과적으로 ζ_vp는 매끄러운 퍼텐셜이 되면서 σ와 ε 사이의 비선형 관계를 제공한다. 특히, ζ₁□ζ₂의 명시적 형태는 Hub­er 함수로, 작은 변형률에서는 순수 점성(크리프) 거동, 큰 변형률에서는 급속 슬립(플라스틱) 거동을 하나의 식으로 연결한다.

세 개 요소(두 점성 + 하나의 소성)를 조합한 경우, 두 가지 설계(그림 4‑좌·우)가 서로 등가임을 증명한다. 좌측 구성은 ζ_vp= (ζ₁□ζ₂)+ζ₃ 형태이며, ζ₃는 추가 점성 항이다. 우측 구성은 eζ₁+eζ₃ 를 먼저 합친 뒤 eζ₂와 인피멀 컨볼루션을 수행해 동일한 ζ_vp를 얻는다. 이때 파라미터 변환 관계 eσ_a=σ_a/(1+D₃/D₂), eD₂=D₂+D₃, eD₃=D₃/(1+D₃/D₂) 가 성립한다. 결과적으로 유효 점도 μ_eff(ε)=min{σ_a|ε|,D₂}+D₃ 로, 변형률이 작을 때는 소성 한계가 지배하고, 크게 되면 점성 항이 우세한다.

다음으로 경험적으로 널리 사용되는 조화 평균식 μ_eff⁻¹=∑μ_i⁻¹ 를 검토한다. 볼록 해석에 기반한 직렬·병렬 결합에서 얻은 μ_eff와는 일반적으로 다르며, 특히 비선형 점성(μ_i가 ε_i에 의존)인 경우 차이가 두드러진다. 논문은 이를 예시(Example 1)와 식(15)·식(3)의 차이로 명확히 보여준다.

마지막으로 비선형 파워‑법칙 점성(Norton‑Hoff) ζ(ε)=n/(n+1)D|ε|^{1+1/n} 을 도입한다. 이 퍼텐셜은 n>1이면 전단‑증가(전단‑두꺼워짐), n<1이면 전단‑감소(전단‑희석) 거동을 나타낸다. 이러한 비선형 퍼텐셜을 기존 소성 요소와 인피멀 컨볼루션으로 결합하면, 전통적 선형 모델을 포함하는 보다 일반적인 점탄소성 프레임워크가 완성된다. 전체적으로 논문은 볼록 해석을 통해 점성·소성·비선형 점성 요소들을 수학적으로 일관된 하나의 퍼텐셜로 통합하고, 기존 경험적 모델과의 차이를 정량적으로 제시함으로써 지구물리·재료공학 분야에서의 모델링 정확성을 크게 향상시킬 수 있음을 입증한다.