3. 柔性铰链的刚度计算与仿真

任务定义

设计一种典型的单轴缺口型柔性铰链（圆弧缺口型），给定几何参数与材料属性，推导其绕弱轴的转动刚度理论公式，并通过有限元静力仿真验证理论计算的正确性。本案例以平面内弯曲（绕Z轴）为主要变形模式，忽略轴向与剪切变形影响。

几何与材料参数

理论刚度公式

对于圆弧缺口型柔性铰链，绕Z轴的转动刚度经典计算公式（Paros-Weisbord简化式）为：

K_θ = (E·b·t²) / (12·R) · f(s)

其中无量纲系数 f(s) 为：

f(s) = 1 / [ (2s + s²) / ((1+s)(2s+s²)^0.5) · arctan(((2s+s²)^0.5)/(1-s)) - 1 ]

式中 s = t/(2R)。代入本案例参数：s = 1/(2×5) = 0.1，计算得 f(0.1) ≈ 0.754。

因此理论转动刚度：

K_θ = (200000 N/mm² × 10 mm × 1² mm²) / (12 × 5 mm) × 0.754
    = 200000 × 10 / 60 × 0.754
    ≈ 25133 N·mm/rad

有限元仿真验证

采用ANSYS APDL建立参数化模型，使用平面应力单元（PLANE182）模拟铰链截面，施加单位弯矩并测量转角，从而反算转动刚度。

! 柔性铰链刚度仿真 - 圆弧缺口型
/PREP7
ET,1,PLANE182
KEYOPT,1,3,0    ! 平面应力

! 材料属性
MP,EX,1,200000
MP,PRXY,1,0.3

! 几何参数
R=5.0
T=1.0
B=10.0
L=20.0        ! 铰链总长（两侧各加10mm刚性段）

! 创建截面轮廓
K,1,-L/2,0
K,2,-R,0
K,3,-R,T/2
K,4,0,T/2+R
K,5,R,T/2
K,6,R,0
K,7,R,-T/2
K,8,0,-T/2-R
K,9,-R,-T/2
K,10,-L/2,-T/2

! 圆弧缺口
LARC,3,4,2,R
LARC,4,5,6,R
LARC,7,8,9,R
LARC,8,9,2,R

! 直线段
L,1,2
L,2,3
L,5,6
L,6,7
L,9,10
L,10,1

! 生成面
AL,ALL

! 网格划分
ESIZE,0.3
AMESH,ALL

! 约束：左端面全约束
NSEL,S,LOC,X,-L/2
D,ALL,ALL,0
ALLSEL

! 在右端面施加单位弯矩（等效为线性分布轴向力）
NSEL,S,LOC,X,L/2
*GET,NCOUNT,NODE,,COUNT
FSUM_LOAD = 1.0    ! 单位弯矩 N·mm

! 计算每个节点的力：F = M·y / I_equiv
! 对矩形截面 I = b·h³/12，此处h为截面高度（含缺口影响，近似用等效高度）
! 简化：在右端面节点施加绕Z轴的力偶
! 取端面中心为原点，施加 ±F 于上下边缘节点

NSEL,R,LOC,Y,GT,0
*GET,NTOP,NODE,,COUNT
F,ALL,FX,0.5/(NTOP*5)   ! 上表面拉力
ALLSEL

NSEL,S,LOC,X,L/2
NSEL,R,LOC,Y,LT,0
*GET,NBOT,NODE,,COUNT
F,ALL,FX,-0.5/(NBOT*5)  ! 下表面压力
ALLSEL

! 求解
/SOLU
SOLVE
FINISH

! 后处理：测量右端面转角
/POST1
NSEL,S,LOC,X,L/2
*GET,UY_TOP,NODE,0,UY
*GET,UY_BOT,NODE,0,UY
! 转角 = (UY_TOP - UY_BOT) / 截面高度
! 截面高度取铰链厚度方向总高 = T + 2R = 11mm
THETA = (UY_TOP - UY_BOT) / 11
K_SIM = 1.0 / THETA
*STATUS,THETA
*STATUS,K_SIM

验证方法

运行上述APDL脚本后，查看输出变量：

• THETA：右端面在单位弯矩作用下的转角（rad）
• K_SIM：仿真得到的转动刚度（N·mm/rad）

判定标准：仿真刚度值与理论值25133 N·mm/rad的相对误差在±10%以内即认为理论公式与仿真互相验证通过。网格密度（ESIZE）可加密至0.2mm进行收敛性检查。

结果解读

理论公式与仿真结果的偏差主要来源于：Paros公式的简化假设（忽略剪切变形、假设圆弧轮廓精确）、有限元网格离散误差以及端部加载方式的近似处理。若误差超出10%，应优先细化网格并检查加载力偶的等效方式。