Hashin破壊基準
理论与物理
Hashin准则是什么
老师，Hashin准则和Tsai-Wu有什么区别？
最大的区别是区分破坏模式。Tsai-Wu是“破坏/未破坏”的二值判定，而Hashin准则（1980）能识别4种独立的破坏模式。
4种破坏模式
Hashin准则的4种模式：
1. 纤维拉伸破坏（Fiber Tension）
2. 纤维压缩破坏（Fiber Compression）
3. 基体拉伸破坏（Matrix Tension）
4. 基体压缩破坏（Matrix Compression）
区分纤维和基体的破坏，并且拉伸和压缩也分开。用4个公式来判定4种模式呢。
是的。这就是Hashin准则的最大优点。知道了破坏模式，就能明确应该提高哪种强度。如果是纤维断裂就增加纤维量，如果是基体开裂就更换基体树脂…对策变得具体。
对应渐进损伤
Hashin准则能用于渐进损伤分析吗？
可以。当各模式的破坏指标超过1时，降低该模式对应的刚度：
- 纤维破坏 → 降低 $E_1$
- 基体开裂 → 降低 $E_2, G_{12}$
- 两者组合 → 降低全部刚度
Abaqus的内置Hashin损伤，将破坏判定+刚度降低+能量耗散一体化实现。是复合材料渐进损伤分析事实上的标准。
Tsai-Wu没有这个功能呢。
Tsai-Wu用一个公式混合了所有模式，因此无法决定“应该降低哪种刚度”。能够进行模式分离的Hashin准则是渐进损伤的基础。
Hashin准则的局限
Hashin准则也有局限：
| 局限 | 说明 |
|---|---|
| 纤维压缩模式过于简单 | 实际由扭结带（纤维局部屈曲）主导 |
| 基体压缩不考虑断裂面角度 | Puck准则或LaRC准则更准确 |
| 无法处理层间剥离 | 需要单独的CZM（内聚力模型） |
Hashin也还有改进的余地呢。
Hashin准则是“模式分离的先驱”，但各模式的物理描述仍较粗糙。LaRC05准则（NASA兰利研究所，2005年）是改善了Hashin局限的最新准则。
总结
我来整理一下Hashin准则。
要点：
- 区分4种破坏模式 — 纤维拉伸/压缩、基体拉伸/压缩
- 对应渐进损伤 — 可实现按模式降低刚度
- Abaqus标准内置 — 内置Hashin损伤
- 比Tsai-Wu物理上更准确 — 但计算更复杂
- 也有局限 — 纤维压缩的扭结带、基体压缩的断裂面角度
Tsai-Wu用于筛选，Hashin用于详细评估，这样的分工明确了。
没错。在设计初期用Tsai-Wu进行概算，在详细设计时用Hashin（或Puck/LaRC）进行精密评估。分阶段的方法是实用的。
Hashin发现CFRP破坏的双轴相互作用
Hashin破坏准则由Zvi Hashin（以色列理工学院）在1973年至1980年间开发。由于传统的最大应力准则在多轴应力场下的复合材料破坏不准确，因此分别定义了纤维断裂、基体断裂、纤维压缩、基体压缩这4种模式，每种模式都有独立的相互作用公式。之后在WWFE（世界范围失效练习，2002-2004年）中，它被评价为CFRP设计中精度最高的准则。
各项的物理意义
- 惯性项（质量项）：$\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统一 |
数值解法与实现
Abaqus的Hashin损伤模型
请告诉我Abaqus的Hashin损伤模型的设置。
Abaqus的内置Hashin损伤分两个阶段设置。
Step 1: 损伤起始准则（Damage Initiation）
```
*DAMAGE INITIATION, CRITERION=HASHIN
X_t, X_c, Y_t, Y_c, S_L, S_T, alpha
```
7个参数：纤维拉伸/压缩强度、基体拉伸/压缩强度、纵向/横向剪切强度、$\alpha$（剪切应力贡献系数）。
Step 2: 损伤演化法则（Damage Evolution）
```
*DAMAGE EVOLUTION, TYPE=ENERGY
G_ft, G_fc, G_mt, G_mc
```
4个断裂能：纤维拉伸/压缩、基体拉伸/压缩的断裂能（$G_c$）。
需要断裂能呢。只有强度值不够吗？
不同求解器的Hashin实现
| 功能 | Abaqus | Nastran | Ansys |
|---|---|---|---|
| Hashin破坏判定 | ○（标准） | △（USDFLD） | ○（ACP Post） |
| 渐进损伤 | ○（DAMAGE EVOLUTION） | △（SOL 400 + 用户） | △（ACP + APDL） |
| 断裂能法 | ○ | — | 有限 |
| 单元删除 | ○（STATUS） | ○（PARAM,ERODEL） | ○（EKILL） |
Abaqus的功能压倒性地丰富呢。
基于Hashin的渐进损伤分析，Abaqus是事实上的标准。论文中说“Hashin damage”几乎都指Abaqus的实现。在Nastran或Ansys中，通过用户子程序或脚本也能实现同等功能，但比较费事。
网格依赖性与正则化
请讲解一下基于断裂能的正则化。
当损伤局部化（集中在一个单元）时，能量耗散会依赖于单元尺寸。使用特征长度（characteristic length）通过单元尺寸对断裂能进行归一化：
$L_c$ 是单元的特征长度（大致为单元尺寸）。这样，即使改变网格，耗散能量也能保持恒定。
单元尺寸越小，破坏应变越大；越大则越小。能量得以守恒。
理解得很完美。但是，如果单元极端大（特征长度 > $2G_c/\sigma_c^2 \cdot E$），会发生回弹，导致数值不稳定。需要注意网格尺寸的上限。
总结
我来整理一下Hashin准则的数值方法。
要点：
相关主题
用語集Hashin破壊基準 — CAE用語解説结构分析Tsai-Wu破壊基準结构分析Mohr-Coulomb破壊基準用語集Tsai-Wu破壊基準 — CAE用語解説结构分析Tsai-Wu破壊基準用語集Tsai-Hill破壊基準 — CAE用語解説この記事の評価ご回答ありがとうございます！参考に
なったもっと
詳しく誤りを
報告