内聚力模型分析及有限元子程序开发.pdf

A th e sis su b m itte d to Z h en g zh o u U n iv ersity fo r th e d eg re e o f M a ster T h e A n a ly sis o f C o h e siv e Z o n e M o d e l an d U ser-d efin ed S ub ro u tin e D ev e lo P m en t in F in ite E lem en t M eth o d B y L iu g an g H u an g S uP erv iso r :P ro f. Ju n Z h an g E n g in e erin g T h erni o P h y sie s S eh o o l o f C h em ic a l E n g in e erin g an d E n erg y M a y 2 0 10 了 摘 要 摘 要 内聚力模 型是弹 塑性 断裂 力 学 中一个被 广 泛 应 用 与研 究 的计 算模 型 。内聚力 模 型避 免 了线 弹 性 断裂 力 学 中 的裂 纹 尖 端 应 力 奇 异 性 ，计 算 了开 裂 过 程 中应 力 以及 断裂 能 。 内聚 力模 型提 出在 裂 纹 尖 端 存 在 微 小 的 内聚 力 区域 ，该 区域 内应 力 为 开裂 位 移值 的函数 ， 即为 张 力位 移 关 系 。在 对 内聚 力模 型不 断研 究 中 ，发 展 出 了不 同张 力 位 移 关 系 ，使 其 能被 广 泛 应 用 于 各 种 韧 性 开 裂 、复合 材 料 界 面 脱 层 以及 粘 接 界 面 开 裂 等 研 究 中 。 指数 内聚力模 型具有 非线性连 续 的张力位 移 关 系 ，考 虑 了开裂 过程 中各 向位 移对 应 力 的祸 合 关系 。研 究证 实 ，只有在 其 参数 q 等 于 1 的条件 下 ，各 向应 力 关 系才 能完 全 祸合 。本 文提 出 了指 数模 型单 向断裂 能计 算 公式 ，在 各 种参 数条 件 下进 行 计 算 测 试 。其 结 果表 明 ，在 复合 开裂 条件 下 ， 当任 一 向应 力 首 先减 小 为零 时 ，其对 应 的单 向断裂 能亦 同时达 到其 临界 的最 大值 。 有 限元 软件 A B A Q U S 中 内聚力 单元仅 包含 双 线 性 内聚 力模 型 ，其他 类 型 内 聚 力模 型应 用 于 A B A Q u S 时 需通 过 子 程 序 编 译 来 实现 。本 文 基 于 内聚 力 单 元 ， 提 出采 用 自定义 材 料 子程序 V U M AT 编 译 内聚 力模 型 。推 导 了张力位 移 关 系在 V U M AT 中 的应 用 方 法 ，并分 别考 虑 线性 与 非 线 性 内聚 力模 型 ，设计 了 内聚 力 模 型 V U M AT 子程 序 的程序 结构 。 编 译 单 向开裂 条 件 下指数 内聚力模 型 V U M AT 子程 序 ，并 与 内聚力 单元包 含 的双 线 性 内聚 力模 型计算 结 果进 行 比较 ，证 明其 能准 确 模 拟 界 面 单 元 的承 载 以 及 开 裂 失 效 过 程 。在 对 指 数 模 型 藕 合 关 系研 究 的基 础 上 ，采 用 单 向断裂 能 作 为 单 元 失效 的判 据 ，编 译 复合 开裂 条 件 下 V U M AT 子 程 序 。将 其 应 用 于 复合 材 料 层 板 脱 层 模 拟 ，与 实验 测 试 值 进 行 比较 ，结 果表 明 ，通 过 设置 合 理 的模