ABAQUS中Cohesive单元建模方法 - 图文
更新时间:2024-01-30 02:01:01 阅读量: 教育文库 文档下载
复合材料模型建模与分析
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模型的界面单元的刚度可以通过一个简单杆的变形公式来理解
??PL (1) AES (2) K其中L为杆长,E为弹性刚度,A为初始截面积,P为载荷。公式(1)又可以写成
??其中S?PA为名义应力,K?EL为材料的刚度。
为了更好的理解K,我们把K?EL写成:
K?EELEL?? (3)
?L1L这里我们用L?来代替1,其中L可以理解为建模厚度,即建模时cohesive interface的几何厚
度;L?为实际厚度,即cohesive interface的真实厚度,这个厚度在cohesive section中定义。
EL可以理解为几何刚度,即模型中cohesive interface所具有的刚度;
EL为cohesive L?interface的真实刚度。当L?为1时,计算界面刚度就采用几何刚度EL,当L?为0.001时,计算时界面刚度变为1000EL。举个小例子,如果界面的实际厚度为0.01,而在建模时就是按照这个厚度建立的,在定义material-section时又specify这层的厚度为0.01,实际上就等于把界面刚度提高了2个数量级,模拟结果当然是不对的,这时定义section时应采用默认厚度1。
ABAQUS在cohesive建模中使用了很“人性化”的设计,实际问题中界面可能很薄,有的只有0.001mm,甚至更小。有些问题cohesive单元的interface还可能是0厚度(比如crack问题),而相对来说整体模型也许很大,如果不引入这两个厚度,我们就要在很大的模型中去创建这个很小的界面这是一个很麻烦的事情。引入这两个厚度,在建模时我们就可以用有限的厚度来代替这个很小的界面厚度,只要在section中定义这个L?就好了。(注:以上大部分内容来自仿真论坛:再议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问题)问题时使用。
0使用图2所示的双线本构模型,其中:tn、ts0及tt0分别代表纯Ⅰ型、纯Ⅱ型或纯Ⅲ破0坏的最大名义应力,?n、?s0,?t0代表相应的最大名义应变,当定义界面单元的初始厚度为
1时,则名义应变等于与之相对应的相对位移?n,?s及?t。
Quade Damage为二次名义应变准则:当名义应变比的平方和等于1时,损伤开始。
??n???s???t??0???0???0??1 ??n???s???t?Maxe Damage为最大名义应变准则:当任何一个名义应变的比值达到1时,损伤开始。
222??n?s?t?max?0,,0??1 0?s?t???nQuads Damage为二次名义应力准则:当各个方向的名义应变比的平方和等于1时,损
伤开始。
?tn??t??tt???0??1 ?0???s0??tn??ts??tt?Maxs Damage为最大名义应力准则:当任何一个名义应力比值达到1时,损伤开始。 max?222?tn0?tn,tstt?,0??1 0tstt?
图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为前文提到的L?,默认值为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,?12?0.25,?13?0.25,?23?0.45;
界面强度T=15MPa,临界能量释放率GIC=0.268 KJ/m2。悬臂梁一端固支,一端施加位移载
荷。
(a)侧视图
(b)俯视图 图9. DCB几何模型
Abaqus和实验[1]得到的力位移曲线如图10所示,从图中可以看出,数值模拟的力位移曲线与实验得到的力位移曲线吻合的很好,数值模拟得到的最大力为65.8N,而实验得到的最大力为62.52N,数值模拟结果略高于实验结果。由此,我们可以得到有限元软件ABAQUS中的cohesive单元可以有效的模拟复合材料层合板的分层。计算得到的变形过程的应力及位移云图如图11、12所示。
7060504030201000246810位移(mm)12141618力(N)
(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.
正在阅读:
ABAQUS中Cohesive单元建模方法 - 图文01-30
MBTI ENFP(外向 直觉 情感 知觉) 公关型 完整分析报告01-16
本科生毕业论文《珠海格力电器股份有限公司财务报告分析》11-11
安徽省黄山市2016-2017学年高二上学期期末质量检测语文试题含答05-17
英语先进教师介绍10-26
《回忆鲁迅先生》教案11-28
沟通技巧之十大常见又很棘手的问题09-04
102液压举升钢管4-M10打孔及攻丝螺纹夹具设计05-25
- exercise2
- 铅锌矿详查地质设计 - 图文
- 厨余垃圾、餐厨垃圾堆肥系统设计方案
- 陈明珠开题报告
- 化工原理精选例题
- 政府形象宣传册营销案例
- 小学一至三年级语文阅读专项练习题
- 2014.民诉 期末考试 复习题
- 巅峰智业 - 做好顶层设计对建设城市的重要意义
- (三起)冀教版三年级英语上册Unit4 Lesson24练习题及答案
- 2017年实心轮胎现状及发展趋势分析(目录)
- 基于GIS的农用地定级技术研究定稿
- 2017-2022年中国医疗保健市场调查与市场前景预测报告(目录) - 图文
- 作业
- OFDM技术仿真(MATLAB代码) - 图文
- Android工程师笔试题及答案
- 生命密码联合密码
- 空间地上权若干法律问题探究
- 江苏学业水平测试《机械基础》模拟试题
- 选课走班实施方案
- 建模
- 单元
- Cohesive
- 图文
- 方法
- ABAQUS
- 广州市建筑报建审批专业管理暂行规定
- 回采地质说明书
- 引风机变频器操作说明及注意事项
- 保代考试题目7
- 1、《邹忌讽齐王纳谏》
- 2016届辽宁省抚顺市抚顺县中考二模数学试卷(带解析)
- 装配式混凝土结构预制构件加工制作前应进行深化设计
- 氨区防止泄漏及人身伤害措施
- 2015年大外科护理管理制度考试试题(A卷)
- 单片机c51下的自动打铃系统
- 论隋唐时期的中日佛教文化交流
- ArcGIS让你的标注炫起来2(Maplex)
- 赤壁赋理解性默写
- 钢筋识图入门
- (埃及)国别报告
- 教育工作会议精神落实情况报告
- 2015年浙江农林大学《风景园林综合理论》考试大纲招收硕士研究生入学考试试题大纲考研大纲
- 钢筋混凝土(结构设计原理)T型截面梁算例
- 家庭理财规划书 -
- 感觉阈值检测仪的临床依据