库仑摩擦模型
库仑摩擦的理论基础
库仑摩擦是什么
老师，请教我有限元中摩擦模型的基础知识。
库仑摩擦是最基本的摩擦模型。切向摩擦力与法向力成正比：
- $\tau$ — 摩擦力（切向应力）
- $p_n$ — 接触压（法向应力）
- $\mu$ — 摩擦系数
固着与滑动
两种状态：
- 固着（stick） — $|\tau| < \mu p_n$。无相对滑动
- 滑动（slip） — $|\tau| = \mu p_n$。摩擦力达到极限而滑动
固着到滑动的转换是非线性的。
摩擦系数的代表值
| 材料配对 | $\mu$（干燥） |
|---|---|
| 钢-钢 | 0.15～0.3 |
| 钢-铝 | 0.2～0.4 |
| 钢-橡胶 | 0.5～0.8 |
| 钢-聚四氟乙烯 | 0.04～0.1 |
| 混凝土-钢 | 0.3～0.5 |
有限元中的摩擦设置
总结
要点：
- $|\tau| \leq \mu p_n$ — 摩擦力的上限
- 固着/滑动的二元状态 — 接触的附加非线性
- 摩擦系数依赖材料配对 — 使用实测值
- 用罚函数法处理摩擦 — 滑动方向也用罚函数弹簧
库仑的1781年实验
法国军事工程师Charles-Augustin de Coulomb在1781年进行了数百次摩擦实验，系统地研究了木材、金属、石头等材料的组合，推导出"库仑摩擦定律"：摩擦力与垂直荷载成正比，与接触面积无关。更有趣的是，追溯更早的历史发现Leonardo da Vinci在1495年的秘密笔记中就已经图示过这个关系，1967年才为人所知，令科学史震惊。
库仑摩擦的数值计算方法
摩擦的罚函数法
摩擦也用罚函数法处理。切向方向的"固着罚函数"：
$k_t$ 是切向罚函数刚度，$\delta_t$ 是切向弹性滑动。当$|\tau| = \mu p_n$时，过渡到滑动。
"弹性滑动"是指在实际滑动前有微小滑动吗？
这是罚函数法的人工"弹性滑动"。物理上，固着状态滑动为零，但数值上有限的切向罚函数刚度导致微小滑动。$k_t$越大，弹性滑动越小。
静摩擦与动摩擦
LS-DYNA可分别设置静摩擦系数 $\mu_s$和动摩擦系数 $\mu_d$。滑动开始时$\mu_s \to \mu_d$（$\mu_d < \mu_s$），用指数函数平滑过渡。
总结
固着-滑动判断
库仑摩擦的有限元实现需要在每个接触点判断"固着"或"滑动"。1970年代Zienkiewicz等人引入的return mapping方法采用两步法：先计算试应力，再向摩擦圆锥投影，大幅降低每迭代步的计算成本。该方法至今仍是ABAQUS和ANSYS的标准算法。
库仑摩擦的实际应用
摩擦的实务
摩擦关键应用：
- 螺栓接合的夹紧力 — 摩擦传递荷载
- 冲压成形 — 模具与毛坯的摩擦影响变形
- 制动器 — 摩擦力 = 制动力
- 管道支架 — 滑动支架的摩擦力
摩擦系数的不确定性
摩擦系数变异性很大（±30%以上）。高度依赖表面状态（粗糙度、润滑、氧化膜）。
应对方法：
- 摩擦系数的敏感性分析（$\mu_{低}$ 和 $\mu_{高}$ 两种工况）
- 使用实测值（试验结果）
- 文献值仅作参考
实务检查清单
制动器尖叫分析
汽车制动器尖叫（brake squeal）由库仑摩擦系数对速度的细微依赖性引起的负阻尼。2000年代福特和TRW共同用Nastran进行的复特征值分析表明：仅将摩擦系数从0.35改为0.40，不稳定模式数就增加3倍，直接指导垫片形状优化。
库仑摩擦的软件对比
摩擦工具
所有求解器都标准支持库仑摩擦。差异出现在高级摩擦模型（磨损、速度相关）。
| 功能 | Abaqus | Nastran | Ansys | LS-DYNA |
|---|---|---|---|---|
| 库仑摩擦 | ○ | ○ | ○ | ○ |
| 静/动摩擦 | △ | — | ○ | ○（FS/FD） |
| 速度相关摩擦 | *FRICTION, EXPONENTIAL | — | — | ○ |
| 磨损模型 | 用户子程序 | — | — | ○（*MAT_WEAR） |
选型指南
摩擦圆锥实现的演变
MSC Nastran在1970年代后期通过SLID/STIFF参数引入简化摩擦模型，但完整库仑摩擦圆锥直到2006年MD Nastran才实现。相比之下，LS-DYNA早在1993年左右的版本930就在type 3接触中允许μs和μd的两级设置，多年来成为碰撞分析的实际标准。
库仑摩擦的先进研究
摩擦的前沿研究
分子动力学摩擦定律的前沿
2010年代以来，分子动力学（MD）仿真揭示了库仑定律不适用的纳米尺度摩擦机制。2019年Nature Materials发表的MoS₂单层实验观测到摩擦力与接触面积的0.5次方成正比的标度关系，开启了连接宏观库仑模型与微观的多尺度摩擦定律研究的新时代。
库仑摩擦的故障排除
摩擦故障排除
摩擦振荡发散失败
1990年代早期碰撞分析中经常使用静态试验的摩擦系数（μ=0.6）直接用于动力分析，导致滑动面产生高频振荡，计算发散。根本原因是缺少速度相关项，改用动摩擦系数（μk≈0.4），或用Stribeck曲线平滑插值，问题迎刃而解。
价值
详细
错误