复合材料机械连接
理论与物理
复合材料螺栓连接
老师，用螺栓连接复合材料与金属有什么不同？
根本不同。金属在螺栓孔周围会发生屈服并重新分配应力，但复合材料会发生脆性破坏。由于无法期待通过屈服进行载荷再分配，螺栓孔的应力集中会直接导致破坏。
复合材料螺栓连接的破坏模式
四种主要破坏模式：
| 模式 | 特征 | 危险度 |
|---|---|---|
| 挤压破坏 | 螺栓孔周围压溃 | 理想的（渐进式） |
| 净截面拉伸破坏 | 螺栓孔截面处拉伸断裂 | 危险的（突发式） |
| 剪切破坏 | 从螺栓孔到边缘发生剪切开裂 | 危险的（突发式） |
| 劈裂破坏 | 从螺栓孔纵向开裂 | 危险的（突发式） |
挤压破坏是“理想的”？
挤压破坏是孔周围逐渐压溃的过程，因此不会导致突发性崩溃。设计时应确定尺寸使挤压破坏优先发生。净截面拉伸和剪切破坏是突发性的，应予以避免。
设计参数
连接部的尺寸参数：
- $e/d$ — 端距/螺栓直径比。$e/d \geq 3$ 可避免剪切破坏
- $w/d$ — 板宽/螺栓直径比。$w/d \geq 5$ 可避免净截面拉伸破坏
- $t/d$ — 板厚/螺栓直径比。$t/d \leq 1$ 以确保挤压强度
$e/d \geq 3$ 的规则和金属一样吗？
金属 $e/d \geq 2$ 就足够，但复合材料因其脆性需要$e/d \geq 3$。这也取决于铺层结构，如果不是像 $[0/\pm45/90]$ 这样平衡良好的铺层，则需要更大的 $e/d$。
FEM中的建模
级别1（简易）: 用弹簧单元表示螺栓
级别2（中等）: 梁单元＋接触面
级别3（详细）: 用实体单元表示螺栓＋孔＋预紧力＋接触＋PDA
级别3非常复杂啊。
用级别3的FEM模拟挤压破坏的过程（孔的压溃→基体开裂→纤维屈曲→最终破坏）。常用Abaqus的Hashin准则和CZM（内聚力模型）进行。
总结
我来整理一下复合材料螺栓连接的理论。
要点：
- 复合材料是脆性破坏 — 没有屈服带来的载荷再分配
- 四种破坏模式 — 挤压（渐进式）、净截面拉伸/剪切/劈裂（突发式）
- 设计使挤压破坏优先发生 — $e/d \geq 3, w/d \geq 5$
- 用FEM级别3再现挤压破坏 — Hashin + CZM + 接触
- 铺层结构对连接强度影响很大 — 平衡良好的铺层是必须的
复合材料连接部强度的非均匀分布
在复合材料螺栓连接部，端部螺栓会集中承受载荷，产生“载荷不均匀分布（Load sharing imbalance）”。根据Hart-Smith的分析（1980年代），在3螺栓的剪切接头中，端部螺栓承担总载荷的45〜55%。这种不均匀性随层压板的面内弹性模量和连接部的几何形状而变化，通过FEM计算详细的载荷分配可以确定最优配置。
各项的物理意义
- 惯性项（质量项）：$\rho \ddot{u}$，即“质量×加速度”。您是否有过急刹车时身体被向前甩出的经历？那种“被带走的感觉”正是惯性力。物体越重越难启动，一旦启动也越难停止。建筑物在地震中摇晃，也是因为地面突然移动，而建筑物的质量“被落下”。静力分析中此项设为零，这是基于“缓慢施加载荷，加速度可忽略”的假设。对于冲击载荷或振动问题，此项绝对不能省略。
- 刚度项（弹性恢复力）：$Ku$ 或 $\nabla \cdot \sigma$。拉弹簧时会感觉到“想要恢复的力”吧？那就是胡克定律 $F=kx$，也是刚度项的本质。那么提问——铁棒和橡皮筋，用相同的力拉伸，哪个伸得更长？当然是橡皮筋。这种“难以拉伸的程度”就是杨氏模量 $E$，它决定了刚度。常见的误解：“刚度高=强度高”是不对的。刚度是“不易变形的程度”，强度是“不易破坏的程度”，是不同的概念。
- 外力项（载荷项）：体积力 $f_b$（如重力）和表面力 $f_s$（如压力、接触力）。可以这样理解——桥上卡车的重量是“作用在整个内容物上的力”（体积力），轮胎压路面的力是“仅作用在表面的力”（表面力）。风压、水压、螺栓的紧固力…这些都是外力。这里容易犯的错误：弄错载荷方向。本想施加“拉伸”却变成了“压缩”——听起来像笑话，但在3D空间中坐标系发生旋转时确实会发生。
- 阻尼项：瑞利阻尼 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。试着弹一下吉他的弦。声音会一直持续吗？不，会逐渐变小。这是因为振动能量通过空气阻力和弦的内部摩擦转化成了热能。汽车的减震器也是同样原理——特意吸收振动能量以提高乘坐舒适性。如果阻尼为零会怎样？建筑物在地震后会一直摇晃不停。实际上不会这样，因此设置适当的阻尼很重要。
假设条件与适用范围
量纲分析与单位制
| 变量 | SI单位 | 注意点・换算备忘 |
|---|---|---|
| 位移 $u$ | m（米） | 输入mm时，载荷・弹性模量也需统一为MPa/N系 |
| 应力 $\sigma$ | Pa（帕斯卡）= N/m² | MPa = 10⁶ Pa。与屈服应力比较时注意单位制不一致 |
| 应变 $\varepsilon$ | 无量纲（m/m） | 注意工程应变与对数应变的区别（大变形时） |
| 弹性模量 $E$ | Pa | 钢: 约210 GPa，铝: 约70 GPa。注意温度依赖性 |
| 密度 $\rho$ | kg/m³ | mm系中为tonne/mm³（钢为 = 10⁻⁹ tonne/mm³） |
| 力 $F$ | N（牛顿） | mm系用N，m系也用N统一 |
数值解法与实现
螺栓连接的FEM详细模型
请讲解一下复合材料螺栓连接的详细FEM模型。
模型构成
- 板: 实体单元（C3D8I 或 SC8R）。铺层用层定义
- 螺栓: 实体单元 或 刚体（根据分析目的）
- 螺栓-孔间接触: 带摩擦的接触。也考虑间隙
- 板-板间接触: 结合面的摩擦。预紧力产生的夹紧
- 损伤模型: Hashin（面内损伤）+ CZM（层间剥离）
还要考虑间隙吗？
螺栓孔的间隙（相对于螺栓直径的孔间隙）会影响挤压载荷的分布。精密加工孔（过盈配合）和标准孔（0.1〜0.3 mm 间隙）的结果会不同。
网格要求
求解器设置
Abaqus设置示例（挤压破坏模拟）:
```
*STEP, NLGEOM=YES, INC=1000
*STATIC
0.01, 1.0, 1e-10, 0.02
*CONTACT
...
*DAMAGE INITIATION, CRITERION=HASHIN
...
*DAMAGE EVOLUTION, TYPE=ENERGY
...
```
单个螺栓连接部的完整模型可达数十万自由度。计算时间数小时至数天。多螺栓连接则规模更大。
总结
我来整理一下复合材料螺栓连接的数值方法。
要点：
- 实体单元＋接触＋PDA — 详细模型的标准构成
- 间隙的建模 — 影响挤压载荷分布
- 孔周围网格为0.5〜1 mm — 捕捉损伤的局部化
- 计算成本高 — 单个螺栓需数小时
- Abaqus Hashin + CZM是标准 — 再现挤压破坏
CFRP螺栓孔周边应力集中的FEM评估
复合材料螺栓连接部的强度评估使用点应力准则（PSC）或平均应力准则（ASC）。PSC将孔端距离材料特性距离D₀处的应力与材料强度进行比较，预测破坏。D₀对于CFRP T300/5208的标准值约为1.5〜3.5mm。通过实验标定此参数并输入FEM，可以实现螺栓孔强度预测精度在±10%以内。
线性单元（1次单元）
节点间线性插值。计算成本低，但应力精度低。注意剪切自锁（可通过减缩积分或B-bar法缓解）。
2次单元（带中间节点）
可表现曲线变形。应力精度大幅提高，但自由度约增加2〜3倍。推荐：应力评估重要时使用。
完全积分 vs 减缩积分
完全积分：有过约束（自锁）风险。减缩积分：有沙漏模式（零能量模式）风险。根据情况选择。
自适应网格
基于误差指标（如ZZ估计量）的自动细化。有效提高应力集中部位的精度。有h法（单元细分）和p法（增加阶次）。
牛顿-拉弗森法
非线性分析的标准方法。每次迭代更新切线刚度矩阵。在收敛半径内具有二次收敛性，但计算成本高。
修正牛顿-拉弗森法
切线刚度矩阵使用初始值或每隔数次迭代更新。每次迭代成本低，但收敛速度为线性。
相关主题
なった
詳しく
報告