ABAQUS中Cohesive单元建模方法Word文件下载.docx
《ABAQUS中Cohesive单元建模方法Word文件下载.docx》由会员分享,可在线阅读,更多相关《ABAQUS中Cohesive单元建模方法Word文件下载.docx(15页珍藏版)》请在冰点文库上搜索。
一种是基于traction-separation描述;
另一种是基于连续体描述。
其中基于traction-separation描述的方法应用更加广泛。
而在基于traction-separation描述的方法中,最常用的本构模型为图2所示的双线性本构模型。
它给出了材料达到强度极限前的线弹性段和材料达到强度极限后的刚度线性降低软化阶段。
注意图中纵坐标为应力,而横坐标为位移,因此线弹性段的斜率代表的实际是cohesive单元的刚度。
曲线下的面积即为材料断裂时的能量释放率。
因此在定义cohesive的力学性能时,实际就是要确定上述本构模型的具体形状:
包括刚度、极限强度、以与临界断裂能量释放率,或者最终失效时单元的位移。
常用的定义方法是给定上述参数中的前三项,也就确定了cohesive的本构模型。
Cohesive单元可理解为一种准二维单元,可以将它看作被一个厚度隔开的两个面,这两个面分别和其他实体单元连接。
Cohesive单元只考虑面外的力,包括法向的正应力以与XZ,YZ两个方向的剪应力。
下文对cohesive单元的参数进展阐述,并介绍参数的选择方法。
图2.双线性本构模型
Cohesive单元的刚度
基于traction-separation模型的界面单元的刚度可以通过一个简单杆的变形公式来理解
〔1〕
其中L为杆长,E为弹性刚度,A为初始截面积,P为载荷。
公式〔1〕又可以写成
〔2〕
其中
为名义应力,
为材料的刚度。
为了更好的理解K,我们把
写成:
〔3〕
这里我们用
来代替1,其中L可以理解为建模厚度,即建模时cohesiveinterface的几何厚度;
为实际厚度,即cohesiveinterface的厚度,这个厚度在cohesivesection中定义。
可以理解为几何刚度,即模型中cohesiveinterface所具有的刚度;
为cohesiveinterface的刚度。
当
为1时,计算界面刚度就采用几何刚度
,当
为0.001时,计算时界面刚度变为1000
。
举个小例子,如果界面的实际厚度为0.01,而在建模时就是按照这个厚度建立的,在定义material-section时又specify这层的厚度为0.01,实际上就等于把界面刚度提高了2个数量级,模拟结果当然是不对的,这时定义section时应采用默认厚度1。
ABAQUS在cohesive建模中使用了很“人性化〞的设计,实际问题中界面可能很薄,有的只有0.001mm,甚至更小。
有些问题cohesive单元的interface还可能是0厚度〔比如crack问题〕,而相对来说整体模型也许很大,如果不引入这两个厚度,我们就要在很大的模型中去创建这个很小的界面这是一个很麻烦的事情。
引入这两个厚度,在建模时我们就可以用有限的厚度来代替这个很小的界面厚度,只要在section中定义这个
就好了。
〔注:
以上大局部内容来自仿真论坛:
再议cohesive应用中对于一些参数的理解〕
“另外有个我的经验公式:
大体上energy>
0.5*〔damageinitiation〕^2/〔stiffness〕
这个公式不难理解,就是锐角三角形的总面积大于一条侧边下的面积,将traction-separation
law画成图线你就一目了然了。
不过根据不同的法如此,会稍微有些区别的。
〞------以上的
话引自dava的个人空间,这里我想解释下这个不等式,有些新手可能一下还看不明白。
damageinitiation为开始破坏时的应力,即三角形的高;
stiffness为刚度,也就是斜率,即
tanq;
所以侧边三角形的底边为damageinitiation/stiffness,0.5*〔damageinitiation〕^2/
〔stiffness〕即为侧边下的三角形面积。
实际上能量还要大于这个侧边下三角形的面积很多,
因为斜率一般都很大。
定义cohesive的材料时,要填入材料的参数,材料参数是材料固有的特性,与几
何没有关系,所以放心大胆的填入吧。
材料参数是由试验得到的,如果不能做实验〔多数情
况如此〕,就去查国际上相关的文献吧,数据甚至比你自己做试验都要详细,在填入数据时
要注意单位的统一。
再说句,断裂能为单位面积上的能量,如你的单位选取N〔力的单位〕和M〔长度单位〕,那么能量的单位为N/M。
下面举例来说明cohesive单元刚度的设置过程,以ABAQUS6.9为例:
进入property界面,点击Material→Creat,在弹出的EditMaterial对话框中,可以编辑新创建的cohesive材料的名称,然后点击Mechanical→Elasticity→Elastic→Traction,在空格中输入相应的刚度。
图3.cohesive单元刚度的定义
1.2.4损伤准如此
1.2.4.1初始损伤准如此
初始损伤对应于材料开始退化,当应力或应变满足于定义的初始临界损伤准如此,如此此时退化开始。
Abaqus的Damagefortractionseparationlaws中包括:
QuadeDamage、MaxeDamage、QuadsDamage、MaxsDamage、MaxpeDamage、MaxpsDamage六种初始损伤准如此,其中前四种用于一般复合材料分层模拟,后两种主要是在扩展有限元法模拟不连续体〔比如crack问题〕问题时使用。
使用图2所示的双线本构模型,其中:
、
与
分别代表纯Ⅰ型、纯Ⅱ型或纯Ⅲ破坏的最大名义应力,
,
代表相应的最大名义应变,当定义界面单元的初始厚度为1时,如此名义应变等于与之相对应的相对位移
QuadeDamage为二次名义应变准如此:
当名义应变比的平方和等于1时,损伤开始。
MaxeDamage为最大名义应变准如此:
当任何一个名义应变的比值达到1时,损伤开始。
QuadsDamage为二次名义应力准如此:
当各个方向的名义应变比的平方和等于1时,损伤开始。
MaxsDamage为最大名义应力准如此:
当任何一个名义应力比值达到1时,损伤开始。
图4.初始损伤准如此定义
EditMaterial对话框中,点击Mechanical→DamageforTractionSeparationLaws,然后根据自己的需要点击相应的损伤准如此。
其中最常用是QuadsDamage。
1.2.4.2损伤演化规律
选择了初始损伤准如此之后,然后点击Suboptions→DamageEvolution,窗口如图5所示。
其中Type包括Displacement和Energy,Displacement为基于位移的损伤演化规律,而Energy为基于能量的损伤演化规律。
Softening中包括Linear,Exponential与Tabular三种刚度退化方式……DamageEvolution中的所有的选项都是用来确定单元达到强度极限以后的刚度降阶方式。
一般常用:
以能量来控制单元的退化,即Type→Energy;
线性软化模型,即Softening→Linear,Degradation→Maximum;
Mixedmodebehavior→BK,Modemixratio→Energy,并选中Power。
图5.损伤演化规律定义
1.3Cohesive单元界面属性
还是在Property界面中,点击Section→Create,在弹出的EditSection对话框中,选择Other→Cohesive。
图6.定义材料的界面属性
在EditSection对话框中,在material的下拉菜单中选择刚刚创建的cohesive材料,也可以点击右侧的create创建一组新的材料;
Response选择tractionseparation。
Initialthickness为前文提到的
,默认值为1,也可以在specify中指定一个特定的值。
1.4将所创建的界面属性赋予几何实体
点击Assign→Section,然后在视图中选中要赋的几何实体,点击左下角的Done,如此弹出如下窗口,在窗口是Section中下拉选中所创建的Cohesive截面,点击OK,操作完成。
图7.给实体赋截面属性
1.5Cohesive单元网格划分
Cohesive单元网格的划分与其他单元根本一致,但是以下几点不同与其他单元,划分网格时应特别注意。
网格密度,cohesive单元的网格尺寸不能太大,通常需要比拟精细的网格,不然容易引起收敛性问题,甚至无法继续计算。
必须使用sweep〔扫掠〕划分网格的方法,并且扫掠的方向垂直于cohesive面,即沿着cohesive单元的厚度方向。
单元种类的选择
图8.cohesive单元种类选择
在单元库中选择cohesive,可以在Viscosity,specify中指定一粘性系数,来改善收敛性,但是粘性系数的设置不能太大,不然会影响计算结果,我们一般设置为0.001;
Elementdeletion:
用于设置单元的删除情况,一般选yes,即当单元完全失效时被删除;
maxdegradation:
一般设置为1,即当SDEG=1时,认为单元失效。
2.Cohesive单元在复合材料分层分析中的应用
为了验证商用有限元软件ABAQUS中的cohesive单元在复合材料分层计算时的有效性,我们通过其与一实验值的比照验证了其计算的准确性。
一DCB试验件,长150mm,宽20mm,单臂厚度1.98mm,预置55mm长的初始裂纹,如图9所示。
材料属性为E11=150GPa,E22=E33=11GPa,G12=G13=6.0GPa,G23=3.7GPa,
;
cohesive单元的材料属性为K=1×
105MPa/mm,界面强度T=15MPa,临界能量释放率GIC=0.268KJ/m2。
悬臂梁一端固支,一端施加位移载荷。
〔a〕侧视图
〔b〕俯视图
图9.DCB几何模型
Abaqus和实验[1]得到的力位移曲线如图10所示,从图中可以看出,数值模拟的力位移曲线与实验得到的力位移曲线吻合的很好,数值模拟得到的最大力为65.8N,而实验得到的最大力为62.52N,数值模拟结果略高于实验结果。
由此,我们可以得到有限元软件ABAQUS中的cohesive单元可以有效的模拟复合材料层合板的分层。
计算得到的变形过程的应力与位移云图如图11、12所示。
〔a〕abaqus计算值
〔b〕实验值
图10.实验与数值模拟结果
〔a〕t=0.25s〔b〕t=0.5s
〔c〕t=0.75s〔d〕t=1.0s
图11.变形过程中应力云图
〔c〕t=0.75s〔d〕t=1.0s
图12.变形过程中位移云图
[1]Camanho,P.P.,Davila,C.G.,2002.Mixed-modedecohesionfiniteelementsforthesimulationofdelaminationinpositematerials.NASA/TM-2002-211737.