비선형·비국소 탄성으로 보는 상미분 장력 모델의 연속체 해석

초록

이 논문은 세포 수준에서 서로 다른 장력을 갖는 두께가 일정한 상미분 장력 모델을 연속체 이론으로 전개하고, 압축에 의해 발생하는 상피 조직의 버클링을 분석한다. 비선형 및 비국소 탄성 항이 나타나며, 버클링 임계값은 차등 장력에 무관함을 보인다.

상세 분석

본 연구는 기존의 이산 셀 기반 차등 장력 모델을 연속체 수준으로 체계적으로 전이시키는 데 초점을 둔다. 먼저, 등면적을 유지하는 등변 사다리꼴 형태의 세포를 가정하고, 각 면에 부여된 장력 Γₐ, Γ_b, Γ을 통해 단일 세포의 에너지를 식 (1)·(2) 로 정의한다. 비차원화 과정에서 세포 면적 √A_c 로 길이를 스케일링하고, λ = L/√A_c, κ = K/√A_c 등 무차원 파라미터를 도입한다. 여기서 α = Γₐ/Γ, β = Γ_b/Γ 로 차등 장력 비율을 나타내며, δ = α−β 가 차등 장력의 비대칭성을 담당한다. 중요한 점은 λ = O(ℓ₀) 와 φ = O(ℓ₀⁻²) 라는 스케일링 가정이다. 이는 세포가 평탄 상태에서 작은 변형을 허용하면서도, φ가 큰 경우(즉, 급격한 기울기)에도 적용 가능하도록 설계된 것이다.

연속체 한계에서는 φ(s) 를 중선 아크길이 s 의 함수로 두고, 전체 에너지 E = ∫ e(φ) ds 로 적분한다. 여기서 e(φ) 는 식 (6) 에서 도출된 비차원 에너지 밀도이며, sec φ 와 sin φ 항이 포함되어 차등 장력 δ 가 비선형적으로 작용한다. φ와 중선 기울기 ψ 사이의 관계는 복잡한 비선형 연산자를 통해 전개된다. 저자들은 Bernoulli 수를 이용한 무한 급수를 도출하고, 이를 역전시켜 φ를 ψ의 도함수들의 급수로 표현한다(식 13). 이는 기존 연구에서 1차 항만 고려했던 것과 달리 2차 항까지 포함함으로써 비선형 및 비국소 효과를 포착한다.

다음 단계에서는 중선 좌표 (x(s), y(s)) 를 ψ와 연결시키는 미분 방정식(식 15·16)을 얻는다. 여기서 f와 g는 ψ와 그 도함수들의 조합으로, O(ℓ₀⁻²) 수준에서 비선형 및 비국소 항이 나타난다. 특히 g 항은 ψ’’ 와 ψ’⁴ 형태로, 전통적인 Euler‑Elastica와는 달리 고차 미분항이 포함된 비국소 탄성 효과를 의미한다.

버클링 해석을 위해 압축 비율 D와 수평 압축력 μ를 라그랑주 승수로 도입하고, 에너지 라그랑지안을 전개한다(식 20). 여기서 ℵ = ℓ₀ Σ 로 정의된 무차원 파라미터가 셀 수 N = 2ℵ 와 직접 연결된다. 차등 장력 δ는 Δ = δℓ₀² 로 결합되어, 최종 지배 방정식(식 21)에서는 Δ와 Λ(세포 두께 비율)만이 등장한다. 이 방정식은 ψ와 그 1·2·3·4 차 미분항이 모두 포함된 비선형 비국소 편미분 방정식이며, 전통적인 선형 탄성 이론과는 근본적으로 다르다.

버클링 임계값을 찾기 위해 작은 진폭 ε 전개를 수행한다. 1차 항에서 얻은 고유값 μ₀ = z₁−z (식 27)와 ψ₀(σ)=Ψ₀ sinπσ 로부터 기본 모드가 결정된다. 2차와 3차 항에서는 비선형 상호작용과 차등 장력 Δ가 나타나며, 특히 μ₁=0 로부터 초임계(supercritical) 버클링임을 확인한다. 최종적으로 에너지 비교를 통해 버클링 임계 압축 D* = 1−√(1−z) 를 도출한다. 흥미롭게도 D*는 Δ(즉, 차등 장력)와 무관하게 결정된다. 이는 차등 장력이 버클링 형태와 포스트버클링 비선형성에 영향을 주지만, 임계점 자체는 전형적인 Euler‑buckling과 동일한 메커니즘에 의해 좌우된다는 의미다.

전체 해석 과정에서 비국소 항(예: ψ’’·ψ’’’ 등)과 비선형 항(ψ⁴, ψ’⁴ 등)이 동시에 등장함으로써, 연속체 모델이 단순한 탄성 막이 아니라, 세포 수준의 장력 불균형이 매크로 스케일에서 비선형·비국소 탄성 효과로 전이되는 메커니즘을 명확히 보여준다. 이는 기존의 차등 장력 모델이 연속체 한계에서 어떻게 복합적인 물리적 거동을 나타내는지를 수학적으로 정량화한 최초의 시도라 할 수 있다.