ABAQUS中Cohesive单元建模方法


    复合材料模型建模与分析 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.1 Cohesive单元的刚度 基于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应用中对于一些参数的理解) 1.2.2 一个解释 “另外有个我的经验公式:大体上energy > 0.5*(damage initiation)^2/(stiffness) 这个公式不难理解,就是锐角三角形的总面积大于一条侧边下的面积,将traction-separation law 画成图线你就一目了然了。不过根据不同的法则,会稍微有些区别的。”------以上的 话引自dava 的个人空间,这里我想解释下这个不等式,有些新手可能一下还看不明白。 damage initiation为开始破坏时的应力,即三角形的高;stiffness 为刚度,也就是斜率,即 tanq ;所以侧边三角形的底边为damage initiation/stiffness,0.5*(damage initiation)^2/ (stiffness)即为侧边下的三角形面积。实际上能量还要大于这个侧边下三角形的面积很多, 因为斜率一般都很大。 1.2.3 关于材料参数 定义cohesive 的材料时,要填入材料的参数,这些材料参数是材料固有的特性,与几 何没有关系,所以放心大胆的填入吧。材料参数是由试验得到的,如果不能做实验(多数情 况如此),就去查国际上相关的文献吧,数据甚至比你自己做试验都要详细,在填入数据时 要注意单位的统一。再说句,断裂能为单位面积上的能量,如你的单位选取N(力的单位)和M(长度单位),那么能量的单位为N/M。 下面举例来说明cohesive单元刚度的设置过程,以ABAQUS6.9为例: 进入property界面,点击Material→Creat,在弹出的Edit Material对话框中,可以编辑新创建的cohesive材料的名称,然后点击Mechanical→Elasticity→Elastic→Traction,在空格中输入相应的刚度。 图3. cohesive单元刚度的定义 1.2.4 损伤准则 1.2.4.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。 1.2.4.2 损伤演化规律 选择了初始损伤准则之后,然后点击Suboptions→Damage Evolution,窗口如图5所示。其中Type包括Displacement和Energy,Displacement为基于位移的损伤演化规律,而Energy为基于能量的损伤演化规律。Softening中包括Linear,Exponential及Tabular三种刚度退化方式……Damage Evolution中的所有的选项都是用来确定单元达到强度极限以后的刚度降阶方式。一般常用:以能量来控制单元的退化,即Type→Energy;线性软化模型,即Softening→Linear,Degradation→Maximum;Mixed mode behavior→BK,Mode mix ratio→Energy,并选中Power。 图5.损伤演化规律定义 1. 3 Cohesive单元界面属性 还是在Property界面中,点击Section→Create,在弹出的Edit Section对话框中,选择Other→Cohesive。 图6. 定义材料的界面属性 在Edit Section对话框中,在material的下拉菜单中选择刚才创建的cohesive材料,也可以点击右侧的create创建一组新的材料;Response选择traction separation。 Initial thickness为前文提到的,默认值为1,也可以在specify中指定一个特定的值。 1.4 将所创建的界面属性赋予几何实体 点击Assign→Section,然后在视图中选中要赋的几何实体,点击左下角的Done,则弹出如下窗口,在窗口是Section中下拉选中所创建的Cohesive截面,点击OK,操作完成。 图7. 给实体赋截面属性 1.5 Cohesive单元网格划分 Cohesive单元网格的划分与其他单元基本一致,但是以下几点不同与其他单元,划分网格时应特别注意。 ①网格密度,cohesive单元的网格尺寸不能太大,通常需要比较精细的网格,不然容易引起收敛性问题,甚至无法继续计算。 ②必须使用sweep(扫掠)划分网格的方法,并且扫掠的方向垂直于cohesive面,即沿着cohesive单元的厚度方向。 ③单元种类的选择 图8.cohesive单元种类选择 在单元库中选择cohesive,可以在Viscosity,specify中指定一粘性系数,来改善收敛性,但是粘性系数的设置不能太大,不然会影响计算结果,我们一般设置为0.001;Element deletion:用于设置单元的删除情况,一般选yes,即当单元完全失效时被删除;max degradation:一般设置为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×105 MPa/mm,界面强度T=15MPa,临界能量释放率GIC=0.268 KJ/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. 变形过程中应力云图 (a)t=0.25s (b)t=0.5s (c)t=0.75s (d)t=1.0s 图12. 变形过程中位移云图 [1] Camanho,P.P.,Davila,C.G.,2002.Mixed-mode decohesion finite elements for the simulation of delamination in composite materials.NASA/TM-2002-211737. 本文档由香当网(https://www.xiangdang.net)用户上传

    下载文档到电脑,查找使用更方便

    文档的实际排版效果,会与网站的显示效果略有不同!!

    需要 2 香币 [ 分享文档获得香币 ]

    下载文档

    相关文档

    建模在化学教学中的应用

    随着时间的推移,教师进行化学教学的过程中建立模型的这种思维方法已经不止是停留在原有的认知基础上,而是能够让学生开阔视野,当面对不同的化学问题时都能迎刃而解的一种有效方法。本文在经过整理调查问卷的...

    1个月前   
    132    0

    abaqus常用技巧总结

    abaqus的step里有maximum number of increment、initial increment、minimum increment 、maximum increment四个...

    2年前   
    1835    0

    业务流程重组中流程建模方法

    摘要 随着我国加入WTO,国内企业面临的市场竞争将比以前更加激烈。在这种条件下,企业迫切需要提高企业的综合竞争能力。而业务流程重组是企业改善管理、提高效益的有效途径。这一思想虽然已经出现了十...

    10年前   
    26351    0

    英语写作手段的运用论文Use of Cohesive Devices on English Compositions

    英语写作手段的运用论文Use of Cohesive Devices on English Compositions By College Students Zhao Jing   A...

    11年前   
    15215    0

    五一建模c题五一建模c题五一建模c题

    五一建模c题 摘要: 0 引言 自改革开放以来,我国在高等教育方面取得了硕果累累的成绩。在实行高等教育扩招后,我国高校的毛入学率逐年攀升,已进入高等教育大众化阶段,各类院校在数量上很大程...

    2年前   
    946    0

    数学建模

    钢管下料问题   队员:周安伟(15895144205)       余乐洲(15189304964)       蔡佳斌 摘要: 在钢材买卖过程中,零售商的刚才都是固定钢管原材料...

    7年前   
    7291    0

    数学建模论文

    舰艇会和问题数学建模论文姓名:班级:学号:舰艇会和问题摘要:当舰艇执行完任务会合航母时,需要采取合适的航行方向与航母会和,可以用坐标系解决这类问题。现代战争中,航空母舰被视为一个国家海军力量的...

    9个月前   
    263    0

    3dsmax在游戏角色建模中的研究与应用

    随着计算机三维影像技术的不断发展,三维图形技术越来越被人类所掌握和应用,由于三维动画比平面图形更直观,更逼真的特点,因此更能给观赏者以身临其境的感觉,

    7个月前   
    437    0

    学习数学建模的心得

    学习数学建模的心得  通过对专题七的学习,我知道了数学探究与数学建模在中学中学习的重要性,知道了什么是数学建模,数学建模就是把一个具体的实际问题转化为一个数学问题,然后用数学方法去解决它,之后...

    7年前   
    153    0

    投资问题数学建模

    数学模型第一次讨论作业问题:某部门现有资金10万元,五年内有以下投资项目供选择:项目A:从第一年到第四年每年初投资,次年末收回本金且获利15%;项目B:第三年初投资,第五年末收回本金且获利25...

    7个月前   
    307    0

    市场营销 课程单元教学设计-定价方法

    《市场营销》课程单元教学设计——《定价方法》一、教案头本次课标题:定价方法授课班级课时2上课地点教学目标能力(技能)目标知识目标能够确定商品的基本价格1.了解影响产品定价的因素、定价目标2.了...

    10个月前   
    245    0

    房地产数学建模

    房地产问题分析 摘要房地产行业与百姓的生活息息相关。近年来,由于房地产价格的不断攀升,房地产行业已经引起了社会的广泛关注。本文分别就影响房地产价格的因素和未来房地产价格的趋势进行了细致的分析研...

    1个月前   
    95    0

    数学建模问题分析

    1、      给出一个所感兴趣的建模的实际问题:上班高峰车辆拥堵情况 (1)      写出问题的实际背景:**发展迅速,人们生活水平提高,私家车越来越多。上班高峰期车辆拥堵严重,通过调查...

    9年前   
    10955    0

    建模注意事项

    建模注意事项首先我们论文出现的毛病:1. 摘要太简单,不吸引人。这说明想吸引美国的评委我们注意一定要把摘要写好写明白写得评委一看就被吸引住了。摘要应该在介绍我们整个论文大体思路的过程中给评委留...

    8年前   
    174    0

    如何写好建模论文

    如何写好建模论文如何写好建模论文一、 写好数模答卷的重要性1.评定参赛队的成绩好坏、高低,获奖级别, 数模答卷,是唯一依据。2. 答卷是竞赛活动的成绩结晶的书面形式。3. 写好答卷的训练,是科...

    6年前   
    139    0

    “自主互助”教学方法在作文教学中的运用

    “自主互助”教学方法在作文教学中的运用 语文是基础学科,肩负着重大的使命。作为语文教师一定要跟上时代的步伐,改变教学理念,面向全体学生,重视“听、说、读、写”训练,全面提高学生的语文素养...

    2年前   
    531    0

    班级工作中的激励方法

    班级工作中的激励方法 激励就是通常所说的调动积极性。在班级管理中能否建立一个良好的班集体,激励起着非常重要的作用。在现实班级管理中,如果找不到新的激励方法,便只有少数人能受到激励,而大多数...

    8年前   
    83    0

    装修中木门刷过漆后的验收方法

    我们知道在装修验收中,如果要说最有难度的东西,那么非木门莫属了,因为当木门做上油漆之后,根本是看不到内部材料的,所以很难区分出好坏;像下面这位业主也知道,不过他却有自己的办法,只是当他看到这种情...

    3年前   
    542    0

    小学作文教学中如何选择恰当的方法

    小学作文是学生认识水平和文字表达能力、思维能力的具体体现,是字、词、句、篇的综合训练。它是一种综合能力的培养。这种复杂能力的培养必须要选择恰当的方法。

    1年前   
    298    0

    婚礼中要掌声的方法

    婚礼中要掌声的方法掌声在婚礼中是最常见的,有时需要我们去引领,有时说到动情处每个人都会自发的鼓掌,为什么我们要在婚礼上要掌声呢?因为掌声是人类最有震撼的肢体语言,是心与心之间最能体现的交流,也...

    7年前   
    100    0

    文档贡献者

    文***品

    贡献于2021-06-09

    下载需要 2 香币 [香币充值 ]
    亲,您也可以通过 分享原创文档 来获得香币奖励!
    下载文档

    该用户的其他文档