ABAQUS中Cohesive单元建模方法.docVIP

ABAQUS中Cohesive单元建模方法 ABAQUS中Cohesive单元建模方法 PAGE / NUMPAGES ABAQUS中Cohesive单元建模方法 v1.0 可编写可改正 复合资料模型建模与剖析 Cohesive 单元建模方法几何模型 使用内聚力模型（ cohesive zone ）模拟裂纹的产生和扩展，需要在估计产生裂纹的区 域加入 cohesive 层。成立 cohesive 层的方法主要有： 方法一、成立完好的构造 （如图 （1 a）所示），而后在上边切割出一个薄层来模拟 cohesive 单元，用这类方法成立的 cohesive 单元与其余单元公用节点，并以此传达力和位移。 方法二、分别成立 cohesive 层和其余构造零件的实体模型，经过“ tie ”绑定拘束， 使得 cohesive 单元双侧的单元位移和应力协调，如图 1（b）所示。 （a） cohesive 单元与其余单元公用节点 （ b）独立的网格经过“ tie ”绑定 图 1. 建模方法 上述两种方法都能够用来模拟复合资料的分层无效，第一种方法区分网格比较复杂； 第二种方法赋资料属性简单，区分网格也方便，可是装置及“ tie ”很繁琐；所以在实质建 模中我们应依据实质构造选用较简单的方法。 资料属性 应 用 cohesive 单 元 模 拟 复 合 材 料 失 效 ， 包 括 两 种 模 型 ： 一 种 是 基 于 traction-separation 描绘；另一种是鉴于连续体描绘。此中鉴于 traction-separation 描 述的方法应用更为宽泛。 而在鉴于 traction-separation 描绘的方法中， 最常用的本构模型为图 2 所示的 双线性 本构模型。它给出了资料达到强度极限前的线弹性段和资料达到强度极限后的刚度线性降低 融化阶段 。 注企图中纵坐标为应力，而横坐标为位移，所以线弹性段的斜率代表的实质是 cohesive 单元的刚度。 曲线下的面积即为资料断裂时的能量开释率 。所以在定义 cohesive 1 v1.0 可编写可改正 的力学性能时， 实质就是要确立上述本构模型的详细形状：包含刚度、 极限强度、 以及临界 断裂能量开释率， 或许最后无效时单元的位移。 常用的定义方法是给定上述参数中的前三项， 也就确立了 cohesive 的本构模型 。Cohesive 单元可理解为一种准二维单元，能够将它看作 被一个厚度分开的两个面，这两个面分别和其余实体单元连结。 Cohesive 单元只考虑面外 的力，包含法向的正应力以及 XZ， YZ 两个方向的剪应力。 下文对 cohesive 单元的参数进行论述，并介绍参数的选择方法。 图 2. 双线性本构模型 Cohesive 单元的刚度 鉴于 traction-separation 模型的界面单元的刚度能够经过一个简单杆的变形公式来 理解 PL （ 1） AE 此中 L 为杆长， E 为弹性刚度， A 为初始截面积， P 为载荷。公式（ 1）又能够写成 S K （ 2） 此中 S P A 为名义应力， K E L 为资料的刚度。 为了更好的理解 K，我们把 K E L写成： K E E L E L （ 3） L 1 L 这里我们用 L 来取代 1，此中 L 能够理解为建模厚度，即建模时 cohesive interface 的几 何厚度； L 为实质厚度， 即 cohesive interface 的真切厚度， 这个厚度在 cohesive section 中定义。 E L 能够理解为几何刚度， 即模型中 cohesive interface 所拥有的刚度； E L 为 L 2 v1.0 可编写可改正 cohesive interface 的真切刚度。 当 L 为 1 时，计算界面刚度就采纳几何刚度 E L ，当 L 为时，计算时界面刚度变成 1000 E L 。举个小例子，假如界面的实质厚度为，而在建模时 就是依据这个厚度成立的，在定义 material-section 时又 specify 这层的厚度为，实质上 就等于把界面刚度提升了 2 个数目级， 模拟结果自然是不对的， 这时定义 section 时应采纳 默认厚度 1。 ABAQUS在 cohesive 建模中使用了很“人性化”的设计， 实质问题中界面可能很薄， 有 的只有 0.001mm，甚至更小。有些问题 cohesive 单元的 interface 还可能是 0 厚度（比方 crack 问题），而相对来说整体模型或许很大，假如不引入这两个厚度，我们就要在很大的 模型中去创立这个很小的界面这是一个很麻烦的事情。 引入这两个厚度， 在建模时我们便可 以用有限的厚度来