粘接界面弹塑性内聚力模型子程序开发.pdfVIP

2014年 1月 郑 州 大 学 学 报 (工 学 版 ) Jan． 2014 第 35卷 第 1期 JournalofZhengzhouUniversity(EngineeringScience) Vo1．35 No．1 文章编号 ：1671—6833(2014)01—0077—04 粘接界面弹塑性内聚力模型子程序开发 张 军 ，贾 宏 ，田 阳 (1．郑州大学 化工与能源学院，河南 郑州 450001；2．天津大学 化学工程学院，天津 300072) 摘 要 ：基于有限元分析软件 ABAQUS内聚力单元 ，采用 自定义材料(VUMAT)编译二维平面状态下弹 塑性 内聚力模型子程序．用所编写的弹塑性 内聚力模型VUMAT子程序对薄板粘接结构的开裂进行模 拟计算，并与双线性 内聚力模型 、多项式 内聚力模型 的VUMAT子程序模拟结果及实验结果进行对比分 析 ．结果表明内聚力关系对粘接结构开裂过程 的拉力变化影响较大，对不同胶体 的粘接结构应选择相应 的内聚 力模 型． 关键词 ：内聚力模型 ；本构关系；自定义材料 中图分类号 ：TQ340．1 文献标志码 ：A doi：10．3969／j．issn．1671—6833．2014．O1．018 0 引言 内聚力模 型 (CZM)是弹塑性断裂力学 中一 个被广泛应用与研究的计算模型，近年来 ，针对不 同的材料性质或结构形式 ，发展出各种不 同形式 一 的内聚力模型 ，使其广泛应用于研究裂纹尖端 塑性变形 、静力和疲劳条件下的蠕变开裂 ，以及岩 图1 弹塑性张力 一位移关系 Fig．1 Relationship ofelasticand plastic 土材料及混凝土 、金属与复合材料界面开裂、粘接 traction-displacem ent 结构开裂等问题 ． 随着计算技术和计算机 的快速发展 ，内聚力 其控制方程为 模型被越来越多地应用于有限元计算 中．内聚力 681； 模型应用于通用有 限元软件 中时，部分形式 的内 聚力模型已经被包含在有限元软件 中，其他形式 61≤6≤62； (1) 的内聚力模型应用时必须进行有 限元子程序开 62≤8≤6 ； 发．笔者在大型有 限元分析软件 ABAQUS的平 台 上 ，借用 ABAQUS的VUMAT接 口利用 FORTARN 6占． 语言编辑了二维弹塑性 内聚力模型，并通过单元 式中： 是粘接界面的张开应力 ；6是张开位移 ； 测试验证其准确性．最后利用编辑 的弹塑性内聚 … 是粘接界面开裂的最大应力 ；6，6，6分别为 力模型对薄板粘接结构的开裂进行 了模拟分析． 开裂时弹塑性关系的形状参数．临界的断裂能值 (即张力位移曲线下包含的面积)为 1 弹塑性 内聚力模型 VUMAT子程序开发