复合材料低速冲击插件:基于ASTM D7136标准自动化建模(Composite Impact Auto‑Builder)
一、低速冲击建模方法的系统评估与ASTM标准
插件的底层建模策略主要依据“Assessment of numerical modeling approaches for thin composite laminates under low-velocity impact”(DOI: 10.1016/j.tws.2023.111053)。该工作围绕碳纤维/环氧复合材料层合板的低速冲击行为,系统比较了不同损伤起始准则、损伤演化方法和界面模型的预测能力,旨在确定一种高精度的数值建模组合。
图 1. 插件计算效果图:左图层间损伤,右图层内损伤
1.1. 模型框架
插件采用切分发创建层内和层间模型,层间界面使用有限厚度(0.001 mm)内聚力单元进行建模。根据文献结果,界面模型的选择从加载初期即显著影响位移和接触时间,零厚度模型会因忽略界面实际厚度而低估最大位移,有限厚度模型则更能准确复现实验响应。
图 2. 自动建模插件
1.2. 材料系统与参数
文献中的材料为 T700/M21,其层内力学参数、强度值、断裂能以及界面参数均见于文献表 3(亦为插件预置值)。层间损伤演化为二次应力准则与 B‑K 混合模式能量准则。这些参数构成插件中 T700 材料数据库的核心。
1.3. 建模细节说明
文献建立的有限元模型中,层合板尺寸为 150 mm × 100 mm,厚度 4 mm,铺层方案包括对称与非对称序列;冲头质量 2 kg,直径 16 mm,冲击能量 25 J(初速 5 m/s)。边界条件参照ASTM标准设置,即在 125 mm × 75 mm 矩形框内支撑试件,仅约束面内平移自由度,不约束法向。插件的边界建模即复现了这一试验构型。
二、插件功能详解
Composite Impact Auto‑Builder 是基于 Abaqus/Explicit 平台的参数化前处理插件,将上述文献验证的高精度建模流程封装为可重复运行的自动化工具,减少人工操作与参数设定偏差。
图 3. Composite Impact Auto‑Builder 1.1 插件界面
2.1 内嵌文献标定的材料参数库
插件内置两组常见的碳纤维复合材料体系:
- T700:DOI: 10.1016/j.tws.2023.111053
- T300:DOI: 10.1177/07316844221147487
用户在 GUI 界面仅需通过下拉菜单选择材料型号(T700 或 T300),VUMAT 实体本构属性及 Cohesive 层截面参数将自动完成赋值。
2.2 精确的有限厚度层间界面建模
与传统的共节点连接或零厚度内聚力模型不同,插件遵循文献结论采用有限厚度的独立内聚力层来模拟分层行为。其实现方式为“逐层剥离”算法:程序读取铺层总数后,在上、下相邻铺层之间自动插入独立的 Cohesive 层几何体,每层厚度为 0.001 mm。
同时,在层合板内每一相邻实体层与 Cohesive 层界面之间,自动建立 Surface‑to‑Surface Penalty 面面接触对(法向硬接触、切向罚摩擦系数 0.3),实现层间正应力与剪应力的真实传递。该设置还原了文献中有限厚度模型对最大中心位移和接触时间更为准确的预测能力。
图 4. 低速冲击建模-模型图
2.3 局部网格控制与铺层自动补齐
网格划分:程序自动识别冲击中心区域(50×50 mm),并在此范围内执行网格加密。用户可独立设置中心加密区与外围粗化区的单元尺寸,兼顾计算精度与效率。所有实体层采用 C3D8R 减缩积分单元并激活单元删除,内聚力层采用 COH3D8 单元,冲头则使用离散刚体单元 R3D4。网格划分基于扫掠技术(Advancing Front)生成。
铺层逻辑:支持非对称铺层序列输入,用户通过逗号分隔输入各层角度。若输入角度数目少于设定的总层数,系统将自动以 0° 铺层补齐所缺层信息,避免因输入遗漏导致模型铺层角度不完整。
图 5. 低速冲击模型-网格划分
2.4 标准试验边界条件的复现
插件内置符合 ASTM D7136 规范的边界模型。一键生成的刚体承载支座包含 125×75mm 的矩形窗口,有限元模型中通过约束支座参考点实现固定。层合板四边的约束条件设置为非完全固支:约束面内位移 U1、U2 以及三个转动自由度 UR1、UR2、UR3,但释放法向位移 U3,从而还原靶板在冲击载荷下的实际弯曲变形形态。
2.5 接触、载荷与分析步自动化
插件自动定义冲头与上表面、各层间界面以及下表面与支座之间的接触关系,统一使用罚函数摩擦系数 0.3,法向硬接触。冲击速度通过预定义场赋予冲头(初始速度沿法向负方向,默认 4430 mm/s,对应约 10 J 能量示例,用户可调)。分析步采用显式动力学,时间周期默认 0.01 s,场输出包含应力 S、应变 E、位移 U、损伤变量 SDEG 和 DMICRT、状态变量 SDV 及 STATUS,历史输出请求接触面法向力 CFN3,便于后处理中快读提取力‑时间/位移曲线。
图 6. 后处理模块可单独显示每层损伤
2.6 版本授权,插件分为两个发布版本
- 基础版(公开发布):面向基础科研与算法验证,核心建模功能与材料库无任何限制,但最高支持生成 4 层的复合材料层合板。在层数 ≤4 时,可直接运行生成模型。
- 专业版(需 License 授权):当用户设定层数 >4 时,程序将调用设备特征码,并提示缺少 license.key 或密钥无效。用户需将机器码发送到微信公众号,有条件获取授权文件,放置于插件目录即可解除层数限制。
基础版获取方式:点赞+推荐+关注后,公众号内回复关键词“冲击插件”。
工程师必备
- 项目客服
- 培训客服
- 平台客服
TOP




















