论坛汇总第一部分
paty之间需要进行asembly的操作吧
还要进行INTERRACTION部分的操作. 在ASSEMBLY中只是位置在一起, 两部分还是没有什么联系, 楼上的，在Interaction里面具体是怎么操作的呢？说具体的好不好啊？ 用INTERACTION中Embedded region 好像可以，可参考帮助文件。 四十一、使用abaqus进行三维地下厂房开挖模拟，有兴趣者请留言！
我正在使用abaqus进行三维地下厂房的开挖模拟，使用joint merterial，drucker-prager模型，模拟开挖以及喷射混凝土变形，有兴趣者请留言
先谢谢了！使用joint material主要是模拟岩石材料的节理（由赤平投影定），使用drucker-prager模型进行弹塑性分析，使用model change模型开挖和混凝土喷射！cafe兄，changqing兄对我建议向你请教：1、能否使用pre-tension 模拟预应力锚杆，对于pre-tension 俺不熟悉，你能否举个简单例子？2、对于断层的模拟有什么较好的方法？谢谢！
使用pre-tension可以模拟预应力锚杆,你可以参考Examples中的manifold这个例子。对于断层的模拟又有很多模式。这要看你的具体要求,另外你是否考虑多条断裂,是否考虑地下水的耦合。这对计算设计有很大影响。个人建议你最好先使用M-C准则算算看。
呵呵，谢谢cafe兄指点！需要模拟的断层不需要考虑水的作用，是一水平断层，没有其它小断层。断层中有方解石脉充填，是否需要用接触来模拟？有什么更好的模拟？请cafe兄指点！
的断层可以使用Interface单元,方解石脉充填带如果不厚的话,可以引入弹簧垫或者直接导入goodman单元来模拟,如果比较后,可以试试节理单元或者使用弱化材料模式。总之,你这个问题的关键不在断层,在充填带的处理上,你可以多做几个模式计算后看结果那个跟实验验证结果接近。 ABAQUS中没有goodman单元,但是可以定制出来。 四十二、在abaqus中如何将几个part merge在一起?
在abaqus中如何将几个part merge在一起呢？ part merge是想做啥？
assembly里, creat instance, （instance就象是part的照片）再把instance之间结合就可以得到你的model了。 将几个part merge在一起，是为了模拟一个体的不同部分具有不同的材料性质！
据我看，还要设成几个part,然后在不同part之间运用*tie，*equation等保证边界变形统一，不能象在ansys里那样有共同边界就自然变形一致了！！
四十三、[讨论]地下隧洞开挖及衬砌计算的一个小疑问 Abaqus Example
1.1.10 Stress-free element reactivation是一个简单的开挖地下隧洞的例子，根据文中叙述： set NLGEOM=no（缺省）时，开挖衬砌后衬砌单元内应力很小； set NLGEOM=yes时，衬砌单元刚加上时有kpa级的应力；但是两种情况最后的计算结果(完全释放开挖荷载)却相差不多！那么实际计算时NLGEOM参数应该如何设置呢？
另外，用marc的单元死活模拟开挖计算了一个实际工程的例子，岩石和混凝土均为弹性材料刚激活的衬砌单元中总是会出现10kpa级的应力，虽然不影响结果，但是同样的模型用abaqus做，刚激活衬砌单元里的应力却很小很小，跟想象的实际情况差不多吧实在想不通为什么会这样？请有经验的网友们帮我分析分析吧。
我想会不会是这样：你得衬砌单元是在开挖之后加上去得吗？若是，那当然应力很小了
用ABAQUS模拟应力释放,没做过,一般有限元程序这点很难处理好.不知ABAQUS可以通过应力释放模拟及时支护还是适时支护,因为,没用过程,也不知是否可以直接加初始应力,若谁知,请讲
刚激活的衬砌单元中总是会出现10kpa级的应力，虽然不影响结果，但是同样的模型用abaqus做，刚激活衬砌单元里的应力却很小很小，跟想象的实际情况差不多吧实在想不通为什么会这样？这个是对的，因为做试验也是很小，但是如果想取较高的可靠度，那么就应该再继续增加材料的强度哦 四十四、abaqus的接触问题
在一个接触问题中，.msg给出一堆诸如以下的警告：
***warning:Contact surface assembly_surf-4 and assembly_surf-3 node number 1225 instance part-7-1 is overclosed by 0.33824 which is too severe--cutback will result.You may want to change the value of hcrit(0.33138) on the *contact pair option.
- 36 -
论坛汇总第一部分
***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来模拟的。 我也是不懂，想了好久也没有好的解决办法，上来集思广益一下。这里高人多。其实我个人觉得先张和后张预应力除了钢筋应力前后变化的不同外，还在预应力的传递上有所不同。广义的预应力应用主要是后张的原理。 比如对如下钢索施加预应力，假定钢筋用杆单元模拟，先张就要保证杆单元（钢索）和块单元（混凝土板）的结点重合以保证预应力的传递；后张的话（应该也是现实中实现的），用于传递预应力的锚具的模拟一直只让我头疼，没有想到较好的、简单的方法实现，所以问问各位的高见。 此主题相关图片如下：
- 39 -