论坛汇总第一部分

在PROPERTY模快的菜单中本身有这项功能，只需直接把任意面定义成蒙皮即可，这跟定义壳截面属性一样，实际上程序内部就是利用实体表面的节点来生成壳单元。 请大侠能否具体谈谈梁和体如何进行多点约束？

SKIN 可以解决\壳、体耦合关系\是不是指壳、体之间无相对滑动,而壳、体之间设置接触面,是否也可以使用蒙皮技术?如果不能,那么壳、体之间设置接触面,自由度是否需要偶合?

壳、体接触问题不用这样解决，只能通过SHELL PART和SOLID PART的接触面来定义，不需要自由度进行藕合 请问楼上的朋友，壳，体之间的耦合用“skin”，如何给“skin”赋材料和截面特性啊？我每次想赋时都无法选中“skin” 与定义壳截面特性一样的，在定义skin时会提示你选择定好的截面的

比如说我建了一个长方体，想把长方体的某个面建成“skin”，但在定义材料特性时只能选择体，而不能选择面，所以就没办法给“skin”赋材料。

你可能对ABAQUS的PROPERTY模快了解的不够，应该没有问题，

如果你对体的property已经有定义,边缘的材料特性也是体的一部分.也就是说你已经定义过嘞,可以简单用1,2等容易数据对skin进行模拟

其实说白了就是你想办法使梁和体在的约束节点处的位移相同就行了，好像和壳差不多，梁体协调，我也没用过。 不是这个意思，我是指体是一个材料，而壳又是另外的材料，如果是这样的，那就没法给体的一个面（设置为“skin”）赋另外的材料了。

如果是不同的维数的part,可以用tie命令找到你需要偶合的部分,用此命令就可以了

在ABAQUS中，看来只能通过不同的PART间的偶合关系来实现，问题是在划分单元是得让PART在交接处的节点基本重合，不知ABAQUS中在划网格时能否很好地控制这一点。 如果几何比较规范，应该自动在边界处协调的

做隧道开挖问题,我们常常用壳单元（三维问题）或梁单元（二维问题）来模拟初期支护结构,有时梁和壳单元与他们相连的实体单元可以处理为完全耦合关系，壳和体有shell to solid couple的选项可以用，我想在壳和体之间设置接触面，不用shell to solid couple选项，可是壳和体之间设置接触面就不收敛，不设接触面收敛，怎么回师？下面是接触定义： ** INTERACTION PROPERTIES **

*Surface Interaction, name=IntProp-1 1.

*Friction, slip tolerance=0.005 0.2,

*Surface Behavior, pressure-overclosure=HARD **

** INTERACTIONS **

** Interaction: Int-1

*Contact Pair,interaction=IntProp-1,ADJUST=0.0 Surf-shell1,Surf-rock1 msg文件里的错误是

***ERROR: TOO MANY ATTEMPTS MADE FOR THIS INCREMENT: ANALYSIS TERMINATED 谁知道怎么回师？ 谢谢！

二十二、〖测试〗〖原创〗关于初始应力计算的一点疑问

此题可以在ANSYS帮助文件的“2.6.13.2. A Sample Initial Stress Problem Using ISFILE”中找到，我这里对初应力文件不是很了解，所以举例想请各位指点一下~ 附件里面的test0.mac表示仅有初应力的工况；test1.mac表示加速度向下作用时，并有初应力的工况；test2.mac表示加速度向上作用时，并有初应力的工况；test.ist是在自重作用下产生的初应力文件，也是程序运行所必需的！ 从位移计算结果可以看出在单纯的初应力作用下（即test0.mac）结构是向上变形的，向上位移是0.379m；在test1.mac作用下位移为0；在test2作用下位

- 21 -

论坛汇总第一部分

移为0.758m。 但是从应力计算结果看出，test0的结果应力最后为零；test1和test2的结果应力却只受外在加速度的影响，而不受初始应力影响。 综合以上的位移和应力分析：位移的结果符合线性叠加原理；而单元和节点的应力却不是这样的，初应力在计算之后马上就释放掉了。问题是有没有办法在计算时，能将初应力和加速度作用下的结果应力一起保存在计算结果中，而不将它释放呢？

二 十三、〖求助〗blast，断裂问题

先看图

取了半个模型，实际计算的时候取了1/4 模型如下图

此主题相关图片如下：

1、在指定seam的时候该选择那个？1线 ？2线？还是其他？我

选择哪个都是不可以的。是不是seam只是针对体内部裂纹的？ 2、crack front说的是什么？ 裂纹拓展面么？图中选择的是那个？

1、seam只是针对体内部裂纹，没错。形态上看3线，最合理。seam是在画网格后才真正出现的界面。获得方法是在实体建好后切一刀。然后指定那一刀为seam。而且，seam不能是将整个实体分为两半的线/面。seam是指实体中在初始太前已有的裂纹上下面，如果裂纹穿透整个实体，显然不合理，意味实体已经完全断裂，没有必要再算了。 2线将整个实体分开了，这样的seam，没有必要分析了1线也就算个裂纹表面，没必要用seam定义用CAE是不能定义动态的扩展裂纹的，裂纹扩展需要直接在inp里自己写

crack front 是指 crack tip裂纹端吧。如果3线是seam，那么s线最上端就是crack front 呵呵，看来我没有把问题说清楚

1、 3下面的线是指引线，不是上面的裂纹。3指示的是交叉点，而不是线。给出这个目的就是为了搞清楚crack front 是不是它。

2、2线是对称面，实际上裂纹是不存在的，裂纹的就是1线构成的那个空缺。2本来是想用来指定为拓展面的。 3、我试过将3定义为裂纹，也就是其中的crack，而se'am是没有的，所以没有指定。但是划分网格的时候，没有得到畸变单元，你知道这个是什么缘故么？

- 22 -

论坛汇总第一部分

如果没定义什么seam，显然3点是应力集中区，也就是crack tip CAE里面不能定义裂纹扩展

在edit crack子窗口下。在那个singularity里面修改midside node parameter就能移动节点到1/4处（设为0.25）在degenerate element。。。。。里面控制裂纹端单元当然，如果选no degeneracy肯定不是奇异的^^b 如果选择了别的，那么在mesh时候就要注意了。挨着tip的那一圈单元要mesh为三角形的。如果没有这么做的话，系统不会自动做出想要的单元

1、 crack的输出不知道怎么选？在output中找不到crack 2、生成的畸变网格不规则，有点乱。 在step里面

最下面两个扭，那个是控制history的。里面能控制crack的输出能输出J\\K\\G\\T 2、 ……这个是分网的技术了，要根据实际情况看

个人的经验是围着crack tip用圆环划出个小区域，区域边均匀设置mesh用的seed。再在小区内mesh，圆域内mesh一层三角形单元

是呀 我就是在history和上面的那个东西里选择了但是 两个都是好用的（对于整个模型） history可以针对整个模型也可以针对个别set surface crack 什么的都可以在step的历史设置里面设定~~毋庸置疑 问题就是在计算的时候，呵呵 它就说没有crack的history

不是很明白……程序会给出没有crack的history的提示？这个我没有遇见过……一般确实有history、fieldout不能输出情况，往往因为单元或材料性质造成的但是简单的crack输出是建立在单元间关系以及应力应变上的，应该不存在这类问题吧~没见过。

返回step，想给crack添加上，呵呵。没有crack了……

这个是操作的失误造成的吧。在CAE里设定好了，存了盘怎么会丢~当然如果建模时侯有错误，可能会令某些属性丢失~不过这样程序应该会给点提示

在采用icem生成好极漂亮的网格导入abaqus后，很多事情都不能办比如说：指定畸变单元不好办，把裂尖的网格指定为wedge网格等等我想我是不是该用hypermesh了，用它事先都设置好，就比较好办了 画网格就是仁者见仁智者见智了顺便提，以前在simwe的hypermesh版查询过奇异单元问题，发现这类问题很少提及。有人提过但是，没人能回答，包括版主。当然我也人为hypermesh能够实现奇异单元，但是自己不会用，不能打保票 cafe斑竹一直喜欢用femap，见过他提到，femap能够实现单元奇异的问题

其实简单的算算是没什么意思的只有复杂的才能现实fem的用处但是复杂的就是要用专业的网格生成工具了 比如icemcfd h'y'p'er'mesh等等接口还是有问题看来必须试试hypermesh了 二十四、〖求助〗关于接触面切向约束问题

在作一个课题时，由于考虑接触面上局部结构的特殊行为，要求不管接触缝面是处于张开还是闭合状态都希望不发生相对的切向滑移。 本来想通过*friction,user中添加Lagrange约束来实现，但该约束仅在接触面闭合时有效；现在写了个uinter子程序来控制切向刚度，但计算时收敛性实在不怎么样，因此特来问问大家有什么更好的思路和办法。 先谢了！

看了看uinter的介绍，似乎该子程序的确可以较好的实现dege兄的目的。但正如UMAT中的DDSDDE一样，DDSDDR似乎对收敛速度影响极大，特别是接触问题涉及到的Discontinuous Iteration，建议dege兄研究一下ABAQUS关于接触问题的算法，或许会有些帮助。

或者利用MPC，对主从面的节点进行约束？？这样可能会导致过约束问题，dege兄可以尝试尝试，有什么结果也可以和大家讨论一下，谢谢！

看看connector element中的rigid beam。这个不会引起过刚度的问题

wjy兄所说的关于DDSDDR的定义确实值得注意，我也正在不断更改和尝试中。做一些小的算例程序都没问题，可是我算的是一大体积混凝土结构在地震动力下非线性响应，似乎还有些困难。 通过mpc给主从面上节点对

- 23 -

论坛汇总第一部分

加切向约束不行，有overconstraint的问题。 xiaohongchen兄对connector elements似乎比较熟悉，能否再给简要谈谈你的rigid beam element想法，你的建议我也会好好读读，谢谢。 别的兄弟有什么好的新建议也望不吝赐教。

二十五〖求助〗接触面上粘聚力c的模拟

最近帮人做一些岩土方面的问题，碰到关于如何考虑沿接触面（如断层、节理面，结构薄弱面等）上的粘聚力c的问题，查看了仿真、傲雪abaqus版块上的精华区和相关文章，没有找到合适的办法。特来请教大家。 个人理解：在初始和未发生张开、滑移阶段，粘聚力作为初始强度发挥作用，而当接触面由完好进入发生张开、滑移阶段，则不予考虑。对于上述力学模型，如何采取一个简单而有效的力学模型进行数值模拟，希望大家给一些提示和帮助。 先谢谢大家。

是否可以这样来考虑，节理层之间用一些弹簧单元来模拟，并设定这种弹簧单元的一个实效应力，该应力也即即节理层展开的最大受力，不知道这样可否。

edodo兄的建议我也考虑了，但是感觉这么一个力学模型对于简单结构在单调荷载作用下的力学行为似乎还可以，复杂应力路径下这么一个弹簧系统模型（包括法向、切向）似乎就不太合适了，且弹簧的等效刚度也不好确定。 除了上述弹簧模型方法，个人还想了这么几类方法： （1）用户定义接触面力学本构（基于uinter）；

(2) 用户定义弱面、节理单元（基于uel）,可能的话与接触面共同使用； （3）基于材料层次去考虑，即将弱面模拟为一层特殊力学行为的材料，

这里也不妨考虑程序本身自带的一些材料本构，如joint material。 请大家帮忙看看。若有别的更好的办法，也请不吝指出。 谢谢！ 用弹簧模拟

复杂问题用最简单办法实现是技术人员应有的本领

此主题相关图片如下：

- 24 -