abaqus6.11一个从初学到精通粘弹性的分析的经验积累

问题积累（待续）

1.abaqus如何调整图例的大小，就是云图左上角那个图框，字太小了看不清！！

直接设置图例的字体大小就可以：工具栏viewport>viewport annotation options>legend(选项卡)>text(选项)>set font(按钮)>size，修改size选项中的数字，就可以修改图例大小了。

2.cohesive element ＡＢＡＱＵＳ

在6.11使用cohesive element，定义cohesive材料属性的时候主要步骤：

1.定义一个材料的名字，比如cohesive，不要去定义任何属性（弹性，弹塑性等等）。

2.打开工具栏model--edit keywords，在inp中手动添加材料的各种属性。

PS: 定义section的时候选cohesive，element control选sweep，element type选cohesive，这些是使用cohesive element的基本步骤。 zero thickness的cohesive section设定abaqus所谓的

zero-thickness，其实就是定义cohesive section的initial thickness=1.0。你可以在定义section的时候定义（specify），也可以用系统默认的thickness（也是1.0），这样有关cohesive element 的计算当中，就有displacement（位移）=strain（应变）*thickness ( 1.0 )=strain的数值。我们知道从1914年Ingless和1921年Griffith提出断裂力学开始，一直到60年代都停留在线弹性断裂力

学(LEFM)的层次。后来由於发现在裂纹尖端进入塑性区后用LEFM仍然无法解决stress singularity的问题。1960年由Barenblatt 和Dugdale率先提出了nonlinear/plastic fracture mechnics的概念，在裂纹前端引入了plastic zone,这也就是我们现在用的 cohesive fracture mechnics的前身。当时这个概念还没引起学术界的轰动。直到1966年Rice发现J-integral及随后发现在LEFM中J-integral 是等于energy release rate的关系(赞大牛Rice一个，在读Ph.D 期间就推导出对断裂力学有重要意义的J-integral了。。。，实令我辈汗颜)。随后在工程中发现了越来越多的LEFM无法解释的问题。cohesive fracture mechnics开始引起更多的关注。在研究以混凝土为代表的quassi-brittle material时，cohesive fracture mechnics提供了非常好的结果，所以在70年代到90年代，cohesive fracture mechnics被大量应用于混凝土研究中。目前比较常用的方法主要是fictitious crack approach和effective-elastic crack approach或是称为equivalent-elastic crack approach. 其中fictitious crack approach只考虑了Dugdale-Barenblatt energy mechanism而effective-elastic crack approach只考虑了基於LEFM 的Griffith-Irwin energy dissipation mechanism，但作了一些修正。简单来说，我读了一下ABAQUS cohesive element的理论帮助，个人觉得ABAQUS的cohesive element采用的是广泛应用于混凝土的类似fictitious crack的方法。只考虑了Dugdale-Barenblatt energy mechanism。这其中softening law 的影响是非常重要的。

但ABAUQS似乎只提供了linear 或者exponential 的softening law，复杂的本构关系还需要另想办法。至於基於Griffith-Irwin energy dissipation mechanism的J-integral值可以在LEFM分析中单独算。(ABAQUS用的是Suo Zhigang 和Hutchinson在1990一篇论文中提出的方法) 目前cohesive fracture mechnics已经被应用于各种材料。不过在使用到纳米或者更小数量级的研究中碰到了不少问题，可能需要结合位错和分子动力学的一些理论，我现在的研究中也碰到类似的一些问题，希望和大家一起讨论。做裂纹ABAQUS有几种常见方法。最简单的是用debond命令, 定义*FRACTURE CRITERION, TYPE=XXX, 参数。。。** *DEBOND, SLAVE=XXX, MASTER=XXX, time increment=XX0,1, ……......time,0要想看到开裂特别注意需要在指定的开裂路径上定义一个*Nset，然后在*INITIAL CONDITIONS, TYPE=CONTACT中定义master, slave, 及指定的Nset这种方法用途其实较为有限。另一种方法，在interaction模块，special, 定义crack seam, 网格最好细化，用collapse element模拟singularity. 这种方法可以计算J积分，应力强度因子等常用的断裂力学参数.裂尖及奇异性定义:在interaction-special,先定义crack, 定义好裂尖及方向, 然后在singularity选择：midside node parameter: 输入0.25, 然后选Collapsed element side, duplicate nodes，8节点单元对应(1/r)+(1/r^1/2)奇异性。这里midside node parameter 选0.25对应裂尖collapse成1/4节点单元。如果midside nodes不移动到1/4处, 则对应(1/r)奇异性, 适合perfect plasticity的情