Hashin破坏准则
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的Built-in Hashin Damage将破坏判定、刚度降低、能量耗散整合为一体。是复合材渐进损伤分析的业界标准。
Tsai-Wu没有这个功能吧。
Tsai-Wu把所有模式混在一个公式中，所以「应该降低哪个刚度」无法确定。能进行模式分离的Hashin准则是渐进损伤分析的基础。
Hashin准则的局限
Hashin准则也有局限：
| 局限 | 说明 |
|---|---|
| 纤维压缩模式过于简化 | 实际上由纤维局部座屈（扭结带）主导 |
| 基体压缩未考虑破坏面角度 | Puck准则或LaRC准则更精确 |
| 无法处理层间剥离 | 需要另外使用CZM（内聚区模型） |
Hashin准则也还有改进空间。
Hashin准则是「模式分离的先驱」，但各模式的物理描述还很粗糙。LaRC05准则（NASA兰利研究所，2005年）改进了Hashin的局限，是最新的准则。
总结
总结Hashin准则。
要点：
- 区分4个破坏模式 — 纤维拉伸/压缩、基体拉伸/压缩
- 对应渐进损伤 — 可进行模式别的刚度降低
- Abaqus标准实现 — Built-in Hashin Damage
- 比Tsai-Wu物理准确 — 但计算更复杂
- 也有局限 — 纤维压缩的扭结带、基体压缩的破坏面角度
明确了Tsai-Wu做初筛，Hashin做详细评估的使用区别。
正是。设计初期用Tsai-Wu做粗算，详细设计用Hashin（或Puck/LaRC）做精密评估。分段式方法在实务上最有效。
Hashin发现了CFRP破坏的二轴相互作用
Hashin破坏准则由Zvi Hashin（以色列理工学院）于1973～1980年开发。传统的最大应力准则在多轴应力场下CFRP破坏预测不准确，所以他独立定义了纤维断裂、基体断裂、纤维压缩、基体压缩4种模式，各自用相互作用式。之后WWFE（世界范围破坏习题，2002～2004年）评估中，该准则被评为CFRP设计中精度最高的准则。
Hashin破坏准则的数值计算方法
Abaqus的Hashin损伤模型
请讲一下Abaqus的Hashin损伤模型设置。
Abaqus的Built-in Hashin Damage分两个阶段设置。
第1步：损伤初始准则（Damage Initiation）
```
*DAMAGE INITIATION, CRITERION=HASHIN
X_t, X_c, Y_t, Y_c, S_L, S_T, alpha
```
7个参数：纤维拉伸/压缩强度、基体拉伸/压缩强度、纵/横剪切强度、$\alpha$（剪切应力贡献系数）。
第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可以通过用户子程序或脚本做到类似的功能，但工作量很大。
网格依赖性与正则化
请讲讲破坏能量的正则化。
损伤局部化（集中在一个单元）时，能量耗散会依赖单元大小。用特征长度来对破坏能量进行单元大小的正规化：
$L_c$ 是单元的特征长度（约为单元大小）。这样无论网格如何变化，耗散能量都保持一致。
单元越小，破坏应变越大；单元越大，应变越小。能量得以保存。
完全理解。但单元太大（特征长度 > $2G_c/\sigma_c^2 \cdot E$）会产生snapback，数值不稳定。注意网格大小上限。
总结
整理Hashin准则的数值方法。
要点：
- Abaqus的*DAMAGE INITIATION + *DAMAGE EVOLUTION是标准 — Hashin损伤的业界事实标准
- 7个材料强度＋4个破坏能量 — 共11个参数
- 破坏能量的正则化 — 消除网格依赖性
- 注意单元大小上限 — 太大易产生snapback
- Nastran和Ansys需要用户子程序 — Abaqus最便利
Hashin准则的4模式计算步骤
Hashin准则计算4个模式：①纤维拉伸(σ11/XT)²+(τ12/Xt)²≥1出现纤维断裂②纤维压缩(σ11/XC)²≥1出现kinking③基体拉伸(σ22/YT)²+(τ12/S12)²≥1出现横向裂纹④基体压缩计算，评估全部4种，首次指标超过1的模式就是初期破损模式。CFRP代表值为XT=2000MPa·XC=1500MPa·YT=50MPa·S12=80MPa。
Hashin破坏准则的实务应用
Hashin准则的实务应用
Hashin准则在实务中怎么用？
主要用于渐进损伤分析（PDA：Progressive Damage Analysis）的基础。
PDA的应用场景
| 场景 | 目的 |
|---|---|
| 飞机的CAI（冲击后压缩） | 预测冲击后的压缩强度 |
| OHC/OHT（开孔压缩/拉伸） | 孔板的拉伸/压缩强度 |
| 缺口效应 | 凹槽处的破坏扩展 |
| 复合材螺栓接合 | 承载破坏的仿真 |
CAI（冲击后压缩）是飞机认证最重要的试验。
正是。复合材结构的设计容限常由CAI强度决定。Hashin＋CZM（层间剥离）的PDA来预测CAI强度，减少试验数。「积木块方法」是航空航天的标准工作流。
分析工作流程
2. 损伤状态转移 — 冲击后的损伤变量转移到Abaqus/Standard
3. 压缩分析（隐式） — 在含损伤的状态下施加压缩荷载，评估残余压缩强度
在显式和隐式间切换。
冲击是短时间动态现象用显式，压缩是准静态用隐式。Abaqus的*IMPORT可进行显式→隐式的状态转移。
获取材料参数
破坏能量 $G_c$ 怎样测量？
| 参数 | 试验方法 | 标准 |
|---|---|---|
| $G_{Ic}$（模式I层间） | DCB（双悬臂梁） | ASTM D5528 |
| $G_{IIc}$（模式II层间） | ENF（端部缺口弯曲） | ASTM D7905 |
| $G_{ft}$（纤维拉伸） | CT试验或推定 | 直接测量困难 |
| $G_{mc}$（基体压缩） | 紧凑压缩 | 研究阶段 |
纤维破坏能量无法直接测？
纤维破坏的能量非常大（$G_{ft} \approx 50 \sim 100$ kJ/m²），没有标准试验法。多用文献值或反解法（用FEM结果与试验对照）估算。
实务检查清单
给个Hashin准则的检查清单。
破坏能量设置最重要，也最难。
强度值可通过标准试验比较容易获得，但破坏能量试验困难或成本高。需对参数不确定性做敏感性分析。
CFRP翼面板破坏包络线设计
飞机CFRP翼面多轴荷载下的安全余度用「破坏包络线」管理。用FEM+Hashin准则在σ11-σ22荷载空间上描绘安全区，确保全部飞行荷载情况都在包络线内。Bombardier C-Series（现A220）的CFRP翼皮用该设计实现轻量化，在全荷载情况下都保证了ALCF（容限设计荷载比）0.85以下。
Hashin破坏准则的软件比较
Hashin准则工具
比较下能用Hashin准则的工具。
Abaqus和LS-DYNA很强。
LS-DYNA的MAT58是Chang-Chang准则（Hashin的变形版）为基础的渐进损伤模型。碰撞时复合材破坏广泛应用。
选择指南
Hashin准则Abaqus完成度最高。
Abaqus的Hashin Damage实现是论文最多引用的复合材损伤模型。用其他求解器做一样的分析需要很大的开发工作量。
Abaqus HashinDamage材料模型
Abaqus的*DAMAGE INITIATION, CRITERION=HASHIN能自动用4种模式评估Hashin准则，初期破坏后进行材料劣化（刚度降低）的渐进分析。DAMAGE EVOLUTION卡设定能量基的刚度低减率。GKN Aerospace用Abaqus+Hashin做CFRP机体部件认证分析，Boeing 787翼形成型面板的静破坏试验精度预测达±8%。
Hashin破坏准则的前沿研究
超越Hashin准则的破坏理论
有改进Hashin准则局限的研究吗？
Hashin准则的各模式在物理上更准确的次世代准则正被开发。
LaRC05准则（NASA兰利）
NASA兰利研究所开发的LaRC05（Langley Research Center Criteria，2005）用物理基础改进了Hashin的4模式：
- 纤维压缩 — 扭结带理论（纤维微观座屈）。比Hashin的简单压缩准则更准确
- 基体压缩 — 基于Mohr-Coulomb的破坏面理论。与Puck准则同等精度
- 纤维拉伸 — 与Hashin同样
- 基体拉伸 — 与Hashin同样但改进了横剪的处理
扭结带是什么？
纤维受压时，纤维初始缺陷（微小的错排）被放大，产生局部座屈。这就是扭结带。压缩强度不是纤维本身的强度，而是纤维的稳定性（初始缺陷和树脂剪切刚度）决定的。
Hashin的 $\sigma_1 / X_c$ 表现不了这种物理。
扭结带的压缩强度为：$X_c \approx G_{12} / (1 + \phi_0 / \gamma_c)$。树脂剪切刚度 $G_{12}$ 和初始纤维错排 $\phi_0$ 决定。Hashin准则看不到这种物理。
Puck准则
Puck准则（德国，2002年）考虑基体破坏的破坏面角度。基体压缩破坏会沿板厚方向的倾斜面发生，破坏面角度成为变量，要找最危险的面。
WWFE（世界范围破坏习题）的对比中，Puck准则最接近试验结果。但计算成本高（破坏面角度优化）。
基于相场法的复合材破坏
把相场法应用于复合材，用连续场追踪纤维断裂和基体裂纹扩展的研究激增。不预先假设Hashin式的模式分离，破坏路径自然决定。
相场法不用预先假设破坏模式？
理想上是那样。但复合材各向异性融入相场法的定式还在研究阶段。
总结
整理超越Hashin的破坏理论。
Hashin准则是复合材破坏的「出发点」，更精密预测要推进到LaRC05/Puck/相场法。
Hashin准则的3D扩展：三维复合材
Hashin准则的2D版无法处理面外应力成分σ33·τ23·τ13，对厚积层板和3D纤维强化材不足。Puck（2001年）和Cuntze（2004年）提出了三维破坏准则，特别改进了受斜冲击荷载的3D CFRP（碳纤维织物、3D编织物）破坏预测精度。宇宙机分厚CFRP结合部设计从2010年代开始采用这些3D扩展准则。
Hashin破坏准则的故障排除
Hashin损伤分析的故障
Hashin损伤分析常见故障。
渐进损伤分析是非线性分析中最难收敛的领域。
收敛困难
损伤开始后分析不收敛。
损伤导致刚度急剧下降，平衡路径变不稳定。对策：
1. 粘性正则化 — *DAMAGE STABILIZATION 添加微小粘性。Abaqus推荐值约 $10^{-5}$
2. 减小增量 — 最小增量设为 $10^{-10}$ 左右
3. 稳定化法 — *STATIC, STABILIZE
4. 切换显式法 — 回避收敛问题
粘性正则化值怎么决定？
越小越物理准确但越难收敛。越大越易收敛但结果不准。从 $10^{-5}$ 开始，验证结果是否与粘性值无关（$10^{-4}$ 和 $10^{-6}$ 结果一样吗）。
网格依赖性
网格变化损伤图样也变化。
检查破坏能量的正则化是否正确运作。尤其：
- 特征长度是否正确 — Abaqus自动计算，但歪斜单元会不准
- 单元大小是否适应破坏过程区 — 基体裂纹的过程区约0.5～2mm。单元应为该尺度
0.5～2mm的单元，网格很细啊。
渐进损伤分析对网格要求苛刻。OHC试验孔周围需0.5mm单元。计算成本大，但精度需要这样。
破坏能量不足
破坏能量对不对怎么验证？
总结
整理Hashin损伤的故障对应。
渐进损伤分析是「FEM高手的领域」。
正是。材料特性理解、非线性分析技能、试验数据对接…全部齐备才能进行有意义的渐进损伤分析。
Hashin准则预测破坏早于试验的问题
Hashin准则比试验First Ply Failure（FPF）早20～30%预测破坏时，常是基体拉伸强度YT或面内剪切强度S12设置过低。特别是湿热吸收后强度下降（YT最多30%、S12约20%）过度应用到设计值时，预测会保守。应该以干燥试验值为基准，环境影响用单独的低减系数来处理，才是正确方法。
相关主题
更详细
错误