Helical Tube - 2026-04-16

날짜: 2026-04-16
작성자: 최원재
프로젝트: Helical Tube

---

목적

진행 내용

버클링(buckling) = 압축을 받는 가늘고 긴 부재가 어느 한계 하중을 넘으면 갑자기 옆으로 휘어져 축력을 더 이상 전달하지 못하는 현상입니다. 케이블을 밀어 넣을 때 이것이 핵심 제약이 됩니다.

1. Euler 의 고전 버클링 양 끝이 단순지지된 길이 $\ell$ 의 기둥이 압축하중 $P$ 를 받을 때, 임계하중:

$$P_{cr}=\frac{{\pi }^{2}EI}{{\ell }^2}$$

$EI$ = 굽힘 강성 (E=영률, I=단면 2차모멘트) $\ell$ = 자유 좌굴 길이(buckling half-wavelength) $P>P_{cr}$ 이면 부재가 반파장 $\ell$ 로 휘어진다. 핵심 통찰: $\ell$ 이 짧을수록 $P_{cr}$ 이 커진다 ($1/{\ell }^2$). 즉 짧게 지지되면 잘 안 꺾인다.

2. 케이블이 튜브 안에 있을 때 — "Self-contact buckling" 케이블이 아무 제약 없는 자유공간이라면 $\ell$ = 전체 길이. 하지만 튜브 안에서는 벽과의 간극 $\delta$ 가 지지 역할을 합니다:

         tube wall
   

   ───────────────────────── ┌─ ┐ ┌─ ┐ ← 케이블이 사인파로 휘어짐 ──┘ └─────┘ └─── 반파장 λ ┌─ ┐ ┌─ ┐ (압축 P) ───────────────────────── tube wall ← λ → 케이블이 사인파로 휘어져 튜브 벽에 닿기 시작하면, 반파장 $\lambda$ 가 자연스럽게 정해집니다. 이 $\lambda$ 가 유효 지지 길이가 되어 $P_{cr}$ 을 결정.

에너지 균형: 굽힘에 저장되는 탄성에너지와 간극 $\delta$ 만큼 옆으로 이동하는 데 필요한 기하적 일이 같아지는 조건에서

$${\lambda }^4=\frac{{\pi }^{4}EI\cdot \delta \cdot \rho }{P}$$

유도 개요 (회전자와 계수 생략):

사인파 변형 $y(s)=\delta \sin (\pi s/\lambda )$ 굽힘에너지 밀도 ∝ $EI/{\lambda }^2$ 경로 곡률 보정 ∝ $1/\rho$ 압축일 ∝ $P\delta /\lambda$ 이를 pull_vs_push.py 72–75줄에 구현:

def P_crit(P): lam = (np.pi4 * EI_c * delta_cable * rho / P)0.25 return np.pi2 * EI_c / lam2 3. 왜 $P_{crit}$ 자체가 $P$ 에 의존하는가 (재귀 구조) 위 식을 합치면:

$$\lambda ={\left(\frac{{\pi }^{4}EI,\delta ,\rho }{P}\right)}^{1/4},P_{cr}=\frac{{\pi }^{2}EI}{{\lambda }^2}$$

$\lambda$ 를 소거하면:

$$P_{cr}={\pi }^{2}EI\cdot {\left(\frac{P}{{\pi }^{4}EI,\delta ,\rho }\right)}^{1/2}=\sqrt{\frac{EI,P}{\delta ,\rho }}$$

흥미로운 결과: $P_{cr}$ 이 $\sqrt{P}$ 에 비례. 즉 압축력이 커질수록 한계도 같이 올라가지만 더 느리게 (제곱근으로). 어느 지점에서 $P>P_{cr}$ 가 되면 좌굴 시작.

좌굴 조건 $P=P_{cr}$ 를 쓰면:

$$P_{buckle}=\frac{EI}{\delta ,\rho }$$

이것이 "자기접촉 좌굴"의 특성 압축력입니다.

4. 본 문제에서의 수치 (PEEK 케이블)

EI = 1.44 N·m² (PEEK, OD 9.5 mm) δ = 2.5 mm (= 12 mm - 9.5 mm, 직경 간극) ρ = 1.652 m (헬릭스 곡률반경)

$$P_{buckle}\approx \frac{EI}{\delta \rho }=\frac{1.44}{0.0025\times 1.652}\approx 349\text{N}$$

케이블에 349 N 이상의 압축력을 주면 튜브 안에서 좌굴. pull_vs_push.py 의 PUSH 해석은 P(0) = 6563 N → 349 N을 압도적으로 초과 → 튜브 초입부터 즉시 좌굴. 출력: Buckling first at s = 0.00 m (시작점부터 이미 불가). 5. 버클링 이후에는 어떻게 되나 — Post-buckling 한 번 좌굴되면 케이블은:

사인 좌굴(sinusoidal buckling) — 튜브 벽과 한쪽만 닿는 얕은 웨이브. $P$ 가 더 커지면 헬리컬 좌굴(helical buckling) — 튜브 내벽을 나선형으로 감싸며 붙어버림. 이후에는 마찰이 급격히 증가 — 새로 생긴 벽 접촉점마다 수직항력 발생. 결국 lock-up — 어떤 압축력을 더해도 전진 못 함. 석유 시추·카테터 문헌 (Dawson-Paslay, Wu-Juvkam-Wold, He-Kyllingstad) 에서 이 단계를 정식으로 모델링합니다.

핵심: 현재 모델은 1단계 (사인 좌굴 시작) 까지만 본다. 이후 post-buckling 은 훨씬 더 복잡하며, 실제로는 좌굴이 시작되자마자 PUSH 가 실패한다고 보면 됩니다.

6. 왜 PULL 은 이 제약이 없나 PULL 에서 케이블은 인장 상태. 인장된 줄은 절대 좌굴하지 않습니다 (인장이 교란을 복원시키는 방향). 대신 인장 파단(breaking strength) 만 체크하면 됩니다.

항목 PULL PUSH 케이블 상태 인장 압축 실패 모드 파단 (7089 N) 좌굴 (349 N) 한계 스케일링 ${\sigma }_y\cdot A$ $EI/(\delta \rho )$ 한계 완화 방법 고강도 재료, 큰 단면적 큰 EI / 작은 δ / 큰 ρ 7. 설계 관점의 함의 $P_{buckle}=EI/(\delta \rho )$ 식에서 PUSH 한계를 높이려면:

변수 효과 구체적 방법 EI ↑ 강성 높은 케이블 금속 와이어, 내장 push-rod, PEEK 대신 탄소섬유 복합재 δ ↓ 튜브-케이블 간극 축소 케이블 지름을 튜브 지름에 근접시키기 (단, 조립/삽입 난이도 ↑) ρ ↑ 완만한 곡률 튜브 설계 단계에서 조정 (R 증가, pitch 감소) 예를 들어 PEEK 케이블을 외경 11 mm (δ=1 mm) 로 키우면 $P_{buckle}=1.44/(0.001\times 1.652)=872$ N 으로 ≈2.5배 개선. EI 를 두 배로 키우면 추가 ×2.

그래도 6563 N 은 못 버티므로, 본 헬리컬에서는 μ 를 낮춰 PUSH 힘 자체를 줄이는 것이 유일한 해결책입니다.

TODO

• [ ]

참고