自己做的一个不同材料坝体接触计算的例子

caoyanhui

这个是自己以前做的一个小例子，为了学习接触计算练手的。
里面的模型参数是根据经验大概拟定的，命令流里面有些细节可能不准确。
希望对初学者有所裨益就行了

shaoyu1979

请问你是用什么做前处理,命令流理的这么清楚,我想不会使用adina-aui做的吧.

另外我发现你的命令流里没有几何信息,而在adina-aui中如果没有几何信息加边界条件和载荷条件包括接触,是很难的,你是如何实现的.

shaoyu1979

有限元网格也非常漂亮,如果你能把如何建立模型的经验详细介绍一下就更好了.建议斑竹加分

caoyanhui

用ANSYS做前处理，然后用APDL命令流方式写成ADINA的输入文件格式，再导入到ADINA。
至于加边界条件和荷载条件可以通过定义结点组（NODESET）和单元面组（ELFACESET）来实现，具体可以参见ADINA的APPLY-LOAD命令的说明。
给一个简单的实例：
1.定义结点组
NODESET 1
ENTRIES   NODE
.
.
.
DATAEND

2.定义单元面组
ELFACESET NAME=1
ENTRIES FACE    EL    GROUP
1  719  1
-1  708  2
1  194  3
-1  95  4
DATAEND

3.施加边界条件
fixity ux
'x-TRANSLATION'
dataend
fixboundary node-set
                1     ux
dataend

4.加外水压力
LOAD PRESSURE NAME=2 MAGNITUDE=9810
APPLY-LOAD
ENTRIES NAME   LTYPE    LNAME  STYPE        SNAME NCUR IDIRN
         2    PRESSURE    2   ELEMENT-FACE     1     2   0
         3    PRESSURE    2   ELEMENT-FACE     2     3   0
         4    PRESSURE    2   ELEMENT-FACE     3     4   0
         5    PRESSURE    2   ELEMENT-FACE     4     5   0
         6    PRESSURE    2   ELEMENT-FACE     5     6   0
         7    PRESSURE    2   ELEMENT-FACE     6     7   0
         8    PRESSURE    2   ELEMENT-FACE     7     8   0
         9    PRESSURE    2   ELEMENT-FACE     8     9   0
        10    PRESSURE    2   ELEMENT-FACE     9    10   0
        11    PRESSURE    2   ELEMENT-FACE    10    11   0
        12    PRESSURE    2   ELEMENT-FACE    11    12   0
        13    PRESSURE    2   ELEMENT-FACE    12    13   0
        14    PRESSURE    2   ELEMENT-FACE    13    14   0
        15    PRESSURE    2   ELEMENT-FACE    14    15   0
        16    PRESSURE    2   ELEMENT-FACE    15    16   0
        17    PRESSURE    2   ELEMENT-FACE    16    17   0
        18    PRESSURE    2   ELEMENT-FACE    17    18   0
        19    PRESSURE    2   ELEMENT-FACE    18    19   0
        20    PRESSURE    2   ELEMENT-FACE    19    20   0
        21    PRESSURE    2   ELEMENT-FACE    20    21   0
        22    PRESSURE    2   ELEMENT-FACE    21    22   0
        23    PRESSURE    2   ELEMENT-FACE    22    23   0
        24    PRESSURE    2   ELEMENT-FACE    23    24   0
DATAEND
如果各位有更好的方法，欢迎交流。

ahali2002

希望谁能提供一个分层填筑（逐级加载）的例子。

ahali2002

我做了一个验证E-B模型的简单例子，模拟土石坝分层填筑，本来想用单元生死来实现，可是计算结果不理想，特别是竖向沉降，不知道是单元生死设置时有问题还是E-B本构有问题，请大家诊断。

shaoyu1979

2.定义单元面组
ELFACESET NAME=1
ENTRIES FACE EL GROUP
1  719  1
-1  708  2
1  194  3
-1  95  4
DATAEND

这个是不是很麻烦，因为单元要一个一个选，而且，你还不知道单元的某个面编号。

caoyanhui

答上楼，这个确实有点麻烦，我是通过以下方法定义面组的：
一是通过在ansys中模型的一一对应来建立组合（该模型本就是通过ansys导入的），
二是通过在adina中define element face set(菜单中有)在模型上点选想要的面，则在条栏中自然就出现了选中的面号、单元号、所属group号，因为adina前处理中选择很不方便，这里点选的时候要比较小心，很容易选错的。

csloy

谢谢曹老大，正在学习中。土坝加固时，在原有土坝的基础上做一80cm厚的防渗墙，请问防渗墙和坝体材料的接触如何做？可以用分层加载和接触单元也采用生死来模拟吗？