脆性材料的失效判定方法
在强度仿真中判定脆性材料的失效，需要根据其“无塑性变形、拉应力主导断裂”的特性，选择合适的失效判据和仿真方法。
以下是关键步骤和常用判据的详细说明：
01
—
脆性材料的失效特点
脆性材料（如硅、玻璃、陶瓷等）的失效通常表现为：
无塑性变形：
断裂前几乎无屈服或延展。
拉应力敏感：
断裂主要由拉应力引发，抗压强度远高于抗拉强度。
缺陷敏感：
强度受内部微裂纹、孔隙等缺陷影响显著（需考虑统计分布，如Weibull理论）。
突然断裂：
失效过程无预警，表现为瞬时脆性断裂。
02
—
常用失效判据
仿真实现方法
最大主拉应力准则（Rankine准则）
原理：
当最大主拉应力（σ₁）超过材料抗拉强度（σ_UTS）时失效。
适用场景：
简单应力状态、初步设计分析。
ANSYS中的实现：
在ANSYS后处理中提取Maximum Principal Stress，手动比较σ₁与σ_UTS。
优点：简单直观，计算量小。
缺点：忽略多轴应力效应和缺陷影响。
Mohr-Coulomb准则
原理：
综合考虑拉应力和剪切应力，适用于拉-剪复合应力状态。
公式：
其中，τ为剪应力，c为内聚力，ϕ为内摩擦角（脆性材料通常ϕ≈0)。
适用场景：
存在压-剪复合应力的结构（如陶瓷轴承）。
ANSYS中实现步骤：
在workbench中定义新材料，并在左侧geomechnical 中找到Mohr-coulomb 。
定义如下5个参数：
有关这5个参数的详细定义，我将在后续的文章中具体解释
Weibull统计失效理论
原理：基于缺陷分布的统计模型，预测失效概率。
公式：
其中，m为Weibull模数，σ0为特征强度，V0为参考体积。
仿真实现：
ANSYS PDS：
将σ_UTS设为随机变量，进行蒙特卡洛模拟。
自定义脚本
导出单元应力/体积数据，计算积分和失效概率。
适用场景：高可靠性要求（如航天器陶瓷隔热瓦）。
断裂力学方法（裂纹扩展分析）
原理：基于应力强度因子（SIF）或能量释放率（G）判定裂纹扩展。
判据：
KI≥KIC（I型裂纹）
G≥GIC（临界能量释放率）
仿真实现：
ANSYS中使用Fracture Tool计算SIF，或插入内聚力单元（CZM）。
适用场景：
已知初始缺陷或裂纹位置（如硅片边缘微裂纹）。
连续损伤力学（CDM）
原理：
通过损伤变量（D）描述材料刚度退化，D=1时完全失效。
仿真实现：
ANSYS APDL中使用TB,DMG定义损伤演化律
适用场景：
模拟损伤累积过程（如陶瓷热震失效）。
03
—
仿真关键步骤
ANSYS为例
材料属性设置：
弹性参数：
EX（弹性模量）、PRXY（泊松比）。强度参数：
σ_UTS（抗拉强度）、K_{IC}（断裂韧性）或Weibull参数。若用CZM：定义法向/切向刚度（
Kn,Ks）和临界断裂能（G_IC）。几何与网格：
在应力集中区域（如孔洞、尖角）加密网格。
裂纹路径或界面处使用内聚力单元（如
INTER205）。载荷与边界条件：
施加工况载荷（力、温度场等）。
包含残余应力可通过
Initial Stress或热膨胀模拟。求解与后处理：
提取最大主应力、SIF或损伤变量。
对比判据阈值（如σ₁ ≥ σ_UTS）。
04
—
注意事项
参数标定：
脆性材料的强度（如σ_UTS、Weibull参数）必须通过实验标定，避免直接使用文献值。
尺寸效应：
小尺寸结构可能表现出更高强度（缺陷少），需调整参数。
残余应力：
制造工艺（如烧结、沉积）引入的残余应力需在仿真中考虑。
网格敏感性：
峰值应力对网格密度敏感，需进行收敛性分析。
05
—
总结建议
快速评估：
优先使用最大主应力准则 + 安全系数。
可靠性分析：
采用Weibull统计方法。
裂纹问题：
选择断裂力学或内聚力模型。
实验验证：
任何仿真结果均需通过断裂实验验证（如三点弯曲、双轴测试）。
