复合材料模型建模与分析
1．Cohesive单元建模方法
1。1 几何模型
使用内聚力模型（cohesive zone）模拟裂纹的产生和扩展,需要在预计产生裂纹的区域加入cohesive层.建立cohesive层的方法主要有：
方法一、建立完整的结构（如图1（a）所示），然后在上面切割出一个薄层来模拟cohesive单元,用这种方法建立的cohesive单元与其他单元公用节点，并以此传递力和位移。
方法二、分别建立cohesive层和其他结构部件的实体模型，通过“tie”绑定约束，使得cohesive单元两侧的单元位移和应力协调，如图1（b）所示。
（a）cohesive单元与其他单元公用节点 （b）独立的网格通过“tie”绑定
图1.建模方法
上述两种方法都可以用来模拟复合材料的分层失效,第一种方法划分网格比较复杂;第二种方法赋材料属性简单，划分网格也方便，但是装配及“tie”很繁琐；因此在实际建模中我们应根据实际结构选取较简单的方法.
1.2 材料属性
应用cohesive单元模拟复合材料失效，包括两种模型：一种是基于traction-separation描述;另一种是基于连续体描述。其中基于traction-separation描述的方法应用更加广泛。
而在基于traction-separation描述的方法中,最常用的本构模型为图2所示的双线性本构模型。它给出了材料达到强度极限前的线弹性段和材料达到强度极限后的刚度线性降低软化阶段。 注意图中纵坐标为应力，而横坐标为位移,因此线弹性段的斜率代表的实际是cohesive单元的刚度。曲线下的面积即为材料断裂时的能量释放率。因此在定义cohesive的力学性能时，实际就是要确定上述本构模型的具体形状：包括刚度、极限强度、以及临界断裂能量释放率，
或者最终失效时单元的位移。常用的定义方法是给定上述参数中的前三项，也就确定了cohesive的本构模型.Cohesive单元可理解为一种准二维单元，可以将它看作被一个厚度隔开的两个面,这两个面分别和其他实体单元连接。Cohesive单元只考虑面外的力，包括法向的正应力以及XZ,YZ两个方向的剪应力.
下文对cohesive单元的参数进行阐述，并介绍参数的选择方法.
图2。 双线性本构模型
1。2。1cohesive单元的刚度
基于traction—separation模型的界面单元的刚度可以通过一个简单杆的变形公式来理解
（1）
其中L为杆长,E为弹性刚度，A为初始截面积，P为载荷。公式（1）又可以写成
（2）
其中为名义应力，为材料的刚度.
为了更好的理解K,我们把写成：
（3）
这里我们用来代替1,其中L可以理解为建模厚度，即建模时cohesive interface的几何厚度；为实际厚度，即cohesive interface的真实厚度，这个厚度在cohesive section中定义。可以理解为几何刚度，即模型中cohesive interface所具有的刚度;为cohesive interface的真实刚度。当为1时，计算界面刚度就采用几何刚度,当为0.001时,计算时界面刚度变为1000.举个小例子，如果界面的实际厚度为0。01，而在建模时就是按照这个厚度建立的，在定义material—section时又specify这层的厚度为0。01,实际上就等于把界面刚度提高了2个数量级，模拟结果当然是不对的，这时定义section时应采用默认厚度1.
ABAQUS在cohesive建模中使用了很“人性化”的设计,实际问题中界面可能很薄，有的只有0。001mm,甚至更小。有些问题cohesive单元的interface还可能是0厚度(比如crack问题），而相对来说整体模型也许很大,如果不引入这两个厚度,我们就要在很大的模型中去创建这个很小的界面这是一个很麻烦的事情。引入这两个厚度，在建模时我们就可以用有限的厚度来代替这个很小的界面厚度，只要在section中定义这个就好了。（注:以上大部分内容来自仿真论坛：再议cohesive应用中对于一些参数的理解）
下面举例来说明cohesive单元刚度的设置过程，以ABAQUS6。9为例：
进入property界面，点击Material→Creat，在弹出的Edit Material对话框中，可以编辑新创建的cohesive材料的名称,然后点击Mechanical→Elasticity→Elastic→Traction,在空格中输入相应的刚度。
图3. cohesive单元刚度的定义
1。2。2 损伤准则
建模方法1。2.2.1 初始损伤准则
初始损伤对应于材料开始退化,当应力或应变满足于定义的初始临界损伤准则,则此时退化开始.Abaqus的Damage for traction separation laws 中包括：Quade Damage、Maxe Damage、Quads Damage、Maxs Damage、Maxpe Damage、Maxps Damage六种初始损伤准则，其中前四种用于一般复合材料分层模拟，后两种主要是在扩展有限元法模拟不连续体（比如crack问题)问题时使用。
使用图2所示的双线本构模型，其中：、及分别代表纯Ⅰ型、纯Ⅱ型或纯Ⅲ破坏的最大名义应力,、，代表相应的最大名义应变,当定义界面单元的初始厚度为1时，则名义应变等于与之相对应的相对位移,及。
Quade Damage为二次名义应变准则：当名义应变比的平方和等于1时，损伤开始。
Maxe Damage为最大名义应变准则:当任何一个名义应变的比值达到1时,损伤开始。
Quads Damage为二次名义应力准则：当各个方向的名义应变比的平方和等于1时，损伤开始。
Maxs Damage为最大名义应力准则：当任何一个名义应力比值达到1时,损伤开始。
图4. 初始损伤准则定义
Edit Material对话框中，点击Mechanical→Damage for Traction Separation Laws，然后根据自己的需要点击相应的损伤准则。其中最常用是Quads Damage。
版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。
发表评论