제4부 기계 설계·구조·재료

4.9 순응 기구: 노치·판스프링 유연 힌지·의사강체 모델(PRBM)·정밀 스테이지

4.9.1 이 장에서 배우는 것 / 선수 지식 / 특히 중요한 독자

이 장은 부재의 탄성 변형으로 운동을 전달하는 기구, 곧 순응 기구(compliant mechanism)의 설계를 다룬다. 회전 조인트를 얇은 부위의 굽힘으로 대체한 부위를 유연 힌지(flexure hinge)라 하며, 이 장은 원형 노치 힌지와 판스프링 힌지의 회전 강성과 항복 회전각, 대변형 외팔보를 하나의 강체 링크와 비틀림 스프링으로 근사하는 의사강체 모델(PRBM)의 계수, 평행사변형 판스프링 안내를 사용하는 정밀 스테이지의 강성·행정·기생 운동·고유 진동수를 계산으로 확인한다. 강체 링크 기구의 해석은 2.18, 말단 강성과 순응 행렬은 2.12, 구조 진동과 공진은 4.8, 정확 구속·키네마틱 커플링·열변형은 4.6, 피로 설계는 4.5, 재료 선정은 4.4, 직렬 탄성·가변 강성 액추에이터는 5.18, 압전·보이스 코일 같은 정밀 액추에이터는 5.22, 연속체·소프트 로봇의 동역학은 3.10, 소프트 로봇은 25.3이 맡는다.

배우는 것:

• 유연 힌지의 강성 — 오일러-베르누이 적분, 포물선 근사의 닫힌 식 $K_\theta\approx2Ebt^{5/2}/(9\pi\sqrt R)$, 판스프링 힌지 $Ebt^3/(12\ell)$, 2차원 유한요소와의 차이.
• 항복 회전각과 행정 — $\sigma=6M/(bt^2)$에서 $\theta_y=\sigma b t^2/(6K_tK_\theta)$, 노치 힌지는 같은 두께의 판스프링 힌지보다 강성이 큰 만큼 회전 범위가 좁다.
• 의사강체 모델 — 특성 반지름 계수 $\gamma$, 스프링 강성 계수 $K_\Theta$, 매개변수 각 계수 $c_\theta$, 탄성곡선 정확해에서의 추출, 선형 극한 $5/6$·$5/2$·$5/4$, 문헌 관용값 $0.85$·$2.65$·$1.24$의 오차.
• 정밀 스테이지 — 고정-안내 판스프링의 $12EI/L^3$ 강성, 응력으로 제한되는 행정 $\sigma L^2/(3Et)$, 기생 수축 $3y^2/(5L)$과 대변형 보정, 유효 질량 $13/35$.

선수 지식: 3.1(강체 동역학), 2.12(정역학과 강성), 4.8(진동과 구조 동특성의 $k/m$ 관계), 재료역학의 보 이론.

특히 중요한 독자: 힌지 없는 소형 기구·그리퍼·정밀 위치 결정 스테이지를 설계하는 기계 설계자, 정밀 구동계(압전·보이스 코일)와 기구를 맞추는 정밀 공학자, 수동 순응 요소로 접촉 힘을 제한하려는 로봇 설계자(13.10).

출처 방침. 기준 시점은 2026-11-05이다. 직접 읽은 자료: G. Chen, B. Xiong, X. Huang, Precision Engineering 35, 505–511(2011)의 서론과 1R PRBM 정의(특성 반지름 계수·스프링 강성·매개변수 각 계수, 힘 하중의 관용값 $0.85$·$2.65$·$1.24$, 순수 모멘트 하중의 $0.7346$·$2.0643$·$1.5164$; 저자 사본 PDF). 서지만 확인: Howell & Midha(1994·1995, J. Mech. Des. 116·117), Howell·Midha·Norton(1996), Smith의 Flexures(2000), Lobontiu의 Compliant Mechanisms(2002), Awtar & Slocum(2007)·Awtar·Slocum·Sevincer(2007), Schotborgh 외(2005)(Crossref 등재 형태). 읽지 못한 것: 위 서지들의 본문, Howell의 Compliant Mechanisms(2001) 교재 본문(Crossref에 도서 서지는 없다), Paros & Weisbord(1965)의 원 문헌(Crossref에 없다). 노치 힌지의 포물선 근사식 $2Ebt^{5/2}/(9\pi\sqrt R)$은 Paros·Weisbord의 근사식으로 알려져 있으나 이 장은 원문을 읽지 못했으며 식은 이 장이 직접 유도하고 계산으로 확인했다 (확인 필요). 이 장의 수치는 이 장이 정한 가정(알루미늄 $E=70$ GPa, 허용 응력 $250$ MPa)과 직접 실행한 계산에서 나온 것이며 실제 제품의 사양이 아니다.

4.9.2 개념 — 순응 기구는 변형을 기능으로 사용하므로 강성·행정·기생 운동을 함께 설계한다

결론부터 적으면 다음과 같다. 순응 기구는 조인트의 마찰과 백래시를 부재의 탄성 변형으로 대체하므로 정밀 위치 결정과 소형화에 유리하지만(Chen 외 2011의 서론이 같은 장점을 든다), 변형이 곧 응력이므로 행정이 항복 응력으로 제한되고, 설계하려는 자유도 외의 방향에도 작은 강성이 남아 기생 운동이 생긴다. 설계의 계산 대상은 세 가지이다. 구동 방향의 강성과 행정, 기생 운동과 구속 방향의 강성, 그리고 이동부 질량과 합쳐진 고유 진동수이다. 변형이 작으면 선형 보 이론으로 충분하고 대변형에서는 탄성곡선 해나 PRBM이 필요하다. PRBM은 탄성 부재를 강체 링크와 비틀림 스프링의 조합으로 바꿔 강체 기구의 해석 도구를 그대로 사용하게 하는 근사이다.

이 장의 계산 결과를 미리 적으면 다음과 같다. ① 원형 노치 힌지(폭 $10$ mm, 반지름 $R=2$ mm, 알루미늄 가정)의 회전 강성은 $t/R=0.1$에서 오일러-베르누이 적분 $0.634$ N·m/rad이고 포물선 근사 닫힌 식은 $t/R=0.05\sim0.8$에서 적분의 $0.994\sim0.900$배이다. ② 같은 형상의 2차원 평면 응력 유한요소는 두꺼운 양 끝을 포함한 오일러-베르누이 적분의 $0.981\sim0.805$배로 $t/R$이 클수록 보 이론이 강성을 과대평가한다. ③ 같은 최소 두께의 판스프링 힌지($\ell=2R$)보다 노치 힌지의 강성이 $5.4$배(t/R=0.1)이고 항복 회전각($K_t=1$)은 $1.51°$ 대 $8.19°$로 그만큼 작다. ④ PRBM 계수는 선형 극한에서 $\gamma=5/6$, $K_\Theta=5/2$, $c_\theta=5/4$이고(solve_bvp로 확인) 끝 처짐각 $70°$에서 $0.851$·$2.741$·$1.221$로 변한다. 문헌 관용값 $0.85$·$2.65$·$1.24$는 끝점 오차를 $\theta_0=10\sim70°$에서 $1.06\,\%$ 이하로 주지만 기울기 오차는 $-6.18\,\%$(10°)~$+2.89\,\%$(70°)이다. ⑤ 판스프링 2개 병렬의 스테이지($L=40$ mm, $t=0.5$ mm, $b=10$ mm)는 강성 $2{,}734$ N/m, 허용 응력 $250$ MPa에서 선형 행정 $3.81$ mm이다. ⑥ 안내 변위 $y/L=0.2$에서 기생 수축의 정확해는 선형식의 $1.0141$배, 구동 강성은 $1.0432$배이다. ⑦ 이동부 $0.1$ kg의 1차 고유 진동수는 $26.27$ Hz이고 유한요소와 $2.5\times10^{-7}$ 안에서 같다. ⑧ scikit-fem·SciPy와 독립 대조했다.

4.9.3 원리

4.9.3.1 유연 힌지의 설계 변수

유연 힌지는 한 방향의 회전을 허용하고 다른 방향의 이동·회전은 구속하도록 단면이 얇아진 부위이다. 결론은 다음과 같다. 원형 노치 힌지는 최소 두께 $t$, 노치 반지름 $R$, 폭 $b$로 정해지고, 회전 강성은 최소 두께의 $5/2$제곱에 비례하며 반지름의 제곱근에 반비례한다. 두께 분포 $h(x)=t+2\big(R-\sqrt{R^2-x^2}\big)$($\lvert x\rvert\le R$)에서 오일러-베르누이 보 이론(단면이 평면을 유지하고 전단 변형을 무시)의 회전 컴플라이언스는 다음이다.

$$\frac{\theta}{M}=\frac{12}{Eb}\int_{-R}^{R}\frac{dx}{h(x)^3},\qquad K_\theta=\frac{M}{\theta}.$$

노치 중앙 근처에서 $h\approx t+x^2/R$이므로 $\int_{-\infty}^{\infty}dx/(t+x^2/R)^3=\tfrac{3\pi}{8}\sqrt R\,t^{-5/2}$이고 닫힌 식이 나온다.

$$K_\theta\approx\frac{2Ebt^{5/2}}{9\pi\sqrt R}.$$

이 식은 $t/R$이 작을 때만 맞는다. 표 ①에서 적분 대비 비가 $t/R=0.05$에서 $0.994$, $0.2$에서 $0.974$, $0.8$에서 $0.900$이다. 판스프링 힌지는 길이 $\ell$의 균일한 얇은 보이며 $K_\theta=Ebt^3/(12\ell)$이다. 힌지 최소 단면의 공칭 굽힘 응력은 $\sigma=6M/(bt^2)$이므로 응력 집중 계수 $K_t$를 곱한 응력이 허용 응력 $\sigma_{allow}$에 이르는 회전각은 $\theta_y=\sigma_{allow}bt^2/(6K_tK_\theta)$이다. 판스프링 힌지는 $K_t=1$이고 $\theta_y=2\sigma_{allow}\ell/(Et)$이다. 같은 $b$·$t$·$\sigma_{allow}$에서 $\theta_y\propto1/K_\theta$이므로 강성이 큰 힌지는 회전 범위가 좁다. 노치 힌지의 $K_t$는 이 장이 계산하지 않았으며 표의 노치 값은 $K_t=1$의 상한이다 (확인 필요). 오일러-베르누이 보 이론은 전단 변형과 응력의 2차원 분포를 무시하므로 풀이 예제에서 2차원 유한요소와 비교했다. 유연 힌지의 회전 중심은 고정되지 않고 하중에 따라 이동하며(회전 정밀도), 이는 4.6의 정확 구속 논의와 연결된다.

4.9.3.2 의사강체 모델(PRBM)

끝단에 하중이 작용하는 외팔보를 PRBM으로 바꾸면 다음과 같다. 고정단에서 $(1-\gamma)l$ 위치에 회전 조인트(특성 피벗)를 두고 그 끝까지 길이 $\gamma l$의 강체 링크를 달며, 조인트에 비틀림 스프링 $k=\gamma K_\Theta EI/l$을 건다. $\gamma$는 특성 반지름 계수, $K_\Theta$는 스프링 강성 계수이다. 의사강체각 $\Theta$와 실제 끝 기울기 $\theta_0$의 비가 매개변수 각 계수 $c_\theta=\theta_0/\Theta$이다(Chen 외 2011). 이 장은 방향이 바뀌지 않는 수직 하중 $F$를 가정한다. 링크 끝의 좌표는 $(1-\gamma(1-\cos\Theta),\ \gamma\sin\Theta)\,l$이고 조인트의 모멘트 평형 $F\gamma l\cos\Theta=k\Theta$에서 하중 매개변수 $\alpha=Fl^2/EI$와 각이 다음 관계를 갖는다.

$$\alpha=\frac{K_\Theta\,\Theta}{\cos\Theta}.$$

PRBM의 정확도는 계수의 값에 달려 있으며, 정확해에서 읽는 방법은 다음과 같다. 탄성곡선의 끝 기울기가 $\theta_0$이면 굽힘 모멘트가 $EI\,d\theta/ds=F(x_{tip}-x)$이므로 $(d\theta/ds)^2=2(F/EI)(\sin\theta_0-\sin\theta)$이고, 보 길이 $l$ 대비 끝점 좌표와 하중은 다음과 같은 적분이다.

$$\alpha=\frac{I^2}{2},\quad \frac al=\frac{J_c}{I},\quad \frac bl=\frac{J_s}{I},\quad (I,J_c,J_s)=\int_0^{\theta_0}\frac{(1,\cos\theta,\sin\theta)\,d\theta}{\sqrt{\sin\theta_0-\sin\theta}}.$$

적분 끝점의 특이성은 $\theta=\theta_0-u^2$ 치환으로 없앤다(구현은 el_int). 끝점 좌표가 PRBM의 식과 같도록 $\tan(\Theta/2)=(1-a/l)/(b/l)$, $\gamma=(b/l)/\sin\Theta$, $K_\Theta=\alpha\cos\Theta/\Theta$, $c_\theta=\theta_0/\Theta$로 각 하중에서 계수를 읽는다. 작은 하중의 선형 극한은 선형 보 이론에서 정해진다. 끝 처짐 $b/l=\alpha/3$, 끝 기울기 $\theta_0=\alpha/2$, 끝 수축 $1-a/l=\tfrac35(b/l)^2$이므로 $\gamma=5/6$, $K_\Theta=3\gamma=5/2$, $c_\theta=3\gamma/2=5/4$이다. 표 ②는 이 값에서 하중과 함께 계수가 변함을 보인다. 문헌의 관용값 $\gamma=0.85$, $K_\Theta=2.65$, $c_\theta=1.24$는 이 장의 계산에서 대변형 구간($\theta_0$ 약 $45\sim70°$)의 값에 해당한다(Chen 외 2011은 이 값을 힘 하중에 "단순화를 위해 관용적으로 선택"한다고 적으며 선택 근거는 읽은 범위에 없다). 하중 방향이 달라지는 경우 계수가 방향각의 함수가 된다는 점은 같은 문헌이 기술한다. 이 장은 수직 하중의 경우만 계산했다.

4.9.3.3 고정-안내 판스프링과 평행사변형 스테이지

양 끝에서 기울기가 $0$으로 유지되는 판스프링(고정-안내)은 평행사변형 안내의 기본 요소이다. 선형 이론에서 길이 $L$·폭 $b$·두께 $t$의 판스프링 하나의 횡방향 강성은 $12EI/L^3=Ebt^3/L^3$이고, $n$개를 병렬로 놓은 스테이지의 강성은 $k=nEbt^3/L^3$이다. 끝에서의 최대 굽힘 모멘트는 $6EIy/L^2$이므로 판스프링 표면 응력은 $\sigma=3Ety/L^2$이고 허용 응력 $\sigma_{allow}$에서 선형 행정은 $y_{max}=\sigma_{allow}L^2/(3Et)$이다. 길이는 행정을 제곱으로, 두께는 반비례로 좌우한다. 안내 변위 $y$에서 판스프링의 양 끝 사이 거리가 줄어드는 기생 운동은 선형 이론에서 $3y^2/(5L)$이다. 이 장은 대변형의 정확해를 다음 관계로 구한다. 고정-안내 보는 중점에 대해 대칭이고 중점에서 모멘트가 $0$이므로 길이 $L/2$인 외팔보 둘이 같은 끝단 하중을 받는 것과 같다. 반 길이 외팔보의 끝점 $(a/l,\,b/l)$와 하중 $\alpha$에서 안내 변위 $y/L=b/l$, 기생 수축 $\Delta/L=1-a/l$, 구동 강성비 $k_{eff}/k_0=\alpha/(3b/l)$이다. 풀이 예제에서 이 관계는 전체 보의 경계값 문제와 $1.4\times10^{-11}$ 상대 차이에서 같다. 하중이 판스프링에 수직을 유지하는 이상적 경우의 결과이며 액추에이터와의 결합 방식에 따라 강성 증가는 달라진다.

고정-안내 판스프링 두 개를 쌓은 이중 평행사변형은 중간 몸체를 거쳐 기생 수축이 상쇄되는 구성으로 알려져 있다. 이 장은 이 성질을 계산하지 않았으며 Awtar와 Slocum의 연구가 이를 다룬다고 서지로만 확인했다 (확인 필요). 스테이지의 1차 고유 진동수는 이동부 질량 $m$과 판스프링 유효 질량으로 $f=\frac1{2\pi}\sqrt{k/(m+n\cdot\frac{13}{35}m_\ell)}$이다. $13/35$는 정적 처짐 형상 $3\xi^2-2\xi^3$으로 구한 질량 계수 $\int_0^1(3\xi^2-2\xi^3)^2d\xi$이고 풀이 예제에서 에르미트 보 유한요소와 $2.5\times10^{-7}$ 안에서 같다. 이 식은 4.8의 $k/m$ 관계의 한 예이며 고유 진동수가 낮으면 스테이지의 정착 시간이 길어진다.

4.9.3.4 정밀 스테이지의 설계 조합

스테이지의 행정 $y_{max}\propto L^2/t$와 강성 $k\propto bt^3/L^3$은 서로 반대 방향으로 움직인다. 행정을 늘리려고 판스프링을 길게 하면 강성과 고유 진동수가 떨어지고, 두께를 줄이면 행정은 늘지만 강성이 $t^3$으로 줄고 구속 방향의 강성도 같이 줄어 외란에 약해진다. 구동계(압전 액추에이터·보이스 코일, 5.22)는 스테이지 강성과 직렬 또는 병렬로 합쳐지므로 구동계의 강성이 스테이지보다 충분히 커야 행정 손실이 작다. 압전 액추에이터의 변위가 스테이지 행정보다 작으면 변위 증폭 기구(순응 레버)가 필요하며 이 장은 이를 계산하지 않았다 (확인 필요). 열팽창에 의한 위치 변화와 정확 구속 설계는 4.6이 다룬다.

4.9.4 알고리즘

입력: 요구 행정 y_req, 구동 방향 강성 k_req 또는 고유 진동수 f_req, 이동부 질량 m, 재료(E, 허용 응력), 구속 방향 요구
1. 형식 선택: 회전 범위가 작고 소형이면 노치 힌지, 범위가 크면 판스프링 힌지, 직선 안내이면 평행사변형(고정-안내 판스프링)
2. 선형 사이징: k = n E b t^3 / L^3, y_max = sigma L^2 / (3 E t), f = sqrt(k / (m + n 13/35 m_l)) / 2 pi
3. 판정: y_max >= y_req (안전 계수 포함), f >= f_req. 미충족이면 L·t·b·n을 조정(행정은 L^2/t, 강성은 b t^3/L^3)
4. 대변형 점검: y/L 가 0.1 이상이면 기생 수축·강성 증가를 탄성곡선 해 또는 PRBM으로 계산(PRBM 계수는 하중 범위에 맞춘다)
5. 노치 힌지: EB 적분으로 K_theta, 2차원 FE로 t/R 가 큰 경우 보정하고 응력 집중 K_t를 FE로 계산해 항복 회전각을 구한다
6. 검증: 유한요소(보 요소 비선형)로 강성·응력·기생 운동·모드를 확인하고 시작품의 강성과 고유 진동수를 측정한다
출력: 치수, 강성, 행정, 기생 운동, 고유 진동수, 피로 수명 확인 항목(4.5)

4.9.5 풀이 예제 — 노치 힌지, PRBM 계수, 평행사변형 스테이지

결론부터 적으면 다음과 같다. ① 포물선 근사 닫힌 식은 $t/R\le0.2$에서 적분과 $2.6\,\%$ 이내로 같고 오일러-베르누이 적분은 2차원 유한요소 대비 $t/R=0.1$에서 약 $4\,\%$, $t/R=0.4$에서 약 $14\,\%$ 강성을 높게 준다. ② PRBM 계수는 하중에 의존하며 문헌 관용값은 이 장의 계산에서 큰 처짐 구간의 값에 해당한다. ③ 스테이지의 기생 수축과 강성 증가는 $y/L\le0.1$에서 각각 $0.4\,\%$·$1.1\,\%$ 이내이다. ④ PRBM은 $\theta_0=45°$에서 끝점을 문헌 관용값 $0.51\,\%$, 이 장의 범위 평균 $0.02\,\%$ 안에서 맞춘다.

구현. k_notch·k_notch_par·k_leaf(힌지 강성), el_int·el_tip(탄성곡선 구적), prbm_par·prbm_tip(PRBM 계수와 예측), guided(고정-안내 판스프링), stage_k·stage_stroke·stage_f(스테이지). 표 생성·도해·캐시는 스크립트에 있다.

from math import sqrt, pi, sin, cos, atan, radians

def simpson(f, a, b, n=2000):
    h = (b - a) / n
    s = f(a) + f(b)
    for i in range(1, n):
        s += (4 if i % 2 else 2) * f(a + i * h)
    return s * h / 3

def notch_h(x, t, R):
    """원형 노치의 두께 h(x), |x| <= R"""
    return t + 2 * (R - sqrt(R * R - x * x))

def k_notch(E, b, t, R, n=4000):
    """굽힘 컴플라이언스 적분 theta/M = 12/(E b) int dx/h^3"""
    f = lambda x: notch_h(x, t, R) ** -3
    return E * b / (12 * simpson(f, -R, R, n))

def k_notch_par(E, b, t, R):
    """h = t + x^2/R 근사의 닫힌 식"""
    return 2 * E * b * t ** 2.5 / (9 * pi * sqrt(R))

def k_leaf(E, b, t, l):
    """균일 판스프링 힌지 E b t^3 / (12 l)"""
    return E * b * t ** 3 / (12 * l)

def el_int(th0, u0=0.0, n=1000):
    """끝단 하중 외팔보의 탄성곡선. th = th0 - u^2로 치환해
    (sin th0 - sin th)^(-1/2)의 특이성을 없앤다.
    (int g, int cos g, int sin g)를 반환.
    u0 = 0이면 전 구간, u0 > 0이면 th = 0 ~ th0 - u0^2 구간"""
    def g(u, w):
        if u == 0:
            return 2 * w(th0) / sqrt(cos(th0))
        q = 2 * cos(th0 - u * u / 2) * sin(u * u / 2)
        return 2 * u * w(th0 - u * u) / sqrt(q)
    ws = (lambda t: 1.0, cos, sin)
    r = sqrt(th0)
    return tuple(simpson(lambda u, w=w: g(u, w), u0, r, n)
                 for w in ws)

def el_tip(th0):
    """(alpha = F l^2/EI, a/l, b/l): 끝 처짐각 th0의 정확해"""
    I, Jc, Js = el_int(th0)
    return I * I / 2, Jc / I, Js / I

def th0_of(al):
    """alpha에 대한 끝 처짐각 (이분법)"""
    lo, hi = 1e-6, 1.45
    for _ in range(60):
        m = (lo + hi) / 2
        if el_tip(m)[0] < al:
            lo = m
        else:
            hi = m
    return (lo + hi) / 2

def prbm_par(th0):
    """정확해에서 읽은 PRBM 계수 (alpha, T, gamma, K, c)"""
    al, a, b = el_tip(th0)
    T = 2 * atan((1 - a) / b)
    return al, T, b / sin(T), al * cos(T) / T, th0 / T

def prbm_tip(al, gam, K, c):
    """PRBM의 끝점 (x/l, y/l)과 끝 기울기.
    평형 K T / cos T = alpha"""
    lo, hi = 0.0, 1.5
    for _ in range(80):
        m = (lo + hi) / 2
        if K * m / cos(m) < al:
            lo = m
        else:
            hi = m
    T = (lo + hi) / 2
    return 1 - gam * (1 - cos(T)), gam * sin(T), c * T

def guided(th0):
    """고정-안내 판스프링(길이 L): 반 길이 외팔보 둘의 합.
    (y/L, 기생 x/L, k_eff/k0)"""
    al, a, b = el_tip(th0)
    return b, 1 - a, al / (3 * b)

def th0_of_b(y):
    lo, hi = 1e-6, 1.45
    for _ in range(60):
        m = (lo + hi) / 2
        if el_tip(m)[2] < y:
            lo = m
        else:
            hi = m
    return (lo + hi) / 2

def stage_k(E, b, t, L, n=2):
    """판스프링 n개 병렬: 각 12 EI/L^3 = E b t^3/L^3"""
    return n * E * b * t ** 3 / L ** 3

def stage_stroke(sig, E, t, L):
    """끝에서 sigma = 3 E t y / L^2 = 허용 응력이 되는 y"""
    return sig * L * L / (3 * E * t)

def stage_f(k, m, ml, n=2):
    """고유 진동수. 판스프링 유효 질량 13/35 ml"""
    return sqrt(k / (m + n * 13 / 35 * ml)) / (2 * pi)

표 ① 원형 노치 힌지의 회전 강성과 항복 회전각.

표 ①에서 $t/R=0.1$이면 $K_\theta=0.634$ N·m/rad이고 같은 최소 두께의 판스프링 힌지($\ell=2R=4$ mm)는 $0.117$ N·m/rad로 노치가 $5.4$배 강하다. 항복 회전각은 노치 $1.51°$($K_t=1$의 상한), 판스프링 $8.19°$이다. 노치 힌지는 작은 크기에서 큰 강성과 정밀한 회전 중심을 주지만 회전 범위가 작다. 도해 2는 양대수 축에서 닫힌 식이 적분과 $t/R\le0.3$에서 시각적으로 겹치고 판스프링 힌지와는 약 한 자릿수 차이임을 보인다.

표 ② 정확해에서 읽은 PRBM 계수(수직 끝단 하중).

선형 극한에서 $\gamma=0.8333$, $K_\Theta=2.5$, $c_\theta=1.25$이며 $5°$에서 $0.8334$·$2.5007$·$1.2499$로 이 값에 이른다. 하중이 커지면 $\gamma$와 $K_\Theta$는 커지고 $c_\theta$는 작아진다($70°$: $0.8514$·$2.7408$·$1.2209$). 문헌 관용값 $0.85$·$2.65$·$1.24$는 $\gamma$가 $\theta_0\approx65\sim70°$, $K_\Theta$가 $60°$ 부근, $c_\theta$가 $45°$ 부근의 값에 해당한다. 이 장의 $5\sim70°$ 평균은 $0.8386$·$2.5720$·$1.2414$이며 이 값에서 $c_\theta$만 관용값과 같다. 정확해의 계수는 solve_bvp로 읽은 값과 $5.8\times10^{-12}$ 안에서 같다.

표 ③ PRBM의 오차(끝점 위치와 끝 기울기, 같은 하중 $\alpha$).

세 계수 집합은 모두 끝점 위치를 $\theta_0=10\sim70°$에서 $1.8\,\%$ 이하로 맞춘다. 가장 큰 차이는 기울기이다. 문헌 관용값은 $K_\Theta=2.65$가 선형 극한 $2.5$보다 커서 작은 하중에서 기울기를 $6.18\,\%$ 낮게 주고(10°) $70°$에서 $+2.89\,\%$이다. 선형 극한 계수는 작은 하중에서 $0.01\,\%$이지만 $70°$에서 끝점 $1.78\,\%$, 기울기 $+5.99\,\%$이다. 이 장의 범위 평균은 $45°$에서 끝점 $0.02\,\%$·기울기 $-0.11\,\%$이나 $70°$에서 $1.23\,\%$·$+4.17\,\%$이다. 어떤 집합도 전 구간에서 최선이 아니며 하중 범위에 맞추어 선택해야 한다. 이 결론은 이 장이 정한 오차 지표(끝점 위치, 끝 기울기)에 대한 것이다.

표 ④ 고정-안내 판스프링의 대변형 효과(L=40 mm, 정확해).

$y/L=0.05$(2 mm)에서 기생 수축이 $60.05$ $\mu$m이고 선형식 $60.00$ $\mu$m와 $0.09\,\%$ 차이이다. $y/L=0.1$(4 mm)에서 정확해가 선형식의 $1.0035$배, 구동 강성이 $1.0104$배이며 $y/L=0.2$(8 mm)에서 각각 $1.0141$배·$1.0432$배이다. 허용 응력으로 정해지는 선형 행정이 $3.81$ mm($y/L=0.095$)이므로 이 스테이지의 사용 범위에서 선형 이론의 오차는 기생 수축 $0.4\,\%$, 강성 $1.1\,\%$ 이내이다. 대변형 효과는 $y_{max}/L=\sigma L/(3Et)$가 커지는 얇고 긴 판스프링이나 허용 응력이 큰 재료의 설계에서 $y/L$이 $0.1$을 넘을 때 의미를 갖는다.

표 ⑤ 스테이지의 강성·행정·고유 진동수.

강성 $2{,}734$ N/m에서 이동부 $0.1$ kg의 고유 진동수가 $26.27$ Hz, $1$ kg이면 $8.32$ Hz이다. 판스프링 유효 질량($2\times13/35\times0.54$ g $=0.40$ g)은 $0.1$ kg 이동부에 비해 $0.4\,\%$로 작다. 행정 $3.81$ mm와 고유 진동수 $26$ Hz를 함께 만족하는 스테이지는 이동부가 가벼워야 하며, 이동부를 $1$ kg으로 늘리면 고유 진동수가 $8$ Hz로 낮아져 4.8의 정착 시간 식 $\ln(50)/(\zeta\omega_n)$에 따라 감쇠비 $0.02$에서 $2\,\%$ 정착 시간이 $0.1$ kg에서 $1.2$ s, $1$ kg에서 $3.7$ s이다.

라이브러리 독립 대조

대조 결과 이상 징후는 없었다. 2차원 유한요소는 격자를 $120\times10\to600\times20$으로 바꿀 때 $t/R=0.1$의 비가 $0.9647\to0.9621$로 수렴하며 이 장이 보고한 $400\times16$의 값(0.9622)은 수렴값과 $0.0001$ 이내이다. 유한요소와 오일러-베르누이 적분의 차이는 노치 힌지의 전단 변형과 응력 분포, 그리고 두꺼운 양 끝의 단부 효과가 함께 만든 것이며 이 장은 이 세 가지를 분리하지 않았다 (확인 필요).

한계. 선형 등방 탄성 재료와 평면 응력만 다루었다. 노치 힌지의 응력 집중 계수 $K_t$와 피로 수명, 노치 형상(타원·모서리 둥글림)의 비교, 힌지의 회전 중심 이동과 축 방향·전단 컴플라이언스, 판스프링의 측면 좌굴과 하중 강성화, 열팽창, 구동계(압전·변위 증폭 레버)의 결합, 이중 평행사변형의 기생 운동 상쇄는 계산하지 않았다 (확인 필요). PRBM은 수직 하중이 방향을 유지하는 끝단 하중의 경우만 다루었고 순수 모멘트나 하중 방향이 바뀌는 경우는 문헌 값(Chen 외 2011의 $0.7346$·$2.0643$·$1.5164$)을 인용하는 데 그친다. 재료 물성과 허용 응력은 이 장의 가정이다.

견해: 유연 힌지 설계의 초기 사이징에는 닫힌 식($K_\theta\approx2Ebt^{5/2}/(9\pi\sqrt R)$, $k=nEbt^3/L^3$, $y_{max}=\sigma L^2/(3Et)$)을 사용하고, $t/R$이 $0.2$ 이상이거나 $y/L$이 $0.1$ 이상인 설계는 2차원 유한요소와 탄성곡선 해로 보정하는 것이 맞다고 본다. 근거는 이 장의 계산에서 닫힌 식이 적분에서 $t/R=0.2$에 $2.6\,\%$, $0.8$에 $10\,\%$ 벗어나고 적분이 2차원 유한요소에서 $t/R=0.4$에 $14\,\%$ 벗어나며, 대변형 효과가 $y/L=0.1$에서 $1\,\%$ 수준이라는 점이다. 이 순서가 실제 시작품의 측정 강성과 일치하는지는 이 장이 검증하지 않았다.

4.9.6 설계 트레이드오프 — 해석 방법의 선택

유연 힌지와 순응 기구의 해석은 다음 중에서 고른다.

• ① 닫힌 식(선형 보 이론). 노치 힌지의 적분과 포물선 근사, 판스프링의 $12EI/L^3$처럼 즉시 계산하며 설계 변수의 영향을 해석적으로 본다. 변형이 크거나 $t/R$이 큰 경우와 복합 형상에서는 오차가 크다(표 ①, ④).
• ② PRBM. 대변형 외팔보와 순응 링크를 강체 링크와 스프링으로 바꿔 기구학·정역학 도구(2.18)를 그대로 사용한다. 계수가 하중에 의존하고 하중 방향이 바뀌면 계수가 변하므로 적용 범위를 정해야 한다(표 ③).
• ③ 유한요소·비선형 보 해석. 형상·응력 집중·대변형·모드를 모두 반영하나 모델링과 계산 비용이 들고 설계 변수의 영향을 해석적으로 보기 어렵다(4.14).

견해: 순응 기구의 설계는 ①로 치수를 정하고 ②로 기구학 개념을 확인한 뒤 ③으로 응력·모드·강성을 확정하며 시작품의 측정으로 보정하는 순서가 합리적이라고 본다. 근거는 이 장의 계산에서 닫힌 식과 PRBM이 각각 오차 한계(표 ①·③)를 가지고 있고 유한요소가 이를 $4\sim14\,\%$ 범위에서 보정함을 보였다는 점이다. 세 방법을 같은 시작품에 적용해 설계 시간과 정확도를 비교한 증거는 이 장에 없다.

4.9.7 응용 축 비교

이 표는 이 책의 정리이며 제품별 요구는 다르다. 응용별 수치는 이 장이 확인하지 않았다.

4.9.8 관련 규격·프로토콜

표준 번호 외에 판·조항은 이 장이 단정하지 않는다. 규격 체계 전반은 27.1이다.

4.9.9 관련 규제·인증 요구와의 연결

순응 기구를 직접 규정하는 로봇 규제는 이 장이 확인하지 못했다 (확인 필요). 간접 연결은 다음과 같다. 수동 순응 요소(탄성 힌지·스프링)로 접촉 힘과 충격을 제한하는 설계는 충돌 안전의 대책이 될 수 있으나(26.9) 이 장은 힘 한계를 계산하지 않았다. 유연 힌지의 피로 파손은 기능 안전 요건의 가정(4.5)을 깨뜨릴 수 있으므로 위험성 평가(26.4)에서 수명과 파손 모드를 고려한다. 힌지 두께의 공차가 강성을 크게 바꾸므로 도면의 공차 표기(4.7)가 성능 증빙의 일부이다.

4.9.10 한계·유의점·미해결

• 선형 등방·평면 응력. 재료 비선형, 이방성(적층 복합재), 소성, 크리프는 다루지 않았다.
• 응력과 피로. 노치의 응력 집중 계수와 피로 수명은 계산하지 않았다. 표의 항복 회전각은 $K_t=1$의 상한이며 실제 범위는 이보다 좁다.
• PRBM의 범위. 수직 끝단 하중에 한정되며 하중 방향이 바뀌는 순응 기구, 순수 모멘트, 변곡점이 있는 변형은 다른 계수 또는 다중 링크 PRBM이 필요하다(Chen 외 2011의 3R PRBM은 읽은 범위의 소개).
• 스테이지의 결합 효과. 구동계 결합, 변위 증폭, 이중 평행사변형, 열 변형, 측면 좌굴은 계산하지 않았다.
• 회전 중심 정밀도. 유연 힌지의 회전 중심 이동과 축 방향 컴플라이언스는 계산하지 않았다.
• 미해결. 노치 형상(타원·모서리 둥글림·하이브리드)의 비교, 설계 변수에 대한 응력 집중 계수의 닫힌 식, 순응 기구의 최적 설계(위상 최적화)와 제조 공차의 영향을 이 장은 다루지 않았다(4.4의 위상 최적화).

4.9.11 체크포인트·연구 질문·출처

체크포인트:

1. $h\approx t+x^2/R$에서 $\int dx/h^3=\tfrac{3\pi}{8}\sqrt R\,t^{-5/2}$을 증명하고 $K_\theta\approx2Ebt^{5/2}/(9\pi\sqrt R)$를 유도하라.
2. 선형 극한에서 $b/l=\alpha/3$, $\theta_0=\alpha/2$, $1-a/l=\tfrac35(b/l)^2$로부터 $\gamma=5/6$, $K_\Theta=5/2$, $c_\theta=5/4$를 유도하라.
3. PRBM의 평형식 $\alpha=K_\Theta\Theta/\cos\Theta$를 하중 방향이 불변이라는 가정과 함께 유도하라.
4. 고정-안내 판스프링이 반 길이 외팔보 둘과 같은 이유를 중점의 모멘트로 설명하고 $k_{eff}/k_0=\alpha/(3b/l)$를 유도하라.
5. $L=60$ mm, $t=0.4$ mm, $b=12$ mm인 알루미늄 스테이지(판스프링 2개)의 강성·선형 행정·1차 고유 진동수(이동부 $0.2$ kg)를 이 장의 식으로 구하라.
6. 노치 힌지의 항복 회전각이 판스프링 힌지의 $1/5.4$인 이유를 $\theta_y\propto1/K_\theta$로 설명하고, 같은 회전 범위를 요구할 때 노치 힌지의 $t$를 어떻게 바꿔야 하는지 논하라.

연구 질문:

• 노치 힌지의 응력 집중 계수 $K_t$가 $t/R$의 함수로 어떻게 변하며, 이를 포함하면 이 장의 항복 회전각 표가 어떻게 달라지는가(이 장은 계산하지 않았다)?
• 하중 방향이 변하는 순응 기구에서 PRBM의 계수를 하중 이력의 함수로 어떻게 정하고 정확도를 어떻게 보증하는가?
• 압전 구동과 변위 증폭 레버를 포함한 정밀 스테이지의 행정·강성·고유 진동수 설계 공간은 어떤 모양이며 열 변형이 이를 어떻게 제한하는가?

1차 출처(조회 2026-11-05).

• G. Chen, B. Xiong, X. Huang, "Finding the optimal characteristic parameters for 3R pseudo-rigid-body model using an improved particle swarm optimizer," Precision Engineering 35, 505–511, 2011, doi:10.1016/j.precisioneng.2011.02.006 — 서론의 1R PRBM 정의와 계수 값을 읽었다.
• L. L. Howell, A. Midha, "A method for the design of compliant mechanisms with small-length flexural pivots," J. Mech. Des. 116, 280–290, 1994, doi:10.1115/1.2919359 — 서지만 확인. Howell·Midha, "Parametric deflection approximations for end-loaded, large-deflection beams in compliant mechanisms," J. Mech. Des. 117, 156–165, 1995, doi:10.1115/1.2826101 — 서지만 확인. Howell·Midha·Norton, "Evaluation of equivalent spring stiffness for use in a pseudo-rigid-body model of large-deflection compliant mechanisms," J. Mech. Des. 118, 126–131, 1996, doi:10.1115/1.2826843 — 서지만 확인.
• S. T. Smith, Flexures: Elements of Elastic Mechanisms, Gordon & Breach/CRC Press, 2000, doi:10.1201/9781482282962 — 서지만 확인. N. Lobontiu, Compliant Mechanisms: Design of Flexure Hinges, CRC Press, 2002, doi:10.1201/9781420040272 — 서지만 확인.
• S. Awtar, A. H. Slocum, "Constraint-based design of parallel kinematic XY flexure mechanisms," J. Mech. Des. 129, 816–830, 2007, doi:10.1115/1.2735342 — 서지만 확인. Awtar·Slocum·Sevincer, "Characteristics of beam-based flexure modules," J. Mech. Des. 129, 625–639, 2007, doi:10.1115/1.2717231 — 서지만 확인. Schotborgh 외, "Dimensionless design graphs for flexure elements and a comparison between three flexure elements," Precision Engineering 29, 41–47, 2005, doi:10.1016/j.precisioneng.2004.04.003 — 서지만 확인.
• L. L. Howell, Compliant Mechanisms, Wiley, 2001 — 서지를 Crossref에서 확인하지 못했다. Paros·Weisbord(1965)의 원 문헌도 Crossref에서 확인하지 못했다 (확인 필요).
• scikit-fem 12.0.2(MeshQuad·ElementQuad2·ElementLineHermite), SciPy 1.18.1(integrate.quad·integrate.solve_bvp·linalg.eigh), numpy 2.5.4.

다음 장. 4.10는 정적 균형과 중력 보상 기구(스프링·카운터웨이트·밸런서)를 다룬다. 이 장의 스프링 강성과 행정 계산은 중력 보상 스프링의 설계로 이어진다.