abaqus问题解析

秦郁文 2016年1月4日浏览：6815 评论：1 收藏：9

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的情况.网格划分:裂尖网格划分有一些技巧需要注意，partition后先处理最外面的正方形，先在对角线和边上布点，记住要点constraint, 然后选第三个选项do not allow the number of elements to change不准seed变化，密度可以自己调整. 最里面靠近圆的正方形可以只在对角线上布点. 也可以进一步分割内圆及在圆周上布点. 里面裂尖周围的内圆选free mesh, element type选cps6或者cpe6，外面四边形选sweep mesh, element type选cps8或者cpe8, 记住把quad下那个缩减积分的勾去掉。

补充一下这种方法的几个常见问题，见不少朋友问过。主要是对断裂力学的理解问题，顺便一起回答一下。1．为什么我设置理想弹塑性(epp)分析的时候得到的xx,yy方向或者最大应力值Sxx, Syy会超过材料的屈服强度Sy呢, 这分析结果可能吗？这是因为在ABAQUS中对应等于材料的屈服强度的是von Mises等效应力Se=Sy，因此在平面应变的条件下，xx方向的应力Sxx=Sy*pi/SRQT(3)>Sy, 而Syy=Sy*(2+pi)/SRQT(3), 大概是3倍的屈服应力。所以得到大于材料的屈服强度的xx及yy方向应力是正常的。2．为什么设置collapse element的时候对弹性分析在中间就一个点而要把单元边上的中点移到1/4处，但弹塑性分析却要在中间设置一圈点并且保持单元边上的中点位置不变呢？这个其实不是随便定的，在有限元中分析裂纹时，对弹性分析需要模拟裂尖1/SQRT(r)的奇异性，这样在把单元边上的中点移到1/4处后计算出来的等参单元拉格郎日型函数对应的u field正好包含1/ SQRT(r)项，事实上这一方法在断裂力学的数值模拟发展史上是很巧妙的一个发现，至今仍然被广泛采用。至于理想弹塑性分析需要模拟裂尖1/r的奇异性, 这样大家都知道在把单元边上的点放在到1/2处后计算出来的正常的等参单元拉格郎日型函数对应的u field包含1/ r项, 可以模拟弹塑性分析需要的裂尖1/r的奇异性。所以在看似动手点几下就能实现的分析模式后面有很清楚漂亮的理论作支持。还有就是比较新的cohesive element单元。需要定义damage initiation和evolution的准则， softening准则目前只有linear和exponential，但对一般材料也够用了。然后通过设置后处理display group可以看到裂纹扩展情况。裂纹扩展不是ABAQUS的强项，目前比较方便的只能用cohesive element，我做过几个模型效果还可以，但对应的参数需要一定的实验数据支持，否则做出来了也不知道对不对。或者使用python控制seam尺寸,然后移动partition和网格，比较麻烦。

3.分析步

分析步主要有初始分析步和后续分析步，每个分析步可以用来描述一个分析过程，例如在后续分析步中施加不同荷载，在初始分析步中施加边界条件等。

4.增量步
增量步是在分析步里面根据模型计算收敛情况设置的，简单模型可以设置较少的增量步，并可使初始增量为1；复杂模型设置多一点增量步，并减少初始增量值。超过设置的允许增量步数，则计算停止。

补充；静力分析时，分析总时间(Time period)，是个相对概念，没有单位的。尽管静力分析，你的荷载可能是一下子定义施加上了，但是aba计算时，数值积分是一步一步迭代计算的，需要N个迭代步，每个迭代步都需要求解联立耦合的方程组，程序计算需要一定的时间，这个时间就是时间增量increments，假若一共迭代N步，这N个increments之和就是分析总时间(Time period)。

一般情况下，程序一开始按照你设定初始时间增量步进行计算，
1）若是在这个初始时间增量步下不能收敛，就会减小时间增量直到减小到你设定的最小时间增量为止，若还是不能收敛，则程序就会给你报错。
2）若是按照初始时间增量步可以顺利计算，程序会自动调大时间增量，一边加快计算进度，时间增量步调大的上限就是你设定的最大增量（Maximum number of increments）。

故而，最大增量步，初始时间增量，最小时间增量，最大时间增量等你的预先设置，在很大程度上影响着程序能不能顺利计算和计算速度的快慢。

5接触问题

接触问题，如果一开始并没接触上，是有间隙的，那么第一步就是要接触上。可以用一个小载荷，但用一个位移，比初始间隙稍微大一点点，或许更好。然后在后继步中，去掉这个位移，并加上载荷。甚至可以加一个大一点的位移，用比较小的增量步，多输出几个增量步的结果，查看接触力的增加。在快要到达所应加的力的地方，做一个重启动，去掉位移，把最后的力加上去。加位移要比加力更容易收敛。

6.abaqus中的收敛

要取得一个精确的结果需要以下这些：

a,网格的收敛

b,瞬态问题的精确的时间整合

c,非线性求解过程适当的收敛

另外，一个精确的结果还需要在创建一个合适的有限元模型时好的工程评价，包括材料，载荷，边界条件以及求解过程。

补充：如果不收敛，一般使计算陷于死循环，为了避免这种情况，程序里都加入是否收敛的判断.
有限元法也是近似解法，收敛是作为计算结束，结果无限接近真实结果的一种判定标准。

7.在abaqus配置完成后，点击运行，如果出现Server Start Failed. The Server May Already Be Running!解决方法是先点击start server右边的stop server，然后再点击start server，如果配置没错的话，这里会Server Start successful的。

8 如何显示梁截面（如何显示三维梁模型）

a）无论是运算还是默认显示，ABA中的梁都是一条线，很多人想看梁截面（一般一个星期有人问一次）。

显示梁截面：view->assembly display option->render beam profiles，自己调节系数

http://forum.simwe.com/viewthread.php?tid=835478&page=1#pid1531086

b）后处理到底能不能显示梁截面？

在deformed shape和undeformed shape都能用上面的方法显示梁截面，在应力云图（contour）不能显示。

c）也经常有人问起如何显示壳单元的厚度

http://forum.simwe.com/thread-865887-1-1.html

9 怎么在局部坐标系下建立参考点

在前处理中，已经建立了局部坐标系，如何在局部坐标系中建立参考点？

这个有点麻烦，看看konadoul图文并茂的示例吧。

http://forum.simwe.com/viewthread.php?tid=863389&highlight=%D7%F8%B1%EA%CF%B5