비틀림 완전 가이드: Tρ/J부터 얇은벽 단면·뒤틀림과 FEA 검증까지

토크에서 전단응력과 비틀림각을 구하는 원리부터 원형축의 재료 효율, 비원형 단면의 뒤틀림, Bredt 전단흐름, 개단면의 구속 뒤틀림과 FEA 검증까지 단계별로 연결합니다.

모터축의 지름을 정하고, 박스 거더의 벽 두께를 배치하고, 채널 보가 예상보다 크게 회전한 이유를 찾는 문제는 모두 비틀림으로 연결됩니다. 그러나 현장에서 자주 쓰는 τ=Tρ/J는 원형축에만 그대로 적용할 수 있습니다. 직사각형 봉은 단면이 뒤틀리고, 얇은 폐단면에서는 재료가 단면 중심보다 벽의 연결 경로를 따라 토크를 전달하며, 개단면은 끝단의 뒤틀림 구속에 따라 전혀 다른 축응력까지 만듭니다.

이 글은 비틀림을 하나의 공식으로 축약하지 않습니다. 외부 토크 T에서 내부 전단응력, 비틀림률 θ′, 단면의 비틀림상수 J, 얇은벽 전단흐름 q, 구속 뒤틀림과 FEA 결과까지 물리적으로 이어지는 계산 사슬을 만듭니다. 모든 수치 예제는 N–mm–MPa 단위계로 다시 계산했으며, 마지막에는 같은 값을 재현하는 Python 코드와 모델 검증 절차를 제공합니다.

학습 결과: 이 글을 마치면 원형축과 비원형축에 적용할 공식을 구분하고, 같은 재료량에서 중공축의 효율을 비교하며, 폐단면과 개단면의 비틀림강성을 계산하고, 자유 뒤틀림과 구속 뒤틀림에 맞는 보·셸·솔리드 모델을 선택해 반력 토크와 회전으로 검증할 수 있습니다.

단계핵심 질문설계 판단

1 · 기초토크는 어떻게 전단응력과 회전을 만드는가?T, τ, θ, GJ의 단위와 역할을 구분합니다.

2 · 원형축왜 τ=Tρ/J가 성립하는가?실축과 중공축을 강도·강성·재료량으로 비교합니다.

3 · 비원형왜 Ip=Ix+Iy를 쓰면 틀리는가?Saint-Venant 뒤틀림과 비틀림상수 J를 구분합니다.

4 · 얇은벽폐단면의 전단흐름은 어떻게 순환하는가?Bredt–Batho 식으로 q, τ, J를 계산합니다.

5 · 전문가구속 뒤틀림은 왜 축응력을 만드는가?warping beam, shell, solid의 적용범위를 판단합니다.

앞선 보의 전단응력 가이드에서 전단중심을 벗어난 횡하중은 T=Ve의 토크를 만든다고 보았습니다. 이번 글은 그 토크가 단면 내부에서 어떻게 전달되고 회전과 응력으로 나타나는지를 이어서 다룹니다.

1단계 — 토크, 전단응력, 회전을 한 식으로 섞지 않는다

토크는 힘이 아니라 힘의 모멘트입니다. 축 방향을 x라 하면 단면의 전단응력 성분 τxy, τxz가 만드는 모멘트의 합이 내부 토크가 됩니다.

T = ∫A (y τxz − z τxy) dA

이 식은 원형, 사각형, 얇은벽, FEA 단면절단 모두에 적용되는 평형식입니다. 반면 τ=Tρ/J는 그 평형을 원형축의 특별한 변형학과 결합한 결과입니다. 먼저 물리량을 구분해야 합니다.

물리량의미N–mm 단위계검산 포인트

T단면이 전달하는 내부 토크N·mm반력 토크의 합과 일치

τ재료점의 전단응력N/mm² = MPa응력 적분의 모멘트가 T

θ두 단면 사이의 상대 회전raddeg로 바꿀 때 180/π

θ′=dθ/dx단위길이당 비틀림각rad/mm균일축이면 θ/L

GJSaint-Venant 비틀림강성N·mm²T/θ′와 일치

비틀림각은 rad로 계산해야 합니다. 2°를 식에 그대로 넣으면 0.0349 rad 대신 2를 사용하게 되어 강성을 57.3배 잘못 평가합니다. 또한 재료의 전단탄성계수 G와 단면의 기하학적 특성 J는 별개입니다. 재료를 바꾸면 G, 단면을 바꾸면 J가 변합니다.

2단계 — 원형축에서 Tρ/J와 TL/GJ를 유도한다

원형축은 비틀린 뒤에도 단면이 평면을 유지하고 반지름선이 직선으로 남는 특별한 형상입니다. 길이 dx 사이의 두 단면이 dθ만큼 회전하면 반지름 ρ의 점이 이동한 원호길이는 ρdθ이고 전단변형률은

γ(ρ) = ρ dθ/dx = ρ θ′

입니다. 선형탄성에서 τ=Gγ이므로 τ=Gρθ′이고, 이를 토크 평형에 넣으면

T = ∫A τρ dA = Gθ′ ∫A ρ²dA = GJθ′

가 됩니다. 원형축에서 J=∫ρ²dA=Ip이므로

τ(ρ)=Tρ/Jθ′=T/(GJ)θ=TL/(GJ) — T, G, J가 길이에 따라 일정할 때

입니다. Illinois의 Mechanics Reference 비틀림 자료는 같은 변형률–응력–평형 순서와 함께 실축 J=πD⁴/32, 중공축 J=π(Do⁴−Di⁴)/32를 정리합니다. 최대응력은 외주 ρ=Do/2에서 생깁니다.

같은 재료량으로 비교해야 중공축의 효율이 보인다

외경이 같은 실축과 중공축을 비교하면 중공축의 J가 작습니다. 이것만 보고 “중공축은 약하다”고 결론 내리면 비교 기준이 잘못된 것입니다. 비틀림에서 재료는 중심보다 바깥에 있을수록 ρ²의 가중치를 받습니다. 따라서 같은 단면적, 즉 같은 길이·재료·질량으로 비교해야 형상 효율을 볼 수 있습니다.

외경 Do=60 mm, 내경 Di=40 mm인 중공축의 면적과 같은 실축 지름은 Ds=√(Do²−Di²)=44.721 mm입니다. 두 축에 T=2.4 kN·m, L=1,200 mm, G=79.3 GPa를 적용합니다.

항목중공축 60/40등면적 실축 44.721해석

단면적1,570.796 mm²1,570.796 mm²재료량 동일

J1.0210×10⁶ mm⁴0.3927×10⁶ mm⁴중공축 2.60배

τmax70.518 MPa136.658 MPa중공축 48.4% 낮음

θ2.038°5.299°중공축 회전 61.5% 감소

이 예제는 탄성 Saint-Venant 비틀림만 비교합니다. 실제 동력축에서는 토크 변동에 의한 피로, 키홈과 스플라인의 응력집중, 임계회전수, 베어링 간 정렬, 얇은 벽의 국부불안정까지 설계기준과 시험으로 확인해야 합니다.

3단계 — 비원형 단면에는 극관성모멘트를 대입하지 않는다

Ip=Ix+Iy는 면적의 기하학적 적분값이고 모든 단면에서 계산할 수 있습니다. 그러나 비틀림상수 J와 같아지는 것은 원형·원환형 단면뿐입니다. 비원형 단면에서는 각 점의 변위가 단순한 ρ θ가 아니고 단면이 축 방향으로 뒤틀립니다. 이 warping 때문에 전단응력은 중심에서의 거리만으로 정해지지 않습니다.

MIT 16.20의 Basic Torsion Theory Handout는 Saint-Venant 문제를 Prandtl 응력함수의 Poisson 방정식과 자유경계조건으로 표현합니다. 즉, 비원형 J는 단순 면적합이 아니라 단면 형상과 경계조건을 푼 결과입니다.

폭 a=120 mm, 두께 b=20 mm인 직사각형의 비틀림상수는 a≥b일 때 다음 근사식으로 계산할 수 있습니다.

J ≈ ab³[1/3 − 0.21(b/a){1 − b⁴/(12a⁴)}]

계산하면 J=0.2864×10⁶ mm⁴입니다. 반면 Ip=ab(a²+b²)/12=2.960×10⁶ mm⁴로 10.335배 큽니다. θ=TL/(GIp)를 사용하면 회전이 실제 Saint-Venant 해의 약 1/10로 계산되어 단면이 10.34배 강한 것처럼 보입니다. 미 공군 응력 매뉴얼을 재구성한 Engineering Library 비틀림 절도 직사각형의 응력과 회전에 별도 형상계수를 사용하며, 얇은 직사각형의 극한에서 J≈bt³/3가 됨을 보여 줍니다.

비원형 응력은 Tρ/J로 구할 수 없다

직사각형의 최대전단응력은 긴 변 중앙에서 생기고 모서리에서는 자유경계조건 때문에 0입니다. 원형 공식이라면 중심에서 먼 모서리가 최대가 되어야 하므로 물리적으로 모순입니다. 직사각형·타원·임의 실단면은 검증된 폐형식 해, 단면 라이브러리, 2D warping 단면해석 또는 3D 모델을 사용해야 합니다.

4단계 — 얇은 폐단면은 전단흐름으로 토크를 전달한다

벽 두께 t가 단면 치수보다 충분히 작고 한 개의 폐셀을 이루면 벽 두께 방향 응력변화를 무시하고 중간선을 따라 전단흐름 q=τt를 정의할 수 있습니다. 중간선이 감싸는 면적을 Am이라 하겠습니다.

폐단면의 작은 벽 조각에 대한 평형은 한 셀을 도는 q가 일정해야 함을 보여 줍니다. 이 흐름이 만드는 토크는 Bredt 식으로

T = 2Am qq = T/(2Am)τ(s)=q/t(s)

입니다. 중요한 점은 두께가 바뀌어도 단일 셀의 q는 같고, 응력 τ만 두께에 반비례한다는 것입니다. 호환조건을 포함하면 Saint-Venant 비틀림상수는

J = 4Am² / ∮[ds/t(s)]

가 됩니다. 이 식은 벽의 중간선을 따라 적분합니다. 외곽면적이나 내측면적을 임의로 섞으면 얇은 벽에서도 오차가 생깁니다. MIT 16.20의 Unit 12: Torsion of Thin Closed Sections와 Handout 3은 q=τt, Bredt 식, J=4A²/∮(ds/t)를 같은 순서로 전개합니다.

두께가 다른 직사각형 박스 예제

중간선 치수 B=120 mm, H=80 mm인 박스에서 위·아래 벽은 2 mm, 좌·우 벽은 4 mm라고 합시다. T=1.0 kN·m이면 Am=9,600 mm²이고

q = 1,000,000/(2×9,600) = 52.083 N/mm

입니다. 따라서 2 mm 벽의 응력은 26.042 MPa, 4 mm 벽은 13.021 MPa입니다. 또한 ∮ds/t=2(120/2)+2(80/4)=160이므로 J=2.304×10⁶ mm⁴입니다.

5단계 — 폐단면과 개단면의 차이는 두께보다 위상에 있다

같은 판재를 사용해도 벽이 닫힌 경로를 만드는지, 슬릿으로 끊어져 있는지에 따라 비틀림강성이 수천 배 달라질 수 있습니다. 얇은 개단면을 길이 bi, 두께 ti인 판 조각의 합으로 보면

Jopen ≈ Σ(bi ti³/3)

입니다. 두께가 세제곱으로 들어가므로 얇은 개단면은 매우 유연합니다. 중간선 120×80 mm, 일정 두께 2 mm인 직사각형 박스를 비교해 보겠습니다.

형상사용식J동일 G, L에서 회전

완전히 닫힌 박스4Am²/∮(ds/t)1,843,200 mm⁴1.0

한 줄 슬릿이 있는 개단면Σ(bti³/3)1,066.667 mm⁴1,728배

T=1.0 kN·m, G=79.3 GPa에서 폐단면의 비틀림률은 0.006842 rad/m, 즉 0.392°/m입니다. 같은 선형 소변형 식을 슬릿 단면에 적용하면 11.82 rad/m가 되어 이미 작은 회전 가정을 심각하게 벗어납니다. 이 결과는 슬릿 단면의 실제 회전을 그대로 예측하려는 값이 아니라, 동일 토크를 지지하는 형상으로 부적합하다는 조기 경고입니다.

다중 셀은 평형만으로 풀리지 않는다

두 개 이상의 폐셀이 있으면 각 셀의 순환 전단흐름 q1, q2, ...가 미지수입니다. 전체 토크 평형 T=2ΣAm,i qi만으로는 식이 부족합니다. 공유 벽에는 인접 셀 흐름의 차이가 작용하고, 모든 셀의 비틀림률이 같다는 호환조건을 추가해야 합니다. 실무 계산서는 다음 세 줄을 분리해서 남기는 편이 안전합니다.

각 셀의 방향을 먼저 정하고 순환 qi의 양의 방향을 표시합니다.

전체 토크 평형과 공유 벽의 실제 흐름 qi−qj를 적습니다.

각 셀에 대해 ∮(q/t)ds = 2AmGθ′ 형태의 동일 비틀림률 조건을 세웁니다.

곡률이 크거나 두께가 단면 치수에 비해 두껍고, 코너 반경·보강재·컷아웃이 중요한 경우에는 Bredt 얇은벽 가정의 오차를 셸 또는 단면해석으로 확인해야 합니다.

6단계 — 자유 뒤틀림과 구속 뒤틀림을 구분한다

Saint-Venant 비틀림은 단면이 축 방향으로 자유롭게 뒤틀릴 수 있다고 봅니다. 균일 토크를 받는 부재의 중간 구간에서는 이 가정이 잘 맞을 수 있습니다. 그러나 I형·채널·앵글 같은 개단면을 두꺼운 엔드플레이트에 용접하거나 다이어프램으로 고정하면 끝단에서 warping 변위가 막힙니다.

뒤틀림이 구속되면 단면의 일부 섬유는 축 방향으로 늘어나고 다른 섬유는 줄어들어 추가 축응력과 bimoment가 생깁니다. 선형 Vlasov 보 이론에서는 흔히

T = GJ θ′ − ECw θ‴B = −ECw θ″

처럼 Saint-Venant 토크와 warping 토크를 분리합니다. 여기서 Cw는 warping constant이며 부호는 좌표와 sectorial-coordinate 규약에 따라 달라질 수 있습니다. 식보다 중요한 것은 J만으로는 구속된 개단면의 축방향 warping 응력을 복원할 수 없다는 점입니다.

Abaqus의 Torsional response—warping 설명은 얇은 개단면이 자유 비틀림에서는 매우 유연하지만 warping이 구속되면 축응력이 발생한다고 구분합니다. 일반 보요소는 자유 뒤틀림을 가정하므로, 개단면의 구속 warping을 모델링하려면 warping 자유도 7을 가진 B31OS/B32OS 계열 또는 충분한 셸/솔리드 모델이 필요합니다.

상황우선 모델반드시 확인할 출력

원형축, 단순 토크 전달일반 beam 또는 solidsection torque, θ, τmax

폐단면, 전체 회전 중심beam + 손계산 BredtGJ, q, 반력 토크

폐단면의 컷아웃·접합부·좌굴shell벽 전단흐름, 국부응력, 변형모드

개단면, warping 자유일반 beam 가능J와 회전, 적용범위

개단면, warping 구속warping beam 또는 shellbimoment, warping 축응력, 경계조건

키홈·필렛·접촉·소성solid국부 메시수렴과 구조적 응력

7단계 — FEA는 색상보다 평형과 회전으로 검증한다

비틀림 해석에서 가장 흔한 실패는 최대 von Mises 응력 한 점만 보고 모델을 통과시키는 것입니다. 비틀림의 1차 검증량은 토크와 비틀림률입니다. 하중이 올바르게 전달되지 않거나 보 로컬축이 틀려도 응력 그림은 그럴듯할 수 있지만, 반력 토크와 T/(GJ)는 즉시 모순을 드러냅니다.

Beam 모델

단면 라이브러리의 J가 Ip인지 Saint-Venant J인지 확인합니다.

요소 로컬축의 토크 성분과 적용한 전역 모멘트 방향을 맞춥니다.

하중점에서 떨어진 균일 구간의 θ′를 두 절점 회전 차이/거리로 구해 T/(GJ)와 비교합니다.

개단면 구속 문제라면 warping 자유도와 bimoment 출력을 지원하는 요소인지 확인합니다.

Shell 모델

중간면 형상과 실제 두께가 Bredt 계산의 Am, t 정의와 일치하는지 확인합니다.

폐단면을 가르는 seam, 중복 절점, 연결 누락이 없는지 확인합니다. 작은 불연속도 개단면으로 바꿀 수 있습니다.

단일 셀의 Saint-Venant 영역에서 q=τt가 둘레를 따라 거의 일정한지, 두께가 다른 곳에서 τ가 역비로 변하는지 봅니다.

코너·용접·컷아웃의 피크는 공칭 Bredt 응력과 분리해 메시수렴과 설계기준의 구조적 응력 방식으로 평가합니다.

Solid 모델

끝면에 한 점 집중모멘트를 직접 걸기보다 분포 coupling 또는 물리적 접촉으로 토크를 전달합니다.

관심 단면에서 ∫(yτxz−zτxy)dA를 계산해 입력 토크를 복원합니다.

고정단 바로 옆은 warping 구속과 Saint-Venant 끝단효과가 섞이므로 균일 구간 검증점과 구분합니다.

응력집중부는 요소 크기를 줄여도 발산하는 특이점인지, 유한한 필렛 응력인지 구분합니다.

메시 수렴표에는 세 값을 함께 남긴다

메시반력 토크 오차구간 비틀림률 오차구조적 응력 변화

coarse평형 오류 확인손계산 T/GJ와 비교기준값 기록

medium0.5% 이내 권장차이 감소 확인coarse 대비 변화율

fine거의 동일목표 오차 이내medium 대비 안정화

허용오차는 해석 목적과 품질계획에서 정해야 하며 위 숫자는 보편 규정값이 아닙니다. 중요한 것은 최대응력 하나가 아니라 평형, 전체 회전, 설계에 사용할 구조적 응력이 각각 안정되는지 기록하는 것입니다.

8단계 — 계산을 Python으로 재현한다

다음 코드는 외부 패키지 없이 이 글의 원형축, 직사각형, 폐단면·개단면 결과를 재현합니다. 단위는 N, mm, MPa입니다.

from math import pi, sqrt, degrees T = 2.4e6 # N·mm L = 1200.0 # mm G = 79_300.0 # MPa = N/mm^2 # Equal-area hollow and solid circular shafts Do, Di = 60.0, 40.0 Ds = sqrt(Do**2 - Di**2) Jh = pi * (Do**4 - Di**4) / 32 Js = pi * Ds**4 / 32 tau_h = T * (Do / 2) / Jh tau_s = T * (Ds / 2) / Js theta_h = T * L / (G * Jh) theta_s = T * L / (G * Js) assert abs(Jh / Js - 2.6) < 1e-12 assert abs(tau_h - 70.517882) < 1e-6 # Noncircular rectangle: torsion constant, not polar area moment a, b = 120.0, 20.0 J_rect = a*b**3 * (1/3 - 0.21*(b/a)*(1 - b**4/(12*a**4))) Ip_rect = a*b*(a**2 + b**2)/12 # Closed box versus the same walls with a slit B, H, t = 120.0, 80.0, 2.0 Am = B*H perimeter = 2*(B + H) J_closed = 4*Am**2/(perimeter/t) J_open = perimeter*t**3/3 print(f"equal-area solid diameter = {Ds:.3f} mm") print(f"hollow J / solid J = {Jh/Js:.3f}") print(f"hollow: tau={tau_h:.3f} MPa, theta={degrees(theta_h):.3f} deg") print(f"solid : tau={tau_s:.3f} MPa, theta={degrees(theta_s):.3f} deg") print(f"rectangle Ip/J = {Ip_rect/J_rect:.3f}") print(f"closed/open J ratio = {J_closed/J_open:.0f}")

예상 출력은 Ds=44.721 mm, Jh/Js=2.600, 중공축 τ=70.518 MPa와 θ=2.038°, 실축 τ=136.658 MPa와 θ=5.299°, 직사각형 Ip/J=10.335, 폐단면/개단면 J 비 1,728입니다.

9단계 — 실무 의사결정 순서

하중경로: 외부 토크의 발생원, 전달부, 반력 위치를 자유물체도로 표시합니다.

설계기준: 허용응력뿐 아니라 허용 비틀림각, 피로수명, 기능 정렬 요구를 정합니다.

단면 분류: 원형/비원형, 실단면/얇은벽, 폐단면/개단면, 단일셀/다중셀을 구분합니다.

경계조건: warping이 자유로운지, 엔드플레이트·다이어프램·접합부가 구속하는지 표시합니다.

1차 계산: 원형축, Bredt, 개단면 Σbt³/3 중 맞는 모델로 J, 응력, 회전을 계산합니다.

상세모델: 컷아웃, 용접, 키홈, 국부좌굴, 구속 warping이 지배할 때만 shell/solid로 확장합니다.

검증: 반력 토크 → 단면 토크 → 비틀림률 → 응력장 → 메시수렴 순서로 닫습니다.

10단계 — 제출 전 체크리스트

T를 N·m와 N·mm 사이에서 1,000배 잘못 변환하지 않았는가?

비원형 단면에 Ip를 J처럼 사용하지 않았는가?

중공축 효율을 비교할 때 같은 외경인지 같은 재료량인지 기준을 명시했는가?

얇은 폐단면의 면적과 적분 경로를 벽 중간선으로 정의했는가?

단일 셀의 q 연속성과 두께별 τ=q/t를 확인했는가?

다중 셀에 토크 평형뿐 아니라 동일 비틀림률 호환조건을 넣었는가?

개단면의 warping 자유/구속을 지지 상세와 일치시켰는가?

FEA 반력 토크와 단면 응력 적분이 입력 토크를 복원하는가?

고정단 피크와 균일 Saint-Venant 영역의 공칭응력을 분리했는가?

강도뿐 아니라 서비스 비틀림각, 피로, 좌굴과 제작 불연속을 검토했는가?

참고자료와 적용범위

이 글의 기본식과 얇은벽 유도는 University of Illinois Torsion, MIT OpenCourseWare 16.20 Handout 3 및 Unit 12, 미 공군 응력 매뉴얼 기반 Beam Torsion, Dassault Systèmes Abaqus 문서의 warping 설명과 beam element 선택 지침을 대조해 작성했습니다. 자료 확인일은 2026-08-17입니다.

수록 식은 균질 선형탄성, 작은 변형, 일정하거나 명시된 단면, Saint-Venant 영역 또는 얇은벽 가정을 전제로 한 교육·1차 검증용입니다. 실제 설계의 허용값, 피로, 소성, 좌굴, 용접·볼트 상세, 복합재 적층, 동적 토크와 안전계수는 해당 산업의 최신 설계기준과 검증된 해석절차를 적용해야 합니다.