论坛汇总第一部分

***warning:Convergence judged unlikely.Increment will be attempted again with a time increment of 6.25E-02.大家都知道，接触问题往往会遇到这种情况。但这段话我不太理解，恳请大家解释一下。另外，hcrit的值在哪里修改？我的inp里*contact pair 中只设置了ajust=0,并没有给出hcrit=0.33138,而警告让我修改该值，这是怎么回事？ 那是因为用的是默认值，你自己加上去就好了。

我在*contact pair option里直接加上了hcrit的值，结果还是同样出错。我把hcrit的值改大，还是不行，再改大，也是不行。这种情况下应该怎么办，总不能不断改hcrit的值吧，如果hcrit的值太大，即使收敛，结果估计也不对，应该怎么办呢？ 请大虾指教。

说真的，我并不是非常了解情况。你可以去参考手册analysis user's manual，我见里面有写，但是记不住。我碰到问题就查，曾经碰到过增大一点就可以的情况。

我把模型简化了，减少了接触面，不收敛。改大hcrit值，收敛了，可后处理显示明显有穿透发生。到底该怎么办，望大虾指点。

解决这个问题一方面需要把单元划分得更细些，使得复杂表面网格划分后不至于变化太大；另一方面需要修改允许的接触容差，通常系统默认hcrit值为最小单元尺寸长度的一半，然而在求解过程中，由于有限元方法通常在单元边上是线性插值，变形后单元边界仍需保持为直线，这很容易导致求解过程中超过设定的容差值，因此需要改写hcrit值。另外，做接触问题用explicit就没有这么高的要求，你可以试试，当然是要以机时为代价的。

多谢大侠指点，确实很有收获。另外问一下，hcrit的具体含义是什么？是否是允许穿透的距离？

Hcrit是每一个increment计算中允许的salve node的穿透距离，当一个increment中salve node的穿透距离大于或者等于hcrit值， abaqus就选用更小的increment步长，此时得到较小的穿透距离，然后迭代计算。Hcrit不是一般意义上的允许穿透距离，因此增大Hcrit的值，只不过为了减小increments，减少迭代而已。一般来说最后的穿透应该和网格的划分密度有关，和hcrit关系不大吧，再有，由于用了一阶单元，单元边是直线，单元边穿透是必然的，只要salve nodes 没有穿透，说明模型没有问题的。因此，1细画salve surface的单元网格，2,减小初始的increment步长。

完全同意，不过在进行这个方法前检验边界条件是否正确，接触对是否正常，因为有些时候结构根本就不再稳定，所以上述2种方法也不好用了。

也就是说100个increments和10000个increments最终算出的结果是一样的（都收敛的情况下）？ 四十五、请教，关于连续配筋混凝土的建模

最近在做些CRCP的计算，其中遇到几个问题希望那位大侠给予指点

1.在ABAQUS中如何模拟CRCP的细微裂缝，其裂缝存在一定的抗剪切能力2.我应用的是CPE8单元，看一些例子时，钢筋直接以*REBAR命令加入，而我建立了模型后可以算，但是有警告信息，说是要用*EMBEDED ELEMENT命令，请问如何应用*EMBEDED命令在平面应变模型中加入加筋？

第二个问题，abaqus只是给你一个建议而已，*rebar的功能也可以用*embeded element 实现。请自己查看用户手册 四十六、讨论：锚杆在abaqus中如何实现 用truss单元。提供一篇文章参考

TiebackWallsinSandNumericalSimulationandDesign.pdf (302.99k) 对于锚杆，我最近尝试用skin来近似模拟，效果还凑合

我做土钉模拟听说flac模拟土更好些我看有把土钉用bar单元土钉和土间用摩擦做的说模拟的效果还可以但是处理连接和一些参数得细心了

我刚看了，好像还是用beam单元做得，可是在三维里做模型还难做啊，beam单元没法用embedded element。 6.5版本现在可以了

用Beam单元模拟锚杆不合理,最好用杆单元!

在三维有限元分析中,一般很难用某种单元模拟锚杆,更何况模拟的效果也不好,因此,还是寻求别的办法,若要模拟的好,必须采用非连续法

烦请具体谈谈非连续法.的确，前面朋友提供的论文很有参考价值，而我仔细研究了这个实例，发现这种思想还是挺好的。就是在加预应力的时候，我遇到了难题。不知作者是如何把预应力加上去的。

ABAQUS中不同类型的单元无法在CAE中直接耦合在一起，不知道你是如果处理得，是否采用了其他的前处理工具来处理，我现在也在做锚索和锚杆的问题，由于模型一直建立不起来，所以进行不下去了，希望能向你学习。

- 37 -

论坛汇总第一部分

你好，谈不上学习，最近因为忙着其他事，锚杆的模拟也就放在一边。我认为处理预应力锚杆在inp里应该是可以用杆单元或弹簧单元来作，关键是加力得方式问题。目前，我还没找到捷径 我觉得你可以试试REBAR单元

rebar单元我试过，我也问过北京那边的办事处专业人员，它可以模拟锚杆，但后处理中显示不好，只是等效形式，而没有单个锚杆的应力或轴力图

用ＡＢＡＱＵＳ来模拟隧道开挖，其模型好建吗？你是用什么来建的？我也是搞隧道开挖方面的，我现在用ＡＮＳＹＳ在做，感觉里面提供的材料模型不太好！ＡＢＡＱＵＳ里的材料模型对模拟岩土类材料要好很多．

- 38 -

论坛汇总第一部分

四十七、Abaqus\\CAE加载技巧

cae确实只能实现有限的功能，这些功能限制可以在手册里查到（表格形式给出）。 载荷可以以力或者位移方式给出，在step里根据计算过程一步步加载。位移加载节约资源，比较容易收敛。承载力可以用求\的方法， 在output requirement 里选上这个选项即可，reaction force 可能以node force 给出，只需要将某个面上所有节点反力加起来，就是该面的反作用力。 abaqus/cae里的后处理，可以实现将你所想要的node set 任意编辑选出，组成显示。这时候，如果你想将挑选出的node set 编组，到input 文件里添上一句指令，node set， followed with 相应的nodes 编号就可以了。

有些时候问题的物理性质决定只能采用力加载，比如两个主应力方向的力比例加载，在非线性情况下位移之间不存在比例关系；在比如简支梁在几个集中力作用下的挠度非线性分析，此时可以通过Riks弧长法来处理说到加力，顺便问一下斑竹，先张预应力和后张预应力怎样加以区别并模拟？ 先张预应力：

*INITIAL CONDITIONS, TYPE=STRESS, REBAR

注意要先进行一次预分析使之达到静力平衡先。后张预应力：

在进行静力预分析时加上*PRESTRESS HOLD使得张拉预应力保持不变。

请问楼主，在包含摩擦接触对的问题中，你的方法有效吗?比如说研究被动桩，桩与土之间采用摩擦接触对，对桩顶采用水平位移与竖向位移成比例加载方式，怎么得到水平承载力？桩必须是刚体吗？好像以前有人说过，对于刚体有办法！

另外一个问题，我大部分用的预应力钢筋并不能用rebar 来模拟。比如张膜结构的钢索。我一般用杆单元，加 *INITIAL CONDITIONS, TYPE=STRESS这样前张和后张是否也可以区分？

没做过，*PRESTRESS HOLD命令需要给定rebar的名称，因此，不用rebar参数可能 很难将先张后张区分，可能得想其他办法。

不过我的理解是，Rebar其实也是杆单元，只是两端节点的坐标根据母体的坐标内插确定，而不需要单独建模，因此，杆单元能模拟的，Rebar应该也可以做吧，钢索可以通过定义材料*No Compression来模拟的。 我也是不懂，想了好久也没有好的解决办法，上来集思广益一下。这里高人多。其实我个人觉得先张和后张预应力除了钢筋应力前后变化的不同外，还在预应力的传递上有所不同。广义的预应力应用主要是后张的原理。 比如对如下钢索施加预应力，假定钢筋用杆单元模拟，先张就要保证杆单元（钢索）和块单元（混凝土板）的结点重合以保证预应力的传递；后张的话（应该也是现实中实现的），用于传递预应力的锚具的模拟一直只让我头疼，没有想到较好的、简单的方法实现，所以问问各位的高见。 此主题相关图片如下：

浙大BBS：abaqus分析技巧

采用abaqus的cae进行力学问题的分析，其对模型的处理存在很多的技巧，对abaqus的一些分析技巧进行一些概述，希望对大家有所帮助 1.abaqus的多图层绘图

abaqus的cae默认一个视区仅仅绘出一个图形，譬如contor图，变形图，x－y曲线图等，其实在abaqus里面存在一个类似于origin里面的图层的概念，对于每个当前视区里面的图形都可以建立一个图层，并且可以将多个图层合并在一个图形里面，称之为Overlay Plot，譬如你可以在同一副图中，左边绘出contor图，右边绘出x－y图等等，并且在abaqus里面的操作也是很简单的。

1.首先进入可视化模块，当然要先打开你的模型数据文件（。odb）

- 39 -

论坛汇总第一部分

2.第一步要先创建好你的图形，譬如变形图等等

3.进入view里面的overlay plot，点击creat，创建一个图层，现在在viewport layer里出现了你创建的图层了

4.注意你创建的图层，可以看到在visible 下面有个选择的标记，表示在视区里面你的图层是否可见，和autocad里面是一样，取消则不可见current表示是否是当前图层，有些操作只能对当前图层操作有效，同cad name是你建立图层的名称，其他的属性值和你的模型数据库及图形的类型有关，一般不能改动的。 5.重复2-4步就可以创建多个图层了

6.创建好之后就可以选择plot/apply，则在视区显示出所有的可见的图层

1.什么是子结构

子结构也叫超单元的（两者还是有点区别的，文后会谈到），子结构并不是abaqus里面的新东东，而是有限元里面的一个概念，所谓子结构就是将一组单元组合为一个单元（称为超单元），注意是一个单元，这个单元和你用的其他任何一种类型的单元一样使用。 2.为什么要用子结构

使用子结构并不是为了好玩，凡是建过大型有限元模型的兄弟们都可能碰到过计算一个问题要花几个小时，一两天甚至由于单元太多无法求解的情况，子结构正是针对这类问题的一种解决方法，所以子结构肯定是对一个大型的有限元模型的，譬如在求解非线性问题的时候，因为对于一个非线性问题，系统往往经过多次迭代，每次这个系统的刚度矩阵都会被重新计算，而一般来说一个大型问题往往有很大一部分的变形是很小的，把这部分作为一个子结构，其刚度矩阵仅要计算一次，大大节约了计算时间。

3.那些情况可以使用子结构

前面提到的非线性问题，包括了很小变形的或者线弹性部分可以使用子结构，特别是当模型中有很多相同的部分时，提到的最多的一个例子就是桌子的四条腿，四条腿作为子结构（因为基本时弹性变形）可以包括了很多的实体单元，可以大大提高效率再一个就是问题确实太大，只有采用子结构将问题分成很多块，计算出结果后再次采用子结构分块计算，一直到能对每块单独计算为止。

4.abaqus中子结构的特点及要注意的问题

子结构是一组单元的集合，但是在子结构中仅仅只有你指定的那些节点的自由度会保留下来而其他节点的自由度都被消除了，其他的节点均是通过线性插值的方式获得求解；

子结构是通过你指定的节点与其他的单元建立联系的；在abaqus的6.4版本中只有子结构这个概念没有超单元了，其区别就是子结构可以求得单元内部准确的解而超单元不行；当你定义子结构的时候不要包含太多的单元，因为单元的刚度矩阵集成的时候会花掉太多的时间，可以用更多的含有较少单元的子结构代替。 5.在abaqus中子结构的用法

一般包括如下部分，首先你要利用substructure generate和retained nodal dofs定义一个子结构，然后你可以定义子结构的内部荷载，边界条件一旦你已经定义了子结构以后你就可以象利用一般的单元一样使用子结构了，譬如输出请求等等。因为内容太多只能另外的文章再述了，大家也可以直接看看手册，要注意的是，abaqus cae是不支持子结构的。

如何在不同的分析步改变材料的参数 我所了解的大概有三种不同的方法：

1.最强大的当然是采用umat的方式，不过需要有深厚的有限元基础，一般人不推荐使用 2.采用场变量，不过功能相对简单

3.采用abaqus的import命令将前面分析的结果传递到新的分析之中 这里介绍下第二种方法 *什么是场变量

所谓场变量，我的理解就是一个环境变量，它建立了一个与材料参数之间的中介，虽然不能直接指定材料参数在不同的分析步具有不同的值，但是通过场变量，间接的达到了目的。 *怎样使用场变量

其实场变量用的较多的实在热力学和流体力学的分析种，这里介绍的仅仅是在固体力学中的用法

- 40 -